skip to main content


Title: A Macroscopic View of Reynolds Scaling and Stretch Effects in Spherical Turbulent Premixed Flames
The burning rate in a spherically turbulent premixed flame is explored using direct numerical simulations, and a model of ordinary differential equations is proposed. The numerical dataset, from a previous work, is obtained from direct numerical simulations of confined spherical flames in isotropic turbulence over a range of Reynolds numbers. We begin the derivation of the model with an equation for the burning rate for the domain under consideration, and using a thin flame assumption and a two-fluid approach, we find the normalized turbulent burning rate to be controlled by the increase in flame surface area due to turbulent wrinkling, and correction factor which is observed to be consistently less than unity. A Reynolds scaling hypothesis for the flame turbulent wrinkling from a previous work using the same numerical dataset is used to model the term controlling the increase in flame surface area. The correction factor is hypothesized to reflect flame stretch effects, and hence this factor is modeled using Markstein theory applied to global averaged quantities. The ordinary differential equations are rewritten to reflect easily observable quantities such as the chamber pressure and mean flame radius, and then expressed in dimensionless form to assess dependence on various dimensionless parameters. The model predictions are found to be in good agreement with the numerical data within expected variances. Additionally, Markstein theory is found to be sufficient in describing the effects of flame stretch in the turbulent premixed flames under consideration.  more » « less
Award ID(s):
1805921
NSF-PAR ID:
10388098
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
AIAA Scitech 2022 Forum. Presented at AIAA Scitech 2022 Forum, San Diego CA & Virtual, Jan 3-7, 2022
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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 ${Re}_{\lambda }^{1.13}$ , 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. 
    more » « less
  2. The objective of this work is to provide physical insight into the mechanisms governing flame–turbulence interactions and explore the impact of the ubiquitous Darrieus–Landau instability on the propagation. It is based on the hydrodynamic theory of premixed flames that considers the flame thickness much smaller than all other fluidynamical length scales. In this asymptotic limit, the flame is thus confined to a surface whilst the diffusion and reaction processes occurring inside the flame zone are accounted for by two parameters: the unburned-to-burned density ratio and the Markstein length. The robust model, which is free of phenomenology and turbulence modelling assumptions, makes transparent the mutual interactions between the flame and the fluid flow, and permits examining trends in flame and flow characteristics while varying the turbulence intensity and mixture properties. It is used in this study to examine the morphological changes of the flame surface that result from the intertwined effects of the turbulence and instability, as demonstrated by the local displacement and curvature of the flame front, the extent of wrinkling and folding of the flame surface, and the overall flame brush thickness. It also provides a direct evaluation of the turbulent flame speed and its dependence on the mean flame curvature and on the hydrodynamic strain that it experiences. Also discussed are the effects of the flame on the flow by examining the various mechanisms of enstrophy and scalar gradient production/destruction, the degree of anisotropy created in the burned gas, and the restructuring of the vortical motion beyond the flame. 
    more » « less
  3. This work describes and tests the calibration process of the chemical-diffusive model (CDM) for the simulation of non-premixed diffusion flames. The CDM is an alternative, simplified approach for incorporating the effects of combustion in a fluid simulation, based on the ideas of regulating the rate of energy release such that the properties of combustion waves (e.g. flames and detonations) are reproduced. Past implementations of the CDM have considered single-stoichiometry fuel-air mixtures or mixtures with variable stoichiom- etry but with premixed modes of combustion. In this work, the CDM is tested and shown to work for non-premixed, low-Mach-number flames (i.e., diffusion flames) by incorporat- ing it into a numerical model which solves the reactive and compressible Navier-Stokes equations with the barely implicit correction (BIC) algorithm, which removes the acoustic limit on the integration time-step size. Simulations of one-dimensional premixed laminar flames reproduce the required premixed laminar flame speed, thickness, and temperature. A two-dimensional, steady-state, laminar coflow diffusion flame is computed, and the result demonstrates the ability of the algorithm to compute a non-premixed flame. Lastly, a two- dimensional simulation of two opposing jets of fuel and air show that the CDM approach can compute the structure of a counter-flow diffusion flame. 
    more » « less
  4. An optical parametric oscillator (OPO) is developed and characterized for the simultaneous generation of ultraviolet (UV) and near-UV nanosecond laser pulses for the single-shot Rayleigh scattering and planar laser-induced-fluorescence (PLIF) imaging of methylidyne (CH) and nitric oxide (NO) in turbulent flames. The OPO is pumped by a multichannel, 8-pulse Nd:YAG laser cluster that produces up to 225 mJ/pulse at 355 nm with pulse spacing of 100 µs. The pulsed OPO has a conversion efficiency of 9.6% to the signal wavelength of∼<#comment/>430nmwhen pumped by the multimode laser. Second harmonic conversion of the signal, with 3.8% efficiency, is used for the electronic excitation of the A-X (1,0) band of NO at∼<#comment/>215nm, while the residual signal at 430 nm is used for direct excitation of the A-X (0,0) band of the CH radical and elastic Rayleigh scattering. The section of the OPO signal wavelength for simultaneous CH and NO PLIF imaging is performed with consideration of the pulse energy, interference from the reactant and product species, and the fluorescence signal intensity. The excitation wavelengths of 430.7 nm and 215.35 nm are studied in a laminar, premixedCH4−<#comment/>H2−<#comment/>NH3–air flame. Single-shot CH and NO PLIF and Rayleigh scatter imaging is demonstrated in a turbulentCH4−<#comment/>H2−<#comment/>NH3diffusion flame using a high-speed intensified CMOS camera. Analysis of the complementary Rayleigh scattering and CH and NO PLIF enables identification and quantification of the high-temperature flame layers, the combustion product zones, and the fuel-jet core. Considerations for extension to simultaneous, 10-kHz-rate acquisition are discussed.

     
    more » « less
  5. Taylor and Francis (Ed.)
    A new computational methodology, termed ‘PeleLM-FDF’ is developed and utilised for high fidelity large eddy simulation (LES) of complex turbulent combustion systems. This methodology is constructed via a hybrid scheme combining the Eulerian PeleLM base flow solver with the Lagrangian Monte Carlo simulator of the filtered density func- tion (FDF) for the subgrid scale reactive scalars. The resulting methodology is capable of simulating some of the most intricate physics of complex turbulence-combustion interactions. This is demonstrated by LES of a non-premixed CO/H2 temporally evolv- ing jet flame. The chemistry is modelled via a skeletal kinetics model, and the results are appraised via a posteriori comparisons against direct numerical simulation (DNS) data of the same flame. Excellent agreements are observed for the time evolution of various statistics of the thermo-chemical quantities, including the manifolds of the multi-scalar mixing. The new methodology is capable of capturing the complex phe- nomena of flame-extinction and re-ignition at a 1/512 of the computational cost of the DNS. The high fidelity and the computational affordability of the new PeleLM-FDF solver warrants its consideration for LES of practical turbulent combustion systems. 
    more » « less