This study presents the numerical analysis of steady laminar flow heat transfer in a horizontal rectangular duct with differential heating on the vertical walls. Three heating configurations: one uniform wall temperature (CS1) and two linearly varying wall temperature cases (CS2 and CS3) are analysed. The study considers the combined effects of natural convection, forced convection and radiation heat transfer on the overall heat transfer characteristics. Air, which is assumed to be a non-participating medium, is chosen as the working fluid. A computational fluid dynamics solver is used to solve a set of governing equations for a range of parameters.For chosen duct aspect ratios, the numerical model simulates the flow and heat transfer for two main effects: buoyancy and radiation heat transfer. Buoyancy effect is represented by Grashof number, which is varied from 2,000 to 1,000,000. The effect of radiation heat transfer is examined by choosing different wall surface emissivity values. The weak and strong radiation effect is represented by the emissivity values of 0.05 and 0.85 respectively. Three duct aspect ratios are considered - 0.5, 1 and 2. The heat transfer characteristics of all the above heating configurations - CS1, CS2, and CS3 are analysed and compared. The numerical results show that, for all heating configurations and duct aspect ratios, the overall heat transfer rate is enhanced when the buoyancy effect increases. Since buoyancy effect induces natural circulation, this circulation is therefore the main mechanism that enhances heat transfer. Radiation heat transfer is found to significantly influence convection heat transfer in high Grashof numbers.

Plumes are formed when a continuous buoyant forcing is supplied at a localized source. Buoyancy can be created by either a heat flux, a compositional difference between the fluid coming from the source and its surroundings, or a combination of both. In this thesis, two types of laminar plumes with different forcing mechanisms were investigated: forced plumes and autocatalytic plumes. The forced plumes were compositionally buoyant and were injected with inertial forcing into a fluid filled tank. The autocatalytic plumes were produced without mechanical forcing by buoyancy that was entirely the consequence of a nonlinear chemical reaction -- the iodate-arsenous acid (IAA) reaction. This reaction propagates as a reacting front and produces buoyancy by its exothermicity, and by the compositional difference between the reactant and product. Both the forced and autocatalytic plumes were examined in starting and steady states. The starting, or transient, state of the plume occurs when it initially rises through a fluid and develops a plume head on top of a trailing conduit. The steady state emerges after the plume head has risen to the top of a fluid filled tank leaving only a persistent conduit. Plume behaviour was studied through experimentation, simulation, and by using simple theoretical analysis. We performed the first ever study of plumes as they crossed over the transition between buoyancy-driven to momentum-driven flow. Regardless of the driving mechanism, forced plumes were found to exhibit a single power law relationship that explains their ascent velocity. However, the morphology of the plume heads was found to depend on the dominating driving mechanism. Confined heads were produced by buoyancy-driven plumes, and dispersed heads by momentum-driven plumes. Autocatalytic plumes were found to have rich dynamics that are a consequence of the interplay between fluid flow and chemical reaction. These plumes produced accelerating heads that detached from the conduit, forming free vortex rings. A second-generation head would then develop at the point of detachment. The detachment process for plumes was sensitively dependent on small fluctuations in their initial formation. In some cases, head detachment could occur multiple times for a single experimental run, thereby producing several generations of autocatalytic vortex rings. Head detachment was reproduced and studied using autocatalytic plume simulations. Autocatalytic flame balls, a phenomenon closely related to autocatalytic plumes, were also simulated. Flame balls were found to have three dynamical regimes. Below a critical radius, the smallest flame balls experienced front death. Above this radius, they formed elongating, reacting tails. The largest flame balls formed filamentary tails unable to sustain a reaction.

Plumes are formed when a continuous buoyant forcing is supplied at a localized source. Buoyancy can be created by either a heat flux, a compositional difference between the fluid coming from the source and its surroundings, or a combination of both. In this thesis, two types of laminar plumes with different forcing mechanisms were investigated: forced plumes and autocatalytic plumes. The forced plumes were compositionally buoyant and were injected with inertial forcing into a fluid filled tank. The autocatalytic plumes were produced without mechanical forcing by buoyancy that was entirely the consequence of a nonlinear chemical reaction -- the iodate-arsenous acid (IAA) reaction. This reaction propagates as a reacting front and produces buoyancy by its exothermicity, and by the compositional difference between the reactant and product. Both the forced and autocatalytic plumes were examined in starting and steady states. The starting, or transient, state of the plume occurs when it initially rises through a fluid and develops a plume head on top of a trailing conduit. The steady state emerges after the plume head has risen to the top of a fluid filled tank leaving only a persistent conduit. Plume behaviour was studied through experimentation, simulation, and by using simple theoretical analysis. We performed the first ever study of plumes as they crossed over the transition between buoyancy-driven to momentum-driven flow. Regardless of the driving mechanism, forced plumes were found to exhibit a single power law relationship that explains their ascent velocity. However, the morphology of the plume heads was found to depend on the dominating driving mechanism. Confined heads were produced by buoyancy-driven plumes, and dispersed heads by momentum-driven plumes. Autocatalytic plumes were found to have rich dynamics that are a consequence of the interplay between fluid flow and chemical reaction. These plumes produced accelerating heads that detached from the conduit, forming free vortex rings. A second-generation head would then develop at the point of detachment. The detachment process for plumes was sensitively dependent on small fluctuations in their initial formation. In some cases, head detachment could occur multiple times for a single experimental run, thereby producing several generations of autocatalytic vortex rings. Head detachment was reproduced and studied using autocatalytic plume simulations. Autocatalytic flame balls, a phenomenon closely related to autocatalytic plumes, were also simulated. Flame balls were found to have three dynamical regimes. Below a critical radius, the smallest flame balls experienced front death. Above this radius, they formed elongating, reacting tails. The largest flame balls formed filamentary tails unable to sustain a reaction.

This study investigates the hydrodynamic variations in the coastal region of Nan-bin, Hualien influence by the topography and water stratification. The data used were collected from four cruises of field observations in 2008/4/29¡B2008/9/5~6¡B2009/3/17~18 and 2009/7/21. Instruments applied include sb-ADCP and CTD. Parameters recorded include flow velocities, water temperature and salinity and tidal elevations. The collected data are analyzed through a variety of time series analysis technique, including buoyancy frequency EOF analysis, potential vorticity and kinetic energy. The results show that (1) the flow field and vorticity in the south side of Hualien Harbor indicate there exist a counterclockwise eddy during flood. The flow pattern reverse to be a clockwise eddy during ebb. The current speed and vorticity is smaller in the surface layer, while the bottom current speed and vorticity is much larger. (2) Currents in transects of along slope, along canyon and along shore all reveal two layers flow separated around 25m, with upper layer flowing eastward and lower layer westward. The current velocity reaches to a maximum at 55m, and the current directions were dominated by the orientation of bottom topography. The baroclinic kinetic energy increase to 0.15m2/s2 along the down slope transect, while the up slope kinetic energy reduced to 0.05m2/s2. (3) The CTD data suggest that the upper layer is warm-fresh water due to river outputs, while the lower layer is dominated by cold-salty open ocean water. (4) The density interface at 25 m is confirmed by buoyance frequency analysis. The EOF analysis of density distribution shows eigenvalue of first mode has maximum at 53m, while the second mode has peaks at 25m and 65m, which match well with the vertical of distribution of flow structures.

The primary objectives of this research were to understand the potential for gravity-stable surfactant floods for enhanced oil recovery without the need for mobility control agents and to optimize the performance of such floods. Surfactants are added to injected water to mobilize the residual oil and increase the oil production. Surfactants reduce the interfacial tension (IFT) between oil and water. This reduction in IFT reduces the capillary pressure and thus the residual oil saturation, which then results in an increase in the water relative permeability. The mobility of the surfactant solution is then greater than the mobility of the oil bank it is displacing. This unfavorable mobility ratio can lead to hydrodynamic instabilities (fingering). The presence of these instabilities results in low reservoir sweep efficiency. Fingering can be prevented by increasing the viscosity of the surfactant solution or by using gravity to stabilize the displacement below a critical velocity. The former can be accomplished by using mobility control agents such as polymer or foam. The latter is called gravity-stable surfactant flooding, which is the subject of this study. Gravity-stable surfactant flooding is an attractive alternative to surfactant polymer flooding under certain favorable reservoir conditions. However, a gravity-stable flood requires a low velocity less than the critical velocity. Classical stability theory predicts the critical velocity needed to stabilize a miscible flood by gravity forces. This theory was tested for surfactant floods with ultralow interfacial tension and found to over-estimate the critical velocity compared to both laboratory displacement experiments and fine-grid simulations. Predictions using classical stability theory for miscible floods were not accurate because this theory did not take into account the specific physics of surfactant flooding. Stability criteria for gravity-stable surfactant flooding were developed and validated by comparison with both experiments and fine-grid numerical simulations. The effects of vertical permeability, oil viscosity and heterogeneity were investigated. Reasonable values of critical velocity require a high vertical permeability without any continuous barriers to vertical flow in the reservoir. This capability to predict when and under what reservoir conditions a gravity-stable surfactant flood can be performed at a reasonable velocity is highly significant. Numerical simulations were also used to show how gravity-stable surfactant flooding can be optimized to increase critical velocity, which shortens the project life and improves the economics of the process. The critical velocity for a stable surfactant flood is a function of the microemulsion viscosity and it turns out there is an optimum value that can be used to significantly increase the velocity and maintain stability. For example, the salinity gradient can be optimized to gradually decrease the microemulsion viscosity. Another alternative is to inject a polymer drive following the surfactant solution, but using polymer complicates the process and adds to its cost without significant benefit in most gravity-stable surfactant floods. A systematic approach was introduced to make decisions on using polymer in applications based on stability criteria and cost. Also, the effect of an aquifer on gravity-stable surfactant floods was investigated as part of a field-scale study and strategies were developed to minimize its effect on the process. This study has provided new insights into the design of an optimized gravity-stable surfactant flood. The results of the numerical simulations show the potential for high oil recovery from gravity-stable surfactant floods using horizontal wells. Application of gravity-stable surfactant floods reduces the cost and complexity of the process. The widespread use of horizontal wells has greatly increased the attractiveness and potential for conducting surfactant floods in a gravity-stable mode. This research has provided the necessary criteria and tools needed to determine when gravity-stable surfactant flooding is an attractive alternative to conventional surfactant-polymer flooding. / text

This thesis reports a new set of experimental data and presents an in-depth analysis of the flow physics of a jet stream produced by a large quarter-round corner-mounted displacement diffuser. The air velocity, temperature and turbulence intensity inside the displacement ventilation (DV) jet have been thoroughly analyzed and compared with the reported findings of previous studies and model predictions. Furthermore, thermal comfort has been analyzed using the measured data following the ASHRAE standard. This thesis also aims at establishing an accurate numerical approach for simulating the heat and fluid flow in a room ventilated by a DV system. The supply boundary condition has been thoroughly investigated, which includes tests of the conventional box and momentum modeling methods, and proposal of a more accurate modeling approach. In addition, the predictive accuracy of the standard k-ϵ, RNG k-ϵ, SST k-ω and RSM turbulence models has been examined against the experimental data.

This thesis reports a new set of experimental data and presents an in-depth analysis of the flow physics of a jet stream produced by a large quarter-round corner-mounted displacement diffuser. The air velocity, temperature and turbulence intensity inside the displacement ventilation (DV) jet have been thoroughly analyzed and compared with the reported findings of previous studies and model predictions. Furthermore, thermal comfort has been analyzed using the measured data following the ASHRAE standard. This thesis also aims at establishing an accurate numerical approach for simulating the heat and fluid flow in a room ventilated by a DV system. The supply boundary condition has been thoroughly investigated, which includes tests of the conventional box and momentum modeling methods, and proposal of a more accurate modeling approach. In addition, the predictive accuracy of the standard k-ϵ, RNG k-ϵ, SST k-ω and RSM turbulence models has been examined against the experimental data.

Buoyant plumes in cold fresh water are of interest because of the possibility of buoyancy reversal due to the nonlinear relation between temperature and density in water. Thus an initially rising plume may become a fountain. This project aims to mathematically model such plumes and fountains using numerical simulation by the means of a commercial software, Comsol Multiphysics. Both turbulent and lam- inar cases were investigated in different geometries, and with the assumption that density is a quadratic function of temperature. The turbulent flow cases as con- sidered here in this thesis are relevant to practical applications such as industrial discharge in cold lakes: whereas, the laminar flow case relates to laboratory experi- ments which are typically at scales too small for the flow to be turbulent. Previous investigation on warm discharge placed more attention on the biological implications of the spread along the lake bed, and not interested in analysing the dynamics of such flow, which turns out to be our focus. Furthermore, investigations on buoyant plumes that become negatively buoyant at later time (fountain flow) as considered previously, are based on the assumption that density is a linear function of tem- perature: where entrainment always reduces buoyancy. Whereas, the consideration of the temperature of maximum density is crucial and realistic in many practical situations, especially the power station warm discharge. Mixing is then bound to produce a mixture that is denser than both the discharge and the ambient water if receiving water is less than Tm: where this situation differs from plumes with linear mixing properties. Therefore, our focus is to better fathom the behaviour of warm discharge so as to give a detailed description of the flow, and also to observe buoyancy reversal whenever water that is denser than both the discharge and the receiving water is produced. The simulations were carried out for Prandtl number Pr = 7 & 11.4 and over the ranges of Froude number 0.1 ≤ Fr ≤ 5 and Reynolds numbers 50 ≤ Re ≤ 106, with source temperatures that are assumed to be higher than the temperature of maximum density Tm, and the ambient water below the Tm. Our results show some distinct behaviours from those experimental investigations by Bukreev, who also considered warm discharge where water that has temperature above the temperature Tm is initiated into a medium below Tm. The results here also showed some differences from those investigations with the linear dependence relation assumption.

The study of magnetohydrodynamics unifies the fields of fluid mechanics and electrodynamics to describe the interactions between magnetic fields and electrically conducting fluids. Flows described by magnetohydrodynamics form a significant aspect in a wide range of engineering applications, from the liquid metal blankets designed to surround and remove heat from nuclear fusion reactors, to the delivery and guidance of nanoparticles in magnetic targeted drug delivery. The ability to optimize these, and other, processes is increasingly reliant on the accuracy and stability of the numerical models used to predict such flows. This thesis addresses this by providing a detailed assessment on the performance of two electromagnetically extended Reynolds-averaged Navier-Stokes models through computations of a number of electromagnetically influenced simple channel and Rayleigh-Bènard convective flows. The models tested were the low-Re k-ε linear eddy-viscosity model of Launder and Sharma (1974), with electromagnetic modifications as proposed by Kenjereš and Hanjalić (2000), and the low-Re stress-transport model of Hanjalić and Jakirlić (1993), with electromagnetic modifications as proposed by Kenjereš and Hanjalić (2004). First, a one-dimensional fully-developed turbulent channel flow was considered over a range of Reynolds and Hartmann numbers with a magnetic field applied in both wall-normal and streamwise directions. Results showed that contributions from the electromagnetic modifications were modest and, whilst both models inherently captured some of the reduction in mean strain that a wall-normal field imposed, results from the stress-transport model were consistently superior for both magnetic field directions. Then, three-dimensional time-dependent Rayleigh-Bènard convection was considered for two different Prandtl numbers, two different magnetic field directions and over a range of Hartmann numbers. Results revealed that, at sufficiently high magnetic field strengths, a dramatic reorganization of the flow structure is predicted to occur. The vertical magnetic field led to a larger number of thinner, more cylindrical plumes whilst the horizontal magnetic field caused a striking realignment of the roll cells' axes with the magnetic field lines. This was in agreement with both existing numerical simulations and physical intuition. The superior performance of the modified stress-transport model in both flows was attributed to both its ability to provide better representation of stress generation and other processes, and its ability to accommodate the electromagnetic modifications in a more natural, and exact, fashion. The results demonstrate the capabilities of the stress-transport approach in modelling MHD flows that are relevant to industry and offer potential for those wishing to control flow structure or levels of turbulence without recourse to mechanical means.

Une fontaine peut se créer quand la flottabilité d'un rejet vertical s'oppose à sa quantité de mouvement. Ce type d'écoulement connaît beaucoup d'applications que ce soit dans la nature (panaches issus des éruptions volcaniques) ainsi que dans l'industrie du bâtiment (chauffage et refroidissement) ou dans le domaine des risques (rejets accidentel de gaz lourd). Dans cette thèse, nous nous focalisons sur l'étude des fontaines turbulentes miscibles. Dans le premier chapitre nous reformulons le modèle théorique de Morton et al. (1956) pour traiter le cas des fontaines en milieu linéairement stratifié. La résolution de ce modèle permet d'obtenir des relations analytiques pour la hauteur de la fontaine et sa hauteur d'étalement. Ce modèle est, par la suite, étendu au cas des panaches et des jets turbulents en milieu linéairement stratifié. Dans le second chapitre, nous proposons un modèle théorique permettant d'étudier une fontaine turbulente miscible en régime établi. Pour calibrer ce modèle, des simulations numériques aux grandes échelles (LES) sont utilisées pour obtenir une estimation des valeurs des constantes associées aux phénomènes d'échanges turbulents entre les parties ascendante et descendante de la fontaine. L'objectif du dernier chapitre est d'apporter, à partir d'expérimentations en laboratoire, des informations quantitatives sur l'influence de forts écarts de masses volumiques dans les écoulements de type fontaine. Les expériences sont réalisées pour des fontaines gazeuses (mélange air/hélium) en régime établi. / A fountain can occur when the buoyancy of a vertically released fluid opposes its momentum. Such flows have many applications in nature (plumes issuing from volcanic eruption), building industry (cooling or heating) or in the area of risk management (accidental release of heavy dangerous gas). In this thesis, we focus on the study of miscible turbulent fountains. In the first chapter, we revisit the theoretical model of Morton et al. (1956) to handle the case of fountains in linearly stratified fluid. The resolution of this model allows us to obtain analytical relations for the fountain height as well as the spreading height of its horizontal layer. This model is subsequently extended to the case of turbulent jets and plumes in linearly stratified fluid. In the second chapter, we propose a theoretical model for the study of a turbulent miscible fountain in a steady state. To calibrate this model, large eddy simulations (LES) are used to obtain an estimate of the values of the constants associated with the additional terms appearing in the equations. The objective of the final chapter is to provide, from laboratory experiments, quantitative information on the influence of strong density differences on the behaviour of a turbulent fountain. These experiments shows that all the classical relations valid for the Boussinesq case can be extended to the non-Boussinesq case by using an appropriate definition of the Froude number.