Development of Unified Flame Surface Density Based Reaction Rate Models for the LES of Turbulent Premixed Flames

Lead Research Organisation: University of Liverpool
Department Name: School of Engineering

Abstract

Abstracts are not currently available in GtR for all funded research. This is normally because the abstract was not required at the time of proposal submission, but may be because it included sensitive information such as personal details.

Publications

10 25 50

 
Description Introduction

----------------------

The aforementioned 3.5 year long project finished in January 2013. The EPSRC final report form in Je-S system has been completed. This short report summarises the major achievements of this project and how the benefits of this project will be utilised in the future.

Background

--------------------

The project was initially conceived to obtain fundamental insight and develop high fidelity Flame Surface Density (FSD) based reaction rate closure for turbulent premixed flames based on a-priori analysis of DNS data. Although FSD closures are relatively well established for Reynolds Averaged Navier Stokes (RANS) simulations, they are relatively rare for LES, and no detailed evaluation of their performance was done so far. In this project, FSD based reaction rate closures for both RANS and LES have been been developed and simultaneously assessed by a-priori analyses of explicitly filtered Direct Numerical Simulation (DNS) data, and a-posteriori evaluations of model performance in actual LES, in a number of configurations for which experimental data is available. The models have been proposed so that they account for fundamental physics of flame propagation and scalar gradient alignment statistics along with the influences of Damkohler, Karlovitz, turbulent Reynolds and global Lewis numbers on these statistics. This ensures that the newly proposed models remain valid for both the corrugated flamelets and thin reaction zones regimes. The performances of these models have been assessed based on a-priori analysis of three-dimensional DNS data. The best models identified based on a-priori analysis have been implemented in actual LES simulations in order to carry out a-posteriori analyses of the models in a number of configurations for which experimental data is available in comparison to LES simulation results. Development of high-fidelity FSD-based reaction rate closure will provide a reliable Computational Fluid Dynamics (CFD) based design tool for clean and cost-effective combustion devices where combustion takes place in the premixed mode (e.g. in Spark Ignition/ Direct Injection engines, Lean Premixed Pre-vaporised (LPP) or industrial gas turbine combustors).



Achievements

----------------------

1. Development of a three-dimensional DNS database of freely propagating statistically planar turbulent premixed flames for different values of root-mean-square turbulent velocity fluctuation, global Lewis number and integral length scales of turbulent velocity fluctuations for a range of different values of turbulent Reynolds number. This provides wealth of important information, which can be used for the purpose of assessing existing FSD-based models and develop new models in the context of RANS and LES. The database developed during the course of this project will continue to be a very useful source of information for gaining fundamental understanding, and for model development, beyond the duration of this project.

2. The performance of the existing algebraic LES-FSD models have been assessed in detail based on a-priori analysis of DNS data for different values of global Lewis number [1,2] and turbulent Reynolds number [3]. The existing models that are the best suited for predicting the FSD in the context of LES have been identified, and a new power-law based algebraic model has been proposed in such as manner so that the effects of Karlovitz, turbulent Reynolds and Lewis numbers on fractal dimension of flame surface and inner cut-off scale obtained from DNS data are explicitly accounted for [1-3].

3. The statistical behaviours of the different terms of the FSD transport equation have been analysed in detail based on explicitly filtered/Reynolds averaged DNS data [4-6]. Models have been proposed for all the unclosed terms for the transport equation of FSD in the context of RANS and LES based on a-priori analysis of DNS data [4-16].

4. The models for the tangential strain rate term of the FSD transport equation have been proposed for both LES and RANS in such a manner that the effects of Lewis number and turbulent Reynolds number on local alignment of scalar gradient with principle strain rates are explicitly accounted for, both in corrugated flamelets and thin reaction zones regime combustion [7,15,16]. The validity of these newly proposed models have been substantiated with respect to the tangential strain rate term extracted from explicitly Reynolds averaged/LES filtered DNS data.

5. The models for the curvature term of the FSD transport equation have been proposed for both LES and RANS in such a manner that the effects of curvature on local flame displacement and the influences of Lewis and turbulent Reynolds numbers on this dependence are explicitly accounted for both in corrugated flamelets and thin reaction zones regime combustion [8-14].

6. The effects of turbulent Reynolds number and global Lewis number on the statistics of flame displacement speed and its components in turbulent premixed flames have been analysed in detail [5,8,10,12,14,17,18]. The methodology of extracting actual three-dimensional displacement speed from two-dimensional experimental measurements has also been assessed based on DNS-based analysis [17].

7. The methodology for obtaining actual three-dimensional Reynolds averaged FSD and the terms of its transport equation from their two-dimensional counterparts, which can be measured experimentally with relative ease, has been devised in the course of this project [19]. The effects of regime of combustion, turbulent Reynolds number and Lewis number on the validity of the aforementioned methodology have been analysed in detail [19].

8. Novel modelling methodologies for turbulent transport of FSD and surface-filtered displacement speed have been devised in this analysis, which accounts for Lewis number effects and local stretch rate dependence of flame propagation.

9. The newly developed models have been implemented in an actual LES code, which has been eventually used for the a-posteriori assessments of the newly proposed models, and for the substantiation of their validity in collaboration with Prof. A. M. Kempf's research group (presently based at University of Duisburg, Germany, and earlier based at Imperial College London until 2011) [20, 21].

In addition, the models and findings have been shared with the Combustion Research Group in Cambridge University (especially with Prof. R. Stewart Cant, and Dr. N. Swaminathan) and Prof. A. M. Kempf's research group (presently based at University of Duisburg, Germany) who are also engaged in computational modelling of turbulent premixed combustion. The progress made during the project was also communicated to Ford, Rolls Royce and Siemens Plcs. and Electricity Supplies Research (ESR) network.



Future work

--------------------

A number of projects are already in motion to exploit the achievements of this project.

1. The fundamental understanding and modelling insight ganined from this project is current contributing to the development of scalar dissipation rate based reaction rate closure in the context of LES in an ongoing EPSRC research project (entitled "Large Eddy Simulation modelling of scalar dissipation rate based reaction rate closure in turbulent premixed flames", ESPRC Ref: EP/1028013/1).

2.The PhD student, Mr. Mohit Katragadda, who worked on this project, is in the process of applying for funding for postdoctoral research and if his application is successful he will be looking into LES modelling of turbulent premixed flames using DNS database in noncanonical configurations (e.g. V-flame configuration) with detailed chemistry.

3. Three papers on curvature and strain rate term modelling of the FSD transport equation are under preparation and will be submitted to reputed journals for publication. A part of this work will be presented in the European Combustion Meeting 2013. A paper has also been submitted to the International Colloquium of Dynamics of Explosions and Reactive Systems 2013 (ICDERS2013).

4.The research team of Principal Investigator (PI) is in the process of assessing the newly developed models using actual LES simulations in different configurations. The PI is also in the process of communicating the research outcomes to the colleagues in industrial sector (e.g. Ford, Rolls Royce, Siemens Plcs.) so that the research outcomes can be assimilated in an industrial set-up relatively rapidly.



References

-------------------

1 M. Katragadda, N. Chakraborty, "A-priori DNS assessment of algebraic LES Flame Surface Density models for non-unity Lewis number flames in the thin reaction zones regime", 33rd International Combustion Symposium, Tsinghua University, Beijing, 1-6th August, 2010.

2. M. Katragadda, N. Chakraborty, R.S. Cant, "A-priori DNS assessment of wrinkling factor based algebraic Flame Surface Density models in the context of Large Eddy Simulations for non-unity Lewis number flames in the thin reaction zones regime." J. Combust. , 794671, 2012.

3. M. Katragadda, N. Chakraborty, R.S. Cant, "Effects of turbulent Reynolds number on the performance of algebraic Flame Surface Density models for Large Eddy Simulation in the thin reaction zones regime: A Direct Numerical Simulation analysis", J. Combust. , 353257, 2012.

4. M. Katragadda, N. Chakraborty, "Effects of Lewis number on Flame Surface Density transport in the context of Large Eddy Simulation", 5th European Combustion Meeting, Cardiff University, 28th June-1st July, 2011.

5. N. Chakraborty, R.S. Cant, "Effects of Lewis number on Flame Surface Density transport in turbulent premixed combustion", Combust. Flame, 158, 1768-1787, 2011.

6. N. Chakraborty, R.S. Cant, "Turbulent Reynolds number dependence of Flame Surface Density transport in the context of Reynolds Averaged Navier Stokes Simulations", Proc. Combust. Instit., 34, 1347-1356, 2013.

7. M. Katragadda, S. P. Malkeson, N. Chakraborty, "Modelling of the tangential strain rate term of the Flame Surface Density transport equation in the context of Reynolds Averaged Navier Stokes Simulation", Proc. Combust. Instit., 33, 1429-1437, 2011.

8. M. Katragadda, S.P. Malkeson, N. Chakraborty, "Modelling of the curvature term of the Flame Surface Density transport equation for Reynolds Averaged Navier Stokes simulations", 7th Mediterranean Combustion Symposium, Chia Laguna, Cagliari, Sardinia, Italy, 11th -15th September, 2011.

9. M. Katragadda, N. Chakraborty, "Modelling of the curvature term of the Flame Surface Density Transport Equation in the context of Reynolds Averaged Navier Stokes Simulation", 34th International Combustion Symposium, Warsaw University of Technology, Warsaw, 29th July-3rd August, 2012.

10. M. Katragadda, N. Chakraborty, "Modelling of the curvature term of the Flame Surface Density Transport Equation for Large Eddy Simulations", 34th International Combustion Symposium, Warsaw University of Technology, Warsaw, 29th July-3rd August, 2012.

11.M. Katragadda, S.P. Malkeson, N. Chakraborty, "Statistical behaviour of the curvature term of the FSD transport equation in the context of RANS", 7th International Symposium on Turbulence, Heat and Mass Transfer, Palermo, Sicily, Italy, 24th -27th September 2012.

12. M. Katragadda, N. Chakraborty, "Effects of Lewis number on the curvature term of the Flame Surface Density transport equation for LES", 7th International Symposium on Turbulence, Heat and Mass Transfer, Palermo, Sicily, Italy, 24th -27th September 2012.

13. M. Katragadda, N. Chakraborty, "Modelling of the curvature term of the Flame Surface Density transport equation for Large Eddy Simulations", J. Combust., 915482, 2012.

14.M. Katragadda, N. Chakraborty, "A-priori Direct Numerical Simulation modelling of the curvature term of the Flame Surface Density transport equation for non-unity Lewis number flames in the context of Large Eddy Simulations." Int. J. Chemical Engg., 103727, 2012.

15.M. Katragadda, N. Chakraborty, "A-priori DNS modelling of the strain rate term of the Flame Surface Density Transport Equation in the Context of LES, 6th European Combustion Meeting, Lund University, Sweden, 25th -28th June, 2013 (accepted).

16.M. Katragadda, N. Chakraborty, "Modelling of the strain rate contribution to the FSD transport for non-unity Lewis number flames in LES", to International Colloquium on Dynamics of Explosions and Reactive Systems 2013 (ICDERS2013) to be held in National Central University, Taipei, Taiwan, 28th July-1st August 2013.

17.N. Chakraborty, M. Klein, R.S. Cant, "Effects of turbulent Reynolds number on the displacement speed statistics in the thin reaction zones regime turbulent premixed combustion." J. Combust., 2011, 473679, 2011.

18.N. Chakraborty, G. Hartung, M. Katragadda, C. F. Kaminski, "A numerical comparison of 2D and 3D density-weighted displacement speed statistics and implications for laser based measurements of flame displacement speed", Combust. Flame, 158, 1372-1390, 2011.

19.N. Chakraborty, E.R. Hawkes, "Determination of 3D Flame Surface Density variables from 2D measurements: Validation using Direct Numerical Simulation" Phys. Fluids, 23, 065113, 2011.

20.T. Ma, O. Stein, N. Chakraborty, A. Kempf, "Comparison of FSD models for LES", 33rd International Combustion Symposium, Tsinghua University, Beijing, 1-6th August, 2010.

21.T. Ma, O. Stein, N. Chakraborty, A. Kempf, "A-posteriori testing of Algebraic Flame Surface Density models for LES", Combust. Theor. Modell. (accepted).
Exploitation Route The major beneficiaries of this work are Internal Combustion (IC) engine and gas turbine manufacturers, who are engaged in developing new concepts for low-pollution and high-efficiency engines. The design process of combustion equipment in both sectors depends heavily on the predictive capability of engineering CFD calculations. In the UK, Ford as an IC engine manufacturer, and Rolls Royce and Siemens Plcs. as gas turbine manufacturers, have shown interest in the outcome of this work on RANS and LES modelling of premixed combustion. Moreover, the Electricity Supplies Research (ESR) network was also involved in project meetings and the model developments during the course of this project. However, the benefits are not limited to the above-mentioned UK manufacturers, as all automotive manufacturers have a program on Spark Ignition (SI) engines, and all gas turbine manufacturers are interested in LPP combustors. The relevant industrial partners were made aware of the developments that took place in this research during interactions at the aforementioned international conferences and UK based meetings. The present fundamental findings concerning turbulent stratified combustion will be of immediate value to all the industries mentioned above.



During the course of this project, Newcastle University has generated new expertise, which will help its scientific base and attract industrial interest. Finally, CFD software companies, who will be able to incorporate state-of-the-art combustion models in their codes and thereby increase their market penetration, will also be interested in the project outcome. Finally, the PhD student (Mr. Mohit Katragadda) is one of the direct beneficiaries of this work, who used this opportunity to work in this project to become a fine researcher with the potential to actively contribute to the UK and international turbulent combustion research in an independent manner in the future years. During the course of this project Mr. Katragadda received thorough training on turbulence, combustion theories, fluid dynamics, Computational Fluid Dynamics (CFD) and high performance computing. Mr. Katragadda collaborated with colleagues based at Imperial College London, and presented his work in a number of progress meetings and conferences, thus enabling him to develop a number of important transferable skills, which will help him immensely in his future assignments. 1. High-quality journal publications

2. Dissemination through conference presentations

3. Effective UK and international collaboration

4. Industrial contacts in Rolls Royce Plc., Siemens, Ford, Shell Plcs., MMI Engineering
Sectors Aerospace, Defence and Marine,Energy,Environment,Transport

 
Description The findings of this project gave a thorough insight into the Flame Surface Density based reaction rate closure for LES and its advantages/disadvantages in comparison to other alternative methods. The new closures developed/recommended in this project are robust as they have gone through both a-priori and a-posteriori test. Moreover, the closures are designed in such a manner that they can address differential diffusion effects which can play potentially a major part for future hydrogen-based economy. Apart from enriching the relevant field of research, the research outcomes provided high-fidelity simulation tools for simulating premixed turbulent combustion for Internal Combustion (IC) engine and gas turbine manufacturers (e.g. Ford, Rolls Royce and Siemens), which will contribute to the development of energy-efficient and environment-friendly devices and wealth generation in the future. This will give rise to considerable socio-economic impact. Moreover, this project gave rise to the development of two highly skilled professionals with expertise of turbulence, combustion, Direct Numerical Simulation, Large Eddy Simulation, parallel computing etc.
First Year Of Impact 2009
Sector Aerospace, Defence and Marine,Energy,Environment,Transport
Impact Types Societal,Economic

 
Description Large Eddy Simulation modelling of scalar dissipation rate based reaction rate closure in turbulent premixed flames
Amount £158,407 (GBP)
Funding ID EP/I028013/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 10/2011 
End 03/2015