A joint numerical study of multi-regime turbulent combustion

Benoît Fiorina; Tan Phong Luu; Samuel Dillon; Renaud Mercier; Ping Wang; Lorenzo Angelilli; Pietro Paolo Ciottoli; Francisco E. Hernández–Pérez; Mauro Valorani; Hong G. Im; James C. Massey; Zhiyi Li; Zhi X. Chen; Nedunchezhian Swaminathan; Sebastian Popp; Sandra Hartl; Hendrik Nicolai; Christian Hasse; Andreas Dreizler; David Butz; Dirk Geyer; Adrian Breicher; Kai Zhang; Christophe Duwig; Weijie Zhang; Wang Han; Jeroen van Oijen; Arthur Péquin; Alessandro Parente; Linus Engelmann; Andreas Kempf; Maximilian Hansinger; Michael Pfitzner; Robert S. Barlow

doi:10.1016/j.jaecs.2023.100221

A joint numerical study of multi-regime turbulent combustion

Benoît Fiorina (Corresponding author)
, Tan Phong Luu
, Samuel Dillon
, Renaud Mercier
, Ping Wang
, Lorenzo Angelilli
, Pietro Paolo Ciottoli
, Francisco E. Hernández–Pérez
, Mauro Valorani
, Hong G. Im
, James C. Massey
, Zhiyi Li
, Zhi X. Chen
, Nedunchezhian Swaminathan
, Sebastian Popp
, Sandra Hartl
, Hendrik Nicolai
, Christian Hasse
, Andreas Dreizler
, David Butz

Dirk Geyer, Adrian Breicher, Kai Zhang, Christophe Duwig, Weijie Zhang, Wang Han, Jeroen van Oijen, Arthur Péquin, Alessandro Parente, Linus Engelmann, Andreas Kempf, Maximilian Hansinger, Michael Pfitzner, Robert S. Barlow

Research output: Contribution to journal › Article › Academic › peer-review

22 Citations (Scopus)

107 Downloads (Pure)

Abstract

This article presents a joint numerical study on the Multi Regime Burner configuration. The burner design consists of three concentric inlet streams, which can be operated independently with different equivalence ratios, allowing the operation of stratified flames characterized by different combustion regimes, including premixed, non-premixed, and multi-regime flame zones. Simulations were performed on three LES solvers based on different numerical methods. Combustion kinetics were simplified by using tabulated or reduced chemistry methods. Finally, different turbulent combustion modeling strategies were employed, covering geometrical, statistical, and reactor based approaches. Due to this significant scattering of simulation parameters, a conclusion on specific combustion model performance is impossible. However, with ten numerical groups involved in the numerical simulations, a rough statistical analysis is conducted: the average and the standard deviation of the numerical simulation are computed and compared against experiments. This joint numerical study is therefore a partial illustration of the community's ability to model turbulent combustion. This exercise gives the average performance of current simulations and identifies physical phenomena not well captured today by most modeling strategies. Detailed comparisons between experimental and numerical data along radial profiles taken at different axial positions showed that the temperature field is fairly well captured up to 60 mm from the burner exit. The comparison reveals, however, significant discrepancies regarding CO mass fraction prediction. Three causes may explain this phenomenon. The first reason is the higher sensitivity of carbon monoxide to the simplification of detailed chemistry, especially when multiple combustion regimes are encountered. The second is the bias introduced by artificial thickening, which overestimates the species’ mass production rate. This behavior has been illustrated by manufacturing mean thickened turbulent flame brush from a random displacement of 1-D laminar flame solutions. The last one is the influence of the subgrid-scale flame wrinkling on the filtered chemical flame structure, which may be challenging to model.

Original language	English
Article number	100221
Number of pages	17
Journal	Applications in Energy and Combustion Science
Volume	16
Early online date	10 Nov 2023
DOIs	https://doi.org/10.1016/j.jaecs.2023.100221
Publication status	Published - Dec 2023

Funding

We are grateful for the opportunities provided by the TNF workshop (https://tnfworkshop.org). The simulations performed by EM2C-CNRS was granted access to the HPC resources from IDRIS, TGCC and CINES under the allocation A0052B10253 made by GENCI (Grand Equipement National de Calcul Intensif) and HPC resources from the Mesocentre computing center of CentraleSupelec and Ecole Normale Supérieure Paris-Saclay supported by CNRS and Région Ile-de-France (http://mesocentre.centralesupelec.fr/). Vincent Moureau and Ghislain Lartigue from CORIA lab are acknowledged for providing the YALES2 code. The author from Jiangsu University acknowledges the support by the Natural Science Foundation of China (NSFC) under Grant No. 91741117. KTH simulations were performed on resources provided by the National Academic Infrastructure for Supercomputing in Sweden (NAISS) at NSC and PDC super-computing centers. The authors from TUDa acknowledge funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Projektnummer 325144795. Work from ULB has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme under grant agreement No 714605 and Marie Skłodowska-Curie grant agreement No 801505. A. Péquin acknowledges the financial support of the Fonds National de la Recherche Scientifique (FRS-FNRS). J.C.M. and N.S. acknowledge the financial support from Mitsubishi Heavy Industries, Ltd. Takasago, Japan. Z.L. and N.S. acknowledge the financial support from the EPSRC (grant number EP/S025650/1). The work of Z.X.C. is supported by the NSFC (grant number 52276096) and the Emerging Interdisciplinary-Young Scholars Project, Peking University, the Fundamental Research Funds for the Central Universities. This work used the ARCHER2 UK National Supercomputing Service (https://www.archer2.ac.uk). J.C.M. Z.L. Z.X.C and N.S. are grateful to the EPSRC (grant number EP/R029369/1) and ARCHER2 for financial and computational support as a part of their funding to the UK Consortium on Turbulent Reacting Flows (https://www.ukctrf.com). We are grateful for the opportunities provided by the TNF workshop ( https://tnfworkshop.org ). The simulations performed by EM2C-CNRS was granted access to the HPC resources from IDRIS, TGCC and CINES under the allocation A0052B10253 made by GENCI (Grand Equipement National de Calcul Intensif) and HPC resources from the Mesocentre computing center of CentraleSupelec and Ecole Normale Supérieure Paris-Saclay supported by CNRS and Région Ile-de-France ( http://mesocentre.centralesupelec.fr/ ). Vincent Moureau and Ghislain Lartigue from CORIA lab are acknowledged for providing the YALES2 code. The author from Jiangsu University acknowledges the support by the Natural Science Foundation of China (NSFC) under Grant No. 91741117 . KTH simulations were performed on resources provided by the National Academic Infrastructure for Supercomputing in Sweden (NAISS) at NSC and PDC super-computing centers. The authors from TUDa acknowledge funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - Projektnummer 325144795 . Work from ULB has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreement No 714605 and Marie Skłodowska-Curie grant agreement No 801505 . A. Péquin acknowledges the financial support of the Fonds National de la Recherche Scientifique (FRS-FNRS) . J.C.M. and N.S. acknowledge the financial support from Mitsubishi Heavy Industries, Ltd., Takasago, Japan . Z.L. and N.S. acknowledge the financial support from the EPSRC (grant number EP/S025650/1 ). The work of Z.X.C. is supported by the NSFC (grant number 52276096 ) and the Emerging Interdisciplinary-Young Scholars Project, Peking University , the Fundamental Research Funds for the Central Universities . This work used the ARCHER2 UK National Supercomputing Service ( https://www.archer2.ac.uk ). J.C.M., Z.L., Z.X.C and N.S. are grateful to the EPSRC (grant number EP/R029369/1 ) and ARCHER2 for financial and computational support as a part of their funding to the UK Consortium on Turbulent Reacting Flows ( https://www.ukctrf.com ).

Keywords

Large Eddy Simulation
Model comparison
Multi regime combustion
Pollutant formation
Turbulent combustion modeling

Access to Document

10.1016/j.jaecs.2023.100221

pdfFinal published version, 7.43 MBLicence: CC BY

Cite this

Fiorina, B., Luu, T. P., Dillon, S., Mercier, R., Wang, P., Angelilli, L., Ciottoli, P. P., Hernández–Pérez, F. E., Valorani, M., Im, H. G., Massey, J. C., Li, Z., Chen, Z. X., Swaminathan, N., Popp, S., Hartl, S., Nicolai, H., Hasse, C., Dreizler, A., ... Barlow, R. S. (2023). A joint numerical study of multi-regime turbulent combustion. Applications in Energy and Combustion Science, 16, Article 100221. https://doi.org/10.1016/j.jaecs.2023.100221

@article{6738aa2d42e241599bca3d35a5a883af,

title = "A joint numerical study of multi-regime turbulent combustion",

abstract = "This article presents a joint numerical study on the Multi Regime Burner configuration. The burner design consists of three concentric inlet streams, which can be operated independently with different equivalence ratios, allowing the operation of stratified flames characterized by different combustion regimes, including premixed, non-premixed, and multi-regime flame zones. Simulations were performed on three LES solvers based on different numerical methods. Combustion kinetics were simplified by using tabulated or reduced chemistry methods. Finally, different turbulent combustion modeling strategies were employed, covering geometrical, statistical, and reactor based approaches. Due to this significant scattering of simulation parameters, a conclusion on specific combustion model performance is impossible. However, with ten numerical groups involved in the numerical simulations, a rough statistical analysis is conducted: the average and the standard deviation of the numerical simulation are computed and compared against experiments. This joint numerical study is therefore a partial illustration of the community's ability to model turbulent combustion. This exercise gives the average performance of current simulations and identifies physical phenomena not well captured today by most modeling strategies. Detailed comparisons between experimental and numerical data along radial profiles taken at different axial positions showed that the temperature field is fairly well captured up to 60 mm from the burner exit. The comparison reveals, however, significant discrepancies regarding CO mass fraction prediction. Three causes may explain this phenomenon. The first reason is the higher sensitivity of carbon monoxide to the simplification of detailed chemistry, especially when multiple combustion regimes are encountered. The second is the bias introduced by artificial thickening, which overestimates the species{\textquoteright} mass production rate. This behavior has been illustrated by manufacturing mean thickened turbulent flame brush from a random displacement of 1-D laminar flame solutions. The last one is the influence of the subgrid-scale flame wrinkling on the filtered chemical flame structure, which may be challenging to model.",

keywords = "Large Eddy Simulation, Model comparison, Multi regime combustion, Pollutant formation, Turbulent combustion modeling",

author = "Beno{\^i}t Fiorina and Luu, \{Tan Phong\} and Samuel Dillon and Renaud Mercier and Ping Wang and Lorenzo Angelilli and Ciottoli, \{Pietro Paolo\} and Hern{\'a}ndez{\textendash}P{\'e}rez, \{Francisco E.\} and Mauro Valorani and Im, \{Hong G.\} and Massey, \{James C.\} and Zhiyi Li and Chen, \{Zhi X.\} and Nedunchezhian Swaminathan and Sebastian Popp and Sandra Hartl and Hendrik Nicolai and Christian Hasse and Andreas Dreizler and David Butz and Dirk Geyer and Adrian Breicher and Kai Zhang and Christophe Duwig and Weijie Zhang and Wang Han and \{van Oijen\}, Jeroen and Arthur P{\'e}quin and Alessandro Parente and Linus Engelmann and Andreas Kempf and Maximilian Hansinger and Michael Pfitzner and Barlow, \{Robert S.\}",

year = "2023",

month = dec,

doi = "10.1016/j.jaecs.2023.100221",

language = "English",

volume = "16",

journal = "Applications in Energy and Combustion Science",

issn = "2666-352X",

publisher = "Elsevier",

}

Fiorina, B, Luu, TP, Dillon, S, Mercier, R, Wang, P, Angelilli, L, Ciottoli, PP, Hernández–Pérez, FE, Valorani, M, Im, HG, Massey, JC, Li, Z, Chen, ZX, Swaminathan, N, Popp, S, Hartl, S, Nicolai, H, Hasse, C, Dreizler, A, Butz, D, Geyer, D, Breicher, A, Zhang, K, Duwig, C, Zhang, W, Han, W, van Oijen, J, Péquin, A, Parente, A, Engelmann, L, Kempf, A, Hansinger, M, Pfitzner, M & Barlow, RS 2023, 'A joint numerical study of multi-regime turbulent combustion', Applications in Energy and Combustion Science, vol. 16, 100221. https://doi.org/10.1016/j.jaecs.2023.100221

TY - JOUR

T1 - A joint numerical study of multi-regime turbulent combustion

AU - Fiorina, Benoît

AU - Luu, Tan Phong

AU - Dillon, Samuel

AU - Mercier, Renaud

AU - Wang, Ping

AU - Angelilli, Lorenzo

AU - Ciottoli, Pietro Paolo

AU - Hernández–Pérez, Francisco E.

AU - Valorani, Mauro

AU - Im, Hong G.

AU - Massey, James C.

AU - Li, Zhiyi

AU - Chen, Zhi X.

AU - Swaminathan, Nedunchezhian

AU - Popp, Sebastian

AU - Hartl, Sandra

AU - Nicolai, Hendrik

AU - Hasse, Christian

AU - Dreizler, Andreas

AU - Butz, David

AU - Geyer, Dirk

AU - Breicher, Adrian

AU - Zhang, Kai

AU - Duwig, Christophe

AU - Zhang, Weijie

AU - Han, Wang

AU - van Oijen, Jeroen

AU - Péquin, Arthur

AU - Parente, Alessandro

AU - Engelmann, Linus

AU - Kempf, Andreas

AU - Hansinger, Maximilian

AU - Pfitzner, Michael

AU - Barlow, Robert S.

PY - 2023/12

Y1 - 2023/12

N2 - This article presents a joint numerical study on the Multi Regime Burner configuration. The burner design consists of three concentric inlet streams, which can be operated independently with different equivalence ratios, allowing the operation of stratified flames characterized by different combustion regimes, including premixed, non-premixed, and multi-regime flame zones. Simulations were performed on three LES solvers based on different numerical methods. Combustion kinetics were simplified by using tabulated or reduced chemistry methods. Finally, different turbulent combustion modeling strategies were employed, covering geometrical, statistical, and reactor based approaches. Due to this significant scattering of simulation parameters, a conclusion on specific combustion model performance is impossible. However, with ten numerical groups involved in the numerical simulations, a rough statistical analysis is conducted: the average and the standard deviation of the numerical simulation are computed and compared against experiments. This joint numerical study is therefore a partial illustration of the community's ability to model turbulent combustion. This exercise gives the average performance of current simulations and identifies physical phenomena not well captured today by most modeling strategies. Detailed comparisons between experimental and numerical data along radial profiles taken at different axial positions showed that the temperature field is fairly well captured up to 60 mm from the burner exit. The comparison reveals, however, significant discrepancies regarding CO mass fraction prediction. Three causes may explain this phenomenon. The first reason is the higher sensitivity of carbon monoxide to the simplification of detailed chemistry, especially when multiple combustion regimes are encountered. The second is the bias introduced by artificial thickening, which overestimates the species’ mass production rate. This behavior has been illustrated by manufacturing mean thickened turbulent flame brush from a random displacement of 1-D laminar flame solutions. The last one is the influence of the subgrid-scale flame wrinkling on the filtered chemical flame structure, which may be challenging to model.

AB - This article presents a joint numerical study on the Multi Regime Burner configuration. The burner design consists of three concentric inlet streams, which can be operated independently with different equivalence ratios, allowing the operation of stratified flames characterized by different combustion regimes, including premixed, non-premixed, and multi-regime flame zones. Simulations were performed on three LES solvers based on different numerical methods. Combustion kinetics were simplified by using tabulated or reduced chemistry methods. Finally, different turbulent combustion modeling strategies were employed, covering geometrical, statistical, and reactor based approaches. Due to this significant scattering of simulation parameters, a conclusion on specific combustion model performance is impossible. However, with ten numerical groups involved in the numerical simulations, a rough statistical analysis is conducted: the average and the standard deviation of the numerical simulation are computed and compared against experiments. This joint numerical study is therefore a partial illustration of the community's ability to model turbulent combustion. This exercise gives the average performance of current simulations and identifies physical phenomena not well captured today by most modeling strategies. Detailed comparisons between experimental and numerical data along radial profiles taken at different axial positions showed that the temperature field is fairly well captured up to 60 mm from the burner exit. The comparison reveals, however, significant discrepancies regarding CO mass fraction prediction. Three causes may explain this phenomenon. The first reason is the higher sensitivity of carbon monoxide to the simplification of detailed chemistry, especially when multiple combustion regimes are encountered. The second is the bias introduced by artificial thickening, which overestimates the species’ mass production rate. This behavior has been illustrated by manufacturing mean thickened turbulent flame brush from a random displacement of 1-D laminar flame solutions. The last one is the influence of the subgrid-scale flame wrinkling on the filtered chemical flame structure, which may be challenging to model.

KW - Large Eddy Simulation

KW - Model comparison

KW - Multi regime combustion

KW - Pollutant formation

KW - Turbulent combustion modeling

UR - https://www.scopus.com/pages/publications/85177060699

U2 - 10.1016/j.jaecs.2023.100221

DO - 10.1016/j.jaecs.2023.100221

M3 - Article

AN - SCOPUS:85177060699

SN - 2666-352X

VL - 16

JO - Applications in Energy and Combustion Science

JF - Applications in Energy and Combustion Science

M1 - 100221

ER -

A joint numerical study of multi-regime turbulent combustion

Abstract

Funding

Keywords

Access to Document

Other files and links

Fingerprint

Cite this