Lecture 12. Droplet Combustion Spray Modeling. Moshe Matalon
|
|
- Ami Crawford
- 5 years ago
- Views:
Transcription
1 Lecture 12 Droplet Combustion Spray Modeling Spray combustion: Many practical applications liquid fuel is injected into the combustion chamber resulting in fuel spray. Spray combustion involves many physical processes, including atomization, oplet collision and agglomeration, vaporization, heat and mass transfer, oplet-air and vapor-air mixing, ignition, turbulence, premixed and/or di usion flames, pollutant production and flame extinction. The liquid fuel is atomized by the combined action of aerodynamical shear, strain and surface tension producing a large number of small oplets which increase the overall surface area exposed to the oxidizer and enhance the rates of heat and mass transfer. The atomization of a 1cm diameter oplet of liquid into oplets of 100 µm diameter, for example, produces a million oplets and increases the overall surface area by a factor of University of Illinois at Urbana- Champaign 1
2 The oplets distribution range from few microns to around 500µm, primarily because practical nozzles cannot produce sprays of uniform op size at typical operating conditions. Furthermore, many of the larger oplets produced in the initial liquid disintegration undergo further breakup into smaller oplets. Depending on the oplet distribution, the mode of operation and the state of the ambience, di erent regimes of burning occur, including a di usion flame lies outside the cloud of oplets burning the ambient oxygen with the fuel that originates from the vaporizing oplets oplet burning occurring all throughout the cloud or in a layer on the edge of the cloud, internal group combustion with an internal flame that separates a group of vaporizing oplets from a group of individually burning oplets individual oplet combustion where individual oplets burn with oxygen di using across the resulting cloud of burning oplets oxidizer entrainment atomizer vaporization oplet burning flame Liquid jet liquid S P R A Y atomization oxidizer entrainment cloud burning University of Illinois at Urbana- Champaign 2
3 Droplet Vaporization and Combustion: Assumptions: oplet is a sphere spherical symmetry single component fuel uniform properties inside the oplet Fuel vapor op no fluid motion inside the oplet quiescent ambience no gravity quasi-steady approximation one-step overall chemical reaction (F + O! Products) Overall mass conservation of the oplet d dt 4 3 lr 3 s = ṁ mass flow rate mass per unit time leaving the oplet 4 r 2 s l s dt =4 r2 s g (v s dt ) Since g / l 1 (typically ) the oplet recedes very slowly compared to the gas-phase processes, and s dt g l v a direct consequence of the jump relation [ (v n V I ) ]] = 0 across the interface r = r s, namely l ṙ s = g (v ṙ s ) after using g / l 1. quasi-steady approximation ṁ or g v can be determined from the steady equations in the gas phase, with the oplet history determined a-posteriori from the relation above University of Illinois at Urbana- Champaign 3
4 Conditions for the mass fractions at the liquid-gas interface can be derived by integrating the species equations across the moving liquid-gas interface, r = r s, or using the previously derived jump condition [[ Y i ((v+v i ) n V I )]]=0 where we have assumed that no reaction occurs along the interface Using the approximation ṙ s ( g / l )v g Y i (v ṙ s ) g D i dy i = ly i ṙ s g vy F g vy O g D F dy F = gv g D O dy O =0 Similarly, integrating the energy (expressed in terms of enthalpy) across the moving liquid-gas interface, r = r s, or using the previously derived jump condition [[ h(v n V I )+pv I + q n ]] = 0 g h g (v ṙ s ) g dt dt = l h l ṙ s l a + a where we have assumed that no reaction occurs along the interface and neglected kinetic energy and viscous dissipation in the energy equation (the small Mach number approximation). Then, using the approximation ṙ s ( g / l )v g v (h g h l ) {z } Q v latent heat of vaporization dt a + dt a =0 University of Illinois at Urbana- Champaign 4
5 g vq v dt a + dt a =0 assuming all the heat conducted inwards goes to heating the liquid fuel, dt 4 a 2 4 dt 3 a3 l c pl dt dt = vq v + a 3 dt lc pl dt {z } l often neglected. a Alternatively it may be accounted for by replacing Q v with an e ective value that accounts for internal heating. dt = v Q v In the following, the gas density will be denoted by (with no subscript) and the oplet radius a. The steady conservation equations are d v dy F v dy O (r2 v) =0 d r 2 dy F D F r 2 D O r 2 d v dt 1 c p r 2 r 2 dy O d r 2 dt = F W F! = O W O! = Q c p! a liquid Conditions far away, as r!1 T = P 0 R/W YF YO! = B e E/RT W F W O T = T 1, = 1 Y F =0 Y O = Y O 1 University of Illinois at Urbana- Champaign 5
6 If one assumes that the phase change at the interface occurs very fast, equilibrium conditions prevail and a definite relation exists between the vapor partial pressure and the oplet surface temperature. This relation is of the form py F = k exp Qv R 1 1 T B T known as the Claussius Clapeyron relation. Instead we shall use for simplicity that the surface temperature is constant, and at r = a T = T B Boundary conditions at r = a vy F vy O D F dy F = v D O dy O =0 dt = vq v T = T B We write dimensionless equations using a as a unit of distance, D th /a as a unit of velocity and normalize the density and temperature with respect to their ambient values 1 and T 1. We also use M = ṁ/4 a( /c p ) to denote the dimensionless mass flux. M = r 2 v T =1 M dy F r 2 M dy O r 2 M dt r 2 Le 1 F r 2 Le 1 O r 2 d d 1 d r 2 r 2 dy F r 2 dy O r 2 dt! = D 2 Y F Y O e 0/T =! =! = q! at r = 1: MY F Le 1 dy F F = M MY O Le 1 dy O O =0 as r!1 dt = ML v, T = T s T = =1, Y O = Y O 1 Y F =0 We have a 6th-order system, with 7 BCs; the extra one for the determination of M. University of Illinois at Urbana- Champaign 6
7 The parameters are = O W O / F W F 0 = E/RT 1 q = Q/ F W F c p T 1 D = 2 1a 2 F B ( /c p )W O L v = Q v /c p T 1 T s = T B /T 1 We note in, particular, the dependence on the Damköhler number on the oplet radius D a 2 and on the ambient pressure D P 0. We also assume unity Lewis numbers. Pure vaporization: D =0 Y F =1 e M/r Y O = Y O 1 e M/r M(1 1/r) T = T s L v + L v e M =ln 1+ 1 T s L v In dimensional form, the vaporization rate is ṁ = 4 a c p ln 1+ c p(t 1 T B ) Q v transfer number B v = c p(t 1 T B) Q v = impetus for transfer resistance to transfer Clearly B v > 0, which implies that, for vaporization, T 1 >T B. University of Illinois at Urbana- Champaign 7
8 T 1 Y F =1 T T B Y F 1 r Combustion: We will consider the Burke-Schumann limit D!1and make use of the coupling functions. T + qy F =1+(T s 1+q L v )(1 e M/r ) T + q 1 Y O =(T s L v )(1 e M/r )+(1+q 1 Y O 1 )e M/r from which we can also write Y F 1 Y O =1 (1 + 1 Y O 1 )e M/r The profiles are readily available 8 < 1 (1 + 1 Y O 1 )e M/r (r<r f ) Y F = : 0 (r>r f ) 8 < 0 (r<r f ) Y O = : ( + Y O 1 )e M/r (r >r f ) University of Illinois at Urbana- Champaign 8
9 8 < T = : (T s L v )(1 e M/r )+(1+q 1 Y O 1 )e M/r (r<r f ) 1+(T s 1+q L v )(1 e M/r ) (r>r f ) There are few conditions that remain to be satisfied; these are T = T s at r = 1, and continuity of all variables at r = r f. The first provides and expression for M and the others determine the flame stando distance and the flame temperature. M =ln 1+ 1+q 1 Y O 1 T s L v In dimensional form, the burning rate is ṁ = 4 a c p ln 1+ c p(t 1 T B )+QY O 1 / O W O Q v B c = c p(t 1 T B )+QY O 1 / O W O Q v = impetus for transfer resistance to transfer transfer number The flame stando distance r f = a ln(1 + B c ) ln(1 + F W F Y O 1 / O W O ) Note that r f > 1 implies that c p (T 1 T B ) > (Q v Q/ F W F ) 1 Y O 1 which is not restrictive because typically Q>Q v. The flame temperature T f = T 1 + (Q/ F W F Q v )/c p (T 1 T B ) 1+ O W O / F W F Y O 1 University of Illinois at Urbana- Champaign 9
10 Y F =1 T Y O 1 T B Y O Y F T 1 1 r f r If the activation energy parameter is considered large, the asymptotic analysis provides a more complete picture spanning the entire range of Damköhler number D. M c BS limit T a Burning rate M v Vaporization limit D ext D ign Damköhler number D D ~ a 2, P 0 Extinction may occur as the oplet size diminishes below a critical value; auto-ignition may occur if the ambient pressure is increased significantly. University of Illinois at Urbana- Champaign 10
11 Droplet history: ( use back the notation a = r s and = g ) The vaporization and burning rates were found to have the form ṁ =4 r s ( /c p )ln(1+b) with the appropriate transfer number. ṁ =4 r 2 s g v ) g v = cp ln(1 + B) 1 r s s dt = g l v = k k = 2 ln(1 + B) with 2r s c p l 2r s s dt = k ) r2 s = r 2 s(0) kt d 2 -law d 2 = d 2 (0) Kt K = 8( /c p) l ln(1 + B) d 2 t life = l g d 2 0 8D th ln(1 + B) 1/K t Nuruzzaman et al. PCI, 1971 The d 2 -law prediction is verified experimentally for most of the oplet lifetime with the exception, perhaps, of a short initial period following ignition University of Illinois at Urbana- Champaign 11
12 The behavior of a single oplet in a combustion chamber depends on the Reynolds number Re = g d 0 u u p /µ g,whereuisacharacteristic velocity of the gas and u p the particle velocity. So far, the e ect of the relative velocity of the ambient gas on the oplet has been neglected, practically taking Re = 0. The spherical symmetry assumption cannot be justified when Re 6= 0. A general approach that has been often used to adess the problem of oplet combustion in a convective field has been to correct a-posteriori the results derived for Re = 0 with empirical correlations. Formulas that have been used are of the form K = K (0) 1+CRe 1/2 Pr 1/3 C 0.3 where K is the evaporation constant, K (0) the evaporation constant for a stationary burning oplet, Pr = µ g /D th g is the Prandtl number. Similarly, burning rate corrections for free convection has been proposed, in the form K = K (0) 1+CGr 0.52 C 0.5 where Gr is the Grashof number, representing the ratio of the buoyancy to viscous forces. Gr = 2 gd 3 0 g T µ 2 g T Here T may be taken as the di erence between the flame temperature and the ambient temperature. The assumption of spherical symmetry for suspended oplets can be justified either in micro-gravity conditions (g 1) or in reduced pressure environment (P 0 g 1) University of Illinois at Urbana- Champaign 12
13 Combustion of a solid particle: Similar to a fuel oplet a solid particle may be surrounded by a di usion flame, so that ṁ = 4 r s c p ln (1 + B) if the flame is close enough to the surface, r f 1 ln (1 + B) =ln 1+ 1 Y O r s solid ṁ =4 r s g D O ln 1+ 1 Y O This is the di usion-controlled limit, where the di usion of oxidizer towards the particle controls the burning. For Le = 1, we have /c p = D O. Using ṁ = d dt 4 3 srs 3 = 4 s rs 2 s dt ṁ =4 r s g D O ln 1+ 1 Y O we find in the di usion-controlled limit r 2 s = r 2 s(0) 2 g D O s ln (1 + 1 Y O ) t d 2 = d g D O s ln (1 + 1 Y O ) t t life = s g d 2 0 8D O ln(1+ 1 Y O ) namely a d 2 -law. University of Illinois at Urbana- Champaign 13
14 If burning occurs on the surface (and possibly inside the pores) of the particle, the burning rate is controlled by the chemical kinetics. For an Arrhenius-type heterogenous surface chemical reaction rate, we have ṁ =4 r 2 sk s Y O r s where k s = B s e E/RTs solid kinetic-controlled burning This is the kinetic-controlled limit, where the chemical reaction rate controls the burning. Using ṁ = d 4 dt 3 srs 3 = 4 s r 2 s s dt k s r s = r s (0) Y O t s d = d 0 2k s g s Y O t namely a d 1 -law (linear). t life = s g d 0 2B s Y O e E/RT t life = 8 >< >: s g s g d 2 0 8D O ln(1+ 1 Y O ) d 0 2B s Y O e E/RTs di usion-controlled kinetic-controlled t independent of T t d 2 0 t e E/RT t d 0 The controlling mechanism is determined by a surface Damköhler number D s = t di usion t kinetic = d 0 B s Y O e E/RTs 4D O ln(1+ 1 Y O ) large particles, or higher ambient temperature ) D s large ) di usion-controlled small particles, or low ambient temperature ) D s small ) kinetic-controlled University of Illinois at Urbana- Champaign 14
15 Spray Modeling: The accurate description of spray combustion, say by DNS, requires identification of whether liquid or gas is present at every point in space, to employ appropriate equations of state in each of the separate phases and to determine the evolution of a large number of liquid-gas interfaces with boundary conditions that involve phase transformation. The large number of oplets required to simulate real dispersed flows in turbulent environments that are encountered in practical applications, make such calculations prohibiitve. There are a practically three approaches to spray combustion modeling. Continuum formulation Lagrangian formulation Probabilistic formulation We will not be discussing any of these in detail, but we will try to get a general idea of what each of these formulations involve. Continuum Formulation. Properties of both the particles and the gas are assumed to exist at every spatial point (the average value of a small parcel that includes both phases) regardless of whether that point is actually in the liquid, solid or gas at that instant. The spray contains N distinct chemical species and M di erent class of oplets identified by fuel content, size and/or velocity. Let m i be the mass of species i and ˆm (k) the mass of oplets of class k. The total volume of the mixture will be denoted by V T and will be given by V T = V g + V l,wherev g and V l are the volumes occupied by the gas and the liquid phases, respectively. ' = V g /V T is the volume fraction of the gas, namely the ratio of the volume occupied by the gas to the total volume of the mixture. 1 ' is the void fraction, or the ratio of the volume occupied by the condensed phase from the total volume V T. University of Illinois at Urbana- Champaign 15
16 The bulk density (mass per total volume V T ) of species i is i. NX The bulk density of the gaseous mixture is = i = 1 NX V T It di ers form the gas density g = 1 V g N X unit volume that includes only gas. Similarly, the density of oplets of class k is ˆ (k), so that the bulk density of oplets is MX ˆ = ˆ (k) = 1 MX ˆm (k) V T k=1 which di ers from the liquid density l = 1 V l M X k=1 i=1 k=1 ˆm (k) = V T V l ˆ = 1 1 ' ˆ i=1 m i = V T V g = 1 ' i=1 m i which is the mass per For a dilute spray the condensed phase volume fraction is small (1 that g,but ˆ =(1 ') l 1 ' = O(1) because g / l 1. g / l ˆ is often referred to as the spray density. ' 1) so Statements of mass conservation of the gas and condensed phases + r v = + r ˆ ˆv = Ṁ Here Ṁ is the total mass (per unit volume, per unit time) added to the gas phase due to vaporization and v, ˆv are the gas and particle velocities. This can be determined, for example, from the expressions we derived for the vaporization rate of isolated oplets. ( +ˆ )+r ( v +ˆ ˆv) i.e., the total mass is conserved. Using ˆ =(1 ') l, we can recast the mass conservation of the oplets in the following + r (ˆv') =Ṁ + r l ' is no longer conserved; the volume occupied by the oplets changes as a result of (i)vaporization, and (ii) distance between oplets University of Illinois at Urbana- Champaign 16
17 Similarly, separate statements of momentum and energy conservation are written for the gas and liquid phases with appropriate source/sink terms representing the total momentum and energy carried to the gas from the vaporizing oplets, the total force exerted on the particles from their surroundings and the internal heating of the oplets which, in general, could have a di erent temperature than the bulk gas. Probabilistic Formulation. In this approach the condensed phase is described by a distribution function, or probability density function (pdf) f(r, x, v,e,t) fdx dv de... is the probable number of particles in the radius range about r, in the spatial range dx about the position x, with velocities in the range dv about v and internal energy in the range e about de, etc..., at time t. The time evolution of the distribution function f is + r x v (ėf)+ = and is known as the spray equation. the terms on the l.h.s. of the spray equation represent op-carrier gas interactions, the source term Q on the r.h.s. represents particle-particle interactions (fragmentation, coalescence, etc..). University of Illinois at Urbana- Champaign 17
18 The quantities in the spray equation are determined by the laws governing the dynamics of individual oplets dx dt = v, dv dt = a, dt =ṙ, de dt =ė, etc... The conservation laws for the gas-oplet mixture will include additional source terms involving the integral of f. For example, the mass and momentum equations are Z f For example, the dependence of ṙ on the radius for liquid ops can be specified by adopting the d 2 -law /c p ln (1 + B) ṙ =, + u ru + r ( f u)= rp + f g 4 r 2 l ṙf dvde, Z 4 3 r3 l a +4 r 2 l ṙ(v u) fdvde, where f is the fluid density, namely the mass of the gas mixture per unit volume of space, which due to the presence of the particles di ers from the gas density g that corresponds to the mass per unit volume occupied by the gas. The main di culty in this approach is the high dimensionality of the distribution function f(r, x, v, e,t) which, under the simplest assumptions, contains nine variables. Lagrangian Discrete Particle Formulation. The oplets are considered as individual point-particles and their trajectories are solved for, separately from the flow field which exists only where there are no oplets. The main di culty is probably the inability to track a su ciently large number of particles in order to University of Illinois at Urbana- Champaign 18
Lecture 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 informationTopics in Other Lectures Droplet Groups and Array Instability of Injected Liquid Liquid Fuel-Films
Lecture Topics Transient Droplet Vaporization Convective Vaporization Liquid Circulation Transcritical Thermodynamics Droplet Drag and Motion Spray Computations Turbulence Effects Topics in Other Lectures
More informationFlame Characteristics of Commercial Fuels With Respect to Single Droplet Combustion - For Spray Applications
International Journal of Theoretical and Applied Mechanics. ISSN 973-685 Volume 1, Number 4 (17) pp. 681-697 Research India Publications http://www.ripublication.com Flame Characteristics of Commercial
More informationTheoretical Developments in Group Combustion of Droplets and Sprays
Theoretical Developments in Group Combustion of Droplets and Sprays William A. Sirignano University of California, Irvine Collaborations: Guang Wu, current student; Randall Imaoka, former student, US Navy;
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 informationConvective Heat and Mass Transfer Prof. A.W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay
Convective Heat and Mass Transfer Prof. A.W. Date Department of Mechanical Engineering Indian Institute of Technology, Bombay Module No. # 01 Lecture No. # 32 Stefan Flow Model We are now familiar with
More informationLaminar Premixed Flames: Flame Structure
Laminar Premixed Flames: Flame Structure Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Course Overview Part I: Fundamentals and Laminar Flames Introduction Fundamentals and mass balances of
More informationSummary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer
1. Nusselt number Summary of Dimensionless Numbers of Fluid Mechanics and Heat Transfer Average Nusselt number: convective heat transfer Nu L = conductive heat transfer = hl where L is the characteristic
More informationHeat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay
Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 18 Forced Convection-1 Welcome. We now begin our study of forced convection
More information2σ e s (r,t) = e s (T)exp( rr v ρ l T ) = exp( ) 2σ R v ρ l Tln(e/e s (T)) e s (f H2 O,r,T) = f H2 O
Formulas/Constants, Physics/Oceanography 4510/5510 B Atmospheric Physics II N A = 6.02 10 23 molecules/mole (Avogadro s number) 1 mb = 100 Pa 1 Pa = 1 N/m 2 Γ d = 9.8 o C/km (dry adiabatic lapse rate)
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 informationLecture 7 Detonation Waves
Lecture 7 etonation Waves p strong detonation weak detonation weak deflagration strong deflagration / 0 v =/ University of Illinois at Urbana- Champaign eflagrations produce heat Thermal di usivity th
More informationModeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!
http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More information2 D. Terminal velocity can be solved for by equating Fd and Fg Fg = 1/6πd 3 g ρ LIQ = 1/8 Cd π d 2 ρ air u
ecture 8 Collection Growth of iquid Hydrometeors. I. Terminal velocity The terminal velocity of an object is the maximum speed an object will accelerate to by gravity. At terminal velocity, a balance of
More informationConvective Mass Transfer
Convective Mass Transfer Definition of convective mass transfer: The transport of material between a boundary surface and a moving fluid or between two immiscible moving fluids separated by a mobile interface
More informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationModeling of dispersed phase by Lagrangian approach in Fluent
Lappeenranta University of Technology From the SelectedWorks of Kari Myöhänen 2008 Modeling of dispersed phase by Lagrangian approach in Fluent Kari Myöhänen Available at: https://works.bepress.com/kari_myohanen/5/
More informationLecture 2 Constitutive Relations
Lecture 2 Constitutive Relations Conservation equations Dv = @ + r v =0 @t rp + r + g DY i + r ( Y iv i )= i i =1, 2,...,N { De = pr v + r q Dh + r q Z Y i p = R h = Y i h o i + c p d o W i We already
More informationChapter 10: Boiling and Condensation 1. Based on lecture by Yoav Peles, Mech. Aero. Nuc. Eng., RPI.
Chapter 10: Boiling and Condensation 1 1 Based on lecture by Yoav Peles, Mech. Aero. Nuc. Eng., RPI. Objectives When you finish studying this chapter, you should be able to: Differentiate between evaporation
More information1. Droplet Growth by Condensation
1. Droplet Growth by Condensation It was shown before that a critical size r and saturation ratio S must be exceeded for a small solution droplet to become a cloud droplet. Before the droplet reaches the
More informationOutline. Definition and mechanism Theory of diffusion Molecular diffusion in gases Molecular diffusion in liquid Mass transfer
Diffusion 051333 Unit operation in gro-industry III Department of Biotechnology, Faculty of gro-industry Kasetsart University Lecturer: Kittipong Rattanaporn 1 Outline Definition and mechanism Theory of
More informationLecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames
Lecture 6 Asymptotic Structure for Four-Step Premixed Stoichiometric Methane Flames 6.-1 Previous lecture: Asymptotic description of premixed flames based on an assumed one-step reaction. basic understanding
More informationUNIVERSITA DEGLI STUDI DI PADOVA. High performance computing of turbulent evaporating sprays
UNIVERSITA DEGLI STUDI DI PADOVA Dipartimento di Ingegneria Industriale DII Corso di Laurea Magistrale in Ingegneria Aerospaziale High performance computing of turbulent evaporating sprays Relatore: Prof.
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More informationPrinciples of Convective Heat Transfer
Massoud Kaviany Principles of Convective Heat Transfer Second Edition With 378 Figures Springer Contents Series Preface Preface to the Second Edition Preface to the First Edition Acknowledgments vii ix
More informationMechanical Engineering. Postal Correspondence Course HEAT TRANSFER. GATE, IES & PSUs
Heat Transfer-ME GATE, IES, PSU 1 SAMPLE STUDY MATERIAL Mechanical Engineering ME Postal Correspondence Course HEAT TRANSFER GATE, IES & PSUs Heat Transfer-ME GATE, IES, PSU 2 C O N T E N T 1. INTRODUCTION
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 informationModeling and Simulation of Fabrication, Transport, and Combustion of. Flow Environments
Modeling and Simulation of Fabrication, Transport, and Combustion of Nano-Engineered Energetic Materials (NEEM) in Flow Environments Vigor Yang The Pennsylvania State University Presented at Nano-Engineered
More informationMixing and Evaporation of Liquid Droplets Injected into an Air Stream Flowing at all Speeds
Mixing and Evaporation of Liquid Droplets Injected into an Air Stream Flowing at all Speeds F. Moukalled* and M. Darwish American University of Beirut Faculty of Engineering & Architecture Mechanical Engineering
More informationClouds associated with cold and warm fronts. Whiteman (2000)
Clouds associated with cold and warm fronts Whiteman (2000) Dalton s law of partial pressures! The total pressure exerted by a mixture of gases equals the sum of the partial pressure of the gases! Partial
More informationInteractive Processes in Sprays - NEL
Interactive Processes in Sprays - NEL During the week of the Study Group essentially all attention was focussed on the reactive and diffusive aspects of a burning spray; at this stage the velocity of the
More informationheat transfer process where a liquid undergoes a phase change into a vapor (gas)
Two-Phase: Overview Two-Phase two-phase heat transfer describes phenomena where a change of phase (liquid/gas) occurs during and/or due to the heat transfer process two-phase heat transfer generally considers
More informationConvection. forced convection when the flow is caused by external means, such as by a fan, a pump, or atmospheric winds.
Convection The convection heat transfer mode is comprised of two mechanisms. In addition to energy transfer due to random molecular motion (diffusion), energy is also transferred by the bulk, or macroscopic,
More informationThe effect of momentum flux ratio and turbulence model on the numerical prediction of atomization characteristics of air assisted liquid jets
ILASS Americas, 26 th Annual Conference on Liquid Atomization and Spray Systems, Portland, OR, May 204 The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization
More informationMULTIPHASE FLOW MODELLING
MULTIPHASE FLOW MODELLING 1 Introduction 2 Outline Multiphase Flow Modeling Discrete phase model Eulerian model Mixture model Volume-of-fluid model Reacting Flow Modeling Eddy dissipation model Non-premixed,
More informationSmoldering combustion of incense sticks - experiments and modeling
Smoldering combustion of incense sticks - experiments and modeling H. S. Mukunda*, J. Basani*, H. M. Shravan** and Binoy Philip*, May 30, 2007 Abstract This paper is concerned with the experimental and
More informationConvective Vaporization and Burning of Fuel Droplet Arrays
Convective Vaporization and Burning of Fuel Droplet Arrays Guang Wu and William A. Sirignano Mechanical & Aerospace Engineering University of California, Irvine May 19, 2009 Background and Task In many
More informationMixing and Combustion in Dense Mixtures by William A. Sirignano and Derek Dunn-Rankin
Mixing and Combustion in Dense Mixtures by William A. Sirignano and Derek Dunn-Rankin At very high pressures and densities, what is different and what is similar about the processes of Injection and Atomization,
More informationNUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP REGIME
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00394 NUMERICAL MODELING OF THE GAS-PARTICLE FLUID FLOW AND HEAT TRANSFER IN THE SLIP
More informationRocket Propulsion Prof. K. Ramamurthi Department of Mechanical Engineering Indian Institute of Technology, Madras
Rocket Propulsion Prof. K. Ramamurthi Department of Mechanical Engineering Indian Institute of Technology, Madras Lecture 32 Efficiencies due to Mixture Ratio Distribution and Incomplete Vaporization (Refer
More informationCombustion MATHEMATICAL MODEL FOR TRANSIENT. S. M. Frolov Λ,F.S.Frolov Λ, and B. Basara y
Combustion MATHEMATICAL MODEL FOR TRANSIENT DROPLET VAPORIZATION S. M. Frolov Λ,F.S.Frolov Λ, and B. Basara y Λ N. N. Semenov Institute of Chemical Physics Russian Academy of Sciences Moscow, Russia y
More information3D spray simulation using advanced intra-droplet and interface modeling
3D spray simulation using advanced intra-droplet and interface modeling TU Darmstadt Mechanical Engineering Simulation of reactive Thermo-Fluid Systems Andrea Pati, Christian Hasse Agenda Introduction
More informationCombustion basics... We are discussing gaseous combustion in a mixture of perfect gases containing N species indexed with k=1 to N:
Combustion basics... T. Poinsot poinsot@imft.fr Only the things you should know to understand the following courses Mainly elements of laminar flame theory for premixed and diffusion flames 1 Copyright
More informationSteady-State Molecular Diffusion
Steady-State Molecular Diffusion This part is an application to the general differential equation of mass transfer. The objective is to solve the differential equation of mass transfer under steady state
More informationChemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore
Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. # 26 Problem solving : Heterogeneous reactions Friends, in last few
More informationTravel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis
Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis DAVID KATOSHEVSKI Department of Biotechnology and Environmental Engineering Ben-Gurion niversity of the Negev
More informationConvection Heat Transfer. Introduction
Convection Heat Transfer Reading Problems 12-1 12-8 12-40, 12-49, 12-68, 12-70, 12-87, 12-98 13-1 13-6 13-39, 13-47, 13-59 14-1 14-4 14-18, 14-24, 14-45, 14-82 Introduction Newton s Law of Cooling Controlling
More informationComputer Fluid Dynamics E181107
Computer Fluid Dynamics E181107 2181106 Combustion, multiphase flows Remark: foils with black background could be skipped, they are aimed to the more advanced courses Rudolf Žitný, Ústav procesní a zpracovatelské
More informationNotes 4: Differential Form of the Conservation Equations
Low Speed Aerodynamics Notes 4: Differential Form of the Conservation Equations Deriving Conservation Equations From the Laws of Physics Physical Laws Fluids, being matter, must obey the laws of Physics.
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 3 LAMINAR BOUNDARY LAYER FLOW LAMINAR BOUNDARY LAYER FLOW Boundary
More information7.2 Sublimation. The following assumptions are made in order to solve the problem: Sublimation Over a Flat Plate in a Parallel Flow
7..1 Sublimation Over a Flat Plate in a Parallel Flow The following assumptions are made in order to solve the problem: 1.. 3. The flat plate is very thin and so the thermal resistance along the flat plate
More informationDroplet behaviour in a Ranque-Hilsch vortex tube
Journal of Physics: Conference Series Droplet behaviour in a Ranque-Hilsch vortex tube To cite this article: R Liew et al 2 J. Phys.: Conf. Ser. 38 523 View the article online for updates and enhancements.
More informationDNS of droplet evaporation and combustion in a swirling combustor
Center for Turbulence Research Annual Research Briefs 28 DNS of droplet evaporation and combustion in a swirling combustor By K. Luo, O. Desjardins AND H. Pitsch. Motivation and objective Turbulent multi-phase
More informationChapter 7: Natural Convection
7-1 Introduction 7- The Grashof Number 7-3 Natural Convection over Surfaces 7-4 Natural Convection Inside Enclosures 7-5 Similarity Solution 7-6 Integral Method 7-7 Combined Natural and Forced Convection
More informationSpray evaporation model sensitivities
Center for Turbulence Research Annual Research Briefs 20 23 Spray evaporation model sensitivities By Shashank, E. Knudsen AND H. Pitsch. Motivation and objective The energy density of solid- and liquid-phase
More informationSection 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow
Section 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow Anderson: Chapter 2 pp. 41-54 1 Equation of State: Section 1 Review p = R g T " > R g = R u M w - R u = 8314.4126
More informationFundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres
Fundamental Stellar Parameters Radiative Transfer Stellar Atmospheres Equations of Stellar Structure Basic Principles Equations of Hydrostatic Equilibrium and Mass Conservation Central Pressure, Virial
More informationDiffusional Growth of Liquid Phase Hydrometeros.
Diffusional Growth of Liquid Phase Hydrometeros. I. Diffusional Growth of Liquid Phase Hydrometeors A. Basic concepts of diffusional growth. 1. To understand the diffusional growth of a droplet, we must
More informationParcel Model. Atmospheric Sciences September 30, 2012
Parcel Model Atmospheric Sciences 6150 September 30, 2012 1 Governing Equations for Precipitating Convection For precipitating convection, we have the following set of equations for potential temperature,
More informationof Nebraska - Lincoln
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Combustion Research at University of Nebraska- Lincoln Mechanical & Materials Engineering, Department of March 2006 Raghavan,
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
More informationRogers and Yau Chapter 10: Drop breakup, snow, precip rate, and bulk models
Rogers and Yau Chapter 10: Drop breakup, snow, precip rate, and bulk models One explanation for the negative exponential (M-P) distribution of raindrops is drop breakup. Drop size is limited because increased
More informationDepartment of Mechanical Engineering BM 7103 FUELS AND COMBUSTION QUESTION BANK UNIT-1-FUELS
Department of Mechanical Engineering BM 7103 FUELS AND COMBUSTION QUESTION BANK UNIT-1-FUELS 1. Define the term fuels. 2. What are fossil fuels? Give examples. 3. Define primary fuels. Give examples. 4.
More informationINSTRUCTOR: PM DR MAZLAN ABDUL WAHID
SMJ 4463: HEAT TRANSFER INSTRUCTOR: PM ABDUL WAHID http://www.fkm.utm.my/~mazlan TEXT: Introduction to Heat Transfer by Incropera, DeWitt, Bergman, Lavine 5 th Edition, John Wiley and Sons Chapter 9 Natural
More informationENVIRONMENTAL FLUID MECHANICS
ENVIRONMENTAL FLUID MECHANICS Turbulent Jets http://thayer.dartmouth.edu/~cushman/books/efm/chap9.pdf Benoit Cushman-Roisin Thayer School of Engineering Dartmouth College One fluid intruding into another
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 information(Heat capacity c is also called specific heat) this means that the heat capacity number c for water is 1 calorie/gram-k.
Lecture 23: Ideal Gas Law and The First Law of Thermodynamics 1 (REVIEW) Chapter 17: Heat Transfer Origin of the calorie unit A few hundred years ago when people were investigating heat and temperature
More informationPHYSICAL MECHANISM OF NATURAL CONVECTION
1 NATURAL CONVECTION In this chapter, we consider natural convection, where any fluid motion occurs by natural means such as buoyancy. The fluid motion in forced convection is quite noticeable, since a
More informationModeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet Model with Heat-Loss Effects
9 th U. S. National Combustion Meeting Organized by the Central States Section of the Combustion Institute May 17-20, 2015 Cincinnati, Ohio Modeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet
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 informationThe Clausius-Clapeyron and the Kelvin Equations
PhD Environmental Fluid Mechanics Physics of the Atmosphere University of Trieste International Center for Theoretical Physics The Clausius-Clapeyron and the Kelvin Equations by Dario B. Giaiotti and Fulvio
More informationFluid Mechanics. du dy
FLUID MECHANICS Technical English - I 1 th week Fluid Mechanics FLUID STATICS FLUID DYNAMICS Fluid Statics or Hydrostatics is the study of fluids at rest. The main equation required for this is Newton's
More informationDiffusion and Adsorption in porous media. Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad
Diffusion and Adsorption in porous media Ali Ahmadpour Chemical Eng. Dept. Ferdowsi University of Mashhad Contents Introduction Devices used to Measure Diffusion in Porous Solids Modes of transport in
More informationCombustion. Indian Institute of Science Bangalore
Combustion Indian Institute of Science Bangalore Combustion Applies to a large variety of natural and artificial processes Source of energy for most of the applications today Involves exothermic chemical
More informationFORMULA SHEET. General formulas:
FORMULA SHEET You may use this formula sheet during the Advanced Transport Phenomena course and it should contain all formulas you need during this course. Note that the weeks are numbered from 1.1 to
More informationENGR Heat Transfer II
ENGR 7901 - Heat Transfer II Convective Heat Transfer 1 Introduction In this portion of the course we will examine convection heat transfer principles. We are now interested in how to predict the value
More informationDimensionless Numbers
1 06.10.2017, 09:49 Dimensionless Numbers A. Salih Dept. of Aerospace Engineering IIST, Thiruvananthapuram The nondimensionalization of the governing equations of fluid flow is important for both theoretical
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 informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
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
More informationReacting Gas Mixtures
Reacting Gas Mixtures Reading Problems 15-1 15-7 15-21, 15-32, 15-51, 15-61, 15-74 15-83, 15-91, 15-93, 15-98 Introduction thermodynamic analysis of reactive mixtures is primarily an extension of the principles
More informationME3560 Tentative Schedule Spring 2019
ME3560 Tentative Schedule Spring 2019 Week Number Date Lecture Topics Covered Prior to Lecture Read Section Assignment Prep Problems for Prep Probs. Must be Solved by 1 Monday 1/7/2019 1 Introduction to
More informationLiquids and solids are essentially incompressible substances and the variation of their density with pressure is usually negligible.
Properties of Fluids Intensive properties are those that are independent of the mass of a system i.e. temperature, pressure and density. Extensive properties are those whose values depend on the size of
More informationWell Stirred Reactor Stabilization of flames
Well Stirred Reactor Stabilization of flames Well Stirred Reactor (see books on Combustion ) Stabilization of flames in high speed flows (see books on Combustion ) Stabilization of flames Although the
More informationEvaluation of equilibrium and non-equilibrium evaporation models for many-droplet gas-liquid ow simulations
PERGAMON International Journal of Multiphase Flow 24 (1998) 1025±1055 Evaluation of equilibrium and non-equilibrium evaporation models for many-droplet gas-liquid ow simulations R. S. Miller, K. Harstad,
More informationDNS of Reacting H 2 /Air Laminar Vortex Rings
46th AIAA Aerospace Sciences Meeting and Exhibit 7-10 January 2008, Reno, Nevada AIAA 2008-508 DNS of Reacting H 2 /Air Laminar Vortex Rings Jeff Doom and Krishnan Mahesh University of Minnesota, Minneapolis,
More informationAN INVESTIGATION OF THE MESH DEPENDENCE OF THE STOCHASTIC DISCRETE DROPLET MODEL APPLIED TO DENSE LIQUID SPRAYS. Simone E.
AN INVESTIGATION OF THE MESH DEPENDENCE OF THE STOCHASTIC DISCRETE DROPLET MODEL APPLIED TO DENSE LIQUID SPRAYS By Simone E. Hieber A THESIS Submitted in partial fulfillment of the requirements for the
More informationNumerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure
21 st ICDERS July 23-27, 27 Poitiers, France Numerical Investigation of Ignition Delay in Methane-Air Mixtures using Conditional Moment Closure Ahmad S. El Sayed, Cécile B. Devaud Department of Mechanical
More informationA FLAMELET MODEL FOR DIFFUSION FLAMES IN POROUS MEDIA. Max Akira Endo Kokubun. Fernando Filho Fachini
A FLAMELET MODEL FOR DIFFUSION FLAMES IN POROUS MEDIA Max Akira Endo Kokubun Instituto Nacional de Pesquisas Espaciais, Rod. Presidente Dutra, km 40 - Cachoeira Paulista/SP, max@lcp.inpe.br Fernando Filho
More informationME3560 Tentative Schedule Fall 2018
ME3560 Tentative Schedule Fall 2018 Week Number 1 Wednesday 8/29/2018 1 Date Lecture Topics Covered Introduction to course, syllabus and class policies. Math Review. Differentiation. Prior to Lecture Read
More informationEvaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows
Evaluation of Liquid Fuel Spray Models for Hybrid RANS/LES and DLES Prediction of Turbulent Reactive Flows by Ali Afshar A thesis submitted in conformity with the requirements for the degree of Masters
More informationJET AND DROPLET BREAKUP MODELLING APPROACHES
Journal of KONES Powertrain and Transport, Vol. 22, No. 3 2015 JET AND DROPLET BREAKUP MODELLING APPROACHES Łukasz Jan Kapusta, Piotr Jaworski Warsaw University of Technology, Institute of Heat Engineering
More informationHeat processes. Heat exchange
Heat processes Heat exchange Heat energy transported across a surface from higher temperature side to lower temperature side; it is a macroscopic measure of transported energies of molecular motions Temperature
More informationEcoulements Diphasiques
Ecoulements Diphasiques Marc Massot Professor at Ecole Centrale Paris Laboratoire EM2C - UPR CNRS 288 - Ecole Centrale Paris Visiting Professor, Department of Mechanical Engineering, Center for Turbulence
More informationNUMERICAL PREDICTIONS OF DEPOSTION WITH A PARTICLE CLOUD TRACKING TECHNIQUE
Committed Individuals Solving Challenging Problems NUMERICAL PREDICTIONS OF DEPOSTION WITH A PARTICLE CLOUD TRACKING TECHNIQUE by James R. Valentine Reaction Engineering International Philip J. Smith Department
More informationCONTENTS Real chemistry e ects Scramjet operating envelope Problems
Contents 1 Propulsion Thermodynamics 1-1 1.1 Introduction.................................... 1-1 1.2 Thermodynamic cycles.............................. 1-8 1.2.1 The Carnot cycle.............................
More informationHeat Transfer Convection
Heat ransfer Convection Previous lectures conduction: heat transfer without fluid motion oday (textbook nearly 00 pages) Convection: heat transfer with fluid motion Research methods different Natural Convection
More informationDroplet Behavior within an LPP Ambiance
Copyright 2010 Tech Science Press FDMP, vol.6, no.4, pp.399-408, 2010 Droplet Behavior within an LPP Ambiance M. Chrigui 1,2, L. Schneider 1, A. Zghal 2, A. Sadiki 1 and J. Janicka 1 Abstract: This paper
More informationStep 1. Step 2. g l = g v. dg = 0 We have shown that over a plane surface of water. g v g l = ρ v R v T ln e/e sat. this can be rewritten
The basic question is what makes the existence of a droplet thermodynamically preferable to the existence only of water vapor. We have already derived an expression for the saturation vapor pressure over
More information