Construction of Libraries for Non-Premixed Tabulated Chemistry Combustion Models including Non-Adiabatic Behaviour due to Wall Heat Losses
|
|
- Coral Armstrong
- 5 years ago
- Views:
Transcription
1 Sonderforschungsbereich/Transregio 40 Annual Report Construction of Libraries for Non-Premixed Tabulated Chemistry Combustion Models including Non-Adiabatic Behaviour due to Wall Heat Losses By G. Frank, J. Zips AND M. Pfitzner Institut für Thermodynamik, Universität der Bundeswehr München Werner-Heisenberg-Weg 39, Neubiberg Non- Flamelet Libraries in the literature are often associated with radiative heat losses, which is oftentimes surpassed by the heat losses due to cooled walls. Different methods of generating non- flamelets are compared with special emphasis on their applicability for flamelets subjected to wall heat losses. 1. Introduction Simulating combustion processes is still a very challenging task for Computational Fluid Dynamics (CFD). A combustion process is described by hundreds of elementary reactions, forming a stiff complex set of partial differential equations to be solved numerically. To get a result in a reasonable time, models simplifying the computational complexity are used for combustion processes. Tabulated chemistry models decouple the time-consuming solution of the complex chemistry from the main fluid flow simulation. For this the thermodynamic state (temperature T and composition Y i ) is parametrized by a set of certain parameters and saved in a table. Transport equations for those parameters are solved during runtime in the flow simulation and the thermodynamic state is then retrieved via table look-up. One of the most often used tabulated chemistry models is the flamelet concept, based on the theory by Peters [1]. Here, the turbulent flame is assumed to consist of a collection of laminar flames, called flamelets. Each thermodynamic state is defined by the mixture fraction Z and the scalar dissipation rate χ. Turbulence-chemistry interaction is oftentimes modelled using assumed probability density functions (PDFs). The thermodynamic state then additionally depends on the variance of the mixture fraction Z 2. Most tabulated chemistry models for non-premixed combustion are based on laminar flamelets, which are generated by solving a one-dimensional laminar counterflow diffusion flame (see Fig.1). Different numerical tools exist with the capability of using detailed chemical reaction mechanisms. In this work the FlameMaster code by Pitsch [2] was used, which solves the flamelet equations in mixture fraction space. For this, a coordinate transformation of the species- and energy conservation equations is done from physical space into mixture fraction space, leading to following equation: ρ Φ i t = ρ χ 2 2 Φ i Z 2 +ω i (1.1)
2 194 G. Frank, J. Zips & M. Pfitzner S-Curve 2000 T max [K] χ st [1/s] FIGURE 1. Schematic counterflow diffusion flame configuration. FIGURE 2. S-Curve of a counterflow diffusion flame. χ = 2D Z 2, (1.2) with Φ i being a reactive scalar (Y i, T ) and ω i being the source term of the scalar. For those equations the Lewis number of each species is assumed to be unity. For details of the calculation of the diffusion coefficient D and the source terms for species mass fractions and temperature the reader is referred to [2]. For steady flamelets, the unsteady term on the left hand side of equation (1.1) vanishes. Those equations can have three solutions, which are shown by the classical S-Curve Fig. 2: An unburnt solution, where no chemical reactions occur (lower branch of the curve), a fully reacted solution (top branch of the curve), and an unstable burning solution (middle branch of the curve). The classic flamelet model by Peters [1] only uses the top branch of the S-Curve with the parameters Z, χ. A Flamelet-Progress Variable (FPV) combustion model can also represent the two other branches of the S-Curve. Combustion processes near cooled walls are subjected to heat loss. Thus the reaction rates of the elementary reactions drop, slow reaction with lower activation energy become more prominent and recombination processes can occur, consequently changing the thermodynamic states compared to the flame. Including heat loss effects into tabulated chemistry models has been studied by many researchers (e.g. [3, 4]) in the scope of radiation heat losses. As additional parameter for the tabulation of the thermodynamic states the enthalpy of the mixture, or a derivative of it, is used. A Direct Numerical Simulation (DNS) by Wang et al. [5] investigating flame-wall interaction validates the use of the enthalpy for wall heat losses. Another near wall effect noticed by [6] and [7] is that burned products can be trapped between the wall and the flame, affecting diffusion and mixing in the near wall region. In the scope of premixed flames, Fiorina et al. [8] extended the intrinsic low-dimensional manifold (ILDM) by an enthalpy coordinate. The ILDM library was built using isothermal burner-stabilized flames. Their model has been verified by two dimensional CFD calculations of premixed and partially premixed burner stabilized flames. In the scope of this work, we will concentrate on the construction of flamelet libraries for non-premixed flames subjected to wall heat losses. First, we will discuss a chemistry table built with the assumption of chemical equilibrium and with the chemistry code Cantera [9]. Afterwards, methods to include non- effects into one dimensional counterflow diffusion flames are presented. The method of frozen chemistry is often found in commercial solvers. Flamelets with a convective sink term as proposed by Lee et al. [10]
3 1 Non-Adiabatic Libraries for Tabulated Chemistry Models hnorm [-] T/T max [-]: FIGURE 3. Equilibrium chemistry look-up table for temperature at a specific pressure parametrized by mixture fraction Z and normalized enthalpy h norm. The temperature ranges from T min = 270 K to T max = 3500 K. are discussed as the third option. The last two methods both change the boundary conditions of the flamelets and represent the methods proposed by Hossain et al. [3] for radiative heat losses and by Wu and Ihme [11] for wall heat losses. Those five methods are of increasing complexity, their advantages, disadvantages and their applicability are discussed. 2. Chemical Equilibrium A very quick and easy way to generate a table is to assume chemical equilibrium. Chemistry is assumed to be infinitely fast, i.e.τ chem. Consequently, the equilibrium state can be calculated by imposing an enthalpy level and a constant operating pressure. In order to compute the thermochemical equilibrium, it is common to use the Gibb s free energy given by G = H TS, (2.1) with H being the enthalpy, T the absolute temperature and S the entropy, respectively. The resulting database can be parametrized in mixture fraction and enthalpy space for a specific pressure level. Turbulence-chemistry interaction can be modeled using a presumed PDF approach which increases the look-up table by one dimension. Fig. 3 shows a generic temperature table in mixture fraction and normalized enthalpy space. The equilibrium approach has two main advantages: first, the table generation is fast and simple to implement. Second, the table generation allows to represent all possible thermochemical states and no interpolation during the look-up process is necessary. However, the assumption of chemical equilibrium is very strict and mostly not fulfilled in a real flame. Obviously the effect of strain on a flame is neglected and close to walls, heat losses cause non-equilibrium effects that cannot be reproduced.
4 196 G. Frank, J. Zips & M. Pfitzner χ st = 0.3 s T [K] FIGURE 4. S-Curve for the frozen chemistry. FIGURE 5. Temperature profiles for frozen chemistry. 3. Counterflow Diffusion Flames Most tabulated chemistry models for non-premixed combustion are based on counterflow diffusion flames. In the following section four methods to include heat loss effects in counterflow diffusion flames are discussed. The reference is a counter flow diffusion flame with methane as fuel and air as oxidizer. The boundary temperatures were chosen as T = 300K. Extinction scalar dissipation rate for this flame is χ 0.3s 1. The detailed chemistry mechanism GRI-3.0 [12] was used Frozen Chemistry The frozen chemistry approach can be found often in commercial solvers (e.g. FLU- ENT [13]). Here, the flamelets are calculated solving steady, counter flow diffusion flames. The table is parameterized in terms of mixture fraction, mixture fraction variance and scalar dissipation rate, for which transport equations are solved during runtime. The composition of the mixture is read from the tabulated flamelets. An additional transport equation for enthalpy is solved to account for heat losses and the temperature is then calculated from the composition and enthalpy. Physically, it can be interpreted as a cooling down of burned products while omitting any recombination effects that may occur. The range of tabulated states for this model is shown in Fig. 4 as grey area. All states below the stable branch of the S-Curve can be calculated. Typical temperature profiles are plotted in Fig. 5. The main disadvantage of this method is, that the influence of the heat loss on the chemistry and the reaction rates is omitted and the expected influence of the radical formation near cooled walls cannot be mapped in the table Mechanisms of Heat Loss Flame-wall interaction can be characterised by the angle between flame front and wall, with the two extrema being either the flame front perpendicular to the wall (head on configuration Fig. 6) or the flame front parallel to the wall (side wall configuration Fig. 7). The heat loss mechanisms for both cases will be investigated in the following sections Flamelets with convective sink term When we assume the very thin flame front being perpendicular to the wall (Fig. 6), then the heat loss vector is perpendicular to the mixture fraction gradient of the flame as well. A one dimensional flamelet at a certain distance from the wall is therefore subjected
5 Non-Adiabatic Libraries for Tabulated Chemistry Models 197 FIGURE 6. Head-on configuration. FIGURE 7. Side wall configuration χ st = 0.3 s T [K] FIGURE 8. S-Curve for flamelets with a convective sink term FIGURE 9. Temperature profiles for flamelets with a convective sink term. to a heat loss of q = α C (T wall T) (3.1) with α C being the convective heat transfer coefficient. Please note, that a configuration as shown in Fig.6 would in reality lead to a triple flame [15]. For this problem it should be seen as simplified thought model. The convective heat transfer coefficient is antiproportional to the convective timescale α C 1/τ conv. Lee et al. [10] propose using a Nusselt correlation to estimate this heat transfer coefficient. With this heat sink term in the unsteady flamelet equations (1.1) the physical temperature dependence of the reaction rates of the elementary reactions is considered. In Fig. 8 a comparison of the S-curve and the S-curve with a convective heat loss is shown. It can be seen that at high scalar dissipation rates the heat loss source terms have little effect on the flame. Since the scalar dissipation rate is antiproportional to the flamelet lifetime, the heat loss is a process on a slow timescale that does little effect on short living flames. Additionally, the flame thickness reduces with higher dissipation rate, and thinner flames can exist nearer to the wall, which was found in [5]. The range of tabulated states with this method is indicated by the grey area in Fig. 8. The area between the stable branch with heat loss and the unburned state needs to be interpolated in the table. For a low scalar dissipation rate of of χ st = 0.3s 1, where the flame extinguishes after a while, the temperature curves are shown in Fig. 9. The highest enthalpy losses
6 198 G. Frank, J. Zips & M. Pfitzner T max [K] χ st [1/s] FIGURE 10. S-Curve for flamelets with adjusted boundaries. T [K] χ st = 0.3 s FIGURE 11. Temperature profiles for flamelets with adjusted boundaries. ( h = c p T ) occur at stoichiometric mixture fraction, where the temperature difference to the wall is the highest. After extinction of the flame, the hot combustion products are cooled down to the boundary temperature. In the instationary simulation also the slow recombination of radicals at low temperatures can be tabulated. A drawback of this method is the huge dependency on the unknown convective time scale τ conv. This time scale firstly determines how much of the area below the S-Curve can be tabulated with actual flamelets, and how much data needs to be interpolated. Moreover, this time scale also influences the production rates of species Flamelets with adjusted boundaries When the flame front is parallel to the wall, the heat loss vector is parallel to the mixture fraction gradient as well. Here heat conduction will occur through the flamelet, a process in physical space, which can be hardly transformed into mixture fraction space, since the flame thickness would be an influencing factor. If a diffusion flame brushes on a surface as shown in Fig. 7, the boundary conditions change. In Fig. 7 the rich fuel side is trapped between flame and wall, it could be also possible that the lean side is between flame and wall. For flamelets with radiation heat losses Hossain et al. [3] propose to feed slightly lean fuel and slightly rich air into the counterflow diffusion flame. This reduces the mixture fraction range on both sides as can be seen in Fig. 11. The maximum flame temperature of the reference flamelet can be reduced by around 300 K, which is sufficient for radiation losses, but for wall heat losses a much higher enthalpy losses are expected. This method seems unsuitable for heat loss effects due to cooled walls Flamelets with permeable wall A fairly new model to incorporate wall heat losses has been proposed by Wu and Ihme [11]. In the counterflow diffusion flame an isothermic permeable wall is introduced that subjects the flame to heat losses (see Fig. 12). In mixture fraction space this translates to an additional thermal boundary condition at a certain mixture fraction Z wall. For Z > Z wall the species simply diffuse to the fuel side. Preliminary results for the resulting temperature profiles are shown in Fig. 13. This method is quite promising, because it seems to allow to tabulate all states below the stable branch in the S-Curve while incorporating heat loss effects on the production rates of the species.
7 Non-Adiabatic Libraries for Tabulated Chemistry Models T [K] FIGURE 12. Schematic counterflow diffusion flame with permeable wall. FIGURE 13. Temperature profiles for flamelets with permeable wall. 4. Conclusions Heat loss effects in tabulated chemistry models have been investigated mainly in scope of thermal radiation heat losses. But heat loss also occurs near cooled walls, where it reduces the reaction speed and possibly quenches a flame. When a flame can be represented by its mixture fraction Z and scalar dissipation rate χ, its composition Y i and temperature T changes according to the enthalpy of the system subjected to wall heat losses. It would be desirable, if models based on tabulated chemistry can reproduce these changes. This work presents a short comparison of five methods to build chemistry tables with heat losses. The first section presented a tabulation method based on chemical equilibrium. This method is very quick and does not need any interpolation in the table. It represents the effect of the reduced enthalpy level on chemical kinetics. However, effects due to strain of the flames are not considered. Additionally, reactions near the wall are probably not in chemical equilibrium. In the next section a simple extention of the original flamelet model, here dubbed frozen chemistry model, as used in many commercial solvers is presented. In this method the composition of an counterflow diffusion flame is assumed to be constant and the temperature is calculated according to the total enthalpy in the system. The flamelet table is easily generated and any state below the stable branch of the S- Curve can be represented. The effects of strain are incorporated. But the effects of the heat losses on the production rates of species is omitted, therefore the composition near the wall is predicted incorrectly. A flamelet with a convective sink term is used to generate the table. Here again strain effects are incorporated in the model. Additionally the effect of heat loss on the chemistry is incorporated for flamelets at low strain. Even extinction can be tabulated. But the table cannot be filled completely with this method and states between the stable burning branch with heat loss and the unburned solution needs to be interpolated. Also species composition and production rates are dependent on an unknown timescale τ conv. For radiation heat losses an adjustment to the boundaries of the counterflow diffusion flame can be used to generate different enthalpy levels in the flame. But the enthalpy losses that the flames are subjected to, are not high enough to sufficiently model heat loss near cooled walls. This model is not suited to include wall heat loss effects. A fairly new model introducing a cooled permeable wall into the counterflow diffusion flame have been proposed by Wu and Ihme [11]. Only preliminary results are available until now, but this models is quite promising. The effort for the table generation is slightly
8 200 G. Frank, J. Zips & M. Pfitzner larger, than for the frozen chemistry model, however it should be possible to tabulate all states below the stable burning branch of the S-Curve while still incoporating the effects of heat losses on the composition. Ongoing work in the project will use four of the presented methods on a generic test case of a reacting flow over a wall and evaluate their capability of predicting the correct heat flux to the wall. Acknowledgments Financial support has been provided by the German Research Foundation (Deutsche Forschungsgemeinschaft DFG) in the framework of the Sonderforschungsbereich Transregio 40. References [1] BRAY, K. AND PETERS, N. (1994) Laminar Flamelets in Turbulent Flames. Turbulent Reacting Flows, Academic Press, London [2] PITSCH, H. (1993). Entwicklung eines Programmpaketes zur Berechnung eindimensionaler Flammen am Beispiel einer Gegenstromdiffusionsflamme. Diploma Thesis, RWTH Aachen [3] HOSSAIN, M., JONES, J. AND MALALASEKERA, W. (2001) Modelling of a Bluff- Body Nonpremixed Flame Using a Coupled Radiation/Flamelet Combustion Model Flow. Turbulence and Combustion 67, [4] COELHO, P.J., TEERLING, O.J. AND ROEKARTS D. (2003) Spectral radiative effects and turbulence/radiation interaction in a non-luminous turbulent jet diffusion flame. Combustion and Flame 133, [5] WANG, Y. AND TROUVÉ, A. (2006) Direct numerical simulation of nonpremixed flame-wall interactions. Combustion and Flame 144, [6] DABIREAU, F., CUENOT, B.. VERMOREL, O. AND POINSOT, T. (2003) Interaction of flames of H 2 + O 2 with inert walls. Combustion and Flame 135, [7] FRANK, G., FERRARO, F. AND PFITZNER, M. (2013). RANS Simulations of Chemical Reactions in Cooling Films. 5 th European Conference for Aeronautics and Space Sciences (EUCASS), Munich [8] FIORINA, B. BARON, R., GICQUEL, O., THEVENIN, D., CARPENTIER, S. AND DARABIHA, N. (2003) Modelling non- partially premixed flames using a flame-prolongation of ILDM. Combustion Theory and Modelling 7(3), [9] CANTERA DEVELOPERS. URL: [10] LEE, D.J., THAKUR, S., WRIGHT, J., IHME, M. AND SHYY, W (2001) Characterization of Flow Field Structure and Species Composition in a Shear Coaxial Rocket GH2/GO2 Injector: Modelling of Wall Heat Losses. 47th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit San Diego, California, paper AIAA [11] WU, H. AND IHME, M. (2015) Modeling of Wall Heat Transfer and Flame/Wall Interaction A Flamelet Model with Heat-Loss Effects. 9th U.S. National Combustion Meeting Cincinnati, Ohio, May [12] SMITH, G. ET AL.. URL: [13] ANSYS, INC. (2012). ANSYS FLUENT Theory Guide. Version 14.5.
9 Non-Adiabatic Libraries for Tabulated Chemistry Models 201 [14] MÜLLER H., FERRARO F. AND PFITZNER M. (2013). Implementation of a Steady Laminar Flamelet Model for Non-Premixed Combustion in LES and RANS simulations. 8th International OpenFOAM Workshop, Jeju, Korea [15] WICHMAN, I. S. (1989) On the Quenching of a Diffusion Flame Near a Cold Wall. Combustion Science and Technology 64(4-6),
A comparison between two different Flamelet reduced order manifolds for non-premixed turbulent flames
8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013 A comparison between two different Flamelet
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 informationSteady Laminar Flamelet Modeling for turbulent non-premixed Combustion in LES and RANS Simulations
Sonderforschungsbereich/Transregio 4 Annual Report 213 139 Steady Laminar Flamelet Modeling for turbulent non-premixed Combustion in LES and RANS Simulations By H. Müller, C. A. Niedermeier, M. Pfitzner
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 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 informationTowards regime identification and appropriate chemistry tabulation for computation of autoigniting turbulent reacting flows
Center for Turbulence Research Annual Research Briefs 009 199 Towards regime identification and appropriate chemistry tabulation for computation of autoigniting turbulent reacting flows By M. Kostka, E.
More informationFlame / wall interaction and maximum wall heat fluxes in diffusion burners
Flame / wall interaction and maximum wall heat fluxes in diffusion burners de Lataillade A. 1, Dabireau F. 1, Cuenot B. 1 and Poinsot T. 1 2 June 5, 2002 1 CERFACS 42 Avenue Coriolis 31057 TOULOUSE CEDEX
More informationREDIM reduced modeling of quenching at a cold inert wall with detailed transport and different mechanisms
26 th ICDERS July 3 th August 4 th, 217 Boston, MA, USA REDIM reduced modeling of quenching at a cold inert wall with detailed transport and different mechanisms Christina Strassacker, Viatcheslav Bykov,
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 informationA validation study of the flamelet approach s ability to predict flame structure when fluid mechanics are fully resolved
Center for Turbulence Research Annual Research Briefs 2009 185 A validation study of the flamelet approach s ability to predict flame structure when fluid mechanics are fully resolved By E. Knudsen AND
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 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 informationFlamelet Analysis of Turbulent Combustion
Flamelet Analysis of Turbulent Combustion R.J.M. Bastiaans,2, S.M. Martin, H. Pitsch,J.A.vanOijen 2, and L.P.H. de Goey 2 Center for Turbulence Research, Stanford University, CA 9435, USA 2 Eindhoven University
More informationDARS overview, IISc Bangalore 18/03/2014
www.cd-adapco.com CH2O Temperatur e Air C2H4 Air DARS overview, IISc Bangalore 18/03/2014 Outline Introduction Modeling reactions in CFD CFD to DARS Introduction to DARS DARS capabilities and applications
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 informationInvestigation of ignition dynamics in a H2/air mixing layer with an embedded vortex
Paper # 070LT-0211 The 8th US National Meeting of the Combustion Institute, Park City, UT, May 19-22, 2013 Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex S.K. Menon
More informationLecture 7 Flame Extinction and Flamability Limits
Lecture 7 Flame Extinction and Flamability Limits 7.-1 Lean and rich flammability limits are a function of temperature and pressure of the original mixture. Flammability limits of methane and hydrogen
More informationApplication of FGM to DNS of premixed turbulent spherical flames
Application of FGM to DNS of premixed turbulent spherical flames R.J.M. Bastiaans, G.R.A Groot, J.A. van Oijen and L.P.H. de Goey, Section Combustion Technology, Department of Mechanical Engineering, Eindhoven
More informationAn Unsteady/Flamelet Progress Variable Method for LES of Nonpremixed Turbulent Combustion
43rd AIAA Aerospace Sciences Meeting and Exhibit, -3 Jan 25, Reno, NV An Unsteady/Flamelet Progress Variable Method for LES of Nonpremixed Turbulent Combustion Heinz Pitsch and Matthias Ihme Stanford University,
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 informationBudget analysis and model-assessment of the flamelet-formulation: Application to a reacting jet-in-cross-flow
Center for Turbulence Research Proceedings of the Summer Program 212 397 Budget analysis and model-assessment of the flamelet-formulation: Application to a reacting jet-in-cross-flow By W. L. Chan, Y.
More informationA Priori Testing of Flamelet Generated Manifolds for Turbulent Partially Premixed Methane/Air Flames
DOI 10.1007/s10494-009-9223-1 A Priori Testing of Flamelet Generated Manifolds for Turbulent Partially Premixed Methane/Air Flames W. J. S. Ramaekers J. A. van Oijen L. P. H. de Goey Received: 7 September
More informationNumerical Simulation of Hydrogen Gas Turbines using Flamelet Generated Manifolds technique on Open FOAM
Numerical Simulation of Hydrogen Gas Turbines using Flamelet Generated Manifolds technique on Open FOAM Alessio Fancello (M.Sc.) Department of Mechanical Engineering Combustion Technology Technische Universiteit
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 informationDNS and LES of Turbulent Combustion
Computational Fluid Dynamics In Chemical Reaction Engineering IV June 19-24, 2005 Barga, Italy DNS and LES of Turbulent Combustion Luc Vervisch INSA de Rouen, IUF, CORIA-CNRS Pascale Domingo, Julien Réveillon
More informationAsymptotic Structure of Rich Methane-Air Flames
Asymptotic Structure of Rich Methane-Air Flames K. SESHADRI* Center for Energy and Combustion Research, Department of Mechanical and Aerospace Engineering, University of California at San Diego, La Jolla,
More informationA Priori Model for the Effective Lewis Numbers in Premixed Turbulent Flames
Paper # 070LT-0267 Topic: Turbulent Flames 8 th US National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013.
More informationLARGE-EDDY SIMULATION OF PARTIALLY PREMIXED TURBULENT COMBUSTION
LARGE-EDDY SIMULATION OF PARTIALLY PREMIXED TURBULENT COMBUSTION Heinz Pitsch Mechanical Engineering Department Stanford University Stanford, CA 94305, USA h.pitsch@stanford.edu ABSTRACT The development
More informationModeling Turbulent Combustion
Modeling Turbulent Combustion CEFRC Combustion Summer School 2014 Prof. Dr.-Ing. Heinz Pitsch Copyright 2014 by Heinz Pitsch. This material is not to be sold, reproduced or distributed without prior written
More informationD. VEYNANTE. Introduction à la Combustion Turbulente. Dimanche 30 Mai 2010, 09h00 10h30
D. VEYNANTE Introduction à la Combustion Turbulente Dimanche 30 Mai 2010, 09h00 10h30 Introduction to turbulent combustion D. Veynante Laboratoire E.M2.C. CNRS - Ecole Centrale Paris Châtenay-Malabry France
More informationLarge-eddy simulation of an industrial furnace with a cross-flow-jet combustion system
Center for Turbulence Research Annual Research Briefs 2007 231 Large-eddy simulation of an industrial furnace with a cross-flow-jet combustion system By L. Wang AND H. Pitsch 1. Motivation and objectives
More informationUnsteady Flamelet Modeling of Soot Formation in Turbulent Diffusion Flames
Unsteady Flamelet Modeling of Soot Formation in Turbulent Diffusion Flames H. Pitsch Department of Applied Mechanics and Engineering Science Center for Energy and Combustion Research, University of California
More informationConsistent turbulence modeling in a hybrid LES/RANS PDF method for non-premixed flames
Consistent turbulence modeling in a hybrid LES/RANS PDF method for non-premixed flames F. Ferraro, Y. Ge and M. Pfitzner Institut für Thermodynamik, Fakultät für Luft- und Raumfahrrttechnik Universität
More 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 informationCoupling tabulated chemistry with large-eddy simulation of turbulent reactive flows
Center for Turbulence Research Proceedings of the Summer Program 2008 237 Coupling tabulated chemistry with large-eddy simulation of turbulent reactive flows By R. Vicquelin, B. Fiorina, N. Darabiha, D.
More informationSection 14.1: Description of the Equilibrium Mixture Fraction/PDF Model. Section 14.2: Modeling Approaches for Non-Premixed Equilibrium Chemistry
Chapter 14. Combustion Modeling Non-Premixed In non-premixed combustion, fuel and oxidizer enter the reaction zone in distinct streams. This is in contrast to premixed systems, in which reactants are mixed
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 informationCharacterization of Flow Field Structure and Species Composition in a Shear Coaxial Rocket GH2/GO2 Injector: Modeling of Wall Heat Losses. D. J.
47th AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit 31 July - 03 August 2011 San Diego California AIAA 2011-6125 Characterization of Flow Field Structure and Species Composition in a Shear Coaxial
More informationExtinction and reignition in a diffusion flame: a direct numerical simulation study
J. Fluid Mech. (24), vol. 518, pp. 231 259. c 24 Cambridge University Press DOI: 1.117/S22112414 Printed in the United Kingdom 231 Extinction and reignition in a diffusion flame: a direct numerical simulation
More informationOverview of Turbulent Reacting Flows
Overview of Turbulent Reacting Flows Outline Various Applications Overview of available reacting flow models LES Latest additions Example Cases Summary Reacting Flows Applications in STAR-CCM+ Ever-Expanding
More informationInvestigation of ignition dynamics in a H2/air mixing layer with an embedded vortex
Paper # 070LT-0211 The 8th US National Meeting of the Combustion Institute, Park City, UT, May 19-22, 2013 Investigation of ignition dynamics in a H2/air mixing layer with an embedded vortex S.K. Menon
More informationLecture 10 Turbulent Combustion: The State of the Art
Lecture 10 Turbulent Combustion: The State of the Art 10.-1 Engineering applications are typically turbulent turbulence models These models use systematic mathematical derivations based on the Navier-
More informationA priori Tabulation of Turbulent Flame Speeds via a Combination of a Stochastic Mixing Model and Flamelet Generated Manifolds 5
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany HEIKO SCHMIDT 1 MICHAEL OEVERMANN 2 ROB J.M. BASTIAANS 3 ALAN R. KERSTEIN 4 A priori Tabulation of Turbulent
More informationCFD Simulation of high pressure real gas flows
CFD Simulation of high pressure real gas flows on the progress form art to physics :o) Maria Magdalena Poschner Prof. Dr. rer. nat. Pfitzner UNIVERSITÄT DER BUNDESWEHR MÜNCHEN Fakultät für Luft und Raumfahrt
More informationRepresenting Turbulence/Chemistry interaction with strained planar premixed
Representing Turbulence/Chemistry interaction with strained planar premixed flames CES Seminar Andreas Pesch Matr. Nr. 300884 Institut für technische Verbrennung Prof. H. Pitsch Rheinisch-Westfälische
More informationIMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE FLUCTUATIONS
Proceedings of the Combustion Institute, Volume 9, 00/pp. 1971 1978 IMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE
More informationAsymptotic Analysis of the Structure of Moderately Rich Methane-Air Flames
Asymptotic Analysis of the Structure of Moderately Rich Methane-Air Flames K. SESHADRI,* X. S. BAI,** H. PITSCH, and N. PETERS Institut für Technische Mechanik, RWTH Aachen, D-52056 Aachen, Federal Republic
More informationA G-equation formulation for large-eddy simulation of premixed turbulent combustion
Center for Turbulence Research Annual Research Briefs 2002 3 A G-equation formulation for large-eddy simulation of premixed turbulent combustion By H. Pitsch 1. Motivation and objectives Premixed turbulent
More 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 informationExamination of the effect of differential molecular diffusion in DNS of turbulent non-premixed flames
Examination of the effect of differential molecular diffusion in DNS of turbulent non-premixed flames Chao Han a, David O. Lignell b, Evatt R. Hawkes c, Jacqueline H. Chen d, Haifeng Wang a, a School of
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 informationNumerical Investigation of Flow and Combustion in a Single-Element GCH 4 /GOX Rocket Combustor
Numerical Investigation of Flow and Combustion in a Single-Element GCH 4 /GOX Rocket Combustor C. Roth, O. Haidn, A. Chemnitz, T. Sattelmayer Technische Universität München, 80333 Munich, Germany Y. Daimon
More informationPublished in: Proceedings of the Fluent Benelux User Group Meeting, 6-7 October 2005, Wavre, Belgium
The application of Flamelet Generated Manifolds in partailly-premixed flames Ramaekers, W.J.S.; Albrecht, B.A.; van Oijen, J.A.; de Goey, L.P.H.; Eggels, R.L.G.M. Published in: Proceedings of the Fluent
More informationLES of an auto-igniting C 2 H 4 flame DNS
Center for Turbulence Research Annual Research Briefs 2011 237 LES of an auto-igniting C 2 H 4 flame DNS By E. Knudsen, E. S. Richardson, J. H. Chen AND H. Pitsch 1. Motivation and objectives Large eddy
More informationLES of the Sandia Flame D Using an FPV Combustion Model
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 82 (2015 ) 402 409 ATI 2015-70th Conference of the ATI Engineering Association LES of the Sandia Flame D Using an FPV Combustion
More 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 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 informationLarge-Eddy Simulation of Turbulent Combustion
Annu. Rev. Fluid Mech. 2006. 38:453 82 The Annual Review of Fluid Mechanics is online at fluid.annualreviews.org doi: 10.1146/annurev.fluid. 38.050304.092133 Copyright c 2006 by Annual Reviews. All rights
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 informationModeling laminar partially premixed flames with complex diffusion in OpenFOAM R
Modeling laminar partially premixed flames with complex diffusion in OpenFOAM R Danny Messig, Jens Keller, Bernd Meyer, Christian Hasse Institute of Energy Process Engineering and Chemical Engineering
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 informationNumerical Investigation of Flow and Combustion in a Single-Element GCH4/GOX Rocket Combustor: Aspects of Turbulence Modeling
Numerical Investigation of Flow and Combustion in a Single-Element GCH4/GOX Rocket Combustor: Aspects of Turbulence Modeling A. Chemnitz, T. Sattelmayer, C. Roth, O. Haidn Technische Universität München,
More informationEffects of Damköhler number on flame extinction and reignition in turbulent nonpremixed flames using DNS
Effects of Damköhler number on flame extinction and reignition in turbulent nonpremixed flames using DNS David O. Lignell a,, Jacqueline H. Chen b, Hans A. Schmutz a a Chemical Engineering Department,
More informationThe Effect of Mixture Fraction on Edge Flame Propagation Speed
8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 213 The Effect of Mixture Fraction on Edge Flame
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 informationA first investigation on using a species reaction mechanism for flame propagation and soot emissions in CFD of SI engines
A first investigation on using a 1000+ species reaction mechanism for flame propagation and soot emissions in CFD of SI engines F.A. Tap *, D. Goryntsev, C. Meijer, A. Starikov Dacolt International BV
More informationSelf-similar behavior of chemistry tabulation in laminar and turbulent multi-fuel injection combustion systems
Center for Turbulence Research Proceedings of the Summer Program 26 39 Self-similar behavior of chemistry tabulation in laminar and turbulent multi-fuel injection combustion systems By A. Naudin, B. Fiorina,
More informationDevelopment of One-Step Chemistry Models for Flame and Ignition Simulation
Development of One-Step Chemistry Models for Flame and Ignition Simulation S.P.M. Bane, J.L. Ziegler, and J.E. Shepherd Graduate Aerospace Laboratories California Institute of Technology Pasadena, CA 91125
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 informationHIGH PRESSURE METHANE-OXYGEN COMBUSTION KINETIC ANALYSIS
HIGH PRESSURE METHANE-OXYGEN COMBUSTION KINETIC ANALYSIS G. Saccone*, P. Natale*, F. Battista* g.saccone@cira.it p.natale@cira.it f.battista@cira.it *CIRA Italian Aerospace Research Centre, Capua Italy,
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 informationTOPICAL PROBLEMS OF FLUID MECHANICS 97
TOPICAL PROBLEMS OF FLUID MECHANICS 97 DOI: http://dx.doi.org/10.14311/tpfm.2016.014 DESIGN OF COMBUSTION CHAMBER FOR FLAME FRONT VISUALISATION AND FIRST NUMERICAL SIMULATION J. Kouba, J. Novotný, J. Nožička
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 informationLarge-Eddy Simulation of Supercritical Mixing Layers
Sonderforschungsbereich/Transregio 4 Annual Report 211 13 Large-Eddy Simulation of Supercritical Mixing Layers By C.A. Niedermeier, M. Jarczyk, S. Hickel, N.A. Adams AND M. Pfitzner Institute of Aerodynamics
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 informationCoupling of ChemApp and OpenFOAM
Coupling of ChemApp and OpenFOAM Messig, Danny; Rehm, Markus; Meyer, Bernd 19th May 2009 Contents 1 Introduction 2 Software OpenFOAM ChemApp Cantera 3 Coupling of chemistry packages and OpenFOAM 4 Testcases
More informationLEWIS NUMBER EFFECTS ON THE STRUCTURE AND EXTINCTION OF DIFFUSION FLAMES DUE TO STRAIN. A. Liñan
LEWIS NUMBER EFFECTS ON THE STRUCTURE AND EXTINCTION OF DIFFUSION FLAMES DUE TO STRAIN A. Liñan Escuela Técnica Superior de Ingenieros Aeronáuticos Universidad Politécnica de Madrid SPAIN In turbulent
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 informationModeling flame brush thickness in premixed turbulent combustion
Center for Turbulence Research Proceedings of the Summer Program 2006 299 Modeling flame brush thickness in premixed turbulent combustion By E. Knudsen, O. Kurenkov, S. Kim, M. Oberlack AND H. Pitsch Turbulent
More informationInteractions between oxygen permeation and homogeneous-phase fuel conversion on the sweep side of an ion transport membrane
Interactions between oxygen permeation and homogeneous-phase fuel conversion on the sweep side of an ion transport membrane The MIT Faculty has made this article openly available. Please share how this
More informationThe development of advanced combustion systems is mainly controlled by the objective to increase fuel
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 29, Orlando, Florida AIAA 29-239 Large-Eddy Simulation of a Turbulent Lifted Flame in a Vitiated
More informationPDF Modeling and Simulation of Premixed Turbulent Combustion
Monte Carlo Methods Appl. Vol. No. (), pp. 43 DOI 5 / MCMA.7. c de Gruyter PDF Modeling and Simulation of Premixed Turbulent Combustion Michael Stöllinger and Stefan Heinz Abstract. The use of probability
More 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 informationFlamelet modelling of non-premixed turbulent combustion with local extinction and re-ignition
INSTITUTE OF PHYSICS PUBLISHING Combust. Theory Modelling 7 (2003) 317 332 COMBUSTION THEORY AND MODELLING PII: S1364-7830(03)38343-3 Flamelet modelling of non-premixed turbulent combustion with local
More informationOverview of Reacting Flow
Overview of Reacting Flow Outline Various Applications Overview of available reacting flow models Latest additions Example Cases Summary Reacting Flows Applications in STAR-CCM+ Chemical Process Industry
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 informationTriple flame: Inherent asymmetries and pentasectional character
Combustion Theory and Modelling, 2014 Vol. 18, No. 3, 454 473, http://dx.doi.org/10.1080/13647830.2014.923116 Triple flame: Inherent asymmetries and pentasectional character Albert Jordà Juanós and William
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 informationCombustion Reaction Model Generation using Principal Component Analysis
Paper # 7F-6 Topic: Turbulent Combustion 7 Fall Meeting of the Western States Section of the Combustion Institute Sandia National Laboratories, Livermore, CA October 6 & 7, 7. Combustion Reaction Model
More informationSubgrid-scale mixing of mixture fraction, temperature, and species mass fractions in turbulent partially premixed flames
Available online at www.sciencedirect.com Proceedings of the Combustion Institute 34 (2013) 1231 1239 Proceedings of the Combustion Institute www.elsevier.com/locate/proci Subgrid-scale mixing of mixture
More informationEffects of radiative heat loss on the extinction of counterflow premixed H 2 air flames
Combust. Theory Modelling 4 (2000) 459 475. Printed in the UK PII: S1364-7830(00)09647-9 Effects of radiative heat loss on the extinction of counterflow premixed H 2 air flames Hongsheng Guo, Yiguang Ju
More informationThe role of diffusion at shear layers in irregular detonations
The role of diffusion at shear layers in irregular detonations Marco Arienti 1 Joseph E. Shepherd 2 1 United Technologies Research Center, 411 Silver Lane, East Hartford, CT 06108 2 California Institute
More informationSmall scales, many species and the manifold challenges of turbulent combustion
Available online at www.sciencedirect.com Proceedings of the Combustion Institute 34 (2013) 1 31 Proceedings of the Combustion Institute www.elsevier.com/locate/proci Small scales, many species and the
More informationModelling Detailed-Chemistry Effects on Turbulent Diffusion Flames using a Parallel Solution-Adaptive Scheme
Modelling Detailed-Chemistry Effects on Turbulent Diffusion Flames using a Parallel Solution-Adaptive Scheme by Pradeep Kumar Jha A thesis submitted in conformity with the requirements for the degree of
More informationThermal NO Predictions in Glass Furnaces: A Subgrid Scale Validation Study
Feb 12 th 2004 Thermal NO Predictions in Glass Furnaces: A Subgrid Scale Validation Study Padmabhushana R. Desam & Prof. Philip J. Smith CRSIM, University of Utah Salt lake city, UT-84112 18 th Annual
More informationLarge-eddy simulation of supersonic reacting flows
Sonderforschungsbereich/Transregio 4 Proceedings of the Summer Program 23 2 Large-eddy simulation of supersonic reacting flows By G. Ribert, L. Bouheraoua AND P. Domingo CORIA - UMR664 CNRS, INSA and Université
More informationDevelopment of Reduced Mechanisms for Numerical Modelling of Turbulent Combustion
Worshop on Numerical Aspects of Reduction in Chemical Kinetics CERMICS-ENPC Cite Descartes - Champus sur Marne, France, September 2nd, 1997 Abstract Development of Reduced Mechanisms for Numerical Modelling
More informationModeling Flame Structure in Wildland Fires Using the One- Dimensional Turbulence Model
Modeling Flame Structure in Wildland Fires Using the One- Dimensional Turbulence Model David O. Lignell Chemical Engineering Dept., Brigham Young University, 350 CB, Provo, UT, USA, davidlignell@byu.edu
More informationCOUNTERFLOW COMBUSTION MODELING OF CH 4 /AIR AND CH 4 /O 2 INCLUDING DETAILED CHEMISTRY
Paper ID ILASS08-A046 ILASS08-11-4 ILASS 2008 Sep. 8-10, 2008, Como Lake, Italy COUNTERFLOW COMBUSTION MODELING OF CH 4 /AIR AND CH 4 /O 2 INCLUDING DETAILED CHEMISTRY D. Urzica, E. Gutheil Interdisziplinäres
More informationAME 513. " Lecture 8 Premixed flames I: Propagation rates
AME 53 Principles of Combustion " Lecture 8 Premixed flames I: Propagation rates Outline" Rankine-Hugoniot relations Hugoniot curves Rayleigh lines Families of solutions Detonations Chapman-Jouget Others
More information