The onset of summer stratification in temperate lakes and reservoirs forces a decoupling of the hypolimnion from the epilimnion that is sustained by strong density gradients in the metalimnion. These strong gradients act as a barrier to the vertical transport of mass and scalars leading to bottom anoxia and subsequent nutrient release from the sediments. The stratification is intermittently overcome by turbulent mixing events that redistribute mass, heat, dissolved parameters and particulates in the vertical. The redistribution of ecological parameters then exerts some control over the ecological response of the lake. This dissertation is focused on the physics of deep vertical mixing that occurs beneath the well-mixed surface layer in stratified lakes and reservoirs. The overall aim is to improve the ability of numerical models to reproduce deep vertical mixing, thus providing better tools for water quality prediction and management. In the first part of this research the framework of a one-dimensional mixed-layer hydrodynamic model was used to construct a pseudo two-dimensional model that computes vertical fluxes generated by deep mixing processes. The parameterisations developed for the model were based on the relationship found between lake-wide vertical buoyancy flux and the first-order internal wave response of the lake to surface wind forcing. The ability of the model to reproduce the observed thermal structure in a range of lakes and reservoirs was greatly improved by incorporating an explicit turbulent benthic boundary layer routine. Although laterally-integrated models reproduce the net effect of turbulent mixing in a vertical sense, they fail to resolve the transient distribution of turbulent mixing events triggered by local flow properties defined at far smaller scales. Importantly, the distribution of events may promote tertiary motions and ecological niches. In the second part of the study a large body of microstructure data collected in Lake Kinneret, Israel, was used to show that the nature of turbulent mixing events varied considerably between the epilimnion, metalimnion, hypolimnion and benthic boundary layer, yet the turbulent scales of the events and the buoyancy flux they produced collapsed into functions of the local gradient Richardson number. It was found that the most intense events in the metalimnion were triggered by high-frequency waves generated near the surface that grew and imparted a strain on the metalimnion density field, which led to secondary instabilities with low gradient Richardson numbers. The microstructure observations suggest that the local gradient Richardson number could be used to parameterise vertical mixing in coarse-grid numerical models of lakes and reservoirs. However, any effort to incorporate such parameterisations becomes meaningless without measures to reduce numerical diffusion, which often dominates over parameterised physical mixing. As a third part of the research, an explicit filtering tool was developed to negate numerical diffusion in a threedimensional hydrodynamic model. The adaptive filter ensured that temperature gradients in the metalimnion remained within bounds of the measured values and so the computation preserved the spectrum of internal wave motions that trigger diapycnal mixing events in the deeper reaches of a lake. The results showed that the ratio of physical to numerical diffusion is dictated by the character of the dominant internal wave motions.

Hydrocarbon reserves are generally produced through wells drilled into reservoir pay zones. During production, gas liberation from the oil phase occurs due to pressure decline in the wellbore. Thus, we expect multiphase flow in some sections of the wellbore. As a multi-phase/multi-component gas-oil mixture flows from the reservoir to the surface, pressure, temperature, composition, and liquid holdup distributions are interrelated. Modeling these multiphase flow parameters is important to design production strategies such as artificial lift procedures. A wellbore fluid flow model can also be used for pressure transient test analysis and interpretation. Considering heat exchange in the wellbore is important to compute fluid flow parameters accurately. Modeling multiphase fluid flow in the wellbore becomes more complicated due to heat transfer between the wellbore fluids and the surrounding formations. Due to mass, momentum, and energy exchange between the wellbore and the reservoir, the wellbore model should be coupled with a numerical reservoir model to simulate fluid flow accurately. This model should be non-isothermal to consider the effect of temperature. Our research shows that, in some cases, ignoring compositional effects may lead to errors in pressure profile prediction for the wellbore. Nearly all multiphase wellbore simulations are currently performed using the "black oil" approach. The primary objective of this study was to develop a non-isothermal wellbore simulator to model transient fluid flow and temperature and couple the model to a reservoir simulator called General Purpose Adaptive Simulator (GPAS). The coupled wellbore/reservoir simulator can be applied to steady state problems, such as production from, or injection to a reservoir as well as during transient phenomena such as well tests to accurately model wellbore effects. Fluid flow in the wellbore may be modeled either using the blackoil approach or the compositional approach, as required by the complexity of the fluids. The simulation results of the new model were compared with field data for pressure gradients and temperature distribution obtained from wireline conveyed pressure recorder and acoustic fluid level measurements for a gas/oil producer well during a buildup test. The model results are in good agreement with the field data. Our simulator gave us further insights into the wellbore dynamics that occur during transient problems such as phase segregation and counter-current multiphase flow. We show that neglecting these multiphase flow dynamics would lead to unreliable results in well testing analysis.

No description available.

No description available.

The instabilities of shallow shear flows are analyzed to study exchanges processes across shear flows in inland and coastal waters, coastal and ocean currents, and winds across the thermal-and-moisture fronts. These shear flows observed in nature are driven by gravity and governed by the shallow water equations (SWE). A highly accurate, and robust, computational scheme has been developed to solve these SWE. Time integration of the SWE was carried out using the fourth-order Runge-Kutta scheme. A third-order upwind bias finite difference approximation known as QUICK (Quadratic Upstream Interpolation of Convective Kinematics) was employed for the spatial discretization. The numerical oscillations were controlled using flux limiters for Total Variation Diminishing (TVD). Direct numerical simulations (DNS) were conducted for the base flow with the TANH velocity profile, and the base flow in the form of a jet with the SECH velocity profile. The depth across the base flows was selected for the' balance of the driving forces. In the rotating flow simulation, the Coriolis force in the lateral direction was perfectly in balance with the pressure gradient across the shear flow during the simulation. The development of instabilities in the shear flows was considered for a range of convective Froude number, friction number, and Rossby number. The DNS of the SWE has produced linear results that are consistent with classical stability analyses based on the normal mode approach, and new results that had not been determined by the classical method. The formation of eddies, and the generation of shocklets subsequent to the linear instabilities were computed as part of the DNS. Without modelling the small scales, the simulation was able to produce the correct turbulent spreading rate in agreement with the experimental observations. The simulations have identified radiation damping, in addition to friction damping, as a primary factor of influence on the instability of the shear flows admissible to waves. A convective Froude number correlated the energy lost due to radiation damping. The friction number determined the energy lost due to friction. A significant fraction of available energy produced by the shear flow is lost due the radiation of waves at high convective Froude number. This radiation of gravity waves in shallow gravity-stratified shear flow, and its dependence on the convective Froude number, is shown to be analogous to the Mach-number effect in compressible flow. Furthermore, and most significantly, is the discovery from the simulation the crucial role of the radiation damping in the development of shear flows in the rotating earth. Rings and eddies were produced by the rotating-flow simulations in a range of Rossby numbers, as they were observed in the Gulf Stream of the Atlantic, Jet Stream in the atmosphere, and various fronts across currents in coastal waters.

No description available.

Thesis (MScEng (Mathematical Sciences. Applied Mathematics))--University of Stellenbosch, 2006. / An analytical modelling procedure for predicting the streamwise pressure gradient for steady laminar incompressible flow of a Newtonian fluid through homogeneous isotropic granular porous media is introduced. The modelling strategy involves the spatial volume averaging of a statistical representative portion of the porous domain to obtain measurable macroscopic quantities from which macroscopic transport equations can be derived. A simple pore-scale model is introduced to approximate the actual complex granular porous microstructure through rectangular cubic geometry. The sound physical principles on which the modelling procedure is based avoid the need for redundant empirical coefficients. The model is generalized to predict the rheological flow behaviour of non-Newtonian purely viscous power law fluids by introducing the dependence of the apparent viscosity on the shear rate through the wall shear stress. The field of application of the Newtonian model is extended to predict the flow behaviour in fluidized beds by adjusting the Darcy velocity to incorporate the relative velocity of the solid phase. The Newtonian model is furthermore adjusted to predict fluid flow through Fontainebleau sandstone by taking into account the effect of blocked throats at very low porosities. The analytical model as well as the model generalizations for extended applicability is verified through comparison with other analytical and semi-empirical models and a wide range of experimental data from the literature. The accuracy of the predictive analytical model reveals to be highly acceptable for most engineering designs.

It was generally accepted that the rate of formation of cerebrospinal °uid (CSF) is independent of intraventricular pressure [26], until A. Sahar and a host of other scientists challenged this belief. A. Sahar substantiated his belief that the rate of (CSF) formation actually depends on intraventricular pressure, see A. Sahar, 1971 [26]. In this work we show that CSF formation depends on some other factors, including the intraventricular pressure. For the purpose of this study, we used the capillary blood °ow model proposed by K.Boryczko et. al., [5] in which blood °ow in the microvessels was modeled as a two-phase °ow; the solid and the liquid volume phase. CSF is formed from the blood plasma [23] which we assume to be in the liquid volume phase. CSF is a Newtonian °uid [2, 23]. The principles and methods of e®ective area" developed by N. Sauer and R. Maritz [21] for studying the penetration of °uid into permeable walls was used to investigate the ¯ltrate momentum °ux from the intracranial capillary wall through the pia mater and epithelial layer of the choroid plexus into the subarachnoid space. We coupled the dynamic boundary equation with the Navier-Stoke's constitutive equation for incompressible °uid, representing the °uid °ow in the liquid volume phase in the capillary to arrive at our model. / Mathematical sciences / M.Sc.

A numerical model has been developed for the study of stratified tidal circulation in Indian Arm - a representative inlet on the southern coast of British Columbia. Equations for horizontal velocity, salt conservation, continuity, density (calculated as a linear function of salinity), and the hydrostatic approximation govern the dynamics. All equations have been laterally integrated under the assumption of negligible cross-inlet variability. The model is time dependent and includes nonlinear advective terms, horizontal and vertical turbulent diffusion of salt and momentum, and variations in width and depth. Provisions for surface wind stress and a flux of fresh water are also included, although neither was utilized in this study. An explicit finite difference scheme centred in both time and space was used to solve for the horizontal and vertical velocity components, salinity, and surface elevation on a staggered rectangular grid. A backward Euler scheme was used to suppress the computational mode. Tests using a semi-implicit scheme to solve the finite difference equations over realistic topography led to numerical instabilities at modest values of the time step - in spite of the unconditional stability criteria - suggesting that linear stability analysis may give misleading results for strongly nonlinear systems. Surface elevations calculated from tidal harmonic analysis and salinity timeseries derived from linearly interpolated CTD casts were prescribed at the open boundary. Initial and boundary conditions based on observations in Burrard Inlet and Indian Arm during the winter of 1974-75 were used to study the inlet's response to tidal forcing and to simulate the deep-water renewal that occurred during this period. Coefficients for the horizontal turbulent diffusion of momentum and salt were set equal to 10⁶ cm² s⁻¹. Reducing this value by a factor of two was found to have little impact on the solution. A further reduction to 10³ cm² s⁻¹ led to numerical instabilities under conditions of dense water inflow. The side friction term in the momentum balance was tuned to match calculated and observed dissipation rates in Burrard Inlet; leading to good agreement between the observed and calculated barotropic tide. Contour plots of tidal amplitudes and phases for model currents and salinities revealed a standing wave pattern for the K₁ and M₂ internal tides in Indian Arm; thus allowing for the possibility of resonance. A comparison of model results with vertical amplitude and phase profiles from harmonic analysis of Cyclesonde current meter timeseries at two locations in Indian Arm was consistent with this result. A least-squares fit was made of the vertical modal structure in the model to the complex tidal amplitudes. This led to calculations of the kinetic energy contained in each of the modes along the model inlet for the M₂ and K₁ constituents. Most of the energy was found to be contained in the barotropic and first baroclinic modes, with the latter dominating in the deep basin, and the former dominating near the sill. Second mode energy was significant for the K₁ constituent at some locations in Indian Arm. There are clear indications in the model of barotropic tidal energy being radiated into the inlet basin via the internal tide. Simulations of the influx of dense water into Indian Arm yielded exchange rates that are consistent with observed values and suggest the possibility of fine-tuning the model coefficients to allow prediction of future overturning events. / Science, Faculty of / Earth, Ocean and Atmospheric Sciences, Department of / Graduate

Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Nuclear Engineering, 1980. / MICROFICHE COPY AVAILABLE IN ARCHIVES AND SCIENCE / Includes bibliographical references. / by John Edward Kelly. / Ph.D.