The flow characteristics of the subcritical and supercritical flows over a free overfall in a rectangular channel are studied experimentally. A series of experiments were conducted in a tilting flume with a wide range of flow rate. Data collected by several researchers are also included. An empirical relationship, which gives the flow rate as a function of the brink depth, the channel bed slope and the bed roughness are confirmed by using data collected in present study. In addition, the behaviors of the ratio of the brink depth to the critical depth according to several flow parameters are examined. Further, the location of the critical depth in subcritical flows while flow is approaching to the fall is investigated. It is concluded that the location of the critical depth in subcritical flow is a function of the Froude number, channel bed slope and the Manning roughness coefficient. The resemblance or the difference in the occurrence of the profile in sub and supercritical flows examined.

This study concentrates on the modelling of three phase flow and miscible CO2 flooding in pore networks that captures the natural porous medium of a reservoir. That is to say, the network, that is a Matlab code, consists of different sided triangles which are located randomly through the grids. The throats that connect the pores are also created by the model. Hence, the lengths and the radii of the throats are varying. The network used in this research is assumed to be representative of mixed-wet carbonates in 2-D. Mixed wettability arises in real porous media when oil renders surfaces it comes into prolonged contact with oil-wet while water-filled nooks and crannies remain water-wet. The model developed is quasi-static approach to simulate two phase and three phase flows. By this, capillary pressures, relative permeabilities, saturations, flow paths are determined for primary drainage, secondary imbibition, and CO2 injection cases. To calculate the relative permeability, capillary entry pressures are first determined. Then, hydraulic conductances and flow rates of the network for each grid are obtained. Phase areas and saturations are also determined. It is accepted that the displacement mechanism in drainage and CO2 injection is piston-like whereas in imbibition it is either piston-like or snap-off. The results of the model are compared with the experimental data from the literature. Although, the pore size distribution and the contact angle of the model are inconsistent with the experimental data, the agreement of the relative permeabilities is promising. The effect of contact angle in the same network for three phase flow where immiscible CO2 is injected as a third phase at supercritical temperature (32 &deg / C) is investigated. And it is found that, the increase in the intrinsic angles causes decrease in relative permeability values. As another scenario, two phase model is developed in which miscible CO2 &amp / #8211 / water is flooded after the primary drainage of the same 2-D network at supercritical temperature (32 &deg / C). This case is compared with the previous case and the effects of miscibility are investigated such that it causes the relative permeability values to increase. Adsorption is another concern of which its effects are analyzed in a single pore model. The model is compared with the reported experimental data at high temperature and pressures. A reasonable fit is obtained.

ABSTRACT GENERALIZED COUETTE FLOW OF A HERSCHEL - BULKLEY FLUID THROUGH ECCENTRIC ANNULUS - AN APPROXIMATE SOLUTION Seyidoglu, Tijen M.S., Department of Chemical Engineering Supervisor: ismail Tosun Co-Supervisor: Ahmet N. Eraslan January 2006, 134 pages Generalized Couette flow of a Herschel-Bulkley fluid in an eccentric annulus is analyzed by approximating the flow geometry as a slit of variable height. Besides an imposed pressure gradient, one of the plates is considered non-stationary to take into account the axial and/or angular motion of the inner pipe in an eccentric annulus system. Depending on the magnitude and the direction of the applied pressure gradient with respect to the plate velocity, three separate flow cases are studied in which the velocity reaches its maximum value either within the plug flow region or at the moving boundary. Velocity distributions are obtained for each case by solving the equations of continuity and motion. Volumetric flow rate expressions are obtained by integrating the velocity distribution over the cross-sectional area. At a given pressure gradient, the results indicate an increase in volumetric flow rate with an increase in eccentricity ratio. Criteria for each flow type is developed in terms of a dimensionless parameter &amp / #923 / , which takes into account the ratio of the imposed pressure gradient to the plate velocity. Volumetric flow rate expressions for Newtonian, Bingham and power-law fluids are obtained by considering the limiting values of the fluid index and yield stress. The validity of the equations are checked by considering the slit height to be a constant, i.e., flow between parallel plates. Surge/swab pressure calculations are carried out for Herschel-Bulkley, power-law and Bingham fluids and the results are expressed as a function of eccentricity ratio, radius ratio, fluid index and yield stress. The results indicate that when the fluid index and the eccentricity ratio are fixed, a slight increase in the radius ratio causes a tremendous increase in surge/swab pressure, especially for low values of fluid index. On the other hand, displacement of the inner pipe from a concentric position causes a decrease in swab/surge pressure when other parameters are held constant. Comparison with the literature values reveals the fact that flow in an eccentric annulus can be modeled as flow between a slit of variable height as long as the radius ratio is greater than 0.5 and the eccentricity ratio is less than 0.3. The results for other values of radius ratio and eccentricity ratio can be used as initial guess values in carrying out numerical calculations.

In this study, a numerical model is developed to determine shoreline changes due to wind wave induced longshore sediment transport, by solving sediment continuity equation and taking one line theory as a base, in existence of seawalls, groins, T-groins, offshore breakwaters and beach nourishment projects, whose dimensions and locations may be given arbitrarily. The model computes the transformation of deep water wave characteristics up to the surf zone and eventually gives the result of shoreline changes with user-friendly visual outputs. A method of representative wave input as annual average wave characteristics is presented. Compatibility of the currently developed tool is tested by a case study and it is shown that the results, obtained from the model, are in good agreement qualitatively with field measurements. In the scope of this study, input manner of long term annual wave data into model in miscellaneous ways is also discussed.

Earthquakes are among natural disasters with significant damage potential / however it is possible to reduce the losses by taking several remedies. Reduction of seismic losses starts with identifying and estimating the expected damage to some accuracy. Since both the design styles and the construction defects exhibit mostly local properties all over the world, damage estimations should be performed at regional levels. Another important issue in disaster mitigation is to determine a robust measure of ground motion intensity parameters. As of now, well-built correlations between shaking intensity and instrumental ground motion parameters are not yet studied in detail for Turkish data. In the first part of this thesis, regional empirical Damage Probability Matrices (DPMs) are formed for Turkey. As the input data, the detailed damage database of the 17 August 1999 Kocaeli earthquake (Mw=7.4) is used. The damage probability matrices are derived for Sakarya, Bolu and Kocaeli, for both reinforced concrete and masonry buildings. Results are compared with previous similar studies and the differences are discussed. After validation with future data, these DPMs can be used in the calculation of earthquake insurance premiums. In the second part of this thesis, two relationships between the felt-intensity and peak ground motion parameters are generated using linear least-squares regression technique. The first one correlates Modified Mercalli Intensity (MMI) to Peak Ground Acceleration (PGA) whereas the latter one does the same for Peak Ground Velocity (PGV). Old damage reports and isoseismal maps are employed for deriving 92 data pairs of MMI, PGA and PGV used in the regression analyses. These local relationships can be used in the future for ShakeMap applications in rapid response and disaster management activities.

Bolu-Ilgaz region was damaged by several large earthquakes in the last century and the structural damage was substantial especially after the 1944 and 1999 earthquakes. The objective of this study is to build the seismic source characterization model for the rupture zone of 1944 Bolu-Gerede earthquake and perform probabilistic seismic hazard assessment (PSHA) in the region. One of the major improvements over the previous PSHA practices accomplished in this study is the development of advanced seismic source models in terms of source geometry and reoccurrence relations. Geometry of the linear fault segments are determined and incorporated with the help of available fault maps. Composite magnitude distribution model is used to properly represent the characteristic behavior of NAF without an additional background zone. Fault segments, rupture sources, rupture scenarios and fault rupture models are determined using the WG-2003 terminology. The Turkey-Adjusted NGAW1 (G&uuml / lerce et al., 2013) prediction models are employed for the first time on NAF system. The results of the study is presented in terms of hazard curves, deaggregation of the hazard and uniform hazard spectrum for four main locations in the region to provide basis for evaluation of the seismic design of special structures in the area. Hazard maps of the region for rock site conditions and for the proposed site characterization model are provided to allow the user perform site-specific hazard assessment for local site conditions and develop site-specific design spectrum. The results of the study will be useful to manage the future seismic hazard in the region.

Slurry walls have been widely used for more than 25 years to control the migration of contaminants in the subsurface. In the USA, vertical barriers are mostly constructed of soil-bentonite using the slurry trench method of construction. In this method, sodium bentonite is mixed with water to form a viscous slurry that is pumped into a trench during excavation to maintain the trench stability. The stable trench is then backfilled with a mixture of soil and slurry having a consistency of high slump concrete. These barriers have been designed primarily for low permeability, generally less than 10&minus / 9 m/s. Some investigations have pointed toward improved performance using admixtures that would provide low permeability. In this study, Soma thermal power plant fly ash, granulated blast furnace slag, lime, and steel slag are used as admixture to improve the performance of slurry walls. Permeability, compressive strength, slump, compressibility properties of the mixtures were found and checked for the minimum requirements. According to the findings of this study, granulated blast furnace slag (GGBS), fly ash and steel slag can be used at certain percentages and curing periods as additive in cement-bentonite barrier wall construction. Permeability of specimens having fly ash decreases by increasing fly ash content. Mixtures having 50 % of GGBS type I with 5 % of lime and 9% bentonite content gave acceptable results in 28 days of curing time. Specimens including 50 % of GGBS type II with 5 % of lime and 9% bentonite content gave the higher permeability value in 28 days of curing time with respect to GGBS type I. In addition, most of the mixtures prepared by steel slag gave the acceptable permeability values in 28 days of curing period. Unconfined compressive strength of all mixtures increase by increasing curing time. Cc, Cr, Cv, kcon values were found from consolidation test results. Permeability values found from consolidation tests are 10 times to 100 times higher than flexible wall k results for the same effective stress of 150 kPa. Generally, mv values are decreasing with increasing curing time. As mv decreases, D increases.

With continuously increasing urbanization, consumer demands and expansion of water supply systems, determination of efficient pump schedules became a more difficult task. Pumping energy costs constitute a significant part of the operational cost of the water distribution networks. This study aims to provide an effective daily pump schedule by minimizing the energy costs for constant and also for multi tariff of electricity (3 Kademeli Elektrik Tarifesi) in water distribution network. A case study has been performed in an area covering N8.3 and N7 pressure zones which are parts of Ankara water distribution network. Both pressure zones consists of 3 multiple pumps in pump station and one tank having 5000 m3 storage volume each. By using genetic algorithm based software (WaterCAD Darwin Scheduler) least-cost pump scheduling and operation policy for each pump station has been determined while satisfying target hydraulic performance requirements such as minimum and maximum service pressures, final water level of storage tank and maximum velocity in pipeline. 32 different alternative scenarios have been created which include multi tariff energy prices, constant tariff energy price, insulated system condition, uninsulated system condition and different pump combinations. The existing base scenario and alternative scenarios which were prepared by using optimal pump schedules have been compared and the achievements of optimizing pump operation have been analyzed. At the end of the study, a satisfying result has been observed that by using determined optimal pump schedule, minimum % 14 of total energy cost can be saved in existing water supply system.

Self-compacting concrete is an innovative construction material / its priority to normal vibrated concrete is that there is not any vibration requirement. Bond strength of reinforcement is one of the key factors that ensures the usefulness of any reinforced concrete structure. In this study, 6 full-scale concrete beams spliced at the mid-span were tested under two-point symmetrical loading. Test variables were bottom cover, side cover, free spacing between longitudinal reinforcement, lap-splice length and presence of transverse reinforcements within the lap-splice region. Specimen SC_22_44_88_800 had cover dimensions close to the code limits and had 36db lap splice length. This specimen showed flexural failure. Specimen SC_44_44_44_710 had 32db lap splice and cover dimensions greater than code minimums. This specimen showed yielding primarily. With the increasing loading, however, bond failure occurred with side splitting. ACI 408 descriptive equation for normal vibrated concrete predicted bar stresses of the unconfined specimens produced with self-compacting concrete acceptably well. The predicted values were lower than the measured values to be on the safe side. The error varied between 3.4% and 6.5%. All predictions of the ACI408 descriptive equation was higher than the measured bar stresses of the confined specimens produced with SCC. All the calculated values were unsafe. The error varied between 10.6% and 34.5%. Specimen SC_44_22_22_530_T4 with 24db lap splice length had side cover and spacing between bars 63.3% and 56.7% less than the ACI 318 limits. The calculated bar stress was 21.6% higher than the measured value. The main reason of the deviation was inadequate cover dimensions. In specimen SC_44_22_22_530_T6, number transverse reinforcement was increased to 6 stirrups to overcome the small cover and spacing problem. However, increased number of stirrups inside a small side and face cover caused weak plane and measured bar stress decreased.

Water distribution network models are widely used for various purposes such as long-range planning, design, operation and water quality management. Before these models are used for a specific study, they should be calibrated by adjusting model parameters such as pipe roughness values and nodal demands so that models can yield compatible results with site observations (basically, pressure readings). Many methods have been developed to calibrate water distribution networks. In this study, Darwin Calibrator, a computer software that uses genetic algorithm, is used to calibrate N8.3 pressure zone model of Ankara water distribution network / in this case study the network is calibrated on the basis of roughness parameter, Hazen Williams coefficient for the sake of simplicity. It is understood that there are various parameters that contribute to the uncertainties in water distribution network modelling and the calibration process. Besides, computer software&rsquo / s are valuable tools to solve water distribution network problems and to calibrate network models in an accurate and fast way using automated calibration technique. Furthermore, there are many important aspects that should be considered during automated calibration such as pipe roughness grouping. In this study, influence of flow velocity on pipe roughness grouping is examined. Roughness coefficients of pipes have been estimated in the range of 70-140.