TY - JOUR

T1 - Reynolds number scaling of burning rates in spherical turbulent premixed flames

AU - Kulkarni, Tejas

AU - Buttay, Romain

AU - Kasbaoui, M. Houssem

AU - Attili, Antonio

AU - Bisetti, Fabrizio

N1 - Funding Information:
This material is based upon work supported in part by the National Science Foundation (NSF) under grant number 1805921. Numerical simulations were carried out on the 'Shaheen' supercomputer at King Abdullah University of Science and Technology (KAUST) and on the 'Stampede 2' supercomputer at the Texas Advanced Computing Center (TACC) through allocation TG-CTS180002 under the Extreme Science and Engineering Discovery Environment (XSEDE). XSEDE is supported by the NSF under grant number ACI-1548562. Any opinions, findings and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the NSF.
Funding Information:
This material is based upon work supported in part by the National Science Foundation (NSF) under grant number 1805921. Numerical simulations were carried out on the ‘Shaheen’ supercomputer at King Abdullah University of Science and Technology (KAUST) and on the ‘Stampede 2’ supercomputer at the Texas Advanced Computing Center (TACC) through allocation TG-CTS180002 under the Extreme Science and Engineering Discovery Environment (XSEDE). XSEDE is supported by the NSF under grant number ACI-1548562. Any opinions, findings and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the NSF.

PY - 2020

Y1 - 2020

N2 - In the flamelet regime of turbulent premixed combustion the enhancement in the burning rates originates primarily from surface wrinkling. In this work we investigate the Reynolds number dependence of burning rates of spherical turbulent premixed methane/air flames in decaying isotropic turbulence with direct numerical simulations. Several simulations are performed by varying the Reynolds number, while keeping the Karlovitz number the same, and the temporal evolution of the flame surface is compared across cases by combining the probability density function of the radial distance of the flame surface from the origin with the surface density function formalism. Because the mean area of the wrinkled flame surface normalized by the area of a sphere with radius equal to the mean flame radius is proportional to the product of the turbulent flame brush thickness and peak surface density within the brush, the temporal evolution of the brush and peak surface density are investigated separately. The brush thickness is shown to scale with the integral scale of the flow, evolving due to decaying velocity fluctuations and stretch. When normalized by the integral scale, the wrinkling scale defined as the inverse of the peak surface density is shown to scale with Reynolds number across simulations and as turbulence decays. As a result, the area ratio and the burning rate are found to increase as, in agreement with recent experiments on spherical turbulent premixed flames. We observe that the area ratio does not vary with turbulent intensity when holding the Reynolds number constant.

AB - In the flamelet regime of turbulent premixed combustion the enhancement in the burning rates originates primarily from surface wrinkling. In this work we investigate the Reynolds number dependence of burning rates of spherical turbulent premixed methane/air flames in decaying isotropic turbulence with direct numerical simulations. Several simulations are performed by varying the Reynolds number, while keeping the Karlovitz number the same, and the temporal evolution of the flame surface is compared across cases by combining the probability density function of the radial distance of the flame surface from the origin with the surface density function formalism. Because the mean area of the wrinkled flame surface normalized by the area of a sphere with radius equal to the mean flame radius is proportional to the product of the turbulent flame brush thickness and peak surface density within the brush, the temporal evolution of the brush and peak surface density are investigated separately. The brush thickness is shown to scale with the integral scale of the flow, evolving due to decaying velocity fluctuations and stretch. When normalized by the integral scale, the wrinkling scale defined as the inverse of the peak surface density is shown to scale with Reynolds number across simulations and as turbulence decays. As a result, the area ratio and the burning rate are found to increase as, in agreement with recent experiments on spherical turbulent premixed flames. We observe that the area ratio does not vary with turbulent intensity when holding the Reynolds number constant.

KW - isotropic turbulence

KW - turbulent reacting flows

UR - http://www.scopus.com/inward/record.url?scp=85095839096&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85095839096&partnerID=8YFLogxK

U2 - 10.1017/jfm.2020.784

DO - 10.1017/jfm.2020.784

M3 - Article

AN - SCOPUS:85095839096

VL - 906

JO - Journal of Fluid Mechanics

JF - Journal of Fluid Mechanics

SN - 0022-1120

M1 - A2

ER -