A Eulerian PDF scheme for LES of nonpremixed turbulent combustion with second-order accurate mixture fraction

Monte Carlo simulations of joint probability density function (PDF) approaches have been developed in the past largely with Reynolds averaged Navier Stokes (RANS) applications. Current interests are in the extension of PDF approaches to large eddy simulation (LES). As LES resolves accurately the large scales of turbulence in time, the Monte Carlo simulation and the flow field need to be tightly coupled. A tight coupling can be achieved if the consistency between the scalar field solution obtained via finite-volume (FV) methods and that from the stochastic solution of the PDF is ensured. For nonpremixed turbulent flames with two distinct streams, the local reactive mixture is described by the mixture fraction. A Eulerian Monte Carlo method is developed to achieve a second-order accuracy in the instantaneous filtered mixture fraction that is consistent with the corresponding FV. The performances of the proposed scheme are extensively evaluated using a one-dimensional model. Then, the scheme is applied to two cases with LES. The first one is a non-reacting mixing flow of two different fluids. The second case is the Sandia piloted turbulent flame D with a steady state flamelet model. Both results confirm the consistency of the proposed method to the level of filtered mixture fraction.     

A LES/PDF simulator on block-structured meshes

A block-structured mesh large-eddy simulation (LES)/probability density function (PDF) simulator is developed within the OpenFOAM framework for computational modelling of complex turbulent reacting flows. The LES/PDF solver is a hybrid solution methodology consisting of (i) a finite-volume (FV) method for solving the filtered mass and momentum equations (LES solver), and (ii) a Lagrangian particle-based Monte Carlo algorithm (PDF solver) for solving the modelled transport equation of the filtered joint PDF of compositions. Both the LES and the PDF methods are developed and combined to form a hybrid LES/PDF simulator entirely within the OpenFOAM framework. The in situ adaptive tabulation method [S.B. Pope, Computationally efficient implementation of combustion chemistry using in situ adaptive tabulation, Combust. Theory Model. 1 (1997), pp. 41–63; L. Lu, S.R. Lantz, Z. Ren, and B.S. Pope, Computationally efficient implementation of combustion chemistry in parallel PDF calculations, J. Comput. Phys. 228 (2009), pp. 5490–5525] is incorporated into the new LES/PDF solver for efficient computations of combustion chemistry with detailed reaction kinetics. The method is designed to utilise a block-structured mesh and can readily be extended to unstructured grids. The three-stage velocity interpolation method of Zhang and Haworth [A general mass consistency algorithm for hybrid particle/finite-volume PDF methods, J. Comput. Phys. 194 (2004), pp. 156–193] is adapted to interpolate the LES velocity field onto particle locations accurately and to enforce the consistency between LES and PDF fields at the numerical solution level. The hybrid algorithm is fully parallelised using the conventional domain decomposition approach. A detailed examination of the effects of each stage and the overall performance of the velocity interpolation algorithm is performed. Accurate coupling of the LES and PDF solvers is demonstrated using the one-way coupling methodology. Then the fully two-way coupled LES/PDF solver is successfully applied to simulate the Sandia Flame-D, and a turbulent non-swirling premixed flame and a turbulent swirling stratified flame from the Cambridge turbulent stratified flame series [M.S. Sweeney, S. Hochgreb, M.J. Dunn, and R.S. Barlow, The structure of turbulent stratified and premixed methane/air flames I: Non-swirling flows, Combust. Flame 159 (2012), pp. 2896–2911; M.S. Sweeney, S. Hochgreb, M.J. Dunn, and R.S. Barlow, The structure of turbulent stratified and premixed methane/air flames II: Swirling flows, Combust. Flame 159 (2012), pp. 2912–2929]. It is found that the LES/PDF method is very robust and the results are in good agreement with the experimental data for both flames.     

Large eddy simulations of turbulent flames using the filtered density function model

Large eddy simulations (LES) of the Sandia/Sydney swirl burners (SM1 and SMA1) and the Sandia/Darmstadt piloted jet diffusion flame (Flame D) are performed. These flames are part of the database of turbulent reacting flows widely considered as benchmark test cases for validating turbulent-combustion models. In the simulations presented in this paper, the subgrid scale (SGS) closure model adopted for turbulence-chemistry interactions is based on the transport filtered density function (FDF) model. In the FDF model, the transport equation for the joint probability density function (PDF) of scalars is solved. The main advantage of this model is that the filtered reaction rates can be exactly computed. However, the density field, computed directly from the FDF solver and needed in the hydrodynamic equations, is noisy and causes numerical instability. Two numerical approaches that yield a smooth density field are examined. The two methods are based on transport equations for specific sensible enthalpy (h_s) and RT, where R is the gas constant and T is the temperature. Consistency of the two methods is assessed in a bluff-body configuration using Reynolds averaged Navier-Stokes (RANS) methodology in conjunction with the transported PDF method. It is observed that the h_s method is superior to the RT method. Both methods are used in LES of the SM1 burner. In the near-field region, the h_s method produces better predictions of temperature. However, in the far-field region, both methods show deviation from data. Simulations of the SMA1 burner and Flame D are also presented using the h_s method. Some deficiencies are seen in the predictions of the SMA1 burner that may be related to the simple chemical kinetics model and mixing model used in the simulations. Simulations of Flame D show good agreement with data. These results indicate that, while further improvements to the methodology are needed, the LES/FDF method has the potential to accurately predict complex turbulent flames.     

Investigation of lengthscales, scalar dissipation, and flame orientation in a piloted diffusion flame by LES

This work investigates the structure of a diffusion flame in terms of lengthscales, scalar dissipation, and flame orientation by using large eddy simulation. This has been performed for a turbulent, non-premixed, piloted methane/air jet flame (Flame D) at a Reynolds-number of 22,400. A steady flamelet model, which was represented by artificial neural networks, yields species mass fractions, density, and viscosity as a function of the mixture fraction. This will be shown to suffice to simulate such flames. To allow to examine scalar dissipation, a grid of 1.97 × 10⁶ nodes was applied that resolves more than 75% of the turbulent kinetic energy. The accuracy of the results is assessed by varying the grid-resolution and by comparison to experimental data by Barlow, Frank, Karpetis, Schneider (Sandia, Darmstadt), and others. The numerical procedure solves the filtered, incompressible transport equations for mass, momentum, and mixture fraction. For subgrid closure, an eddy viscosity/diffusivity approach is applied, relying on the dynamic Germano model. Artificial turbulent inflow velocities were generated to feature proper one- and two-point statistics. The results obtained for both the one- and two-point statistics were found in good agreement to the experimental data. The PDF of the flame orientation shows the tilting of the flame fronts towards the centerline. Finally, the steady flamelet approach was found to be sufficient for this type of flame unless slowly reacting species are of interest.     

湍流扩散火焰局部熄火和再燃现象的PDF模拟

对一个值班湍流CH4／O2／N2射流扩散火焰(Sandia Flame D)进行了数值模拟研究．所采用的数学物理模型包括双尺度的k—ε湍流模型，标量联合的概率密度函数(PDF)输运方程方法，甲烷氧化的ARM简化化学反应机理(包含16种组分，12步总包反应)和欧几里德最小生成树(EMST)小尺度混合模型．将计算结果和实验数据进行了比较，不仅对于平均量，对于标量的散点分布和条件概率密度分布也是如此．计算结果表明文中采用的模型不仅能够预测宏观的火焰结构，而且预测了湍流燃烧中复杂的局部熄火和再燃过程．     

Eulerian transported probability density function sub-filter model for large-eddy simulations of turbulent combustion

Reactive flow simulations using large-eddy simulations (LES) require modelling of sub-filter fluctuations. Although conserved scalars like mixture fraction can be represented using a beta-function, the reactive scalar probability density function (PDF) does not follow an universal shape. A one-point one-time joint composition PDF transport equation can be used to describe the evolution of the scalar PDF. The high-dimensional nature of this PDF transport equation requires the use of a statistical ensemble of notional particles and is directly coupled to the LES flow solver. However, the large grid sizes used in LES simulations will make such Lagrangian simulations computationally intractable. Here we propose the use of a Eulerian version of the transported-PDF scheme for simulating turbulent reactive flows. The direct quadrature method of moments (DQMOM) uses scalar-type equations with appropriate source terms to evolve the sub-filter PDF in terms of a finite number of delta-functions. Each delta-peak is characterized by a location and weight that are obtained from individual transport equations. To illustrate the feasibility of the scheme, we compare the model against a particle-based Lagrangian scheme and a presumed PDF model for the evolution of the mixture fraction PDF. All these models are applied to an experimental bluff-body flame and the simulated scalar and flow fields are compared with experimental data. The DQMOM model results show good agreement with the experimental data as well as the other sub-filter models used.     

Computationally-efficient and accurate particle PDF simulations of turbulent combustion using coupled pre-partitioned adaptive chemistry and tabulation

LES/PDF methods are known to provide accurate results for challenging turbulent combustion configurations with strong turbulence-chemistry interactions. These methods are generally applicable as they do not make any assumptions on the topology of the underlying flame structure. However, this added generality comes at an increased computational cost. To mitigate this added cost, the majority of the LES/PDF computations performed to date utilize reduced mechanisms. We recently presented a coupled pre-partitioned adaptive chemistry (PPAC) and tabulation (ISAT) methodology (Newale et al., Comb. Th. Mod., 2019), which retains the fidelity of the detailed mechanism, while keeping the computational cost affordable. This methodology was tested in a partially-stirred reactor configuration. In this work, we describe the developments required for a holistic integration of PPAC-ISAT with a LES/PDF framework. We examine the performance of this coupled methodology in two LES/PDF configurations of Sandia flame D. A smaller simulation domain is initially utilized to characterize the efficiency and accuracy of standalone PPAC and coupled PPAC-ISAT in detail. Then, the performance of PPAC-ISAT is examined in a full-scale LES/PDF simulation. We show that the coupled PPAC-ISAT LES/PDF captures the resolved mean and RMS profiles of temperature and major species mass fractions to within 2% and OH to within 5%, with a reduction in the average simulation wall clock time per time step of 39% over an ISAT implementation using the detailed mechanism.     

Experimental study of velocity-scalar filtered joint density function for LES of turbulent combustion

The velocity-scalar filtered joint density function (FJDF) used in large eddy simulation (LES) of turbulent combustion is experimentally studied. Measurements are made in the fully developed region of an axisymmetric turbulent jet using an array consisting of three X-wires and resistance-wire temperature sensors. Filtering in the cross-stream and streamwise directions is realized by using the array and by invoking Taylor’s hypothesis, respectively. The means of the FJDF conditional on the subgrid-scale (SGS) turbulent kinetic energy and the SGS scalar variance at a given location range from close to joint normal to bimodal with the peaks separated in both velocity and scalar spaces, which correspond to qualitatively different mixing regimes. For close to joint normal FJDFs, the SGS fields are well mixed. For bimodal FJDFs, the conditionally filtered scalar diffusion and dissipation strongly depend on the SGS velocity and scalar, consistent with a combination of diffusion layers and plane strain in the SGS fields, which is similar to the counter-flow model for laminar flamelets. The results suggest that in LES, both mixing regimes could potentially be modeled accurately. The velocity field affects the SGS variance and the filtered scalar dissipation rate primarily by changing the degree of nonequilibrium of the SGS scalar and the SGS time scale, respectively. This study further demonstrates the importance of including velocity in mixing models.     

Unstrained and strained flamelets for LES of premixed combustion

The unstrained and strained flamelet closures for filtered reaction rate in large eddy simulation (LES) of premixed flames are studied. The required sub-grid scale (SGS) PDF in these closures is presumed using the Beta function. The relative performances of these closures are assessed by comparing numerical results from large eddy simulations of piloted Bunsen flames of stoichiometric methane–air mixture with experimental measurements. The strained flamelets closure is observed to underestimate the burn rate and thus the reactive scalars mass fractions are under-predicted with an over-prediction of fuel mass fraction compared with the unstrained flamelet closure. The physical reasons for this relative behaviour are discussed. The results of unstrained flamelet closure compare well with experimental data. The SGS variance of the progress variable required for the presumed PDF is obtained by solving its transport equation. An order of magnitude analysis of this equation suggests that the commonly used algebraic model obtained by balancing source and sink in this transport equation does not hold. This algebraic model is shown to underestimate the SGS variance substantially and the implications of this variance model for the filtered reaction rate closures are highlighted.     

湍流扩散燃烧的数值研究-PDF方法和火焰面模型的性能比较

采用标量概率密度函数(PDF)方法、稳态和非稳态火焰面模型三种方法对一个值班湍流CH_4/O_2/N_2射流扩散火焰(Sandia Flame D)进行数值计算,以比较不同燃烧模型的性能。PDF方法通过计算反应标量的PDF输运方程来得到标量分布,而火焰面模型只求解单标量混合物分数的PDF方程,组分和温度分布通过火焰面方程的求解或者火焰面数据库的插值得到。计算结果和实验数据对比表明PDF方法计算结果最好但计算量相当大,稳态火焰面模型则反之。综合而言,非稳态火焰面模型的预测结果相对稳态模型有了非常大的改进,而计算量仍然容易接受,非常适合工程应用。     

Studies of the flow and turbulence fields in a turbulent pulsed jet flame using LES/PDF

A turbulent piloted jet flame subject to a rapid velocity pulse in its fuel jet inflow is proposed as a new benchmark case for the study of turbulent combustion models. In this work, we perform modelling studies of this turbulent pulsed jet flame and focus on the predictions of its flow and turbulence fields. An advanced modelling strategy combining the large eddy simulation (LES) and the probability density function (PDF) methods is employed to model the turbulent pulsed jet flame. Characteristics of the velocity measurements are analysed to produce a time-dependent inflow condition that can be fed into the simulations. The effect of the uncertainty in the inflow turbulence intensity is investigated and is found to be very small. A method of specifying the inflow turbulence boundary condition for the simulations of the pulsed jet flame is assessed. The strategies for validating LES of statistically transient flames are discussed, and a new framework is developed consisting of different averaging strategies and a bootstrap method for constructing confidence intervals. Parametric studies are performed to examine the sensitivity of the predictions of the flow and turbulence fields to model and numerical parameters. A direct comparison of the predicted and measured time series of the axial velocity demonstrates a satisfactory prediction of the flow and turbulence fields of the pulsed jet flame by the employed modelling methods.     

Effects of combined dimension reduction and tabulation on the simulations of a turbulent premixed flame using a large-eddy simulation/probability density function method

A turbulent lean-premixed propane–air flame stabilised by a triangular cylinder as a flame-holder is simulated to assess the accuracy and computational efficiency of combined dimension reduction and tabulation of chemistry. The computational condition matches the Volvo rig experiments. For the reactive simulation, the Lagrangian Large-Eddy Simulation/Probability Density Function (LES/PDF) formulation is used. A novel two-way coupling approach between LES and PDF is applied to obtain resolved density to reduce its statistical fluctuations. Composition mixing is evaluated by the modified Interaction-by-Exchange with the Mean (IEM) model. A baseline case uses In Situ Adaptive Tabulation (ISAT) to calculate chemical reactions efficiently. Its results demonstrate good agreement with the experimental measurements in turbulence statistics, temperature, and minor species mass fractions. For dimension reduction, 11 and 16 represented species are chosen and a variant of Rate Controlled Constrained Equilibrium (RCCE) is applied in conjunction with ISAT to each case. All the quantities in the comparison are indistinguishable from the baseline results using ISAT only. The combined use of RCCE/ISAT reduces the computational time for chemical reaction by more than 50%. However, for the current turbulent premixed flame, chemical reaction takes only a minor portion of the overall computational cost, in contrast to non-premixed flame simulations using LES/PDF, presumably due to the restricted manifold of purely premixed flame in the composition space. Instead, composition mixing is the major contributor to cost reduction since the mean-drift term, which is computationally expensive, is computed for the reduced representation. Overall, a reduction of more than 15% in the computational cost is obtained.     

Large-eddy simulation of MILD combustion using partially stirred reactor approach

Subgrid-scale (SGS) parameterization and method for calculating filtered reaction rate are critical components of an accurate large-eddy simulation (LES) of turbulent flames. In this study, we integrate gradient-type structural SGS models with a partially stirred reactor approach by using detailed chemical kinetics to simulate a turbulent methane/hydrogen jet flame under moderate or intense low-oxygen dilution (MILD) conditions. The study examines two oxygen dilution levels. The framework is assessed through a systematic and comprehensive comparison of temperature, and mass fractions of major and minor species with experimental data and other reference simulation results. Overall, the statistics of the combustion field show excellent agreement with measurements at different axial locations, and a significant improvement compared to some previous simulations. It suggests that the proposed nonlinear LES framework is able to accurately model MILD combustion with reasonable computational cost.     

A consistent LES/filtered-density function formulation for the simulation of turbulent flames with detailed chemistry

A hybrid large-Eddy simulation/filtered-density function (LES–FDF) methodology is formulated for simulating variable density turbulent reactive flows. An indirect feedback mechanism coupled with a consistency measure based on redundant density fields contained in the different solvers is used to construct a robust algorithm. Using this novel scheme, a partially premixed methane/air flame is simulated. To describe transport in composition space, a 16-species reduced chemistry mechanism is used along with the interaction-by-exchange with the mean (IEM) model. For the micro-mixing model, typically a constant ratio of scalar to mechanical time-scale is assumed. This parameter can have substantial variations and can strongly influence the combustion process. Here, a dynamic time-scale model is used to prescribe the mixing time-scale, which eliminates the time-scale ratio as a model constant. Two different flame configurations, namely, Sandia flames D and E are studied. Comparison of simulated radial profiles with experimental data show good agreement for both flames. The LES–FDF simulations accurately predict the increased extinction near the inlet and re-ignition further downstream. The conditional mean profiles show good agreement with experimental data for both flames.     

受限湍流射流扩散火焰的PDF模拟

采用k-ε双方程模型和概率密度函数(PDF)相结合的办法,研究受限条件下的湍流射流扩散火焰,着重考虑受限条件下固体壁面、压力梯度等因素对速度场和标量场求解的影响,并在此基础上对两个不同尺寸的受限湍流燃烧场进行计算,分别研究了受限湍流射流扩散火焰的流场结构、火焰结构和火焰形状.最后给出结果,定性分析,并得出结论.     

Premixed turbulent combustion modeling using tabulated detailed chemistry and PDF

Tabulated chemistry and presumed probability density function (PDF) approaches are combined to perform RANS modeling of premixed turbulent combustion. The chemistry is tabulated from premixed flamelets with three independent parameters: the equivalence ratio of the mixture, the progress of reaction, and the specific enthalpy, to account for heat losses at walls. Mean quantities are estimated from presumed PDFs. This approach is used to numerically predict a turbulent premixed flame diluted by hot burnt products at an equivalence ratio that differs from the main stream of reactants. The investigated flame, subjected to high velocity fluctuations, has a thickened-wrinkled structure. A recently proposed closure for scalar dissipation rate that includes an estimation of the coupling between flame wrinkling and micromixing is retained. Comparisons of simulations with experimental measurements of mean velocity, temperature, and reactants are performed.     

A multivariate quadrature based moment method for LES based modeling of supersonic combustion

The transported probability density function (PDF) approach is a powerful technique for large eddy simulation (LES) based modeling of scramjet combustors. In this approach, a high-dimensional transport equation for the joint composition-enthalpy PDF needs to be solved. Quadrature based approaches provide deterministic Eulerian methods for solving the joint-PDF transport equation. In this work, it is first demonstrated that the numerical errors associated with LES require special care in the development of PDF solution algorithms. The direct quadrature method of moments (DQMOM) is one quadrature-based approach developed for supersonic combustion modeling. This approach is shown to generate inconsistent evolution of the scalar moments. Further, gradient-based source terms that appear in the DQMOM transport equations are severely underpredicted in LES leading to artificial mixing of fuel and oxidizer. To overcome these numerical issues, a semi-discrete quadrature method of moments (SeQMOM) is formulated. The performance of the new technique is compared with the DQMOM approach in canonical flow configurations as well as a three-dimensional supersonic cavity stabilized flame configuration. The SeQMOM approach is shown to predict subfilter statistics accurately compared to the DQMOM approach.     

Direct evaluation of the subgrid scale scalar flux in turbulent premixed flames with conditioned dual-plane stereo PIV

With the dual-plane stereo PIV technique the instantaneous three-dimensional resolved rate-of-strain tensor is directly measured in turbulent premixed flames. Simultaneously, also the instantaneous subgrid scale (SGS) scalar flux is measured with fine resolution, where for the latter term the conditioned particle image velocimetry (CPIV) technique is applied. The subgrid resolution reaches 118 μm, allowing a 9 × 9 resolution of a subgrid filter with width Δ = 1 mm. This combined measurement approach allows the a-priori comparison of models for the SGS scalar flux term with direct measurements which is important for large eddy simulation methods in turbulent premixed flames. Two different flame conditions of a premixed V-shaped turbulent flame are investigated where the turbulence intensity is varied by a factor of nearly three. The instantaneous radial and axial SGS fluxes are compared with the following three models: gradient model with Smagorinsky approach for the turbulent viscosity, Clark model, and extended gradient model with an anisotropy term. None of these models shows a good correlation with the directly measured flux. The anisotropy term alone (being nearly similar to the Clark model) shows, however, a right trend behaviour. An analysis of the data indicates a significant dependency of the experimentally determined SGS flux on the Favre averaged reaction progress (spatially averaged over the SGS area). A relatively simple closure for the SGS flux, which describes the dilatation due to the gasdynamic expansion, and which is a function proportional to , shows a rather good correlation with direct measurement for some of the components. A successful SGS scalar flux model for premixed turbulent flames most likely needs to include at least two different effects.     

The effect of timescale variation in multiple mapping conditioning mixing of PDF calculations for Sandia Flame series D–F

A stochastic implementation of the multiple mapping conditioning (MMC) model has been used for the modelling of turbulence–chemistry interactions in a series of turbulent jet diffusion flames with varying degrees of local extinction (Sandia Flames D–F). The mapping function approximates the cumulative probability distribution of mixture fraction and the corresponding variance can be controlled by a standard implementation of the scalar mixing timescale. The conditional fluctuations are controlled by a minor dissipation timescale, τ_min. The results show a clear dependence of the conditional fluctuations on the choice of the minor timescale, and the appropriate value for turbulent jet flames is similar to values determined in related direct numerical simulation (DNS) studies of homogeneous turbulent reacting flows. The predictions of means and variances of temperature and species mass fractions are very good for all flames, indicating an appropriate modelling of the conditional variances. Further sensitivity studies with respect to particle number density demonstrate a relative insensitivity of the results to the particle number in the numerical solution procedure. Good results can be obtained with as few as 10 particles per cell, allowing for a computationally inexpensive implementation of a Monte Carlo/probability density function (PDF) method.