Lecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
|
|
- Allen Fowler
- 5 years ago
- Views:
Transcription
1 Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of velocities and state variables in space and time. This occurs when the Reynolds number Re = UL/ν 1 L is a characteristic dimension of the vessel U is a characteristic velocity of the flow ν is the kinematic viscosity µ/ρ (cm 2 /s). Re 1 = L/U L 2 /ν = t a overall time to advect a fluid particle a distance L = 1 t d viscous dissipation time over distance L Disturbances in the flow, due to instabilities, may develop large velocity gradients. At Re 1, viscous dissipation prevents the development of these large velocity gradients. But at Re 1, viscous dissipation is unable to smear out the large developing gradients. 2 permission of the owner, M. Matalon. 1
2 Statistical description The aim is to describe the fluctuating velocity and scalar fields in terms of their statistical distributions. 1. Averaging time averaging ū(x) = lim T T 0 u(t, x)dt u(t) ū ensemble averaging compiled from many different realizations where u(t n, x) isthen th realization 1 N u(t, x) = lim u(t n, x) N N n=1 good also if ū depends on time ū t 1 t 2 t n measurements taken at a fixed point t t If we decompose u(t, x) into a mean and a fluctuation u (t, x) then, by definition, u (t, x) =u (t, x) = Probability Density Functions (pdf) The distribution function F u (U) of u is defined by the probability of finding a value of u<u F u (U) =p(u <U) U is the sample space of the random variable u consists of all possible realizations of u. The probability of finding a value of u in an interval U<u<U+ U is given by p(u<u<u+ U) =F u (U + U) F u (U) The probability density function (pdf) of u is defined as Clearly P u (U) = df u(u) du P u (U)dU =1 4 permission of the owner, M. Matalon. 2
3 In turbulent flows, the pdf of any stochastic variable depends, in principle, on the position x and on time t, or P u (U; x, t) is the probability density in U-space and a function of x and t. The simplified notation P (u; x, t) will be used from now on. Once the pdf of a variable is known one can define the average or mean value the mean of a function g(u) ū(x, t) = ḡ(x, t) = σ 2 = up(u; x, t)du g(u) P (u; x, t)du and the variance (a measure of the width of the pdf) (u ū) 2 P (u; x, t))du σ 2 = u 2 2uū +ū 2 = u 2 ū 2 σ 2 = u 2 u =ū + u u 2 ū 2 = ū 2 +2u ū u 2 ū 2 = u 2 mean of the fluctuations squared 5 u(t) ū u bimodal pdf P (u) ū u normal pdf t P (u) 6 permission of the owner, M. Matalon. 3
4 u 2 is the variance - the extent of variation from the mean the root mean square u 2 is the intensity the correlation of the fluctuations of two variables is often of interest uv= (ū + u )( v + v )=ū v + ū v +ū v + u v =ū v + u v correlation R uv R uv indicates the extent of the interdependence of the two variables the correlation coefficient R uv is often normalized such that 1 <R uv < 1 R uv = u v u 2 v 2 R uv 0 R uv ±1 R uv > 0 u and v are poorly correlated u and v are nearly perfectly correlated (or anticorrelated) u and v are correlated (and have the same sign most of the time) 7 Turbulent Scales Turbulent flows are said to be made up of small eddies with a multitude of sizes and vorticities. A fluid eddy is a macroscopic fluid element in which the microscopic elements comprising the eddy behave as a unit. An eddy is identified by a characteristic size and a characteristic velocity. A number of small eddies may be embedded in a large eddy. A characteristic of a turbulent flow is the existence of a wide range of length scales, or eddy sizes. The largest size L, is the macroscale or the flow scale. Eddies of size L are characterized by the mean flow velocity U. The Reynolds number of these eddies Re = UL/ν is large so that the effect of viscosity is negligibly small. The large eddies, however, break up as a result of instabilities, transferring their energy to smaller eddies, characterized by smaller Reynolds numbers. This process continues until the Reynolds number is sufficiently small; now viscosity acts to smear out large velocity gradients or dissipate the available kinetic energy. This is basically the eddy cascade hypothesis, proposed by Kolmogorov for homogeneous isotropic turbulence. 8 permission of the owner, M. Matalon. 4
5 In order to characterize the distribution of eddy length scales at any position in the flow field, we use the correlation of the axial velocity between two points. The correlation of the axial velocity u at time t between two points separated by a distance r, i.e., at x and x + r, is R(x, r,t)=u (x,t) u (x + r,t) it shows how related, or coherent, are the velocities at these two points in space For homogeneous isotropic turbulence the location x is arbitrary and r may be replaced by its absolute value r = r. The normalized correlation f(r, t) =R(r, t)/u 2 (t) defines a length scale (t) = 0 f(r, t)dr called the integral length scale f(r, t) 1 f(0,t) = 1 due to the normalization and goes to zero as r, because very large eddies occur less frequent. r 9 If v is a typical velocity scale, i.e., a measure of the velocity fluctuations around the mean, it can be defined from the kinetic energy k = 1 2 v i v i = 3 2 v2 2k v = 3 v represents the turnover time of the integral scale eddies. The turnover time of these eddies, t = /v is then proportional to the integral time scale. 10 permission of the owner, M. Matalon. 5
6 For very small values of r only very small eddies fit into the distance between x and x + r. The motion of these small eddies is influenced by viscosity which provides an additional dimensional quantity for scaling. We denote the dissipation of kinetic energy, per unit mass, by. [ ]= 1 s (m/s)2 =m 2 /s 3 [ ν ]=m 2 /s Dimensional analysis then yields the Kolmogorov length scale η ν 3 η = 1/4 We can also define a Kolmogorov time scale ν 1/2 t η = f(r, t) 1 and velocity scale v η =(ν) 1/4 η r inertial range 11 According to Kolmogorov s theory the energy transfer from the large eddies at the integral scale (i.e, of length scale ) is equal to the dissipation of energy at the Kolmogorov scale. Since at the integral scale the energy transfer rate = k3/2 = v2 /v = v3 = k3/2 If we define the turbulent Reynolds number as Re T = v /ν, the ratios of the smallest to largest length, time and velocity scales are η = (ν3 /) 1/4 ν 3/4 3/4 = = Re v T t η t = (ν/)1/2 v / ν 1/2 1/2 = = Re v T v η v = (ν)1/4 ν 1/4 1/4 v = = Re v T 12 permission of the owner, M. Matalon. 6
7 Various approaches Direct Numerical Simulation (DNS). Resolves the entire range of turbulent length scales, without recourse to any modeling assumptions. All the turbulence scales must be resolved in the computational mesh, from the smallest Kolmogorov scale up to the integral scale. DNS is being used for fundamental turbulent studies in a limited way. Because of the high computational cost, however, it is not a feasible tool for wide parametric studies and, at the present, is not a useful approach that can be used to simulate practical engineering problems. 13 DNS of a nonpremixed methane-air jet flame with a 4-step chemistry model and 8 species. Computations used 100 millions grid points, 340, 000 processor hours (4 months) to cover a region of approximately 1 cm 3 in physical space. Pantano, JFM 2004 methane oxygen Temperature isocontours water CO H 2 H Formation of holes Species mass fraction at a given instant in time 14 permission of the owner, M. Matalon. 7
8 Reynolds-Averaged Navier-Stokes equations (RANS). In this approach any instantaneous quantity is decomposed into a time or ensemble average and a fluctuating component, and introduced into the governing equations. The equations for the mean quantities, however, contain unknown correlation terms that call for the equations for the fluctuations. These however, contain higher order correlations and the system does not close. The closure problem is particularly critical in combustion problems because of the numerous quantities arising from the chemical source terms, the species and energy equations that need to be modeled, in addition to the usual Reynolds stresses encountered in ordinary fluid mechanics. 15 Consider first the incompressible NS equations (we will be using index notation) vi ρ t + v j v j x j =0 v i x j = σ ij x j Averaging the equations, using that the fluctuations have zero mean, yields v j =0 x j vi ρ t + v v i j = σ ij ρ(vi x j x v j ) j the extra term on the RHS resulting from the fact that Dv/Dt = D v / Dt The mean velocity field must also satisfy the continuity equation. However in the momentum equations there are additional terms or effective stresses, called Reynolds stresses that give the flow an added viscosity. 16 permission of the owner, M. Matalon. 8
9 One can write equations for the fluctuating quantities, but these produce additional unknown terms, of the form v i v j v k. This is known as the closure problem of turbulence. Turbulence modeling is an attempt to write empirical expressions for the Reynolds stress. The level of closure depends on the complexity of the model. The simplest closure approach is to model the additional stresses ρ(vi v j ) in an analogous way to the stress-strain relation for a Newtonian fluid - the turbulentviscosity hypothesis. ρ(vi v j )+ 2 3 ρkδ ij = ρν vi T + v j x j x i deviatoric stress where k = 1 2 v i v i is half the kinetic energy ν T is the turbulent viscosity (also referred to as eddy viscosity) 17 Reynolds Averaged Navier-Stokes equations (RANS) v j x j =0 v i t + v v i j = 1 p + 2 x j ρ x 3 ρkδ ij + j x j ūi ν eff + ū j x j x i with ν eff = ν + ν T the effective viscosity Unlike ν, the turbulence viscosity is not a material property; it is a function of space and time; i.e. ν T = ν T (x,t). 18 permission of the owner, M. Matalon. 9
10 In combustion problems, the continuity equation must include density variations. When averaging the continuity equation ρ t + (ρv j) =0 x j the additional term ρ vj arises. Similarly, triple correlations arise of the form ρ vi v j arise from the momentum equations. It is therefore convenient to introduce a density-weighted average velocity, called the Favre average 1 average, ṽ i = ρv i ρ and the fluctuations about the Favre-average, denoted by v i = v i ṽ i. When averaging the continuity equation, Favre-averaging ensures that ρvi vanishes and the averaged continuity equation retains the same form as the continuity equation itself, namely ρ t + ( ρṽ j) =0 x j 1 Introduced by Favre in Note that the mean of the velocity fluctuation, in the Favre decomposition v i =ṽ i + vi, does not vanish, i.e., v i = 0 rather ρv i =0 More impressive is the simplification obtained in the momentum equation, because ρv i v j = ρṽ i ṽ j + ρ v i v j instead of four terms in the more conventional approach. (Note that it has also the same structure as the conventional average.) There are of course questions as to whether ṽ is readily measurable experimentally, and in extending the Reynolds-stress idea to ρ v i v j (instead of the conventional ρvi v j ), using ρ( v i v j )+ 2 3 ρ kδ ṽi ij = ρν T + ṽ j with k = 1 2v x j x i v i i but the mathematical simplification is quite significant. 20 permission of the owner, M. Matalon. 10
11 There are of course additional terms that need to be modeled, when averaging the combustion equations. These include ρ v i ψ for any scalar ψ, such as mass fraction and temperature, arising from Dψ/Dt and other arising from the highly nonlinear chemical source term ω. To show the difficulties associated with the chemical source term, consider the simple form ω = F (T ), where F (T )=B(T a T ) e E/RT F = ρb(t a T ) 1 T T a T F = F ( T ) 1 T T a T E/R T e 1+ ET exp ET R T R T 2 Because of the large Zel dovich number, even small temperature fluctuations in the reaction zone cause significant fluctuations of the chemical source term. 21 PDF Transport Equation. This approach consists of developing and solving model transport equations for the pdf of a set of variables that describe the hydrodynamics and/or thermochemical state of the mixture. An advantage of this approach is that the pdf of a certain variable contains more information than just the mean and the correlation. Another advantage is that some aspects of turbulent transport are in closed form, which limits the modeling requirements. The main drawback, however, is the large number of independent variables that characterizes the pdf making it computationally costly. 22 permission of the owner, M. Matalon. 11
12 Large Eddy Simulations (LES). LES operates on the governing equations in order to reduce the range of length scales that need to be computed, thus lowering the computational cost. A low-pass filtering is used to remove the small scales, which are replaced by subgrid models. The range of length scale filtered out, and reduction in the computational cost compared to DNS, depend on the computational resources. 23 Regimes of Premixed Turbulent Combustion We assume equal diffusivities for all reactive scalars, and a Schmidt number equal to one D th = D F = D O D; ν/d =1 The flame thickness and flame (residence) time are then given by l f = D/S L t f = D/S 2 L As a result, the turbulent Reynolds number, may be expressed as Re T = v ν Re T = v D = v S L l f The Damköhler number number (flow time/reaction time) may be defined as D = /v = S L l f /S L l f v 24 permission of the owner, M. Matalon. 12
13 The various regimes of turbulent combustion will be presented on a log-log graph depicting v /S L versus /l f. The first characterization is the line Re T =1separating laminar and turbulent combustion. Turbulent combustion corresponds to large Reynolds numbers, i.e., the region to the right, sufficiently removed from the line Re T = 1. Another distinction is between small and large Damköhler numbers. Flames are typically characterized by fast chemistry, or D 1, while the region D 1 corresponds to well stirred reactors where all turbulent scales are smaller than the chemical time scale. v /S L Re T =1 well-stirred reactor D =1 turbulent flames 1 laminar flames /l f 25 Further refinement of this diagram involves the Karlovitz number Ka, which relates the flame scales to the Kolmogorov scales. Two such numbers can be defined, one based on the flame thickness l f and the other on the thickness of the reaction zone l R. They represents the ratio of the residence time in the flame zone, or reaction zone, relative to the smallest turbulent turnover time, and thus whether turbulence penetrate and distort the flame or reaction zone structures. Ka = t f = l2 f /ν t η η 2 /ν = l2 f η 2 or Ka = t f = ν2 /lf 2 t η ν 2 /vη 2 = v2 η SL 2 and using δ = l R /l f Ka δ = l2 R η 2 = δ2 Ka Re T = v SL v 2 v 2 vη = S L l f v = D l f S L v η v η v = (ν)1/4 v = Re 1/4 T S L 2 Re T = D Re 1/2 Ka T Re T = D 2 Ka 2 It should be noted that equalities were set here for scaling purpose only. 26 permission of the owner, M. Matalon. 13
14 The relation Re T = D 2 Ka 2 implies that v = Re S T = D 2 Ka 2 = Ka 2 SL L l f l f v l f l f v 1/3 = Ka 2/3 S L l f and may be used to identify, within the regime of large D, conditions under which the smallest eddies are or are not able to penetrate the flame zone. This depends on whether Ka is smaller/larger than one (say). 1 The boundary Ka = 1 implies v 1/3 = S L l f 27 When Ka < 1 the flame residence time is much shorter than any turbulent time scales and the flame thickness is smaller than the smallest turbulent scale. The flame retains the laminar structure but is wrinkled by the turbulence motions. (a) In the sub-region where v < S L, known as the wrinkled flamelet regime, theturnovervelocityv of even the large eddies is not sufficient to compete with the advancement of the flame front at the laminar speed S L. The only effect of turbulence is to disturb the flame front, resulting in wrinkled flames. (b) In the sub-region where v >S L, referred to as the corrugated flamelet regime, the entire flame is embedded in eddies of the size of the Kolmogorov scale with velocities exceeding S L, which are able to wrinkle the flame and form pockets of fresh and burned gas. When Ka > 1 the flame residence time is on the order of the turbulent times (based on integral scale) but much larger then the Kolmogorov turnover times, and the Kolmogorov scales are smaller than the laminar flame thickness. Here the small eddies can presumably interact and modify the flame internal structure. 28 permission of the owner, M. Matalon. 14
15 v /S L distributed reaction zones Ka δ =1 thin reaction zones Ka =1 (δ =0.1) 1 Re T =1 laminar flames corrugated flamelets wrinkled flames /l f Regime diagram of premixed turbulent combustion on a log-log plot 29 When Ka > 1 the small eddies penetrate the flame zone, but may or may not modify the much thinner reaction zone. The distinction is between what has been referred to as thin reaction zones and broken reaction zones. I the latter, the laminar structure could no longer be identified. The boundary between these two regions depends on the relative thickness of the reaction zone δ. Typically δ 10 1,whichimpliesKa = 100, or v 1/3 = 10 4/3 S L l f 30 permission of the owner, M. Matalon. 15
Introduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationExercises in Combustion Technology
Exercises in Combustion Technology Exercise 4: Turbulent Premixed Flames Turbulent Flow: Task 1: Estimation of Turbulence Quantities Borghi-Peters diagram for premixed combustion Task 2: Derivation of
More informationTurbulent Premixed Combustion
Turbulent Premixed Combustion Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Example: LES of a stationary gas turbine velocity field flame 2 Course Overview Part II: Turbulent Combustion Turbulence
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationLecture 15. The Turbulent Burning Velocity
Lecture 15 The Turbulent Burning Velocity 1 The turbulent burning velocity is defined as the average rate of propagation of the flame through the turbulent premixed gas mixture. In the laminar case, solutions
More informationTURBULENT COMBUSTION NORBERT PETERS. Institut für Technische Mechanik Rheinisch-Westfälische Technische Hochschule Aachen, Germany
TURBULENT COMBUSTION NORBERT PETERS Institut für Technische Mechanik Rheinisch-Westfälische Technische Hochschule Aachen, Germany PUBLISHED BY THE PRESS SYNDICATE OF THE UNIVERSITY OF CAMBRIDGE The Pitt
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationA G-equation formulation for large-eddy simulation of premixed turbulent combustion
Center for Turbulence Research Annual Research Briefs 2002 3 A G-equation formulation for large-eddy simulation of premixed turbulent combustion By H. Pitsch 1. Motivation and objectives Premixed turbulent
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationD. VEYNANTE. Introduction à la Combustion Turbulente. Dimanche 30 Mai 2010, 09h00 10h30
D. VEYNANTE Introduction à la Combustion Turbulente Dimanche 30 Mai 2010, 09h00 10h30 Introduction to turbulent combustion D. Veynante Laboratoire E.M2.C. CNRS - Ecole Centrale Paris Châtenay-Malabry France
More informationMass Transfer in Turbulent Flow
Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.
More informationLecture 10 Turbulent Combustion: The State of the Art
Lecture 10 Turbulent Combustion: The State of the Art 10.-1 Engineering applications are typically turbulent turbulence models These models use systematic mathematical derivations based on the Navier-
More informationEulerian models. 2.1 Basic equations
2 Eulerian models In this chapter we give a short overview of the Eulerian techniques for modelling turbulent flows, transport and chemical reactions. We first present the basic Eulerian equations describing
More information7. Turbulent Premixed Flames
7. Turbulent Premixed Flames Background: - Structure of turbulent premixed flames; 7. Turbulent Premixed Flames 1 AER 1304 ÖLG - Instantaneous flame fronts (left) and turbulent flame brush envelope (right).
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationModeling of turbulence in stirred vessels using large eddy simulation
Modeling of turbulence in stirred vessels using large eddy simulation André Bakker (presenter), Kumar Dhanasekharan, Ahmad Haidari, and Sung-Eun Kim Fluent Inc. Presented at CHISA 2002 August 25-29, Prague,
More information2. Conservation Equations for Turbulent Flows
2. Conservation Equations for Turbulent Flows Coverage of this section: Review of Tensor Notation Review of Navier-Stokes Equations for Incompressible and Compressible Flows Reynolds & Favre Averaging
More informationTutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace
Tutorial School on Fluid Dynamics: Aspects of Turbulence Session I: Refresher Material Instructor: James Wallace Adapted from Publisher: John S. Wiley & Sons 2002 Center for Scientific Computation and
More informationA Simple Turbulence Closure Model
A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: Mean velocity: Turbulent velocity: Gradient operator: Advection operator: V =
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationModeling flame brush thickness in premixed turbulent combustion
Center for Turbulence Research Proceedings of the Summer Program 2006 299 Modeling flame brush thickness in premixed turbulent combustion By E. Knudsen, O. Kurenkov, S. Kim, M. Oberlack AND H. Pitsch Turbulent
More informationLecture 9 Laminar Diffusion Flame Configurations
Lecture 9 Laminar Diffusion Flame Configurations 9.-1 Different Flame Geometries and Single Droplet Burning Solutions for the velocities and the mixture fraction fields for some typical laminar flame configurations.
More informationA Simple Turbulence Closure Model. Atmospheric Sciences 6150
A Simple Turbulence Closure Model Atmospheric Sciences 6150 1 Cartesian Tensor Notation Reynolds decomposition of velocity: V = V + v V = U i + u i Mean velocity: V = Ui + V j + W k =(U, V, W ) U i =(U
More informationSimulations for Enhancing Aerodynamic Designs
Simulations for Enhancing Aerodynamic Designs 2. Governing Equations and Turbulence Models by Dr. KANNAN B T, M.E (Aero), M.B.A (Airline & Airport), PhD (Aerospace Engg), Grad.Ae.S.I, M.I.E, M.I.A.Eng,
More informationChapter 7. Basic Turbulence
Chapter 7 Basic Turbulence The universe is a highly turbulent place, and we must understand turbulence if we want to understand a lot of what s going on. Interstellar turbulence causes the twinkling of
More informationProcess Chemistry Toolbox - Mixing
Process Chemistry Toolbox - Mixing Industrial diffusion flames are turbulent Laminar Turbulent 3 T s of combustion Time Temperature Turbulence Visualization of Laminar and Turbulent flow http://www.youtube.com/watch?v=kqqtob30jws
More informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
More informationPublished in: Proceedings of the Fluent Benelux User Group Meeting, 6-7 October 2005, Wavre, Belgium
The application of Flamelet Generated Manifolds in partailly-premixed flames Ramaekers, W.J.S.; Albrecht, B.A.; van Oijen, J.A.; de Goey, L.P.H.; Eggels, R.L.G.M. Published in: Proceedings of the Fluent
More informationThe Johns Hopkins Turbulence Databases (JHTDB)
The Johns Hopkins Turbulence Databases (JHTDB) HOMOGENEOUS BUOYANCY DRIVEN TURBULENCE DATA SET Data provenance: D. Livescu 1 Database Ingest and Web Services: C. Canada 1, K. Kalin 2, R. Burns 2 & IDIES
More informationMath 575-Lecture Viscous Newtonian fluid and the Navier-Stokes equations
Math 575-Lecture 13 In 1845, tokes extended Newton s original idea to find a constitutive law which relates the Cauchy stress tensor to the velocity gradient, and then derived a system of equations. The
More informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More informationLecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory
Lecture 8 Laminar Diffusion Flames: Diffusion Flamelet Theory 8.-1 Systems, where fuel and oxidizer enter separately into the combustion chamber. Mixing takes place by convection and diffusion. Only where
More informationProject Topic. Simulation of turbulent flow laden with finite-size particles using LBM. Leila Jahanshaloo
Project Topic Simulation of turbulent flow laden with finite-size particles using LBM Leila Jahanshaloo Project Details Turbulent flow modeling Lattice Boltzmann Method All I know about my project Solid-liquid
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationChapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows
Chapter 7 The Time-Dependent Navier-Stokes Equations Turbulent Flows Remark 7.1. Turbulent flows. The usually used model for turbulent incompressible flows are the incompressible Navier Stokes equations
More informationA Priori Model for the Effective Lewis Numbers in Premixed Turbulent Flames
Paper # 070LT-0267 Topic: Turbulent Flames 8 th US National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013.
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationRegularization modeling of turbulent mixing; sweeping the scales
Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationDr.-Ing. Frank Beyrau Content of Lecture
Content of Lecture 1. Phenomenology of Combustion 2. Thermodynamic Fundamentals 3. Chemical Reaction Kinetics 4. Ignition and Ignition Limits 5. Laminar Flame Theory 6. Turbulent Combustion 7. Pollutants
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationScalar gradient and small-scale structure in turbulent premixed combustion
Center for Turbulence Research Annual Research Briefs 6 49 Scalar gradient and small-scale structure in turbulent premixed combustion By S. H. Kim AND H. Pitsch. Motivation and objectives The scalar gradient
More information7. Basics of Turbulent Flow Figure 1.
1 7. Basics of Turbulent Flow Whether a flow is laminar or turbulent depends of the relative importance of fluid friction (viscosity) and flow inertia. The ratio of inertial to viscous forces is the Reynolds
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationReaction Rate Closure for Turbulent Detonation Propagation through CLEM-LES
5 th ICDERS August 7, 05 Leeds, UK through CLEM-LES Brian Maxwell, Matei Radulescu Department of Mechanical Engineering, University of Ottawa 6 Louis Pasteur, Ottawa, KN 6N5, Canada Sam Falle School of
More informationFlow and added small-scale topologies in a turbulent premixed flame
Flow and added small-scale topologies in a turbulent premixed flame L. Cifuentes*, A. Kempf* and C. Dopazo** luis.cifuentes@uni-due.de *University of Duisburg-Essen, Chair of Fluid Dynamics, Duisburg -
More informationTurbulence and its modelling. Outline. Department of Fluid Mechanics, Budapest University of Technology and Economics.
Outline Department of Fluid Mechanics, Budapest University of Technology and Economics October 2009 Outline Outline Definition and Properties of Properties High Re number Disordered, chaotic 3D phenomena
More informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal and Vijaya Rani PI: Sarma L. Rani Department of Mechanical and Aerospace
More informationTurbulence Instability
Turbulence Instability 1) All flows become unstable above a certain Reynolds number. 2) At low Reynolds numbers flows are laminar. 3) For high Reynolds numbers flows are turbulent. 4) The transition occurs
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationInsights into Model Assumptions and Road to Model Validation for Turbulent Combustion
Insights into Model Assumptions and Road to Model Validation for Turbulent Combustion Venke Sankaran AFRL/RQR 2015 AFRL/RQR Basic Research Review UCLA Jan 20, 2015 AFTC PA Release# 15011, 16 Jan 2015 Goals
More informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal PI: Sarma L. Rani Department of Mechanical and Aerospace Engineering The
More informationPDF Modeling and Simulation of Premixed Turbulent Combustion
Monte Carlo Methods Appl. Vol. No. (), pp. 43 DOI 5 / MCMA.7. c de Gruyter PDF Modeling and Simulation of Premixed Turbulent Combustion Michael Stöllinger and Stefan Heinz Abstract. The use of probability
More informationLES Approaches to Combustion
LES Approaches to combustion LES Approaches to combustion LES Approaches to Combustion W P Jones Department of Mechanical Engineering Imperial College London Exhibition Road London SW7 2AZ SIG on Combustion
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationTurbulent Rankine Vortices
Turbulent Rankine Vortices Roger Kingdon April 2008 Turbulent Rankine Vortices Overview of key results in the theory of turbulence Motivation for a fresh perspective on turbulence The Rankine vortex CFD
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationCharacteristics of Linearly-Forced Scalar Mixing in Homogeneous, Isotropic Turbulence
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-1103 Characteristics of Linearly-Forced Scalar Mixing in Homogeneous, Isotropic Turbulence
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationEffects of turbulence and flame instability on flame front evolution
PHYSICS OF FLUIDS 18, 104105 2006 Effects of turbulence and flame instability on flame front evolution Jiao Yuan, Yiguang Ju, and Chung K. Law a Department of Mechanical and Aerospace Engineering, Princeton
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
More informationLecture 4: The Navier-Stokes Equations: Turbulence
Lecture 4: The Navier-Stokes Equations: Turbulence September 23, 2015 1 Goal In this Lecture, we shall present the main ideas behind the simulation of fluid turbulence. We firts discuss the case of the
More informationLecture 2. Turbulent Flow
Lecture 2. Turbulent Flow Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of this turbulent water jet. If L is the size of the largest eddies, only very small
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationConsistent turbulence modeling in a hybrid LES/RANS PDF method for non-premixed flames
Consistent turbulence modeling in a hybrid LES/RANS PDF method for non-premixed flames F. Ferraro, Y. Ge and M. Pfitzner Institut für Thermodynamik, Fakultät für Luft- und Raumfahrrttechnik Universität
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 19 Turbulent Flows Fausto Arpino f.arpino@unicas.it Introduction All the flows encountered in the engineering practice become unstable
More informationRepresenting Turbulence/Chemistry interaction with strained planar premixed
Representing Turbulence/Chemistry interaction with strained planar premixed flames CES Seminar Andreas Pesch Matr. Nr. 300884 Institut für technische Verbrennung Prof. H. Pitsch Rheinisch-Westfälische
More informationPDF modeling and simulation of premixed turbulent combustion
Monte Carlo Methods Appl. Vol. 4 No. 4 (8), pp. 343 377 DOI. / MCMA.8.6 c de Gruyter 8 PDF modeling and simulation of premixed turbulent combustion Michael Stöllinger and Stefan Heinz Abstract. The use
More informationA Ghost-fluid method for large-eddy simulations of premixed combustion in complex geometries
Center for Turbulence Research Annual Research Briefs 2005 269 A Ghost-fluid method for large-eddy simulations of premixed combustion in complex geometries By V. Moureau, P. Minot, C. Bérat AND H. Pitsch
More informationTowards regime identification and appropriate chemistry tabulation for computation of autoigniting turbulent reacting flows
Center for Turbulence Research Annual Research Briefs 009 199 Towards regime identification and appropriate chemistry tabulation for computation of autoigniting turbulent reacting flows By M. Kostka, E.
More informationMulti-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah
Multi-Scale Modeling of Turbulence and Microphysics in Clouds Steven K. Krueger University of Utah 10,000 km Scales of Atmospheric Motion 1000 km 100 km 10 km 1 km 100 m 10 m 1 m 100 mm 10 mm 1 mm Planetary
More informationExtinction and reignition in a diffusion flame: a direct numerical simulation study
J. Fluid Mech. (24), vol. 518, pp. 231 259. c 24 Cambridge University Press DOI: 1.117/S22112414 Printed in the United Kingdom 231 Extinction and reignition in a diffusion flame: a direct numerical simulation
More informationLES of the Sandia Flame D Using an FPV Combustion Model
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 82 (2015 ) 402 409 ATI 2015-70th Conference of the ATI Engineering Association LES of the Sandia Flame D Using an FPV Combustion
More information36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson
36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through
More informationBest Practice Guidelines for Combustion Modeling. Raphael David A. Bacchi, ESSS
Best Practice Guidelines for Combustion Modeling Raphael David A. Bacchi, ESSS PRESENTATION TOPICS Introduction; Combustion Phenomenology; Combustion Modeling; Reaction Mechanism; Radiation; Case Studies;
More informationThe mean shear stress has both viscous and turbulent parts. In simple shear (i.e. U / y the only non-zero mean gradient):
8. TURBULENCE MODELLING 1 SPRING 2019 8.1 Eddy-viscosity models 8.2 Advanced turbulence models 8.3 Wall boundary conditions Summary References Appendix: Derivation of the turbulent kinetic energy equation
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationModeling Turbulent Combustion
Modeling Turbulent Combustion CEFRC Combustion Summer School 2014 Prof. Dr.-Ing. Heinz Pitsch Copyright 2014 by Heinz Pitsch. This material is not to be sold, reproduced or distributed without prior written
More informationFLAME AND EDDY STRUCTURES IN HYDROGEN AIR TURBULENT JET PREMIXED FLAME
FLAME AND EDDY STRUCTURES IN HYDROGEN AIR TURBULENT JET PREMIXED FLAME M. Shimura, K. Yamawaki, Y.-S. Shim, M. Tanahashi and T. Miyauchi Department of Mechanical and Aerospace Engineering Tokyo Institute
More informationJ. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and
J. Szantyr Lecture No. 4 Principles of the Turbulent Flow Theory The phenomenon of two markedly different types of flow, namely laminar and turbulent, was discovered by Osborne Reynolds (184 191) in 1883
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationCorso di Laurea in Ingegneria Aerospaziale. Performance of CFD packages for flow simulations in aerospace applications
Universitá degli Studi di Padova Dipartimento di Ingegneria Industriale DII Department of Industrial Engineering Corso di Laurea in Ingegneria Aerospaziale Performance of CFD packages for flow simulations
More informationA priori Tabulation of Turbulent Flame Speeds via a Combination of a Stochastic Mixing Model and Flamelet Generated Manifolds 5
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany HEIKO SCHMIDT 1 MICHAEL OEVERMANN 2 ROB J.M. BASTIAANS 3 ALAN R. KERSTEIN 4 A priori Tabulation of Turbulent
More informationMax Planck Institut für Plasmaphysik
ASDEX Upgrade Max Planck Institut für Plasmaphysik 2D Fluid Turbulence Florian Merz Seminar on Turbulence, 08.09.05 2D turbulence? strictly speaking, there are no two-dimensional flows in nature approximately
More information1 Mean Flow Equations
ME 543 Averaged Equations 23 January 2018 These notes are meant as a supplement to the text, and as an aid for the class notes, but not as an independent document. 1 Mean Flow Equations 1.1 Mean flow equations
More informationDNS and LES of Turbulent Combustion
Computational Fluid Dynamics In Chemical Reaction Engineering IV June 19-24, 2005 Barga, Italy DNS and LES of Turbulent Combustion Luc Vervisch INSA de Rouen, IUF, CORIA-CNRS Pascale Domingo, Julien Réveillon
More informationSurvey of Turbulent Combustion Models for Large Eddy Simulations of Propulsive Flowfields
AIAA SciTech 5-9 January 2015, Kissimmee, Florida 53rd AIAA Aerospace Sciences Meeting AIAA 2015-1379 Survey of Turbulent Combustion Models for Large Eddy Simulations of Propulsive Flowfields Downloaded
More informationNonequilibrium Dynamics in Astrophysics and Material Science YITP, Kyoto
Nonequilibrium Dynamics in Astrophysics and Material Science 2011-11-02 @ YITP, Kyoto Multi-scale coherent structures and their role in the Richardson cascade of turbulence Susumu Goto (Okayama Univ.)
More informationIntroduction of compressible turbulence
Introduction of compressible turbulence 1 Main topics Derive averaged equations for compressible turbulence Introduce a math. technique to perform averaging in presence of density variation Favre average
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationLarge Eddy Simulation of Bluff Body Stabilized Premixed Flame
44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 2006, Reno, Nevada AIAA 2006-152 Large Eddy Simulation of Bluff Body Stabilized Premixed Flame Ionut Porumbel and Suresh Menon School of Aerospace
More informationIMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE FLUCTUATIONS
Proceedings of the Combustion Institute, Volume 9, 00/pp. 1971 1978 IMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE
More informationNote the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size
L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous
More informationThermal NO Predictions in Glass Furnaces: A Subgrid Scale Validation Study
Feb 12 th 2004 Thermal NO Predictions in Glass Furnaces: A Subgrid Scale Validation Study Padmabhushana R. Desam & Prof. Philip J. Smith CRSIM, University of Utah Salt lake city, UT-84112 18 th Annual
More information4. DETAILED MODELLING OF FA AND DDT
4. DETAILED MODELLING OF FA AND DDT 4.1 Introduction 4.1.1 On the Notion of State of the Art : Research versus Applications The state of the art in numerical turbulent combustion modelling is a notion
More informationRayleigh-Taylor Unstable Flames
Rayleigh-Taylor Unstable Flames Elizabeth P. Hicks 1,2 and Robert Rosner 2 CIERA, Northwestern University 1 University of Chicago 2 CIERA Conference: September 2, 2011 1 Type Ia Supernovae Image: NASA
More informationIsotropic homogeneous 3D turbulence
Chapter 6 Isotropic homogeneous 3D turbulence Turbulence was recognized as a distinct fluid behavior by Leonardo da Vinci more than 500 years ago. It is Leonardo who termed such motions turbolenze, and
More informationmeters, we can re-arrange this expression to give
Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.
More information