On velocity and reactive scalar spectra in turbulent premixed flames

AbstractKinetic energy and reactive scalar spectra in turbulent premixed flames are studied from compressible three-dimensional direct numerical simulations (DNS) of a temporally evolving rectangular slot-jet premixed flame, a statistically one-dimensional configuration. The flames correspond to a lean premixed hydrogen–air mixture at an equivalence ratio of 0.7, preheated to 700 K and at 1 atm, and three DNS are considered with a fixed jet Reynolds number of 10 000 and a jet Damköhler number varying between 0.13 and 0.54. For the study of spectra, motivated by the need to account for density change, which can be locally strong in premixed flames, a new density-weighted definition for two-point velocity/scalar correlations is proposed. The density-weighted two-point correlation tensor retains the essential properties of its constant-density (incompressible) counterpart and recovers the density-weighted Reynolds stress tensor in the limit of zero separation. The density weighting also allows the derivation of balance equations for velocity and scalar spectrum functions in the wavenumber space that illuminate physics unique to combusting flows. Pressure–dilatation correlation is a source of kinetic energy at high wavenumbers and, analogously, reaction rate–scalar fluctuation correlation is a high-wavenumber source of scalar energy. These results are verified by the spectra constructed from the DNS data. The kinetic energy spectra show a distinct inertial range with a $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}-5/3$ scaling followed by a ‘diffusive–reactive’ range at higher wavenumbers. The exponential drop-off in this range shows a distinct inflection in the vicinity of the wavenumber corresponding to a laminar flame thickness, $\delta _L$, and this is attributed to the contribution from the pressure–dilatation term in the energy balance in wavenumber space. Likewise, a clear spike in spectra of major reactant species (hydrogen) arising from the reaction-rate term is observed at wavenumbers close to $\delta _L$. It appears that in the inertial range classical scaling laws for the spectra involving the Kolmogorov scale are applicable, but in the high-wavenumber range where chemical reactions have a strong signature the laminar flame thickness produces a better collapse. It is suggested that a full scaling should perhaps involve the Kolmogorov scale, laminar flame thickness, Damköhler number and Karlovitz number.This is the accepted manuscript. The final version is available from CUP at http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=9317552&fileId=S0022112014003929

Tabulation of turbulent burning rates via a combination of a stochastic mixing model and tabulated chemistry

In this paper we propose a technique for a priori turbulent flame speed tabulation (TFST) for a given parameter space in standard combustion- regime diagrams. It can be used as a subgrid-scale (SGS) model in Large Eddy Simulation (LES). In a first step, stationary laminar flamelets are computed and stored over the progress variable following the ideas of flamelet generated manifolds (FGM). In a second step, the incompressible one-dimensional Navier-Stokes equations supplemented by the equation for the progress variable are solved on a grid that resolves all turbulent scales. Additionally, turbulent transport is implemented via the linear eddy model (LEM). The turbulent flame structures are solved until a statistically stationary mean value of the turbulent flame speed has been reached. The results are stored in a table that could be used by large scale premixed combustion models, e.g. front tracking schemes. Results are compared to an algebraic model and to direct numerical simulations (DNS)

A subgrid-scale model for large-eddy simulation of liquid/gas interfaces based on one-dimensional turbulence

The interface/turbulence interaction between two fluids in a turbulent environment has an important role in many technical processes, e.g. primary liquid atomization in combustion devices. Primary atomization has a significant role in spray formation and its characteristics. The resulting dynamics typically span 4–6 orders of magnitude in length scales, making detailed numerical simulations exceedingly expensive. This motivates the need for modeling approaches based on spatial filtering such as large-eddy simulation (LES). In this paper, a new approach based on One-Dimensional turbulence (ODT) is presented to describe the subgrid interface dynamics. ODT is a stochastic model simulating turbulent flow evolution along a notional one-dimensional line of sight by applying instantaneous maps that represent the effects of individual turbulent eddies on property fields. It provides affordable high resolution of interface creation and property gradients within each phase, which are key for capturing the local behavior as well as overall trends. ODT has previously been shown to reproduce the main features of an experimentally determined regime diagram for primary jet breakup. Here a new approach called VODT is presented which produces a size-conditioned as well as a total time rate of generation of droplets for given flow conditions at an interface. At the LES level, the total droplet generation from VODT is interpreted as a rate of mass conversion of LES-resolved liquid into unresolved droplets. Preliminary results of applying VODT to a cell with a planar-shear-layer are discussed at the end of the paper