NUMERICAL INVESTIGATION OF THE TURBULENCE- COMBUSTION INTERACTION IN NON-PREMIXED CH 4 /AIR FLAMES
|
|
- Millicent Phillips
- 5 years ago
- Views:
Transcription
1 UMERICAL IVESTIGATIO OF THE TURBULECE- COMBUSTIO ITERACTIO I O-PREMIXED CH 4 /AIR FLAMES K. M. Saqr, M. M. Sies, and M. A. Wahid High-Speed Reacting Flow Laboratory, Faculty of Mechanical Engineering, Universiti Tenologi Malaysia, Sudai - Johor, Malaysia haledsaqr@gmail.com Received 9 April 009; accepted 9 July 009 ABSTRACT The interaction between turbulence and combustion in non-premixed flames is numerically analyzed. The reaction rate is expressed as a function of the large-eddy time scale through applying the Spalding eddy dissipation model. The turbulence energy production and dissipation rates are simulated by the turbulence model. The reacting S equations were spatially discretized and solved through a finite volume scheme and a decoupled pressure-velocity approach, respectively. The flame was assumed to be steady-state, two dimensional and axisymmetric. The reported results include the velocity, temperature and turbulent reaction rate along the flame propagation field. It is found that the elevated augmentation of free stream turbulence intensity reduces the reaction zone significantly, hence, induces the flame extinction process. Keywords: CFD, Combustion, non-premixed flames, turbulence modeling, eddy dissipation model 1 ITRODUCTIO The numerical simulation of turbulent non-premixed flames encounters a major difficulty when it comes to calculating the reaction rate. This difficulty comes, in principle, from the interaction between turbulence and reaction regions (Magnussen, 1981). In non-premixed flames, combustion occurs on the interface between fuel and oxidizer. For this reason, the mixing rate, which radically depends on turbulence, controls the chemical reaction (Warnatz, 006). In the wor of Spalding (Spalding, 1970). the reaction rate is controlled by the largeeddy mixing time. This time equals to the ratio between the local TKE and its dissipation rate. Although turbulence is thought of as a promoting factor to several applied combustion processes, free stream turbulence can be an extinction driving factor for non-premixed flames. When the oxidant is introduced to the combustion chamber with certain turbulence intensity, it achieves good mixing with the fuel. evertheless, if such turbulence intensity is exceedingly high, the flame becomes unstable, feeble, and it is susceptible to extinct. This negative effect, theoretically, comes from the competition between mass and heat diffusion.
2 70 K. M. Saqr et al. When the turbulence intensity is very high, the interface between the fuel and oxidizer become very chaotic. This in return, prevents the flame from propagating in a homogenous and stable region. At one point, the reaction is both possible and impossible at very close instants of time. The theoretical study of such interrelationship must involve a robust method to express the reciprocating effects of turbulence and heat/mass diffusion of the reaction. Magnussen et al. proposed several mechanisms to couple the effect of molecular mixing due to turbulence and chemical reaction (Magnussen 1976; Magnussen 1981) and (Magnussen et al. 1979). Their most nown and established model is the Eddy Dissipation Model (EDM). The EDM is, in its essence, a modification of the Spalding eddy-brea-up model. The foundation of such wor is to pass over the influence of chemical inetics by replacing the chemical time scale of a certain one-step reaction by the turbulent time scale ( / ). This is based on the idea that the mixing of reactants and products is the factor that determines the chemical reaction rate, in a turbulent flow (Peters 000) Although EDM was firstly devised for premixed combustion, when it comes to integrating such model in a CFD code, it can be used to model non-premixed flames with a reasonable degree of accuracy. A finite volume solver would trac the mass fraction of fuel and oxidizer in each and every cell. The existence of a nonzero value for the turbulent time scale is strictly related with the mixing process. This close relationship comes from the coupling between the convection-diffusion equation for each specie, the conservation equations, and turbulence model equations. There are several variants of the turbulence model (Cebeci 004). The most well nown and established variants are the Launder and Spalding (Launder and Spalding 197), renormalization group (RG) model (Yahot and Orszag 1986), and realizable model (Shih et al. 1995). The ey differences between the three variants exist in the method of calculating turbulence viscosity, the turbulence diffusion of and, and the generation and destruction terms in the ε equation (Jiyuan et al. 008) and (Salas et al. 1999). In this paper, we use the turbulence model, coupled with the EDM, is used in order to investigate the effect of turbulence intensity in the air stream on the combustion characteristics and flow dynamics of methane-air non-premixed flames. The mathematical treatment of the problem is firstly highlighted, followed by an explanation of the implemented numerical method. The results are discussed in the final section of this paper. MATHEMATICAL MODEL The governing equations of the flow field are the reacting conservation equations in the XY Cartesian coordinates. The following assumptions are made to simplify the mathematical treatment of the problem: 1. The reacting flow is assumed to be steady-state. The effect of gravity (i.e. buoyancy) is neglected 3. Radiation heat transfer is neglected
3 umerical Investigation of the Turbulence-Combustion Interaction Combustion is assumed to be stoichiometric and adiabatic (no heat transfer to surroundings) 5. Thermal energy created by viscous shear is neglected ( Br <<1) The governing equations, after simplification, are given below in their differential forms:.1 Conservation of Mass and Species Total mass conservation: V 0 Mass conservation for species V (1) Y () Where and are the mass diffusion velocity and reaction rate for species, respectively. The diffusive term in the equation is subjected to Y 1 0 (3) And the reaction rate is governed by 1 0 (4). Conservation of Momentum for a Reacting Flow p Y f x (5) x x y xx yx X direction: uv 1 p (6) y x y xy yy Y direction: vv Y f y.3 Conservation of Energy The following version of the energy equation was derived based on (Kuo 1986; Poinsot and Veynante 001; Peters 000): 1 h (7) 1 1 M V E p et h V where p V E h h Y 1 T T1 Cp dt (8)
4 7 K. M. Saqr et al. S V 3 1 T and S V V (9).4 Turbulence Modeling The turbulence model presented by (Spalding 1970) is used to model the turbulent inetic energy and its dissipation rate. Since the flow is treated according to the ideal gas low, only the effect of thermal compressibility is considered to influence the dissipation of turbulent inetic energy (Sarar and Balarishnan 1990). The turbulent inetic energy is obtained from the following equation: (10) t V P Dt t V C1 P C (11) where t C and Dt (1) RT Turbulence intensity is expressed as I 3 (13) u v.5 Combustion Modeling The current model adopts the EDM in order to calculate the effect of turbulent chemical reaction rate. This model taes the minimum of the following rates as the local average reaction rate in equation (Peters 000): Fuel Oxygen AY Fuel Y A Oxygen (14) (15) Products A B Y 1 Product (16) 3 UMERICAL METHOD
5 umerical Investigation of the Turbulence-Combustion Interaction 73 The computational domain represent a horizontal, axisymmetric reactive flow field, with 00 mm length and 0 mm height. Methane is injected from the left side, from a nozzle of 4 mm diameter. Air is injected in the same flow direction of methane, from a port which surrounds the methane inlet port with an net diameter of 16 mm. The flow field geometry is illustrated in figure mm 0 mm Air Methane Air Figure 1: The flow field geometry for the non-premixed flame The computational domain was spatially discretized to 8385 variable density quadrilateral grid cells as in figure. The flow field variables in each cell was assumed to have an average value at the cell center. This first order upwind scheme was used in order to calculate the pressure, velocity, and temperature through the discretization of equations (1) to (), (5) to (7) (Kari and Patanar 1989). A momentum-weighted averaging method was implemented in order to interpolate the cell-center velocity values to the face values (Rhie and Chow 1983). This technique aims to prevent the instability in the pressure calculations. b Figure : A schematic showing a section of the variable density grid. (a) high density (i.e. fine) grid (b) low density (i.e. coarse) grid The Semi-Implicit Method for Pressure-Lined Equations (SIMPLE) was used to couple the pressure and velocity iterative calculations (Patanar and Spalding 197). An implicit relaxation approach was also employed in order to stabilize nonlinear iterative solution convergence of the pressure field. a 4 RESULTS AD DISCUSSIO Six cases were solved in order to show the effects of inlet air turbulence intensity on the flame characteristics. The air to fuel ratio were fixed at 6.00 in all cases. The turbulence intensity varied from 5% to 100%. Table 1 lists the boundary conditions of the studies cases. Table 1. Boundary conditions of the six cases Case Inlet air turbulence intensity Total inlet mass flow a 5% g/s b 0% g/s c 40% g/s d 60% g/s e 80% g/s f 100% g/s
6 74 K. M. Saqr et al. Figure 3 shows the temperature pattern for the six cases. The temperature distribution indicates the reaction zone, which in turns, refers to the structure. It is clear that the first three cases are approximately identical in terms of the flame structure, holder, and temperature gradient. The reason behind such similarity is the dominance of the flame induced TKE over the free stream TKE. Thus, the reaction zone is formed and solely controlled by the turbulence generated from the interception between the fuel and air flow streams in the reaction zone. On the contrary, in the three cases corresponding to 60%, 80% and 100% free stream turbulence, the TKE generated in the flame flow field is less than such of the inlet air. Figure 4 draws a quantitative comparison of the flow field TKE and free stream TKE for the six cases. When the free stream TKE prevails over the flow field TKE, two major changes occur in the reaction zone. Firstly, the location of maximum temperature is withdrawn near to higher turbulence intensity (i.e. TKE), which is near to the air inlet port. Secondly, the interception between the inlet air stream, which contains higher TKE than the average of the flame flow field, creates a region of low turbulence intensity, as in figure 5. (a) (b) (c) (d) (e) (f) Figure 3. Flame temperature patterns of the six cases
7 umerical Investigation of the Turbulence-Combustion Interaction 75 In figure 4, it is clear that the TKE of the inlet air stream increases over the average TKE of the flame flow field starting from case (d). Referring to equation (14), region (I) in figure 4 has high turbulence intensity generated in the combustion zone than that provided in by the free stream turbulence. On the other hand, region (II) has higher turbulence intensity provided by the inlet air than that generated in the combustion zone. The effect of such disagreement is evident in the temperature pattern. Turbulence Kinetic Energy (m.s - ) Air Free Stream TKE Average Flame Induced TKE I II 0 a b c d e f Figure 4. Comparison of air free stream TKE and average TKE of the flame flow field In order to verify the effect of the free stream turbulence on the reaction zone, a third evaluation is introduced for the axial velocity and turbulent reaction rate in the six cases, as in figures 5 and 6. The turbulent reaction rate is numerically observed to be consistent in the first three cases, as well as the axial velocity profile. This steadiness comes from the dependence of the turbulent reaction rate on the flow induced eddies. In the first three cases, the axial velocity, in fact, increases due to combustion. The heat released by combustion is converted to inetic energy increasing both axial velocity and TKE. The reaction zone length, in the first three cases, has a total axial length and maximum rate of 0.15 m and g.mol/m 3.s, respectively. The flame propagation velocity reaches the maximum value of 97 m/s at the outlet. In contrast, in the last three cases, where the inlet air turbulence intensity is higher than the turbulence intensity produced by the strain rate of the flame flow field, the reaction zone is intensified near to the inlet port. Figure 5 compares the turbulent reaction rate and axial velocity profiles along the flame axis. The axial velocity, on the contrary from the previous cases, deteriorates after the reaction zone. The highest turbulent reaction rate exist in a region near to the air inlet port, where the turbulence intensity (i.e. TKE) is much higher than that developed downstream the flame propagation field. Turbulent reaction rate reaches a maximum value of one order of magnitude higher than that in the first three cases. However, the axial reaction zone length decreases drastically to 0.07 m. On the other side, the axial velocity profile taes a route that is strongly lied to the relationship between inlet TKE and the produced TKE. In this sense, the substantially high turbulence dissipation rate, ε, associated with the inlet air stream dissipates the produced TKE after the reaction zone. Case
8 76 K. M. Saqr et al. (a) (b) (c) Figure 5: Flame propagation velocity and turbulent reaction rate along the flame centerline for cases (a-c) (d) (e) (f) Figure 6: Flame propagation velocity and turbulent reaction rate along the flame centerline for cases (d-f)
9 5 COSLUSIO umerical Investigation of the Turbulence-Combustion Interaction 77 The effects of free stream turbulence on the non-premixed methane-air flame were investigated numerically. The implementation of the EDM with the -ε turbulence model revealed important information about the turbulence-combustion interaction. The elevated values of free stream turbulence were found to have negative effects on the combustion process and flame propagation. Such negative effects decrease the reaction zone size, and axial flame propagation velocity. The competition between the free stream turbulence and flow induced turbulent inetic energy, via the EDM perspective, were the theoretical foundation to justify such negative effects. Future research should verify these results by using different numerical techniques such as the probability density function and different turbulence models. Experimental investigation using laser diagnostics would reveal much more information about the turbulence-combustion interaction at these severe conditions. omenclature Br Cp E f h j Brinman number (dimensionless) Specific heat at constant pressure for species (j.g-1.-1) Stored energy per unit mass (j.g-1) Volume force acting on species in direction () Sensible enthalpy of species (j.g-1) Turbulence inetic energy (TKE) (m.s ) M Molecular weight of species (mol) P T Y Turbulence production (g.m-1.s-3) Temperature () Isentropic expansion factor (dimensionless) Mass fraction of species (dimensionless) Turbulent Prandtl number for, 1.3 (dimensionless) Reaction rate of species (s-1) Dynamic viscosity (g.m-1.s 1) Stress tensor as defined in [16] Dt Turbulent dilatation dissipation term C Constant, C Constant, 0.09 C Constant, 1.9 h e p R V 1 Enthalpy of formation for species (j.g-1) Effective conductivity (W.m-1.-1) Pressure (Pa) Universal gas constant (dimensionless) Velocity components in X and Y directions (m.s-1)
10 78 K. M. Saqr et al. t Turbulent Prandtl number for, 1.0 (dimensionless) Macroscopic density (g.m-3) Diffusion velocity of species (m.s-1) Turbulent (eddy) viscosity Differential operator, Kronecer tensorial symbol REFERECES Cebeci T (004). Turbulence models and their application, Springer, pp Jiyuan T, Guan HY, and Chaoqun L (008). Computational fluid dynamics: a practical approach, Butterworth-Heinemann. pp Kari KC and Patanar SV (1989). Pressure-Based Calculation Procedure for Viscous Flows at All Speeds in Arbitrary Configurations. AIAA Journal, 7, pp Kuo KK (1986). Principles of Combustion, John Wily and Sons. ew Yor, USA. pp Launder BE and Spalding DB (197). Lectures in Mathematical Models of Turbulence. Academic Press, London, England. Magnussen BF (1981). On the Structure of Turbulence and a Generalized Eddy Dissipation Concept for Chemical Reaction in Turbulent Flow. 19th AIAA Aerospace Science Meeting, Jan. pp Magnussen BF and Hjertager BH (1976). On mathematical models of turbulent combustion with special emphasis on soot formation and combustion. Magnussen BF, Hjertager BH, Olsen JG, and Bhaduri D (1979). Effect of Turbulent Structure and Local Concentration on Soot Formation and Combustion in CH Diffusion Flames. 17th International Symposium on Combustion (1978). Combustion Institute, Pittsburg, Pennsylvania, pp , Patanar SV and Spalding DB (197). A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows. International Journal of Heat and Mass Transfer, 15, pp Peters (000). Turbulent Combustion, Cambridge University Press. pp Poinsot T and Veynante D (001). Theoretical and umerical Combustion. R.T. Edwards Publishing. PA, USA. pp Rhie CM and Chow WL (1983). umerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation. AIAA Journal, 1, pp
11 umerical Investigation of the Turbulence-Combustion Interaction 79 Salas MD, Hefner J, and Saell L (1999). Modeling Complex Turbulent Flows, Springer. pp Sarar S and Balarishnan L (1990). Application of a Reynolds-Stress Turbulence Model to the Compressible Shear Layer. ICASE Report 90-18, ASA CR Shih TH, Liou WW, Shabbir A, Yang Z, and Zhu J (1995). A ew Eddy-Viscosity Model for High Reynolds umber Turbulent Flows - Model Development and Validation. Computers and Fluids, 4(3), pp Spalding DB (1970). Mixing and chemical reaction in steady confined turbulent flames. International 13th Symposium on Combustion. The Combustion Institute. Warnatz J, Maas M, and Dibble RW (006). Combustion: physical and chemical fundamentals, modeling and simulation, experiments, pollutant formation. Springer, pp Yahot V and Orszag SA (1986). Renormalization Group Analysis of Turbulence: I. Basic Theory. Journal of Scientific Computing, 1(1), pp
A NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER
SCIENTIFIC RESEARCH AND EDUCATION IN THE AIR FORCE-AFASES 2016 A NUMERICAL ANALYSIS OF COMBUSTION PROCESS IN AN AXISYMMETRIC COMBUSTION CHAMBER Alexandru DUMITRACHE*, Florin FRUNZULICA ** *Institute of
More informationEVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS
EVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS A Aroussi, S Kucukgokoglan, S.J.Pickering, M.Menacer School of Mechanical, Materials, Manufacturing Engineering and
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 informationPROPOSAL OF HEAT EXCHANGER IN MICRO COGENERATION UNIT, CONFIGURATION WITH BIOMASS COMBUSTION. Jozef HUŽVÁR, Patrik NEMEC
PROPOSAL OF HEAT EXCHANGER IN MICRO COGENERATION UNIT, CONFIGURATION WITH BIOMASS COMBUSTION Jozef HUŽVÁR, Patri NEMEC Authors: Worplace: Jozef Hužvár, MSc. Eng. Patri Nemec, MSc. Eng. Faculty of Mechanical
More informationENTROPY PRODUCTION MODELING IN CFD OF TURBULENT COMBUSTION FLOW
ENTROPY PRODUCTION MODELING IN CFD OF TURBULENT COMBUSTION FLOW Ivar S. Ertesvåg, Jostein Kolbu Department of Energy and Process Engineering Norwegian University of Science and Technology NO-7491 Trondheim,
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 informationNumerical Simulation of Entropy Generation in Hydrogen Enriched Swirl Stabilized Combustion
Saqr & Wahid CFD Letters Vol. 5(1) 13 www.cfdl.issres.net Vol. 5 (1) March 13 Numerical Simulation of Entropy Generation in Hydrogen Enriched Swirl Stabilized Combustion Khalid M. Saqr 1,* and Mazlan A.
More informationNumerical investigation of swirl flow inside a supersonic nozzle
Advances in Fluid Mechanics IX 131 Numerical investigation of swirl flow inside a supersonic nozzle E. Eslamian, H. Shirvani & A. Shirvani Faculty of Science and Technology, Anglia Ruskin University, UK
More informationADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER
ADVANCED DES SIMULATIONS OF OXY-GAS BURNER LOCATED INTO MODEL OF REAL MELTING CHAMBER Ing. Vojtech Betak Ph.D. Aerospace Research and Test Establishment Department of Engines Prague, Czech Republic Abstract
More informationComparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
Proceedings of 4 th ICCHMT May 17-0, 005, Paris-Cachan, FRANCE 381 Comparison of two equations closure turbulence models for the prediction of heat and mass transfer in a mechanically ventilated enclosure
More informationCOMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE
COMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE L. Velázquez-Araque 1 and J. Nožička 2 1 Division of Thermal fluids, Department of Mechanical Engineering, National University
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 informationDEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
DRAFT Proceedings of ASME IMECE: International Mechanical Engineering Conference & Exposition Chicago, Illinois Nov. 5-10, 2006 IMECE2006-14867 DEVELOPMENT OF CFD MODEL FOR A SWIRL STABILIZED SPRAY COMBUSTOR
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 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 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 informationJJMIE Jordan Journal of Mechanical and Industrial Engineering
JJMIE Jordan Journal of Mechanical and Industrial Engineering Volume Number, June.6 ISSN 995-6665 Pages 99-4 Computational Fluid Dynamics of Plate Fin and Circular Pin Fin Heat Sins Mohammad Saraireh *
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 informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationCFD modeling of combustion
2018-10 CFD modeling of combustion Rixin Yu rixin.yu@energy.lth.se 1 Lecture 8 CFD modeling of combustion 8.a Basic combustion concepts 8.b Governing equations for reacting flow Reference books An introduction
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationNUMERICAL ANALYSIS OF TURBULENT FLAME IN AN ENCLOSED CHAMBER
NUMERICAL ANALYSIS OF TURBULENT FLAME IN AN ENCLOSED CHAMBER Naveen Kumar D 1*, Pradeep R 2 and Bhaktavatsala H R 3 1 Assistant Professor Department of Mechanical Engineering, M S Engineering College,
More informationA Comparative Analysis of Turbulent Pipe Flow Using k And k Models
A Comparative Analysis of Turbulent Pipe Flow Using k And k Models 1 B. K. Menge, 2 M. Kinyanjui, 3 J. K. Sigey 1 Department of Mathematics and Physics. Technical University of Mombasa, P.O BOX 90420-80100,Mombasa,
More informationImpact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions
Journal of Physics: Conference Series PAPER OPEN ACCESS Impact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions To cite this article: A Rosiak and
More informationCFD and Kinetic Analysis of Bluff Body Stabilized Flame
CFD and Kinetic Analysis of Bluff Body Stabilized ame A. Dicorato, E. Covelli, A. Frassoldati, T. Faravelli, E. Ranzi Dipartimento di Chimica, Materiali e Ingegneria Chimica, Politecnico di Milano, ITALY
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 informationVERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION
VERTICAL TURBULENT BUOYANT HELIUM JET CFD MODELING AND VALIDATION Cheng Z, Agranat V.M. and Tchouvelev A.V. A.V.Tchouvelev & Associates, Inc., 659 Spinnaker Circle, Mississauga, Ontario, Canada L5W R Hydrogenics
More informationWall treatments and wall functions
Wall treatments and wall functions A wall treatment is the set of near-wall modelling assumptions for each turbulence model. Three types of wall treatment are provided in FLUENT, although all three might
More informationCFD study of gas mixing efficiency and comparisons with experimental data
17 th European Symposium on Computer Aided Process Engineering ESCAPE17 V. Plesu and P.S. Agachi (Editors) 2007 Elsevier B.V. All rights reserved. 1 CFD study of gas mixing efficiency and comparisons with
More informationMUSCLES. Presented by: Frank Wetze University of Karlsruhe (TH) - EBI / VB month review, 21 September 2004, Karlsruhe
MUSCLES Modelling of UnSteady Combustion in Low Emission Systems G4RD-CT-2002-00644 R&T project within the 5 th Framework program of the European Union: 1 Numerical computations of reacting flow field
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More informationMathematical Investigation and CFD Simulation of Monolith Reactors: Catalytic Combustion of Methane
Excerpt from the Proceedings of the COMSOL Conference 8 Hannover Mathematical Investigation and CFD Simulation of Monolith Reactors: Catalytic Combustion of Methane Maryam Ghadrdan *,, Hamid Mehdizadeh
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 informationNO x Emissions Prediction from a Turbulent Diffusion Flame
NO x Emissions Prediction from a Turbulent Diffusion Flame Akpan, Patrick 1 and Njiofor, Tobenna 2 1 Department of Mechanical Engineering, University of Nigeria, Nsukka 2 Process Systems Engineering, Cranfield
More informationARRHENIUS LAW MODIFICATION FOR TURBULENT COMBUSTION MODELING
ARRHENIUS LAW MODIFICATION FOR TURBULENT COMBUSTION MODELING M. Javadi, M. Moghiman, A. Zamani Department of mechanical engineering, Ferdowsi University Mashhad, Iran, P.O.Box: 91775-1111 Mohammad.Javadi@gmail.com
More informationSoot formation in turbulent non premixed flames
Soot formation in turbulent non premixed flames A. Cuoci 1, A. Frassoldati 1, D. Patriarca 1, T. Faravelli 1, E. Ranzi 1, H. Bockhorn 2 1 Dipartimento di Chimica, Materiali e Ing. Chimica, Politecnico
More informationCalculating equation coefficients
Fluid flow Calculating equation coefficients Construction Conservation Equation Surface Conservation Equation Fluid Conservation Equation needs flow estimation needs radiation and convection estimation
More informationLecture 12. Droplet Combustion Spray Modeling. Moshe Matalon
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
More informationMeysam ATASHAFROOZ, Seyyed Abdolreza GANDJALIKHAN NASSAB, and Amir Babak ANSARI
THERMAL SCIENCE: Year 014, Vol. 18, No., pp. 479-49 479 NUMERICAL INVESTIGATION OF ENTROPY GENERATION IN LAMINAR FORCED CONVECTION FLOW OVER INCLINED BACKWARD AND FORWARD FACING STEPS IN A DUCT UNDER BLEEDING
More informationCFD Analysis for Thermal Behavior of Turbulent Channel Flow of Different Geometry of Bottom Plate
International Journal Of Engineering Research And Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 13, Issue 9 (September 2017), PP.12-19 CFD Analysis for Thermal Behavior of Turbulent
More informationTransactions on Engineering Sciences vol 18, 1998 WIT Press, ISSN
Simulation of natural convection in a reservoir P. Jelmek*, V. Havlik\ R. Cerny\ P. Pfikryl" * Czech Technical University, Faculty of Civil Engineering, Department of Physics, Thdkurova 7, 166 29 Prague
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
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 informationAdvanced near-wall heat transfer modeling for in-cylinder flows
International Multidimensional Engine Modeling User s Group Meeting at the SAE Congress April 20, 2015 Detroit, MI S. Šarić, B. Basara AVL List GmbH Advanced near-wall heat transfer modeling for in-cylinder
More informationHydrodynamic Characteristics of Gradually Expanded Channel Flow
International Journal of Hydraulic Engineering 2013, 2(2): 21-27 DOI: 10.5923/j.ijhe.20130202.01 Hydrodynamic Characteristics of Gradually Expanded Channel Flow Edward Ching-Rue y LUO Department of Civil
More informationObjectives. Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation
Objectives Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation Conservation of Mass Conservation of Mass Mass, like energy, is a conserved
More informationComputational fluid dynamics study of flow depth in an open Venturi channel for Newtonian fluid
Computational fluid dynamics study of flow depth in an open Venturi channel for Newtonian fluid Prasanna Welahettige 1, Bernt Lie 1, Knut Vaagsaether 1 1 Department of Process, Energy and Environmental
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More informationCHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH
82 CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH The coefficient of lift, drag and power for wind turbine rotor is optimized using an iterative approach. The coefficient
More informationNumerical investigation of dual-stage high velocity oxy-fuel (HVOF) thermal spray process: A study on nozzle geometrical parameters
Numerical investigation of dual-stage high velocity oxy-fuel (HVOF) thermal spray process: A study on nozzle geometrical parameters Jafari, H., Emami, S., & Mahmoudi, Y. (2017). Numerical investigation
More informationNUMERICAL AND EXPERIMENTAL INVESTIGATIONS OF AIR FLOW AND TEMPERATURE PATTERNS OF A LOW VELOCITY DIFFUSER
NUMERICAL AND EXPERIMENTAL INVESTIGATIONS OF AIR FLOW AND TEMPERATURE PATTERNS OF A LOW VELOCITY DIFFUSER M Cehlin and B Moshfegh Division of Energy and Mechanical Engineering, Department of Technology,
More informationS. Kadowaki, S.H. Kim AND H. Pitsch. 1. Motivation and objectives
Center for Turbulence Research Annual Research Briefs 2005 325 The dynamics of premixed flames propagating in non-uniform velocity fields: Assessment of the significance of intrinsic instabilities in turbulent
More informationA copy can be downloaded for personal non-commercial research or study, without prior permission or charge
Oni, T. O., and Paul, M. C. (2014) Numerical simulation of turbulent heat transfer and fluid flow in different tube designs. In: World Congress on Engineering 2014, 2-4 July 2014, London, UK. Copyright
More informationNUMERICAL ANALYSIS OF THE THREE-MATERIAL DOWNHOLE FLOW FIELD IN HYDROTHERMAL JET DRILLING
2017 WJTA-IMCA Conference and Expo October 25-27, 2017 New Orleans, Louisiana Paper NUMERICAL ANALYSIS OF THE THREE-MATERIAL DOWNHOLE FLOW FIELD IN HYDROTHERMAL JET DRILLING Xianzhi Song, Zehao Lyu, Haizhu
More informationCombustion and Mixing Analysis of a Scramjet Combustor Using CFD Pradeep Halder 1 Edla Franklin 2 Dr. P. Ravinder Reddy 3
IJSRD - International Journal for Scientific Research & Development Vol. 2, Issue 09, 2014 ISSN (online): 2321-0613 Combustion and Mixing Analysis of a Scramjet Combustor Using CFD Pradeep Halder 1 Edla
More informationXXXVIII Meeting of the Italian Section of the Combustion Institute
Coupling a Helmholtz solver with a Distributed Flame Transfer Function (DFTF) to study combustion instability of a longitudinal combustor equipped with a full-scale burner D. Laera*, S.M. Camporeale* davide.laera@poliba.it
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationNatural Convection in Parabolic Enclosure Heated from Below
www.ccsenet.org/mas Modern Applied Science Vol. 5, No. 3; June 011 Natural Convection in Parabolic Enclosure Heated from Below Dr. Ahmed W. Mustafa (Corresponding auther) University of Tikrit, College
More informationLarge-eddy simulations of turbulent reacting stagnation point flows
Copyright 1997, American Institute of Aeronautics and Astronautics, Inc. AIAA Meeting Papers on Disc, January 1997 A9715437, AIAA Paper 97-0372 Large-eddy simulations of turbulent reacting stagnation point
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 informationFlow Structure Investigations in a "Tornado" Combustor
Flow Structure Investigations in a "Tornado" Combustor Igor Matveev Applied Plasma Technologies, Falls Church, Virginia, 46 Serhiy Serbin National University of Shipbuilding, Mikolayiv, Ukraine, 545 Thomas
More informationCoupled Fluid and Heat Flow Analysis Around NACA Aerofoil Profiles at Various Mach Numbers
International Journal of Engineering Research and Technology. ISSN 0974-3154 Volume 6, Number 2 (2013), pp. 249-258 International Research Publication House http://www.irphouse.com Coupled Fluid and Heat
More informationOn the transient modelling of impinging jets heat transfer. A practical approach
Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,
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 informationTHERMAL ANALYSIS OF SECOND STAGE GAS TURBINE ROTOR BLADE
Polymers Research Journal ISSN: 195-50 Volume 6, Number 01 Nova Science Publishers, Inc. THERMAL ANALYSIS OF SECOND STAGE GAS TURBINE ROTOR BLADE E. Poursaeidi, M. Mohammadi and S. S. Khamesi University
More informationNUMERICAL SIMULATION OF LDI COMBUSTOR WITH DISCRETE-JET SWIRLERS USING RE-STRESS MODEL IN THE KIVA CODE
NUMERICAL SIMULATION OF LDI COMBUSTOR WITH DISCRETE-JET SWIRLERS USING RE-STRESS MODEL IN THE KIVA CODE S. L. Yang, C. Y. Teo, and Y. K. Siow Department of Mechanical Engineering Engineering Mechanics
More informationAdvanced Turbulence Models for Emission Modeling in Gas Combustion
1 Advanced Turbulence Models for Emission Modeling in Gas Combustion Ville Tossavainen, Satu Palonen & Antti Oksanen Tampere University of Technology Funding: Tekes, Metso Power Oy, Andritz Oy, Vattenfall
More informationExploring STAR-CCM+ Capabilities, Enhancements and Practices for Aerospace Combustion. Niveditha Krishnamoorthy CD-adapco
Exploring STAR-CCM+ Capabilities, Enhancements and Practices for Aerospace Combustion Niveditha Krishnamoorthy CD-adapco Outline Overview of modeling capability Applications, Practices and Enhancements
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 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 informationInternational Journal of Research in Advent Technology, Vol.6, No.11, November 2018 E-ISSN: Available online at
Comparative analysis of cylindrical and helical coil counter flow type of heat exchanger used in thermoelectric generator for waste heat recovery using CFD fluent Chanchal Kumar 1, a, Dr. Savita Vyas 2,b
More informationSimplified kinetic schemes for oxy-fuel combustion
Simplified kinetic schemes for oxy-fuel combustion 1 A. Frassoldati, 1 A. Cuoci, 1 T. Faravelli, 1 E. Ranzi C. Candusso, D. Tolazzi 1 Dipartimento di Chimica, Materiali e Ingegneria Chimica, Politecnico
More informationAlgebraic Turbulence-Chemistry Interaction Model
50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 09-12 January 2012, Nashville, Tennessee AIAA 2012-0183 50th AIAA Aerospace Sciences Meeting, January 9-12,
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
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 informationNUMERICAL ANALYSIS OF THE EFFECTS OF STREAMLINING GEOMETRY AND A VECTOR WALL ON THE THERMAL AND FLUID FLOW IN A SRU THERMAL REACTOR.
NUMERICAL ANALYSIS OF THE EFFECTS OF STREAMLINING GEOMETRY AND A VECTOR WALL ON THE THERMAL AND FLUID FLOW IN A SRU THERMAL REACTOR Chun-Lang Yeh Department of Aeronautical Engineering, National Formosa
More informationEXPERIMENTAL SETUP The experimental test-room was representative of a full-sized office space. This was a natural
Seventh International IBPSA Conference Rio de Janeiro, Brazil August 13-15, 2001 ON THE NATURAL DISPLACEMENT VENTILATION FLOW THROUGH A FULL SCALE ENCLOSURE, DRIVEN BY A SOURCE OF BUOYANCY AT FLOOR LEVEL
More informationSimulation of Turbulent Lifted Flames and their Transient Propagation
25 th ICDERS August 2-7th, 2015 Leeds, UK Simulation of Turbulent Lifted Flames and their Transient Propagation S. Ruan, Z. Chen, N. Swaminathan University of Cambridge Cambridge, UK 1 Introduction Turbulent
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationEffect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process
, pp. 844 850 Effect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process Baokuan LI and Fumitaka TSUKIHASHI 1) Department of Thermal Engineering, The School
More informationANSYS Advanced Solutions for Gas Turbine Combustion. Gilles Eggenspieler 2011 ANSYS, Inc.
ANSYS Advanced Solutions for Gas Turbine Combustion Gilles Eggenspieler ANSYS, Inc. 1 Agenda Steady State: New and Existing Capabilities Reduced Order Combustion Models Finite-Rate Chemistry Models Chemistry
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationApplication of Computational Fluid Dynamics for Different Fire Strengths in a Compartment Using Combustion Modelling
Fire Science and Technology Vol.33 No. (4) 35-46 35 Research Paper Application of Computational Fluid Dynamics for Different Fire Strengths in a Compartment Using Combustion Modelling Aravind Kumar. A,
More informationNumerical investigation on the effect of inlet conditions on the oblique shock system in a high-speed wind tunnel
American Journal of Physics and Applications 2013; 1(3): 91-98 Published online November 30, 2013 (http:// www.sciencepublishinggroup.com/j/ajpa) doi: 10.11648/j.ajpa.20130103.16 Numerical investigation
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 informationComputational Fluid Dynamics
Computational Fluid Dynamics A Practical Approach Jiyuan Tu RMIT University, Australia Guan Heng Yeoh Australian Nuclear Science and Technology Organisation Chaoqun Liu University of Texas, Arlington ^fl
More informationTutorial: Premixed Flow in a Conical Chamber using the Finite-Rate Chemistry Model
Tutorial: Premixed Flow in a Conical Chamber using the Finite-Rate Chemistry Model Introduction The purpose of this tutorial is to provide guidelines and recommendations for setting up and solving the
More informationIntroduction Flares: safe burning of waste hydrocarbons Oilfields, refinery, LNG Pollutants: NO x, CO 2, CO, unburned hydrocarbons, greenhouse gases G
School of Process, Environmental and Materials Engineering Computational study of combustion in flares: structure and emission of a jet flame in turbulent cross-flow GERG Academic Network Event Brussels
More informationTurbulent Natural Convection in an Enclosure with Colliding Boundary Layers
Turbulent Natural Convection in an Enclosure with Colliding Boundary Layers Abstract Mutuguta John Wanau 1* 1. School of Pure and Applied Sciences, Murang a University of Technology, P.O box 75-10200,
More informationPredicting NO Formation with Flamelet Generated Manifolds
Predicting NO Formation with Flamelet Generated Manifolds J. A. van Oijen and L. P. H. de Goey Dept. Mechanical Engineering, Technische Universiteit Eindhoven P.O. Box, 6 MB Eindhoven, The Netherlands
More informationINTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---
ILASS-Europe 2002 Zaragoza 9 11 September 2002 INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING--- Z. Tapia, A. Chávez e-mail: ztapia@imp.mx Instituto Mexicano del Petróleo Blvd. Adolfo Ruiz Cortines
More informationNumerical solution of the flow in the gas transfer system of an industrial unit
Technical Journal of Engineering and Applied Sciences Available online at www.tjeas.com 213 TJEAS Journal-213-3-24/3688-3694 ISSN 251-853 213 TJEAS Numerical solution of the flow in the gas transfer system
More informationDetailed Outline, M E 320 Fluid Flow, Spring Semester 2015
Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous
More informationTutorial for the supercritical pressure pipe with STAR-CCM+
Tutorial for the supercritical pressure pipe with STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities with
More informationJoule Heating Effect on the Coupling of Conduction with Magnetohydrodynamic Free Convection Flow from a Vertical Flat Plate
Nonlinear Analysis: Modelling and Control, 27, Vol. 12, No. 3, 37 316 Joule Heating Effect on the Coupling of Conduction with Magnetohydrodynamic Free Convection Flow from a Vertical Flat Plate M. A. Alim
More informationNumerical studies on natural ventilation flow in an enclosure with both buoyancy and wind effects
Numerical studies on natural ventilation flow in an enclosure with both buoyancy and wind effects Ji, Y Title Authors Type URL Numerical studies on natural ventilation flow in an enclosure with both buoyancy
More informationA Unified 3D CFD Model for Jet and Pool Fires
A Unified 3D CFD Model for Jet and Pool Fires Chenthil Kumar K. 1, Anil Kumar K. R. 1 and Amita Tripathi 2 1 Fluidyn Consultancy (P) Ltd, 146, Ring Road, 5, HSR Layout, Bangalore - 560102, India 2 Fluidyn
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 informationActive Control of Separated Cascade Flow
Chapter 5 Active Control of Separated Cascade Flow In this chapter, the possibility of active control using a synthetic jet applied to an unconventional axial stator-rotor arrangement is investigated.
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More information