Flow Dependent Turbulent Schmidt Number in Containment Atmosphere Mixing

Total Page:16

File Type:pdf, Size:1020Kb

Flow Dependent Turbulent Schmidt Number in Containment Atmosphere Mixing Flow Dependent Turbulent Schmidt Number in Containment Atmosphere Mixing Rok Krpan, Iztok Tiselj, Ivo Kljenak Jožef Stefan Institute, Reactor Engineering Division Jamova cesta 39 SI-1000, Ljubljana, Slovenia [email protected], [email protected], [email protected] ABSTRACT In turbulence modelling, constant value of turbulent Schmidt number is usually used. However, when compared to experiments, constant value does not always provide the best results. A new model for prescribing different values of turbulent Schmidt number to different flow regions is proposed and validated against three experiments on containment mixing performed in SPARC and PANDA experimental facilities. 1 INTRODUCTION During a severe accident in a light water nuclear reactor, hydrogen combustion could threaten the integrity of the nuclear power plant containment, which could lead to release of radioactive material into the environment. The study of hydrogen distribution in the containment is thus important to predict the occurrence of regions with high local hydrogen concentrations and flammable mixture in order to effectively install hydrogen mitigation systems in containments. Various experiments are being performed to simulate atmosphere mixing occurring in containment during severe accidents and results are used to validate Computational Fluid Dynamics (CFD) codes in order to simulate phenomena in actual power plants [1]. Turbulent Schmidt number (푆푐푡) and turbulent Prandtl number (푃푟푡) are two non- dimensional numbers used in turbulence modelling to describe the turbulent transport of mass and heat, respectively. These two numbers are also connected to turbulent fluctuations and turbulent transport of momentum. In CFD calculations of the injection of fluid from a nozzle of circular cross-section into a reservoir containing stagnant fluid of similar density, constant values of 푆푐푡 and 푃푟푡 are usually used. Typically, their values are specified with comparison of calculation and experimental results, and are in the range of 0.7 to 1 [2]. However, a constant value does not always provide the best results when compared to experiments. Values of these turbulent numbers may change throughout the flow field [3] and within the boundary layer [4]. Several authors proposed different models for 푃푟푡 and 푆푐푡 values. Yimer et al. [5] proposed a parabolic curve fit for 푆푐푡. This function is based on several turbulent round free jet experiments, which showed that the 푆푐푡 increases monotonously from a value of 0.62 on the jet axis to 0.82 on the jet edge. However, a constant average value of 0.7 was recommended for use in CFD applications involving axisymmetric free-jet flows. Sturgess and McManus [6] suggested a formulation for 푆푐푡 based on the 푘 − 휀 turbulence model. One of the model constants, 퐶휇, is modelled using the ratio of turbulence kinetic energy production to its dissipation rate, instead of keeping its value constant. Keistler et al. [7] proposed a new set of transport equations for enthalpy and convection-diffusion equation for mean mass fraction including variance and dissipation rate in order to determine 푆푐푡 and 푃푟푡. However, this set of 717.1 717.2 equations involve 26 additional constants. Goldberg [8] proposed a method for calculating 푆푐푡 and 푃푟푡 values based on algebraic Reynolds stress model and Reynolds stress anisotropy. This method gives in the bulk flow, a constant 푆푐푡value of 0.7, while in the near wall region a lower 푆푐푡 mean value of 0.34 is obtained. An experiment on containment mixing performed in SPARC experimental facility is used to validate the model proposed in the present work. The interaction of vertical axisymmetric air or steam jet on a horizontal layer of helium-steam or helium-air mixture is simulated with open- source CFD code OpenFOAM. 2 LIMITED VERTICAL JET Vertical jet, injected in the opposite direction of the gravity, considered in the present work is limited in vertical flow (streamwise) direction, at first with a light gas layer and later with the ceiling of the vessel (Figure 1). Consequently, the jet at some point changes direction, flows downstream around the main jet and generates a recirculation flow. As shown in Figure 1, different flow regions can be defined using the mean flow properties: A. main rising jet, B. shear flow region, C. returning jet, D. quiescent environment. Figure 1: Absolute value of vertical velocity in a vertical jet. Let us define a Cartesian coordinate system with origin located in the axis of the injection. Figure 2 left shows vertical velocity radial (x-direction) profile. The velocity is the highest at the jet centre (in the coordinate system origin), meaning its gradient in radial direction is zero. Further from the jet origin, vertical velocity decreases towards the jet edge (gradient is negative), becomes zero on the boundary between main jet and the down flowing returning jet (in the centre of the shear flow region), gets slightly negative in the returning down flowing jet (gradient is zero again in the axis of returning jet), and then the velocity increases (gradient is positive) back to zero in the quiescent environment. Figure 2 right shows vertical velocity vertical profile in the jet’s axis. The velocity decreases (velocity gradient is negative) until the jet reaches layer consisting of light gas or ceiling, where it changes direction. Proceedings of the International Conference Nuclear Energy for New Europe, Portorož, Slovenia, September 7 ̶ 10, 2020 717.3 0.8 3.5 0.7 3 0.6 2.5 0.5 2 0.4 Vertical velocity [m/s] Vertical velocity [m/s] 1.5 0.3 1 0.2 0.1 0.5 0 0 -0.1 -0.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 Radial distance [m] Vertical distance [m] Figure 2: Vertical velocity radial profile (left) and vertical velocity vertical profile (right). In Figure 3 are shown vertical velocity gradients in radial x-direction (휕푢푧/휕푥) (Figure 3 left) and vertical (휕u푧/휕푧) direction (Figure 3 right. 휕푢푧/휕푥 is multiplied with x-coordinate sign to maintain the same sign when crossing the origin. Since the jet is axisymmetric, the gradient in y-direction 휕푢푧/휕푦 is equal to 휕푢푧/휕푥. Regarding the sign of 휕푢푧/휕푥, an area which comprises main jet, shear flow region, and part of returning jet to its axis, can be defined. With the sign of 휕u푧/휕푧 similar area can be defined with details on top of the jet. Vertical artefacts on both sides of the jet, are visual errors generated by the visualization software, while computing the gradients. 휕푢 휕푢 푧 푧 휕푥 휕푧 Figure 3: Vertical velocity radial gradient (left) and vertical velocity vertical gradient (right). 3 PHYSICAL MODEL The containment mixing experiments were simulated with the OpenFOAM CFD code. The atmosphere in the vessel was considered as a compressible mixture of ideal gases. The erosion and mixing process was modelled as a two-component single-phase flow, where single continuity, momentum and total energy equation were solved. Convection-diffusion equation was solved only for a single specie, while the other species mass fraction was calculated using ∑푖=1,2 푌푖 = 1. Proceedings of the International Conference Nuclear Energy for New Europe, Portorož, Slovenia, September 7 ̶ 10, 2020 717.4 Turbulent Schmidt (푆푐푡) and turbulent Prandtl (푃푟푡) numbers are two non-dimensional numbers defined as the ratio of momentum eddy diffusivity (휈푡) to mass eddy diffusivity (퐷푡), and as the ratio of momentum eddy diffusivity to thermal eddy diffusivity (훼푡), respectively. In other words, they describe the ratio of the rates of turbulent transport of momentum (momentum diffusion) to the turbulent transport of mass (turbulent mass diffusion) or heat (turbulent heat diffusion) [9]: 휈푡 휈푡 푆푐푡 = , 푃푟푡 = (1) 퐷푡 훼푡 푆푐푡 takes place in the convection-diffusion equation: 휕휌푌퐻푒 휕 휕 휇푡 휕푌퐻푒 + (휌푢푗푌퐻푒) = ((휌퐷 + ) ) (2) 휕푡 휕푥푗 휕푥푗 푆푐푡 휕푥푗 where 휌, 푌퐻푒, 푡, 푢푗, 퐷, 휇푡 and 푆푐푡 are density, helium mass fraction, time, velocity, diffusion coefficient and eddy viscosity, respectively. On the other hand, 푃푟푡 takes place in the total energy equation: 휕휌ℎ 휕 휕휌퐾 휕 휕푝 + (휌푢푗ℎ) + + (휌푢푗퐾) − = 휕푡 휕푥푗 휕푡 휕푥푗 휕푡 (3) 휕 휇 휇푡 휕ℎ 휕 휇푡 휕푌푖 (( + ) ) − ∑ ( (ℎ푖 ⋅ (휌퐷 + )) ) + 휌푘푢푘 휕푥푗 푃푟 푃푟푡 휕푥푗 휕푥푗 푆푐푡 휕푥푗 푖 where ℎ, 퐾, 푝, 푃푟 and ℎ푖 are enthalpy, kinetic energy, pressure, Prandtl number and gas species enthalpy, respectively. 3.1 Variable Turbulent Schmidt and Turbulent Prandtl Number Model The vertical velocity and its gradients are used to determine the shear flow in the vertical limited jet (i.e. the jet boundary) and to specify the regions with different 푆푐푡 value. The proposed model to prescribe different 푆푐푡 value to different flow regime region is: 1 Sc = Sc + (Sc − Sc ) t t,min 2 t,max t,min (4) 푟푣 휕푢푧 휕푢푧 휕푢푧 |푢| 1 − tanh c푘,1 ( + + ) − − 푐푘,2 max(|푢|) ⏟휕푥 휕푦 휕푧 ⏟max (| 푢 |) velocity gradients normalized velocity { [ ( ) ]} where 푆푐푡,푚푖푛, 푆푐푡,푚푎푥, 푐푘,1, 푐푘,2 and 푟푣 are minimal and maximal 푆푐푡 values specified by the user, two model constants and inlet radius, respectively. Term with velocity gradients prevails in the main and in the returning jet, while term with normalized velocity prevails in the axis of the main jet. The constants 푐푘,1 defines the behaviour of the function regarding the mean flow properties. The constant 푐푘,2 lowers the function argument, where vertical velocity and its gradients are small, but their values are still noticeable. It sets the relative velocity value cut-off and selects the regions to which the function is applied. The values of the constants were set in such manner, that the selected area contains main jet, shear flow region, and part of the returning jet.
Recommended publications
  • Turbulent-Prandtl-Number.Pdf
    Atmospheric Research 216 (2019) 86–105 Contents lists available at ScienceDirect Atmospheric Research journal homepage: www.elsevier.com/locate/atmosres Invited review article Turbulent Prandtl number in the atmospheric boundary layer - where are we T now? ⁎ Dan Li Department of Earth and Environment, Boston University, Boston, MA 02215, USA ARTICLE INFO ABSTRACT Keywords: First-order turbulence closure schemes continue to be work-horse models for weather and climate simulations. Atmospheric boundary layer The turbulent Prandtl number, which represents the dissimilarity between turbulent transport of momentum and Cospectral budget model heat, is a key parameter in such schemes. This paper reviews recent advances in our understanding and modeling Thermal stratification of the turbulent Prandtl number in high-Reynolds number and thermally stratified atmospheric boundary layer Turbulent Prandtl number (ABL) flows. Multiple lines of evidence suggest that there are strong linkages between the mean flowproperties such as the turbulent Prandtl number in the atmospheric surface layer (ASL) and the energy spectra in the inertial subrange governed by the Kolmogorov theory. Such linkages are formalized by a recently developed cospectral budget model, which provides a unifying framework for the turbulent Prandtl number in the ASL. The model demonstrates that the stability-dependence of the turbulent Prandtl number can be essentially captured with only two phenomenological constants. The model further explains the stability- and scale-dependences
    [Show full text]
  • Analysis of Mass Transfer by Jet Impingement and Heat Transfer by Convection
    Analysis of Mass Transfer by Jet Impingement and Study of Heat Transfer in a Trapezoidal Microchannel by Ejiro Stephen Ojada A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Mechanical Engineering Department of Mechanical Engineering University of South Florida Major Professor: Muhammad M. Rahman, Ph.D. Frank Pyrtle, III, Ph.D. Rasim Guldiken, Ph.D. Date of Approval: November 5, 2009 Keywords: Fully-Confined Fluid, Sherwood number, Rotating Disk, Gadolinium, Heat Sink ©Copyright 2009, Ejiro Stephen Ojada i TABLE OF CONTENTS LIST OF FIGURES ii LIST OF SYMBOLS v ABSTRACT viii CHAPTER 1: INTRODUCTION AND LITERATURE REVIEW 1 1.1 Introduction (Mass Transfer by Jet Impingement ) 1 1.2 Literature Review (Mass Transfer by Jet Impingement) 2 1.3 Introduction (Heat Transfer in a Microchannel) 7 1.4 Literature Review (Heat Transfer in a Microchannel) 8 CHAPTER 2: ANALYSIS OF MASS TRANSFER BY JET IMPINGEMENT 15 2.1 Mathematical Model 15 2.2 Numerical Simulation 19 2.3 Results and Discussion 21 CHAPTER 3: ANALYSIS OF HEAT TRANSFER IN A TRAPEZOIDAL MICROCHANNEL 42 3.1 Modeling and Simulation 42 3.2 Results and Discussion 48 CHAPTER 4: CONCLUSION 65 REFERENCES 67 APPENDICES 79 Appendix A: FIDAP Code for Analysis of Mass Transfer by Jet Impingement 80 Appendix B: FIDAP Code for Fluid Flow and Heat Transfer in a Composite Trapezoidal Microchannel 86 i LIST OF FIGURES Figure 2.1a Confined liquid jet impingement between a rotating disk and an impingement plate, two-dimensional schematic. 18 Figure 2.1b Confined liquid jet impingement between a rotating disk and an impingement plate, three-dimensional schematic.
    [Show full text]
  • Turbulence, Entropy and Dynamics
    TURBULENCE, ENTROPY AND DYNAMICS Lecture Notes, UPC 2014 Jose M. Redondo Contents 1 Turbulence 1 1.1 Features ................................................ 2 1.2 Examples of turbulence ........................................ 3 1.3 Heat and momentum transfer ..................................... 4 1.4 Kolmogorov’s theory of 1941 ..................................... 4 1.5 See also ................................................ 6 1.6 References and notes ......................................... 6 1.7 Further reading ............................................ 7 1.7.1 General ............................................ 7 1.7.2 Original scientific research papers and classic monographs .................. 7 1.8 External links ............................................. 7 2 Turbulence modeling 8 2.1 Closure problem ............................................ 8 2.2 Eddy viscosity ............................................. 8 2.3 Prandtl’s mixing-length concept .................................... 8 2.4 Smagorinsky model for the sub-grid scale eddy viscosity ....................... 8 2.5 Spalart–Allmaras, k–ε and k–ω models ................................ 9 2.6 Common models ........................................... 9 2.7 References ............................................... 9 2.7.1 Notes ............................................. 9 2.7.2 Other ............................................. 9 3 Reynolds stress equation model 10 3.1 Production term ............................................ 10 3.2 Pressure-strain interactions
    [Show full text]
  • Huq, P. and E.J. Stewart, 2008. Measurements and Analysis of the Turbulent Schmidt Number in Density
    GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L23604, doi:10.1029/2008GL036056, 2008 Measurements and analysis of the turbulent Schmidt number in density stratified turbulence Pablo Huq1 and Edward J. Stewart2 Received 18 September 2008; revised 27 October 2008; accepted 3 November 2008; published 3 December 2008. [1] Results are reported of measurements of the turbulent where rw is the buoyancy flux, uw is the Reynolds shear Schmidt number Sct in a stably stratified water tunnel stress. The ratio Km/Kr is termed the turbulent Schmidt experiment. Sct values varied with two parameters: number Sct (rather than turbulent Prandtl number) as Richardson number, Ri, and ratio of time scales, T*, of salinity gradients in water are used for stratification in the eddy turnover and eddy advection from the source of experiments. Prescription of a functional dependence of Sct turbulence generation. For large values (T* 10) values on Ri = N2/S2 is a turbulence closure scheme [Kantha and of Sct approach those of neutral stratification (Sct 1). In Clayson, 2000]. Here the buoyancy frequency N is defined 2 contrast for small values (T* 1) values of Sct increase by the gradient of density r as N =(Àg/ro)(dr/dz), and g is with Ri. The variation of Sct with T* explains the large the gravitational acceleration. S =dU/dz is the velocity scatter of Sct values observed in atmospheric and oceanic shear. Analyses of field, lab and numerical (RANS, DNS data sets as a range of values of T* occur at any given Ri. and LES) data sets have led to the consensus that Sct The dependence of Sct values on advective processes increases with Ri [see Esau and Grachev, 2007, and upstream of the measurement location complicates the references therein].
    [Show full text]
  • The Analogy Between Heat and Mass Transfer in Low Temperature Crossflow Evaporation
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by AUT Scholarly Commons The analogy between heat and mass transfer in low temperature crossflow evaporation Reza Enayatollahi*, Roy Jonathan Nates, Timothy Anderson Department of Mechanical Engineering, Auckland University of Technology, Auckland, New Zealand. Corresponding Author: Reza Enayatollahi Email Address: [email protected] Postal Address: WD308, 19 St Paul Street, Auckland CBD, Auckland, New Zealand Phone Number: +64 9 921 9999 x8109 Abstract This study experimentally determines the relationship between the heat and mass transfer, in a crossflow configuration in which a ducted airflow passes through a planar water jet. An initial exploration using the Chilton-Colburn analogy resulted in a coefficient of determination of 0.72. On this basis, a re-examination of the heat and mass transfer processes by Buckingham’s-π theorem and a least square analysis led to the proposal of a new dimensionless number referred to as the Lewis Number of Evaporation. A modified version of the Chilton-Colburn analogy incorporating the Lewis Number of Evaporation was developed leading to a coefficient of determination of 0.96. 1. Introduction Heat and mass transfer devices involving a liquid interacting with a gas flow have a wide range of applications including distillation plants, cooling towers and aeration processes and desiccant drying [1-5]. Many studies have gone through characterising the heat and mass transfer in such configurations [6-9]. The mechanisms of heat and mass transfer are similar and analogical. Therefore, in some special cases where, either the heat or mass transfer data are not reliable or may not be available, the heat and mass transfer analogy can be used to determine the missing or unreliable set of data.
    [Show full text]
  • 1 Pe = Re× Sc Peclet Number
    Passive Scalar Diffusion in the Near Field Region of Turbulent Rectangular Submerged Free Jets Andrea Boghia, Ivan Di Venutab, Fabio Gorib,* a School of Water, Energy and Agrifood, Cranfield University, Cranfield, Bedfordshire MK43 0AL, United Kingdom, b Department of Industrial Engineering, University of Rome “Tor Vergata”, Via del Politecnico 1, 00133 Rome, Italy. * Corresponding Author, [email protected] Abstract Jets are a common way to transfer mass among fluids, or from a fluid to a surface. At moderate Reynolds numbers and low turbulent intensities the jet exhibits a Near Field Region (NFR) several diameters long. Numerical results and a theoretical model are presented for the passive scalar diffusion in the NFR of a submerged free jet. Large Eddy Simulations (LES), in the Reynolds number 5000-40,000 and the Schmidt number range 1-100, are performed obtaining the passive scalar fields. Three mathematical models for the passive scalar diffusion are presented; the first one is valid in the NFR, specifically in the Undisturbed Region of Flow (URF), and the other two, obtained under the hypotheses of Tollmien and Görtler momentum spreadings, are valid in the Potential Core Region (PCR). The last two models employ a turbulent Schmidt number inversely proportional to the mean velocity gradient, conclusion obtained by the LES numerical results. The self-similar solutions of the passive scalar show good agreement with the LES results. The wide range of Reynolds and Schmidt numbers investigated gives generality to the results. Keywords: Submerged Rectangular Free Jet; Turbulent Flow; Near Field Region; Undisturbed Region of Flow; Passive Scalar; Large Eddy Simulation; Self-Similarity.
    [Show full text]
  • Compressible Effects Modelling in Turbulent Cavitating Flows Jean Decaix, Eric Goncalvès Da Silva
    Compressible effects modelling in turbulent cavitating flows Jean Decaix, Eric Goncalvès da Silva To cite this version: Jean Decaix, Eric Goncalvès da Silva. Compressible effects modelling in turbulent cav- itating flows. European Journal of Mechanics - B/Fluids, Elsevier, 2013, 39, pp.11-31. 10.1016/j.euromechflu.2012.12.001. hal-00768757 HAL Id: hal-00768757 https://hal.archives-ouvertes.fr/hal-00768757 Submitted on 12 Feb 2013 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Compressible Effects Modelling in Turbulent Cavitating Flows Jean Decaix∗, Eric Goncalv`es∗ LEGI-Grenoble INP, 1025 rue de la Piscine, 38400 St Martin d’Heres, France Abstract A compressible, multiphase, one-fluid RANS solver has been developed to study turbulent cavitating flows. The interplay between turbulence and cav- itation regarding the unsteadiness and structure of the flow is complex and not well understood. This constitutes a determinant point to accurately sim- ulate the dynamic behaviour of sheet cavities. In the present study, different formulations including compressibility effects on turbulence are investigated. Numerical results are given for two partial cavities on Venturi geometries and comparisons are made with experimental data.
    [Show full text]
  • Stability and Instability of Hydromagnetic Taylor-Couette Flows
    Physics reports Physics reports (2018) 1–110 DRAFT Stability and instability of hydromagnetic Taylor-Couette flows Gunther¨ Rudiger¨ a;∗, Marcus Gellerta, Rainer Hollerbachb, Manfred Schultza, Frank Stefanic aLeibniz-Institut f¨urAstrophysik Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam, Germany bDepartment of Applied Mathematics, University of Leeds, Leeds, LS2 9JT, United Kingdom cHelmholtz-Zentrum Dresden-Rossendorf, Bautzner Landstr. 400, D-01328 Dresden, Germany Abstract Decades ago S. Lundquist, S. Chandrasekhar, P. H. Roberts and R. J. Tayler first posed questions about the stability of Taylor- Couette flows of conducting material under the influence of large-scale magnetic fields. These and many new questions can now be answered numerically where the nonlinear simulations even provide the instability-induced values of several transport coefficients. The cylindrical containers are axially unbounded and penetrated by magnetic background fields with axial and/or azimuthal components. The influence of the magnetic Prandtl number Pm on the onset of the instabilities is shown to be substantial. The potential flow subject to axial fieldspbecomes unstable against axisymmetric perturbations for a certain supercritical value of the averaged Reynolds number Rm = Re · Rm (with Re the Reynolds number of rotation, Rm its magnetic Reynolds number). Rotation profiles as flat as the quasi-Keplerian rotation law scale similarly but only for Pm 1 while for Pm 1 the instability instead sets in for supercritical Rm at an optimal value of the magnetic field. Among the considered instabilities of azimuthal fields, those of the Chandrasekhar-type, where the background field and the background flow have identical radial profiles, are particularly interesting.
    [Show full text]
  • On the Values for the Turbulent Schmidt Number in Environmental Flows
    Università di Napoli Federico II, 1224 On the Values for the Turbulent Schmidt Number in Environmental Flows Carlo Gualtieri University of Napoli Federico II, Napoli, Italy Acknowledgements Athanasios Angeloudis1, Fabian Bombardelli2, Sanjeev Jha3, Thorsten Stoesser 4 1 Imperial College, London, UK 2 University of California at Davis, Davis, USA Lyon, April 27, 2018 3 University of Manitoba, Winnipeg, Canada 4 Cardiff University, Cardiff, UK Outline of seminar Introduction Some remarks on turbulence and turbulence modeling The Turbulent Schmidt Number within the RANS Approach Review of the Literature on the Parameterization of the Turbulent Schmidt Number Water Systems Atmosphere Systems Case Studies Contaminant Dispersion Due to Transverse Turbulent Mixing in a Shallow Water Flow Tracer Transport in a Contact Tank Sediment Transport in Suspension Discussion Concluding remarks On the Values for the Turbulent Schmidt Number in Environmental Flows Università di Napoli Federico II, 1224 Importance of fluids in the environment – 1 Think about it. All living creatures are immersed in a moving fluid or another, in the air of the atmosphere or water of rivers, lakes and sea. http://quagroup.com/hydrologic-cycle/ On the Values for the Turbulent Schmidt Number in Environmental Flows Università di Napoli Federico II, 1224 Importance of fluids in the environment – 2 Flow transport has a number of benefits: It takes things away, and usually to a larger place (rivers to lake and sea, local airshed to broader atmosphere, troposphere to stratosphere, etc.) It creates time delay to allow modifications to take place along the way (ex: biodegradation of organic wastes down a river) Fluid motion at human and larger scales is almost always unstable, and the ensuing turbulence generates mixing , which • reduces concentrations, and • promotes encounters (fluid-fluid, interfaces with other fluid or soil), which in turn promote chemical reactions and transformation.
    [Show full text]
  • Effect of Schmidt Number on the Structure and Propagation of Density Currents
    OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible This is an author’s version published in: http://oatao.univ-toulouse.fr/9409 Official URL: http://doi.org/10.1007/s00162-008-0085-2 To cite this version: Bonometti, Thomas and Balachandar, Sivaramakrishnan Effect of schmidt number on the structure and propagation of density currents. (2008) Theoretical and Computational Fluid Dynamics, 22 (5). 341-361. ISSN 0935-4964 Any correspondence concerning this service should be sent to the repository administrator: [email protected] DOI 10.1007/s00162-008-0085-2 Thomas Bonometti · S. Balachandar Effect of Schmidt number on the structure and propagation of density currents Abstract The results of a numerical study of two- and three-dimensional Boussinesq density currents are described. They are aimed at exploring the role of the Schmidt number on the structure and dynamics of den- sity driven currents. Two complementary approaches are used, namely a spectral method and a finite-volume interface capturing method. They allow for the first time to describe density currents in the whole range of 2 4 Schmidt number 1 ≤ Sc ≤∞and Reynolds number 10 ≤ Re ≤ 10 . The present results confirm that the Schmidt number only weakly influences the structure and dynamics of density currents provided the 4 Reynolds number of the flow is large, say of O(10 ) or more. On the contrary low- to moderate-Re density currents are dependant on Sc as the structure of the mixing region and the front velocities are modified by diffusion effects.
    [Show full text]
  • Mixing in High Schmidt Number Turbulent Jets Differs from Gas-Phase
    Mixing in High Schmidt Number Turbulent Jets Thesis by Paul Lewis Miller In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasadena., California 1991 (submitted April 25, 1991) 01991 Paul Lewis Miller All rights reserved Acknowledgements I wish to acknowledge the following people for their assistance and/or sup- port during the time I have been at Caltech and GALCIT. First and foremost, I thank my advisor, Paul Dimotakis, for his encouragement throughout my graduate studies, and for many suggestions and contributions to this work. His scientific understanding and his personal character have made an indelible impression upon me. In collaboration with Paul, Dan Lang is responsible for a large number of accomplishments while at GALCIT that have greatly facilitated and improved this thesis. These include, among others, development of the HYDRA computer system and libraries, data acquisition support, and electronic design. His efforts have been a tremendous help. I give special thanks to Gene Broadwell, whose ideas concerning turbulent mixing have had a strong influence on this thesis. His friendship over the past several years is especially appreciated. Several of my fellow students have contributed help. In particular I'd like to thank Cliff Frieler, who is always interested in a problem, usually has a solution, and never refuses a request for assistance, and Dave Dowling, for assistance at the onset of the experiments and the indirect influence of his impressive thesis. Many staff members at GALCIT provided important contributions in the course of this work. Herb Gaebler offered assistance with a number of mechani- cal matters, Harry Hamaguchi produced the photographs of the facility and con- ducted the difficult printing of the jet growth-rate pictures, Pave1 Svitek aided in construction of the Phase I1 improvements and with the jet calibration, and Phil Wood created a piece of art with his machining of the new nozzle.
    [Show full text]
  • Characterising the Heat and Mass Transfer Coefficients for a Crossflow
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by AUT Scholarly Commons Characterising the heat and mass transfer coefficients for a crossflow interaction of air and water Reza Enayatollahi*, Roy Jonathan Nates, Timothy Anderson Department of Mechanical Engineering, Auckland University of Technology, Auckland, New Zealand. Corresponding Author: Reza Enayatollahi Email Address: [email protected] Postal Address: WD308, 19 St Paul Street, Auckland CBD, Auckland, New Zealand Phone Number: +64 9 921 9999 x8109 Abstract An experimental study was performed in order to characterise the heat and mass transfer processes, where an air stream passes through a sheet of falling water in a crossflow configuration. To achieve this, the hydrodynamics of a vertical liquid sheet in a ducted gaseous crossflow were studied. Four distinct flow regimes were identified (a stable sheet, a broken sheet, a flapping sheet and a lifted sheet) and mapped using Reynolds and Weber numbers. Subsequently, the Buckingham π theorem and a least squares analyses were employed leading to the proposal of two new dimensionless numbers referred to as the Prandtl Number of Evaporation and the Schmidt Number of Evaporation. These describe the heat and mass transfer in low temperature evaporation processes with crossflow interaction. Using these dimensionless numbers, empirical correlations for Sherwood and Nusselt numbers for the identified flow regimes were experimentally determined. These correlations were in a good agreement with their corresponding experimental values. It was found that flapping sheets have the strongest heat and mass transfer intensities whereas the weakest intensities were seen for the “stable” sheets.
    [Show full text]