A mixing controlled direct chemistry (MCDC) model for diesel engine combustion modelling using large eddy simulation

A mixing controlled direct chemistry (MCDC) combustion model with sub-grid scale (SGS) mixing effects and chemical kinetics has been evaluated for Large Eddy Simulation (LES) of diesel engine combustion. The mixing effect is modelled by a mixing timescale based on mixture fraction variance and sub-grid scalar dissipation rate. The SGS scalar dissipation rate is modelled using a similarity term and a scaling factor from the analysis of Direct Numerical Simulation (DNS) data. The chemical reaction progress is estimated from a kinetic timescale based on local internal energy change rate and equilibrium state internal energy. An optical research engine operating at conventional operating conditions and Low Temperature Combustion (LTC) conditions was used for evaluation of the combustion model. From the simulation results, the effect of SGS scalar mixing is evaluated at different stages of combustion. In the context of LES, the new approach provides improved engine modelling results compared to the Direct Chemistry Solver (DCS) combustion model.     

Modelling of lifted turbulent diffusion flames in a channel mixing layer by the flame hole dynamics

The partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics as an effort to develop a prediction model for the turbulent flame lift off. The essence of the flame hole dynamics is derivation of the random walk mapping, from the flame-edge theory, which governs expansion or contraction of the quenching holes initially created by the local quenching events. The numerical simulation for the flame hole dynamics is carried out in two stages. First, a direct numerical simulation is performed for a constant-density fuel–air channel mixing layer to obtain the background turbulent flow and mixing fields, from which a time series of two-dimensional scalar-dissipation-rate array is extracted. Subsequently, a Lagrangian simulation of the flame hole random walk mapping, projected to the scalar dissipation rate array, yields a temporally evolving turbulent extinction process and its statistics on partial quenching characteristics. In particular, the probability of encountering the reacting state, while conditioned with the instantaneous scalar dissipation rate, is examined to reveal that the conditional probability has a sharp transition across the crossover scalar dissipation rate, at which the flame edge changes its direction of propagation. This statistical characteristic implies that the flame edge propagation instead of the local quenching event is the main mechanism controlling the partial quenching events in turbulent flames. In addition, the conditional probability can be approximated by a heavyside function across the crossover scalar dissipation rate.     

Conditional statistics of inert droplet effects on turbulent combustion in reacting mixing layers

Direct numerical simulation (DNS) of turbulent reacting mixing layers laden with evaporating inert droplets is used to assess the droplet effects in the context of the conditional moment closure (CMC) for multiphase turbulent combustion. The temporally developing mixing layer has an initial Reynolds number of 1000 based on the vorticity thickness with more than 16 million Lagrangian droplets traced. An equivalent mixture fraction incorporating the inert vapour mass fractions clearly demonstrates the effects of vapour dilution on the mixture. Instantaneous fields and conditional statistics, such as the singly conditioned scalar dissipation rate, the gas temperature 〈 T _g|η 〉, conditional variance of the gas temperature 〈 T _g ^”2|η 〉 and conditional covariance between the fuel mass fraction and gas temperature 〈 Y _f ^” T _g ^”|η 〉 show considerable droplet effects. Comparison between the droplet-free and droplet-laden reacting mixing layer cases suggests significant extinction in the latter case. The resulting large conditional fluctuations around the conditional means contradict the basic assumption behind the first-order singly conditioned CMC. More sophisticated CMC approaches, such as doubly conditioned or second-order CMCs are, in principle, better able to incorporate the effects of evaporating droplets, but significant modelling challenges exist. The scalar dissipation rate doubly conditioned on the mixture fraction and a normalized gas temperature 〈 χ | η, ζ 〉 exemplifies the modelling complexity in the CMC of multiphase combustion.     

Scalar dissipation rates in isothermal and reactive turbulent opposed-jets: 1-D-Raman/Rayleigh experiments supported by LES

Isothermal and reactive turbulent opposed flows are presented, which are appropriate to test the applicability and performance of models for turbulence, mixing, chemical reaction, and turbulence-chemistry interaction. Transient flow and scalar fields are measured using laser Doppler velocimetry and one-dimensionally resolved Raman/Rayleigh spectroscopy. Aside of statistical moments of temperature, mean species, and velocity components, scalar dissipation rate across the mixing and reaction layer is determined on a single-shot base. Using large eddy simulation in connection with a steady flamelet model, it is shown how numerical data can serve to estimate the influence of experimental noise upon a measured quantity, such as scalar dissipation. As a key result, it is shown that an increase in scalar rate of dissipation by chemical reactions is caused by a significant increase in the mixture fraction diffusivity, which outweighs the decrease in mixture fraction gradients. In mixture fraction space, local maxima of scalar dissipation rate are found on the rich side, which cannot be correctly reproduced by the steady flamelet model assuming equal species diffusivity. Furthermore, the impact of experimental noise on conditional probability density functions of scalar dissipation rate is shown (exemplary) to lead to erroneous conclusions from experimental data.     

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.     

Conditional dissipation of scalars in homogeneous turbulence: Closure for MMC modelling

While the mean and unconditional variance are to be predicted well by any reasonable turbulent combustion model, these are generally not sufficient for the accurate modelling of complex phenomena such as extinction/reignition. An additional criterion has been recently introduced: accurate modelling of the dissipation timescales associated with fluctuations of scalars about their conditional mean (conditional dissipation timescales). Analysis of Direct Numerical Simulation (DNS) results for a passive scalar shows that the conditional dissipation timescale is of the order of the integral timescale and smaller than the unconditional dissipation timescale. A model is proposed: the conditional dissipation timescale is proportional to the integral timescale. This model is used in Multiple Mapping Conditioning (MMC) modelling for a passive scalar case and a reactive scalar case, comparing to DNS results for both. The results show that this model improves the accuracy of MMC predictions so as to match the DNS results more closely using a relatively-coarse spatial resolution compared to other turbulent combustion models.     

A quantitative method for a priori evaluation of combustion reaction models

In recent years, direct numerical simulations have been used increasingly to evaluate the validity and performance of combustion reaction models. This study presents a new, quantitative method to determine the ideal model performance attainable by a given parameterization of the state variables. Data from direct numerical simulation (DNS) of unsteady CO/H₂–air jet flames is analysed to determine how well various parameterizations represent the data, and how well specific models based on those parameterizations perform. Results show that the equilibrium model performs poorly relative to an ideal model parameterized by the mixture fraction. The steady laminar flamelet model performs quite well relative to an ideal model parameterized by mixture fraction and dissipation rate in some cases. However, at low dissipation rates or at dissipation rates exceeding the steady extinction limit, the steady flamelet model performs poorly. Interestingly, even in many cases where the steady flamelet model fails (particularly at low dissipation rate), the DNS data suggests that the state may be parameterized well by the mixture fraction and dissipation rate. A progress variable based on the CO₂ mass fraction is proposed, together with a new model based on the CO₂ progress variable. This model performs nearly ideally, and demonstrates the ability to capture extinction with remarkable accuracy for the CO/H₂ flames considered.     

Turbulence/flame/wall interactions in non-premixed inclined slot-jet flames impinging at a wall using direct numerical simulation

In the present work, three-dimensional turbulent non-premixed oblique slot-jet flames impinging at a wall were investigated using direct numerical simulation (DNS). Two cases are considered with the Damköhler number (Da) of case A being twice that of case B. A 17 species and 73-step mechanism for methane combustion was employed in the simulations. It was found that flame extinction in case B is more prominent compared to case A. Reignition in the lower branch of combustion for case A occurs when the scalar dissipation rate relaxes, while no reignition occurs in the lower branch for case B due to excessive scalar dissipation rate. A method was proposed to identify the flame quenching edges of turbulent non-premixed flames in wall-bounded flows based on the intersections of mixture fraction and OH mass fraction iso-surfaces. The flame/wall interactions were examined in terms of the quenching distance and the wall heat flux along the quenching edges. There is essentially no flame/wall interaction in case B due to the extinction caused by excessive turbulent mixing. In contrast, significant interactions between flames and the wall are observed in case A. The quenching distance is found to be negatively correlated with wall heat flux as previously reported in turbulent premixed flames. The influence of chemical reactions and wall on flow topologies was identified. The FS/U and FC/U topologies are found near flame edges, and the NNN/U topology appears when reignition occurs. The vortex-dominant topologies, FC/U and FS/S, play an increasingly important role as the jet turbulence develops.     

Conditional velocity statistics in the double scalar mixing layer – A mapping closure approach

In this work we use 3D direct numerical simulations (DNS) to investigate the average velocity conditioned on a conserved scalar in a double scalar mixing layer (DSML). The DSML is a canonical multistream flow designed as a model problem for the extensively studied piloted diffusion flames. The conditional mean velocity appears as an unclosed term in advanced Eulerian models of turbulent non-premixed combustion, like the conditional moment closure and transported probability density function (PDF) methods. Here it accounts for inhomogeneous effects that have been found significant in flames with relatively low Damköhler numbers. Today there are only a few simple models available for the conditional mean velocity and these are discussed with reference to the DNS results. We find that both the linear model of Kutznetzov and the Li and Bilger model are unsuitable for multi stream flows, whereas the gradient diffusion model of Pope shows very close agreement with DNS over the whole range of the DSML. The gradient diffusion model relies on a model for the conserved scalar PDF and here we have used a presumed mapping function PDF, that is known to give an excellent representation of the DNS. A new model for the conditional mean velocity is suggested by arguing that the Gaussian reference field represents the velocity field, a statement that is evidenced by a near perfect agreement with DNS. The model still suffers from an inconsistency with the unconditional flux of conserved scalar variance, though, and a strategy for developing fully consistent models is suggested.     

Influence of advective heat flux on extinction scalar dissipation rate and velocity in negative edge flames

A counterflow flame geometry, which has previously been experimentally shown to produce stable negative edge flames, was studied using numerical simulations. In this geometry, the flame edge is formed off the counterflow centreline owing to a local increase in scalar dissipation rate. Hot products from the stable nonpremixed flame on the centreline flow through the edge at velocities of ～ 1–5 m/s. The size of the counterflow burner and the gas flowrates are varied in the simulations to alter the flame strength and velocity at the flame edge. The advection of products through the edge is shown to extend the flame extinction to higher scalar dissipation rates than required for centreline extinction. For high velocities, the scalar dissipation rate required for flame extinction can be related to the centreline extinction value by considering only the effect of energy addition to the flame edge via advection. However, for lower edge flame velocities, the effects of increased thermal and species diffusion through the edge must also be included. Since the advection at the edge is a product of both the local velocity and temperature gradient, a single correlation between the scalar dissipation rate and the negative edge flame velocity does not exist.     

A priori assessment of the subgrid scale viscous/scalar dissipation closures in compressible turbulence

An appraisal is made of several subgrid scale (SGS) viscous/scalar dissipation closures via a priori analysis of direct numerical simulation data in a temporally evolving compressible mixing layer. The effects of the filter width, the compressibility level and the Schmidt number are studied for several models. Based on the scaling of SGS kinetic energy, a new formulation for SGS viscous dissipation is proposed. This yields the best overall prediction of the SGS viscous dissipation within the inertial subrange. An SGS scalar dissipation model based on the proportionality of the turbulent time scale with the scalar mixing time scale also performs the best for the filter widths in the inertial subrange. Two dynamic methods are implemented for the determination of the model coefficients. The one based on the global equilibrium of dissipation and production is shown to be more satisfactory than the conventional dynamic model.     

A DNS study on the flame structures and flame stabilization mechanism of a laboratory-scale lean premixed jet flame in crossflow

In the present work, direct numerical simulation (DNS) of a laboratory-scale lean premixed reacting jet flame in crossflow was performed to understand the flame structures and the flame stabilization mechanism. In the DNS, an ethylene-air jet with an equivalence ratio of 0.6 was injected into a hot vitiated crossflow. The jet Reynolds number reaches 6161. The DNS results were compared with those of the experiment with a good agreement. It was found that the windward and leeward branches of the flame show significantly different behaviors. The windward flame branch, appearing lifted and discontinuous, is located in the shear layer regions with high temperature, low vorticity and low scalar dissipation rate. The location of the peak heat release rate shifts to a higher mixture fraction with increasing distance from the jet exit. The leeward branch of the flame anchors in the shear layer near the jet exit. The recirculation zone in the wake of the jet facilitates the stabilization of the leeward flame. The chemical explosive mode analysis (CEMA) and species budget analysis were employed to characterize the local combustion mode. Auto-ignition plays a key role in the stabilization of the windward flame where a large range of extinction is also found due to the high strain rate. In contrast, premixed flame propagation is dominant on the leeward side.     

A priori analysis of sub-grid variance of a reactive scalar using DNS data of high Ka flames

Direct numerical simulations (DNS) of low and high Karlovitz number (Ka) flames are analysed to investigate the behaviour of the reactive scalar sub-grid scale (SGS) variance in premixed combustion under a wide range of combustion conditions (regimes). An order of magnitude analysis is performed to assess the importance of various terms in the variance evolution equation and the analysis is validated using the DNS results. This analysis sheds light on the relative behaviour among turbulent transport and production, scalar dissipation and chemical processes involved in the evolution of the SGS variance at different Ka. The common expectation is that the variance equation shifts from a reaction-dissipation balance at low Ka to a production–dissipation balance at high Ka with diminishing reaction contribution. However, in large eddy simulation (LES), a high Ka alone does not make the reaction term negligible, as the relative importance of the reaction term has a concurrent increase with filter size. The filter size can be relatively large compared with the Kolmogorov length scale in practical LES of high Ka flames, and as a consequence a reaction–production–dissipation balance may prevail in the variance equation even in a high Ka configuration, and this possibility is quantified using the DNS analysis in this work. This has implications from modelling perspectives, and therefore two commonly used closures in LES for the SGS scalar dissipation rate are investigated a priori to estimate the importance of the above balance in LES modelling. The results are explained to highlight the interplay among turbulence, chemistry and dissipation processes as a function of Ka.     

聚变堆液态包层准二维磁流体湍流模型研究

为研究聚变堆氚增殖包层中液态金属湍流磁流体动力学(MHD)效应,开发了一种新的准二维单方程 MHD 湍流模型,并进行了相关数值模拟程序的编制及验证。对于矩形管道中的准二维 MHD 湍流流动,三维流 动主要发生在哈德曼层中,中心的主流区呈现出二维流动。为了反映这种特殊的流动特征,新湍流模型在标准 k-ε 模型的基础之上去掉了传统的耗散项,代之以电磁耗散项来模拟湍流 MHD 效应。同时,采用 Bradshaw 假设来对 湍流涡粘系数进行模化。为验证该湍流模型是否合理,编制了相关数值模拟程序,并利用直接数值模拟(DNS)结 果对该程序进行了校正,数值模拟结果与 DNS 结果吻合较好。计算结果表明,该湍流模型可应用于聚变堆液态 包层 MHD 湍流流动的数值模拟。     

Conditional Moment Closure (CMC) applied to autoignition of high pressure methane jets in a shock tube

Autoignition of non-premixed methane–air mixtures is investigated using first-order Conditional Moment Closure (CMC). Turbulent velocity and mixing fields simulations are decoupled from the CMC calculations due to low temperature changes until ignition occurs. The CMC equations are cross-stream averaged and finite differences are applied to discretize the equations. A three-step fractional method is implemented to treat separately the stiff chemical source term. Two detailed chemical kinetics mechanisms are tested as well as two mixing models. The present results show good agreement with published experimental measurements for the magnitude of both ignition delay and kernel location. The slope of the predicted ignition delay is overpredicted and possible sources of discrepancy are identified. Both scalar dissipation rate models produce comparable results due to the turbulent flow homogeneity assumption. Further, ignition always occurs at low scalar dissipation rates, much lower than the flamelet critical value of ignition. Ignition is found to take place in lean mixtures for a value of mixture fraction around 0.02. The conditional species concentrations are in qualitative agreement with previous research. Homogeneous and inhomogeneous CMC calculations are also performed in order to investigate the role of physical transport in the present autoignition study. It is found that spatial transport is small at ignition time. Predicted ignition delays are shown to be sensitive to the chemical kinetics. Reasonable agreement with previous simulations is found. Improved formulations for the mixing model based on non-homogeneous turbulence are expected to have an impact.     

Assessment of dynamic closure for premixed combustion large eddy simulation

Turbulent piloted Bunsen flames of stoichiometric methane–air mixtures are computed using the large eddy simulation (LES) paradigm involving an algebraic closure for the filtered reaction rate. This closure involves the filtered scalar dissipation rate of a reaction progress variable. The model for this dissipation rate involves a parameter β_c representing the flame front curvature effects induced by turbulence, chemical reactions, molecular dissipation, and their interactions at the sub-grid level, suggesting that this parameter may vary with filter width or be a scale-dependent. Thus, it would be ideal to evaluate this parameter dynamically by LES. A procedure for this evaluation is discussed and assessed using direct numerical simulation (DNS) data and LES calculations. The probability density functions of β_c obtained from the DNS and LES calculations are very similar when the turbulent Reynolds number is sufficiently large and when the filter width normalised by the laminar flame thermal thickness is larger than unity. Results obtained using a constant (static) value for this parameter are also used for comparative evaluation. Detailed discussion presented in this paper suggests that the dynamic procedure works well and physical insights and reasonings are provided to explain the observed behaviour.     

Subgrid-scale scalar flux modelling based on optimal estimation theory and machine-learning procedures

New procedures are explored for the development of models in the context of large eddy simulation (LES) of a passive scalar. They rely on the combination of the optimal estimator theory with machine-learning algorithms. The concept of optimal estimator allows to identify the most accurate set of parameters to be used when deriving a model. The model itself can then be defined by training an artificial neural network (ANN) on a database derived from the filtering of direct numerical simulation (DNS) results. This procedure leads to a subgrid scale model displaying good structural performance, which allows to perform LESs very close to the filtered DNS results. However, this first procedure does not control the functional performance so that the model can fail when the flow configuration differs from the training database. Another procedure is then proposed, where the model functional form is imposed and the ANN used only to define the model coefficients. The training step is a bi-objective optimisation in order to control both structural and functional performances. The model derived from this second procedure proves to be more robust. It also provides stable LESs for a turbulent plane jet flow configuration very far from the training database but over-estimates the mixing process in that case.     

Effects of body force on the statistical behaviour and modelling of scalar variance in turbulent premixed flames

The effects of body force/external pressure gradient on the statistical behaviours of the reaction progress variable variance and the terms of its transport equation have been investigated for different turbulence intensities using DNS data of statistically planar flames. Since the extent of flame wrinkling increases with the strengthening of body force promoting unstable stratification, the scalar variance has been found to decrease under strong body force promoting stability. This trend is particularly strong for low turbulence intensities where the probability density function of the reaction progress variable cannot be approximated by a bimodal distribution. Therefore, an algebraic relation for the reaction progress variable variance, derived based on a presumed bimodal probability density function of reaction progress variable, cannot be used for general flow conditions. The contributions of chemical reaction and scalar dissipation rates in the scalar variance transport equation remain leading order source and sink, respectively for all cases irrespective of the strength and direction of the body force. The counter-gradient type transport is found to weaken with increasing body force magnitude when the body force is directed from the heavier unburned gas to the lighter burned gas side of the flame brush, and vice versa. Although a scalar dissipation rate-based reaction rate closure can be utilised to model the reaction rate contribution to the scalar variance transport accurately, the dissipation rate contribution due to the gradient of the Favre-averaged reaction progress variable cannot be ignored and it plays a key role for large magnitudes of body force promoting stable stratification. An algebraic closure of the scalar dissipation rate, originally proposed for high Damköhler number combustion, has been modified for the thin reaction zones regime combustion by incorporating the effects of Froude number. This model has been shown to predict the scalar dissipation rate accurately for all cases considered here.     

Turbulent scalar flux transport in head-on quenching of turbulent premixed flames: a direct numerical simulations approach to assess models for Reynolds averaged Navier Stokes simulations

The statistical behaviour and the modelling of turbulent scalar flux transport have been analysed using a direct numerical simulation (DNS) database of head-on quenching of statistically planar turbulent premixed flames by an isothermal wall. A range of different values of Damköhler, Karlovitz numbers and Lewis numbers has been considered for this analysis. The magnitudes of the turbulent transport and mean velocity gradient terms in the turbulent scalar flux transport equation remain small in comparison to the pressure gradient, molecular dissipation and reaction-velocity fluctuation correlation terms in the turbulent scalar flux transport equation when the flame is away from the wall but the magnitudes of all these terms diminish and assume comparable values during flame quenching before vanishing altogether. It has been found that the existing models for the turbulent transport, pressure gradient, molecular dissipation and reaction-velocity fluctuation correlation terms in the turbulent scalar flux transport equation do not adequately address the respective behaviours extracted from DNS data in the near-wall region during flame quenching. Existing models for transport equation-based closures of turbulent scalar flux have been modified in such a manner that these models provide satisfactory prediction both near to and away from the wall.     

A priori evaluation of large eddy simulation subgrid-scale scalar flux models in isotropic passive-scalar and anisotropic buoyancy-driven homogeneous turbulence

Direct numerical simulation (DNS) of passive (non-buoyant) and active (buoyant) scalar homogeneous turbulence is carried out using a standard pseudo-spectral numerical method. The flow settings simulated include stationary forced and decaying passive-scalar turbulence, as well as decaying anisotropic active-scalar turbulence. The Schmidt number is unity in all cases. The results are compared with, and are found to be in very good agreement with, previous similar DNS studies. The well-validated DNS data are divided into 19 sets, and are employed to study different large eddy simulation (LES) subgrid-scale (SGS) models for the SGS scalar flux. The models examined include three eddy-viscosity-type models (Smagorinsky, Vreman and Sigma with a constant SGS Schmidt number), a Dynamic Structure model and two versions of the Gradient (Gradient and Modulated Gradient) model. The models are investigated with respect to their ability to predict the orientation, and the magnitude, of the SGS scalar flux. Eddy-viscosity models are found to predict the magnitude of the SGS scalar flux accurately, but are poor at predicting the orientation of the SGS scalar flux. The Dynamic Structure and Gradient models are better than eddy-viscosity models at predicting both the magnitude and direction. However, neither of them can be realised in an actual LES, without carrying additional transport equations. Based on these observations, four new models are proposed – combining directions from Dynamic Structure and Gradient models, and magnitudes from Smagorinsky and Vreman eddy-viscosity models. These models are expected to be better than eddy-viscosity and Modulated Gradient models, and this is confirmed by preliminary a posteriori tests.