柴油/空气湍流扩散燃烧的一个小火焰模型

本文将小火焰（flamelet）理论应用于分析柴油/空气湍流扩散燃烧的小火焰结构,以正十二烷同空气的一步反应为基础,建立柴油机燃烧的Flamelet模型,利用数值方法求出了柴油机湍流扩散燃烧的Flamelet结构.并采用假定PDF的方法,选取截尾式高斯分布的概率密度分布函数,将其与Flamelet结构相结合,求得燃烧过程中各参数的时均值,分析得出湍流脉动和非平衡作用对燃烧过程的影响.     

LES of a premixed jet flame DNS using a strained flamelet model

Turbulent premixed flames in the thin and broken reaction zones regimes are difficult to model with Large Eddy Simulation (LES) because turbulence strongly perturbs subfilter scale flame structures. This study addresses the difficulty by proposing a strained flamelet model for LES of high Karlovitz number flames. The proposed model extends a previously developed premixed flamelet approach to account for turbulence’s perturbation of subfilter premixed flame structures. The model describes combustion processes by solving strained premixed flamelets, tabulating the results in terms of a progress variable and a hydrogen radical, and invoking a presumed PDF framework to account for subfilter physics. The model is validated using two dimensional laminar flame studies, and is then tested by performing an LES of a premixed slot-jet direct numerical simulation (DNS). In the premixed regime diagram this slot-jet is found at the edge of the broken reaction zones regime. Comparisons of the DNS, the strained flamelet model LES, and an unstrained flamelet model LES confirm that turbulence perturbs flame structure to leading order effect, and that the use of an unstrained flamelet LES model under-predicts flame height. It is shown that the strained flamelet model captures the physics characterizing interactions of mixing and chemistry in highly turbulent regimes.     

Numerical analysis of gaseous hydrogen/liquid oxygen flamelet at supercritical pressures

Supercritical conditions are typically encountered in high-pressure combustion devices such as liquid propellant rockets and gas turbine engines. Significant real fluid behaviors including steep property variations occur when the fluid mixtures pass through the thermodynamic transcritical regime. The laminar flamelet concept is a robust and reliable approach that correctly accounts for real fluid effects, the large variation in thermophysical properties, and the detailed chemical kinetics for turbulent flames at transcritical and supercritical conditions. In the present study, the flamelet equations in the mixture fraction space are extended to treat the flame field of general fluids over transcritical and supercritical states. Flamelet computations are carried out for gaseous hydrogen and cryogenic liquid oxygen flames under a wide range of thermodynamic conditions. Based on numerical results, the detailed discussions are made for the effects of real fluid, pressure, and differential diffusion on the local flame structure and the characteristics encountered in liquid propellant rocket engines.     

Characteristics of flamelets in spray flames formed in a laminar counterflow

A two-dimensional numerical simulation of a spray flame formed in a laminar counterflow is presented, and the flamelet characteristics are studied in detail. The effects of strain rate, equivalence ratio, and droplet size are examined in terms of mixture fraction and scalar dissipation rate. n-Decane (C₁₀H₂₂) is used as a liquid spray fuel, and the droplet motion is calculated by the Lagrangian method without the parcel model. A one-step global reaction is employed for the combustion reaction model. The results show that there appear large differences in the trends of gaseous temperature and mass fractions of chemical species in the mixture fraction space between the spray flame and the gaseous diffusion flame. The gas temperature in the spray flame is much higher than that in the gaseous diffusion flame. This is due to the much lower scalar dissipation rate and the coexistence of premixed and diffusion-limited combustion in the spray flame. For the spray flames, gas temperature and mass fractions of chemical species are not unique functions of the mixture fraction scalar dissipation rate. This is because the production rate of the mixture fraction, namely evaporation rate of the droplets, in the upstream region is not in proportion to its transport-diffusion rate in the downstream region. The behavior shows marked differences as the strain rate decreases, the equivalence ratio increases, or the droplet size decreases.     

湍流扩散火焰中对基于混合分数的湍流燃烧模型的数值研究

根据条件矩模型(CMC)和小火焰面模型在模型构建上的相似,针对具有不同大小雷诺数和湍流-化学相互作用特性的非预混湍流射流火焰,对这两种模型进行了数值研究和比较.湍流燃烧模型采用Lagrangian型非稳态小火焰模型(LFM)和径向加权积分的CMC模型,而在H2/N2火焰的数值研究中还考虑了稳态小火焰模型的数值模拟结果....     

A filtered tabulated chemistry model for LES of premixed combustion

A new modeling strategy called F-TACLES (Filtered Tabulated Chemistry for Large Eddy Simulation) is developed to introduce tabulated chemistry methods in Large Eddy Simulation (LES) of turbulent premixed combustion. The objective is to recover the correct laminar flame propagation speed of the filtered flame front when subgrid scale turbulence vanishes as LES should tend toward Direct Numerical Simulation (DNS). The filtered flame structure is mapped using 1-D filtered laminar premixed flames. Closure of the filtered progress variable and the energy balance equations are carefully addressed in a fully compressible formulation. The methodology is first applied to 1-D filtered laminar flames, showing the ability of the model to recover the laminar flame speed and the correct chemical structure when the flame wrinkling is completely resolved. The model is then extended to turbulent combustion regimes by including subgrid scale wrinkling effects in the flame front propagation. Finally, preliminary tests of LES in a 3-D turbulent premixed flame are performed.     

采用考虑详细化学反应机理的火焰面模型模拟湍流扩散火焰

采用详细的甲烷氧化化学反应动力学机理(GRI-Mech3．0)对不同拉伸率条件下的拉伸层流扩散火焰面结构进行了数值计算，建立了一个包含一系列拉伸层流火焰面结构的火焰面数据库．将这些层流火焰面结构和美国Sandia国家实验室测得的湍流扩散火焰(FlameD)的平均火焰结构进行了对比，发现层流火焰面所覆盖的范围基本包含了所考虑的湍流火焰中不同位置的平均火焰结构，这表明火焰面模型是合理的．然后，采用火焰面模型对该湍流扩散火焰进行了数值模拟并和实验数据进行了比较，考察了火焰面模型的精确程度和模拟深度.     

An a priori model for the effective species Lewis numbers in premixed turbulent flames

An a priori model for the effective species Lewis numbers in premixed turbulent flames is presented. This a priori   analysis is performed using data from a series of direct numerical simulations (DNS) of lean (?=0.4$?=0.4$) premixed turbulent hydrogen flames, with Karlovitz number ranging from 10 to 1562 (Aspden et al., 2011). The conditional mean profiles of various species mass fraction versus temperature are evaluated from the DNS and compared to unstretched laminar premixed flame profiles. The turbulent flame structure is found to be different from the laminar flame structure. However, the turbulent flame can still be mapped onto a laminar flame with an appropriate change in the Lewis numbers of the different species. A transition from “laminar” Lewis numbers to unity Lewis numbers as the Karlovitz number increases is clearly captured. A model for those effective Lewis numbers with respect to the turbulent Reynolds number is developed. This model is derived from a Reynolds-averaged Navier–Stokes (RANS) formulation of the reactive scalar and temperature balance equations. The dependency of the effective Lewis numbers on the Karlovitz number instead of the Reynolds number is discussed in this paper. Unfortunately, given that the ratio of the integral length to the laminar flame thickness is fixed throughout this series of DNS, a change in the Karlovitz number is equivalent to a change in the Reynolds number. Incorporating these effective Lewis numbers in simulations of turbulent flames would have several impacts. First, the associated laminar flame speed and laminar flame thickness vary by a factor of two through the range of obtained effective Lewis numbers. Second, the turbulent premixed combustion regime diagram changes because a unique pair of laminar flame speed and laminar flame thickness cannot be used, and a dependency on the effective Lewis numbers has to be introduced. Finally, a turbulent flame speed model that takes into account these effective Lewis numbers is proposed.     

Evaluation of the laminar diffusion flamelet model in the calculation of an axisymmetric coflow laminar ethylene-air diffusion flame

A numerical study of an axisymmetric coflow laminar ethylene-air diffusion flame at atmospheric pressure was conducted using detailed chemistry and complex thermal and transport properties and two different methodologies: (1) the direct simulation method of solving the two-dimensional axisymmetric elliptic governing equations, and (2) the steady-state stretched diffusion flamelet model. Soot formation and radiative heat transfer were not taken into account in these calculations, both for simplicity and to avoid the complications associated with the issues of how to incorporate these chemical and physical processes into the flamelet model. The same reaction mechanism and thermal and transport properties were used in the 2D direct simulation and the generation of the flamelet library. The flamelet library was generated from the solutions of counterflow ethylene-air diffusion flames at a series of stretch rates. Results of the 2D direct simulation and the flamelet model are compared in physical space. Although the overall results of the flamelet model are qualitatively similar to those of the direct simulation, significant differences exist between the results of the two methods even for temperature and major species. The direct simulation method predicts that the peak concentrations of CO₂ and H₂O occur in different regions in the flame, while the flamelet model results show that the peak concentrations of CO₂ and H₂O occur in the same region. The flamelet model predicts an overly rapid approach to the equilibrium structure in the downstream region, leading to significantly higher flame temperatures. The main reason for the failure of the flamelet model in the downstream region is due to the neglect of the effects of multidimensional convection and diffusion and the fundamental difference in the chemical structure between a coflow diffusion flame and a counterflow diffusion flame. The findings of this paper are highly relevant to understanding the flamelet model results in the calculations of multidimensional turbulent diffusion flames.     

A filtered-laminar-flame PDF sub-grid scale closure for LES of premixed turbulent flames. Part I: Formalism and application to a bluff-body burner with differential diffusion

A sub-grid scale closure for Large Eddy Simulation (LES) of turbulent combustion based on physical-space filtering of laminar flames is discussed. Applied to an unstructured grid, the combustion LES filter size is not fixed in this novel approach devoted to LES with refined meshes, but calibrated depending on the local level of unresolved scalar fluctuations. The context is premixed or stratified flames, the derived model relies on four balance equations for mixture fraction and its variance, and a progress variable and its variance. The proposed formalism is based on a presumed probability density function (PDF) derived from the filtered flames. Closures for the terms of the equations that are unresolved over LES grids are achieved through the PDF. The method uses flamelet tabulated detailed chemistry and is first applied to the simulation of laminar flames (1D and 2D) over various grids for validation, before simulating a turbulent burner studied experimentally by Sweeney et al. (2012). Since this burner also features differential diffusion effects, the numerical model is modified to account for accumulation of carbon in the recirculation zone behind the bluff-body. A differential diffusion number based on the gradient of residence times is proposed, in an attempt to globally quantify differential diffusion effects in burners.     

Strained flamelets for turbulent premixed flames II: Laboratory flame results

The predictive ability of strained flamelets model for turbulent premixed flames is assessed using Reynolds Averaged Navier Stokes (RANS) calculations of laboratory flames covering a wide range of conditions. Reactant-to-product (RtP) opposed flow laminar flames parametrised using the scalar dissipation rate of reaction progress variable are used as strained flamelets. Two turbulent flames: a rod stabilised V-flame studied by Robin et al. [Combust. Flame 153 (2008) 288-315] and a set of pilot stabilised Bunsen flames studied by Chen et al. [Combust. Flame 107 (1996) 223-244] are calculated using a single set of model parameters. The V-flame corresponds to the corrugated flamelets regime. The strained flamelet model and an unstrained flamelet model yield similar predictions which are in good agreement with experimental measurements for this flame. On the other hand, for the Bunsen flames which are in the thin reaction zones regime, the unstrained flamelet model predicts a smaller flame brush compared to experiment. The predictions of the strained flamelets model allowing for fluid-dynamics stretch induced attenuation of the chemical reaction are in good agreement with the experimental data. This model predictions of major and minor species are also in good agreement with experimental data. The results demonstrate that the strained flamelets model using the scalar dissipation rate can be used across the combustion regimes.     

Three-dimensional temperature gradients in premixed turbulent flamelets via crossed-plane Rayleigh imaging

A new technique for obtaining instantaneous, high-resolution, three-dimensional thermal structure data from turbulent flames, crossed-plane Rayleigh imaging is described and then demonstrated. Quantitative Rayleigh imaging measurements are made simultaneously in two orthogonal, intersecting laser-sheet illumination planes. At points along the line of intersection of the two laser sheets, instantaneous, three-dimensional temperature gradient data are measured. The technique has higher resolution than parallel plane measurement techniques, which have limited resolution in the direction orthogonal to the parallel planes. The technique is used to measure temperature gradient data for a lean, premixed, methane-air turbulent V-flame with an equivalence ratio of Φ=0.7, and normalized turbulence intensity (, where u^′ is the turbulence intensity and is the unstretched laminar flame speed. Measurements are also presented for a laminar V-flame and a laminar Bunsen flame for comparison. Finally, an unstretched laminar flame calculation is made. Quantitative estimates of the experimental uncertainty are presented. The primary source of uncertainty in the data is due to shot noise. Measured temperature gradient data for laminar flames differ from that of the unstretched laminar flame calculation, especially in the oxidation layer. Turbulent flame temperature gradient data indicate that the turbulent V-flame thermal structure is not significantly perturbed from the measured laminar V-flame structure. For the flame studied, the flamelet approximation is valid if the flamelet used is based on the measured laminar flame structure of the V-flame. Isothermal surface orientation data are presented and are close to parallel for most realizations. Isothermal surface density is calculated from the distribution of isothermal surface orientations and from conditional averages of the magnitude of the temperature gradient. Isothermal surface density does not vary for different isothermal surfaces.     

Modeling curvature effects in diffusion flames using a laminar flamelet model

The goal of this paper is to investigate the effects of curvature of mixture fraction iso-surfaces on the transport of species in diffusion flames. A general flamelet formulation is derived mathematically considering both curvature effects and differential diffusion effects. These theoretical results suggest that curvature does not play a role in the transport process irrespective of the flame curvature if species transport is described with a unity Lewis number. On the other hand, a curvature-induced term becomes explicit when differential diffusion effects are considered, and it acts as a convective term in mixture fraction space. It is found that this term needs to be taken into account when the radius of curvature is comparable or smaller than the local flame thickness. For the proper integration of the flamelet equations, the scalar dissipation rate and curvature dependences on mixture fraction are modeled by considering two basic curved one-dimensional flame configurations. The flamelet equations accounting for curvature effects are solved with various prescribed curvature values. Results indicate that the mass fraction profiles of species with very small or large Lewis numbers are shifted significantly in mixture fraction space by the inclusion of curvature. Differential diffusion effects are enhanced by negative curvature values and suppressed by positive curvature values. In cases where flame curvature is not uniform, the curvature-induced convective term generates gradients along mixture fraction iso-surfaces, which enhance tangential diffusion effects. Budget analysis is performed on an axisymmetric laminar coflow diffusion flame to highlight the importance of the curvature-induced convective term compared to other terms in the full flamelet equation. A comparison is made between full chemistry simulation results and those obtained using planar and curved flamelet-based chemistry tabulation methods.     

DNS of spark ignition and edge flame propagation in turbulent droplet-laden mixing layers

A parametric study of forced ignition at the mixing layer between air and air carrying fine monosized fuel droplets is done through one-step chemistry direct numerical simulations to determine the influence of the size and volatility of the droplets, the spark location, the droplet-air mixing layer initial thickness and the turbulence intensity on the ignition success and the subsequent flame propagation. The propagation is analyzed in terms of edge flame displacement speed, which has not been studied before for turbulent edge spray flames. Spark ignition successfully resulted in a tribrachial flame if enough fuel vapour was available at the spark location, which occurred when the local droplet number density was high. Ignition was achieved even when the spark was offset from the spray, on the air side, due to the diffusion of heat from the spark, provided droplets evaporated rapidly. Large kernels were obtained by sparking close to the spray, since fuel was more readily available. At long times after the spark, for all flames studied, the probability density function of the displacement speed was wide, with a mean value in the range 0.55-0.75S_L, with S_L the laminar burning velocity of a stoichiometric gaseous premixed flame. This value is close to the mean displacement speed in turbulent edge flames with gaseous fuel. The displacement speed was negatively correlated with curvature. The detrimental effect of curvature was attenuated with a large initial kernel and by increasing the thickness of the mixing layer. The mixing layer was thicker when evaporation was slow and the turbulence intensity higher. However, high turbulence intensity also distorted the kernel which could lead to high values of curvature. The edge flame reaction component increased when the maximum temperature coincided with the stoichiometric contour. The results are consistent with the limited available experimental evidence and provide insights into the processes associated with ignition of practical spray flames.     

Hydroxyl time-series measurements and simulations for turbulent premixed jet flames in the thickened preheat regime

Quantitative measurements of OH concentration time series are presented for turbulent lean-premixed, methane-air jet flames theoretically in the thickened preheat regime. Picosecond time-resolved laser-induced fluorescence (PITLIF) reveals unique differences between these premixed flames and previous non-premixed jet flames. Time-averaged [OH] measurements are used to identify mean flame structures and to discern how these structures are affected by varying bulk flow velocities and heat release. More importantly, hydroxyl time series are inspected to distinguish among three main regions in these turbulent premixed flames. These regions include the reacting side of the flame brush, the mixing side of the flame brush (radially outside the location of heat release), and above the flame tip. Although the main reaction zone appears to be broadened by its associated high turbulent intensity, a combination of statistical analysis plus flamelet simulations suggests that the primary internal structure responsible for the OH distribution remains constant across the mean flame brush. Therefore, the absolute concentration of OH depends principally on the intermittency of this instantaneous internal structure. Outside the mean flame brush, mixing of OH with co-flow air shifts the distribution of absolute OH concentrations. Distinct autocorrelation functions are found within the three different regions identified for these premixed flames. Across the flame brush, integral time scales are dominated by turbulent convection, as verified by flamelet simulations. Above the flame tip, integral time scales are determined by a competition between turbulent convection and the reaction rate for OH destruction.     

A DNS study of hydrogen/air swirling premixed flames with different equivalence ratios

Hydrogen/air swirling premixed flames with different equivalence ratios are studied using direct numerical simulation. A fourth-order explicit Runge–Kutta method for time integration and an eighth-order central differencing scheme for spatial discretization are used to solve the full Navier–Stokes (N–S) equation system. A 9 species 19-step reduced mechanism for hydrogen/air combustion is adopted. The flames are stabilized with the help of a recirculation zone characterizing a high swirling flow. The vortex structures of the swirling premixed flames are presented. The flame structures are investigated in terms of the flame front curvature and tangential strain rate probability density functions (pdfs). The local flamelet temperature profiles are also extracted randomly along the flame front and compared with the corresponding laminar flame temperature profile. In order to study preferential diffusion effects, direct numerical simulation of two additional freely propagating planar flames in isotropic turbulence is conducted. Preferential diffusion effects observed in the planar flames are suppressed in the swirling flames. Further analysis confirms that the coherent small-scale eddies play important roles in the interactions between turbulence and the flame front. They are able to change the dynamic properties of the flame font and lead to enhanced burning intensity in the flame front with negative curvature for both stoichiometric and fuel-lean flames.     

Effects of radiation on spray flame characteristics and soot formation

Two-dimensional numerical simulations are applied to spray flames formed in a laminar counterflow and the effects of radiation on spray flame characteristics and soot formation are studied. N-Decane (C₁₀H₂₂) is used as the liquid fuel, and the droplet motion is calculated by the Lagrangian method. A single-step global reaction is employed for the combustion reaction model. A kinetically based soot model with a flamelet model is used to predict soot formation. Radiation is taken into account using the discrete ordinate method. The results show that radiation strongly affects the spray flame behavior and soot formation. Without the radiation model, flame temperature and soot volume fraction are greatly overestimated. The soot is formed in the diffusion flame regime, and its radiation emission increases with the increase in the equivalence ratio of the droplet fuel. This trend is in good agreement with that of the luminous flame behavior observed in the experiments.     

Modeling ignition and chemical structure of partially premixed turbulent flames using tabulated chemistry

Correctly reproducing the autoignition and the chemical composition of partially premixed turbulent flames is a challenge for numerical simulations of industrial applications such as diesel engines. A new model DF-PCM (diffusion flame presumed conditional moment) is proposed based on a coupling between the FPI (flame prolongation of ILDM) tabulation method and the PCM (presumed conditional moment) approach. Because the flamelets used to build the table are laminar diffusion flames, DF-PCM cannot be used for industrial applications like Diesel engines due to excessive CPU requirements. Therefore two new models called AI-PCM (autoignition presumed conditional moment) and ADF-PCM (approximated diffusion flames presumed conditional moment) are developed to approximate it. These models differ from DF-PCM because the flamelet libraries used for the table rely on PSR calculations. Comparisons between DF-PCM, AI-PCM, and ADF-PCM are performed for two fuels, n-heptane, representative of diesel fuels, and methane, which does not exhibit a “cool flame” ignition regime. These comparisons show that laminar diffusion flames can be approximated by flamelets based on PSR calculations in terms of autoignition delays and steady state profiles of the progress variable. Moreover, the evolution of the mean progress variable of DF-PCM can be correctly estimated by the approximated models. However, as discussed in this paper, errors are larger for CO and CO₂ mass fractions evolutions. Finally, an improvement to ADF-PCM, taking into account ignition delays, is proposed to better reproduce the ignition of very rich mixtures.