Influence of turbulence–chemistry interaction for n-heptane spray combustion under diesel engine conditions with emphasis on soot formation and oxidation

The influence of the turbulence–chemistry interaction (TCI) for n-heptane sprays under diesel engine conditions has been investigated by means of computational fluid dynamics (CFD) simulations. The conditional moment closure approach, which has been previously validated thoroughly for such flows, and the homogeneous reactor (i.e. no turbulent combustion model) approach have been compared, in view of the recent resurgence of the latter approaches for diesel engine CFD. Experimental data available from a constant-volume combustion chamber have been used for model validation purposes for a broad range of conditions including variations in ambient oxygen (8?21% by vol.), ambient temperature (900 and 1000 K) and ambient density (14.8 and 30 kg/m³). The results from both numerical approaches have been compared to the experimental values of ignition delay (ID), flame lift-off length (LOL), and soot volume fraction distributions. TCI was found to have a weak influence on ignition delay for the conditions simulated, attributed to the low values of the scalar dissipation relative to the critical value above which auto-ignition does not occur. In contrast, the flame LOL was considerably affected, in particular at low oxygen concentrations. Quasi-steady soot formation was similar; however, pronounced differences in soot oxidation behaviour are reported. The differences were further emphasised for a case with short injection duration: in such conditions, TCI was found to play a major role concerning the soot oxidation behaviour because of the importance of soot-oxidiser structure in mixture fraction space. Neglecting TCI leads to a strong over-estimation of soot oxidation after the end of injection. The results suggest that for some engines, and for some phenomena, the neglect of turbulent fluctuations may lead to predictions of acceptable engineering accuracy, but that a proper turbulent combustion model is needed for more reliable results.     

Comparison of well-mixed and multiple representative interactive flamelet approaches for diesel spray combustion modelling

The application of detailed chemistry to the computational fluid dynamics simulation of combustion process in diesel engines has many potentials, including the possibility to predict auto-ignition, diffusion flame structure, stabilisation and soot formation in a wide range of operating conditions, also taking into account the effects of different fuel types. Among the approaches that were proposed over the years, the ones that are mostly used in practical calculations can be divided into two main categories: the first assumes each cell to be a well-stirred reactor, while the second employs the flamelet assumption to describe both auto-ignition and turbulent diffusion flame propagation. Despite the fact that both types of model have been widely validated over the years, a detailed comparison between them appears to be very useful in order to understand better the relevant parameters governing auto-ignition, flame stabilisation and the formation of pollutant emissions. This work is focused on a comparison of two different combustion models that were recently implemented by the authors in an open-source code. The first assumes each cell to be a homogeneous reactor and neglects interaction between turbulence and chemistry, while in the second, multiple laminar flamelets are used to represent the structure of a turbulent diffusion flame. Suitable techniques for online reaction rate tabulation and chemical mechanism reduction are also incorporated, to make the use of bigger mechanisms possible (up to 150 species). The two models are compared and validated by simulating constant-volume diesel combustion in a wide range of operating conditions, including variations of ambient temperature and oxygen concentration. Comparison between the computed and experimental data on flame structure, auto-ignition and flame lift-off enables an understanding of the main relevant differences between the models in the way both auto-ignition and flame stabilisation processes are predicted.     

Large eddy simulation of turbulent premixed flames using level-set G-equation

Level-set G-equation and stationary flamelet chemistry are used in large eddy simulation of a propane/air premixed turbulent flame stabilized by a bluff body. The aim was to study the interaction between the flame front and turbulent eddies, and in particular to examine the effect of sub-grid scale (SGS) eddies on the wrinkling of the flame surface. The results indicated that the two types of turbulence eddies—the resolved large scale eddies and the unresolved SGS eddies—have different effects on the flame. The fluctuation of the flame surface, which is responsible for the broadening of the time averaged mean flame brush by turbulence, depends on the large resolved turbulence eddies. Time averaged mean flow velocity, temperature, and major species concentrations mainly depend on the large scale resolved eddies. The unresolved SGS eddies contribute to the wrinkling at the SGS level and play an important role in the enhancement of the propagation speed of the resolved flame front. In addition, the spatially filtered intermediate species, such as radicals, and the spatially filtered reaction rates strongly depend on the small SGS eddies. The asymptotic behavior of flame wrinkling by the SGS eddies, with respect to the decrease in filter size and grid size, is investigated further using a simplified level-set equation in a model shear flow. It is shown that to minimize the influence of the SGS eddies, fine grid and filter size may have to be used.     

Numerical study of the effects of gravity on soot formation in laminar coflow methane/air diffusion flames under different air stream velocities

Numerical simulations of laminar coflow methane/air diffusion flames at atmospheric pressure and different gravity levels were conducted to gain a better understanding of the effects of gravity on soot formation by using relatively detailed gas-phase chemistry and complex thermal and transport properties coupled with a semi-empirical two-equation soot model. Thermal radiation was calculated using the discrete-ordinates method coupled with a non-grey model for the radiative properties of CO, CO₂, H₂O, and soot. Calculations were conducted for three coflow air velocities of 77.6, 30, and 5 cm/s to investigate how the coflowing air velocity affects the flame structure and soot formation at different levels of gravity. The coflow air velocity has a rather significant effect on the streamwise velocity and the fluid parcel residence time, especially at reduced gravity levels. The flame height and the visible flame height in general increase with decreasing the gravity level. The peak flame temperature decreases with decreasing either the coflow air stream velocity or the gravity level. The peak soot volume fraction of the flame at microgravity can either be greater or less than that of its normal gravity counterpart, depending on the coflow air velocity. At sufficiently high coflow air velocity, the peak soot volume fraction increases with decreasing the gravity level. When the coflow air velocity is low enough, soot formation is greatly suppressed at microgravity and extinguishment occurs in the upper portion of the flame with soot emission from the tip of the flame owing to incomplete oxidation. The numerical results provide further insights into the intimate coupling between flame size, residence time, thermal radiation, and soot formation at reduced gravity level. The importance of thermal radiation heat transfer and coflow air velocity to the flame structure and soot formation at microgravity is demonstrated for the first time.     

Soot formation and temperature field structure in co-flow laminar methane-air diffusion flames at pressures from 10 to 60 atm

The effects of pressure on soot formation and the structure of the temperature field were studied in co-flow methane-air laminar diffusion flames over a wide pressure range, from 10 to 60 atm in a high-pressure combustion chamber. The selected fuel mass flow rate provided diffusion flames in which the soot was completely oxidized within the visible flame envelope and the flame was stable at all pressures considered. The spatially resolved soot volume fraction and soot temperature were measured by spectral soot emission as a function of pressure. The visible (luminous) flame height remained almost unchanged from 10 to 100 atm. Peak soot concentrations showed a strong dependence on pressure at relatively lower pressures; but this dependence got weaker as the pressure is increased. The maximum conversion of the fuel’s carbon to soot, 12.6%, was observed at 60 atm at approximately the mid-height of the flame. Radial temperature gradients within the flame increased with pressure and decreased with flame height above the burner rim. Higher radial temperature gradients near the burner exit at higher pressures mean that the thermal diffusion from the hot regions of the flame towards the flame centerline is enhanced. This leads to higher fuel pyrolysis rates causing accelerated soot nucleation and growth as the pressure increases.     

Compressible turbulent flame speeds of highly turbulent standing flames

In this paper we present the first measurement of turbulent burning velocities of a highly turbulent compressible standing flame induced by shock-driven turbulence in a Turbulent Shock Tube. High-speed schlieren, chemiluminescence, PIV, and dynamic pressure measurements are made to quantify flame–turbulence interaction for high levels of turbulence at elevated temperatures and pressure. Distributions of turbulent velocities, vorticity and turbulent strain are provided for regions ahead and behind the standing flame. The turbulent flame speed is directly measured for the high-Mach standing turbulent flame. From measurements of the flame turbulent speed and turbulent Mach number, transition into a non-linear compressibility regime at turbulent Mach numbers above 0.4 is confirmed, and a possible mechanism for flame generated turbulence and deflagration-to-detonation transition is established.     

热电偶沉积法测量乙烯非预混火焰的试验研究

烟黑容积份额的测量是研究烟黑生成的反应机理的额的薪方法.本文详细描述了采用热电偶沉积法测量烟黑容积份额的理论基础和数据处理过程,并将此方法应用于层流乙烯非预混火焰的测量中.测量结果表明,该火焰中烟黑容积份额的分布同火焰结构和火焰温度都有关.     

The effect of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame

The influence of preferential diffusion on soot formation in a laminar ethylene/air diffusion flame was investigated by numerical simulation using three different transport property calculation methods. One simulation included preferential diffusion and the other two neglected preferential diffusion. The results show that the neglect of preferential diffusion or the use of unity Lewis number for all species results in a significant underprediction of soot volume fraction. The peak soot volume fraction is reduced from 8.0 to 2.0 ppm for the studied flame when preferential diffusion is neglected in the simulation. Detailed examination of numerical results reveals that the underprediction of soot volume fraction in the simulation neglecting preferential diffusion is due to the slower diffusion of some species from main reaction zone to PAH and soot formation layer. The slower diffusion of these species causes lower PAH formation rate and thus results in lower soot inception rate and smaller particle surface area. The smaller surface area further leads to smaller surface growth rate. In addition, the neglect of preferential diffusion also leads to higher OH concentration in the flame, which causes the higher specific soot oxidation rate. The lower inception rate, smaller surface growth rate and higher specific oxidation rate results in the lower soot volume fraction when preferential diffusion is neglected. The finding of the paper implies the importance of preferential diffusion for the modeling of not only laminar but maybe also some turbulent flames.     

The effects of swirling partially premixed flame on scaled kinetic energy transport in a gas turbine-like combustor

The multi-scale interaction between combustion and turbulence is of great importance in modifying the small-scale flame structure and kinetic energy, especially in swirling flames under practical conditions. In the present study, direct numerical simulation of swirling partially premixed flame is conducted within a model combustor under gas turbine conditions. The reactive flow is compared to the corresponding non-reactive one to investigate the influence of combustion on the scaled kinetic energy transport. Kinetic energy spectra demonstrate that the turbulent kinetic energy is reduced in the dissipative subrange while enhanced in the energetic one by the flame. The critical scale is located in the inertial subrange and close to the estimated turbulent flame thickness. Filtering analyses show that the resolved-scale kinetic energy is augmented by the increased large-scale pressure-gradient work in the reactive flow, while the subgrid-scale kinetic energy is attenuated by the enhanced small-scale viscous dissipation. The backscatter prevails in the heat release regions when the filter size is larger than the laminar flame thickness, and this effect decreases with the swirling flow developing downstream. The interaction between the kinetic energy flux and the local dilatation as well as the subgrid-scale pressure-gradient work is also investigated to achieve a comprehensive understanding about the effects of combustion on the backscatter.     

Simulations of sooting turbulent jet flames using a hybrid flamelet/stochastic Eulerian field method

The stochastic Eulerian field method is applied to simulate 12 turbulent C₁?C₃ hydrocarbon jet diffusion flames covering a wide range of Reynolds numbers and fuel sooting propensities. The joint scalar probability density function (PDF) is a function of the mixture fraction, enthalpy defect, scalar dissipation rate and representative soot properties. Soot production is modelled by a semi-empirical acetylene/benzene-based soot model. Spectral gas and soot radiation is modelled using a wide-band correlated-k model. Emission turbulent radiation interactions (TRIs) are taken into account by means of the PDF method, whereas absorption TRIs are modelled using the optically thin fluctuation approximation. Model predictions are found to be in reasonable agreement with experimental data in terms of flame structure, soot quantities and radiative loss. Mean soot volume fractions are predicted within a factor of two of the experiments whereas radiant fractions and peaks of wall radiative fluxes are within 20%. The study also aims to assess approximate radiative models, namely the optically thin approximation (OTA) and grey medium approximation. These approximations affect significantly the radiative loss and should be avoided if accurate predictions of the radiative flux are desired. At atmospheric pressure, the relative errors that they produced on the peaks of temperature and soot volume fraction are within both experimental and model uncertainties. However, these discrepancies are found to increase with pressure, suggesting that spectral models describing properly the self-absorption should be considered at over-atmospheric pressure.     

Transported JPDF modelling and measurements of soot at elevated pressures

Accurate measurements and modelling of soot formation in turbulent flames at elevated pressures form a crucial step towards design methods that can support the development of practical combustion devices. A mass and number density preserving sectional model is here combined with a transported joint-scalar probability density function (JDPF) method that enables a fully coupled scalar space of soot, gas-phase species and enthalpy. The approach is extended to the KAUST turbulent non-premixed ethylene-nitrogen flames at pressures from 1 to 5 bar via an updated global bimolecular (second order) nucleation step from acetylene to pyrene. The latter accounts for pressure-induced density effects with the rate fitted using comparisons with full detailed chemistry up to 20 bar pressure and with experimental data from a WSR/PFR configuration and laminar premixed flames. Soot surface growth is treated via a PAH analogy and soot oxidation is considered via O, OH and O₂ using a Hertz-Knudsen approach. The impact of differential diffusion between soot and gas-phase particles is included by a gradual decline of diffusivity among soot sections. Comparisons with normalised experimental OH-PLIF and PAH-PLIF signals suggest good predictions of the evolution of the flame structure. Good agreement was also found for predicted soot volume statistics at all pressures. The importance of differential diffusion between soot and gas-phase species intensifies with pressure with the impact on PSDs more evident for larger particles which tend to be transported towards the fuel rich centreline leading to reduced soot oxidation.     

Prediction of soot and thermal radiation in a model gas turbine combustor burning kerosene fuel spray at different swirl levels

Combustion of kerosene fuel spray has been numerically simulated in a laboratory scale combustor geometry to predict soot and the effects of thermal radiation at different swirl levels of primary air flow. The two-phase motion in the combustor is simulated using an Eulerian–Lagragian formulation considering the stochastic separated flow model. The Favre-averaged governing equations are solved for the gas phase with the turbulent quantities simulated by realisable k–? model. The injection of the fuel is considered through a pressure swirl atomiser and the combustion is simulated by a laminar flamelet model with detailed kinetics of kerosene combustion. Soot formation in the flame is predicted using an empirical model with the model parameters adjusted for kerosene fuel. Contributions of gas phase and soot towards thermal radiation have been considered to predict the incident heat flux on the combustor wall and fuel injector. Swirl in the primary flow significantly influences the flow and flame structures in the combustor. The stronger recirculation at high swirl draws more air into the flame region, reduces the flame length and peak flame temperature and also brings the soot laden zone closer to the inlet plane. As a result, the radiative heat flux on the peripheral wall decreases at high swirl and also shifts closer to the inlet plane. However, increased swirl increases the combustor wall temperature due to radial spreading of the flame. The high incident radiative heat flux and the high surface temperature make the fuel injector a critical item in the combustor. The injector peak temperature increases with the increase in swirl flow mainly because the flame is located closer to the inlet plane. On the other hand, a more uniform temperature distribution in the exhaust gas can be attained at the combustor exit at high swirl condition.     

Hydrogen enrichment enhances soot formation in swirl-stabilized non-premixed turbulent combustion of ethylene in a model gas turbine combustor

A switch from fossil fuels to hydrogen is currently not feasible mostly due to supply and infrastructure issues. One of the possible approaches, and this is now practiced to a limited extent in industrial gas turbines, is to blend relatively small amounts of hydrogen with fossil fuels curbing the carbon dioxide emissions. However, studies assessing the influence of modest amounts of hydrogen blending with hydrocarbon fuels on soot processes yielded contradictory results. Most of these experimental and numerical studies were performed on laminar diffusion flames and studies on turbulent flames are scarce. One of the confounding factors in assessing the influence of hydrogen is selection of a control experiment in which the fossil fuel is blended with the same amount of an inert diluent. Using helium in the control experiment is preferable because of its similar transport properties and heat capacity to those of hydrogen. Hence, we studied the soot processes in a model gas turbine combustor in which the flame is stabilized by an air swirl. Swirl-stabilized platform ensures that with and without hydrogen/helium dilution, the hydrodynamics of the combustor stays fixed. Base fuel ethylene is supplemented with hydrogen or helium by the same amount to separate the dilution affects and assess the direct chemical interaction of hydrogen related to soot formation. Soot volume fraction and primary soot particle diameters were measured by auto-compensating laser induced-incandescence for all cases. Flow field data obtained using stereoscopic particle image velocimetry is utilized to ascertain the hydrodynamic effects on soot distribution due to addition of lighter species. Soot formation was found to be enhanced by the addition of hydrogen when allowance was made for the dilution effect using the helium doped flame experiments. Possible causes of this observation including the molecular diffusivities of hydrogen and helium, and chemical interaction are discussed.     

Two-dimensional soot distributions in buoyant turbulent fires

Spatially resolved two-dimensional soot volume fractions were measured using laser-induced incandescence in 7.1 cm methane and ethylene turbulent buoyant flames to study the distributions of soot in vertical and horizontal planes, and to provide data for soot model validation. Factors affecting the LII signals were considered including the laser energy profile and the laser attenuation effects. The absolute soot volume fractions were obtained by comparison to existing extinction measurements. The instantaneous soot images were collected to cover the entire flame height. Statistical quantities of soot volume fractions including mean, root mean square, probability density function, and spatial correlation coefficient were calculated at five downstream locations. The results show that instantaneous distributions of soot volume fractions exhibit significant differences compared to the ensemble averages, strong fluctuation around the mean, relatively homogeneous probability density function, and highly anisotropic spatial correlation.     

3-D simulation of soot formation in a direct-injection diesel engine based on a comprehensive chemical mechanism and method of moments

Here, we propose both a comprehensive chemical mechanism and a reduced mechanism for a three-dimensional combustion simulation, describing the formation of polycyclic aromatic hydrocarbons (PAHs), in a direct-injection diesel engine. A soot model based on the reduced mechanism and a method of moments is also presented. The turbulent diffusion flame and PAH formation in the diesel engine were modelled using the reduced mechanism based on the detailed mechanism using a fixed wall temperature as a boundary condition. The spatial distribution of PAH concentrations and the characteristic parameters for soot formation in the engine cylinder were obtained by coupling a detailed chemical kinetic model with the three-dimensional computational fluid dynamic (CFD) model. Comparison of the simulated results with limited experimental data shows that the chemical mechanisms and soot model are realistic and correctly describe the basic physics of diesel combustion but require further development to improve their accuracy.     

On the combined effects of compositional inhomogeneity and ammonia addition to turbulent flames of ethylene

This paper is part of a broader program aimed at investigating the effects of co-firing clean fuels such as ammonia or hydrogen with hydrocarbons. The focus is on soot formation as well as flame stability in turbulent mixed-mode combustion, which is highly relevant in practical combustors. Ammonia substitution for nitrogen results in reduced flame stability, and this is correlated to differences in flame speed and extinction strain rate. While it is known that the addition of ammonia suppresses soot, visual inspection of compositionally inhomogeneous flames of ethylene-ammonia indicates a reduction in ammonia's ability to suppress soot formation. Measurements of soot volume fraction and laser-induced fluorescence in selected UV and visible bands are made along the centreline in selected flames to test this hypothesis. Experimental results are then compared to simulations in laminar diffusion flames, stratified counterflow flames, and partially premixed flames. All results confirm the soot-inhibiting ability of ammonia. Increasing inhomogeneity, leading to higher centreline mixture fractions, enhances soot formation, and the level of enhancement is greater for flames with ammonia than without. Moreover, it is found that partial premixing is ultimately responsible for determining the amount of soot formed as opposed to stratification of fuel mixtures near the pilot.     

Modeling nongray gas-phase and soot radiation in luminous turbulent nonpremixed jet flames

Much progress has been made in radiative heat transfer modeling with respect to treatment of nongray radiation from both gas-phase species and soot particles, while radiation modeling in turbulent flame simulations is still in its infancy. Aiming at reducing this gap, this paper introduces state-of-the-art models of gas-phase and soot radiation to turbulent flame simulations. The full-spectrum k-distribution method (Modest, M.F., 2003, Journal of Quantitative Spectroscopy & Radiative Transfer, 76, 69–83) is implemented into a three-dimensional unstructured CFD code for nongray radiation modeling. The mixture full-spectrum k-distributions including nongray absorbing soot particles are constructed from a narrow-band k-distribution database created for individual gas-phase species, and an efficient scheme is employed for their construction in CFD simulations. A detailed reaction mechanism including NO _x and soot kinetics is used to predict flame structure, and a detailed soot model using a method of moments is employed to determine soot particle size distributions. A spherical-harmonic P₁ approximation is invoked to solve the radiative transfer equation. An oxygen-enriched, turbulent, nonpremixed jet flame is simulated, which features large concentrations of gas-phase radiating species and soot particles. Nongray soot modeling is shown to be of greater importance than nongray gas modeling in sooty flame simulations, with gray soot models producing large errors. The nongray treatment of soot strongly influences flame temperatures in the upstream and the flame-tip region and is essential for accurate predictions of NO. The nongray treatment of gases, however, weakly influences upstream flame temperatures and, therefore, has only a small effect on NO predictions. The effect of nongray soot radiation on flame temperature is also substantial in downstream regions where the soot concentration is small. Limitations of the P₁ approximation are discussed for the jet flame configuration; the P₁ approximation yields large errors in the spatial distribution of the computed radiative heat flux for highly anisotropic radiation fields such as those in flames with localized, near-opaque soot regions.     

Modelling nongrey gas-phase and soot radiation in luminous turbulent nonpremixed jet flames

Much progress has been made in radiative heat transfer modelling with respect to the treatment of nongrey radiation from both gas-phase species and soot particles, while radiation modelling in turbulent flame simulations is still in its infancy. Aiming at reducing this gap, this paper introduces state-of-the-art models of gas-phase and soot radiation to turbulent flame simulations. The full-spectrum k-distribution method (M.F. Modest, 2003, Journal of Quantitative Spectroscopy & Radiative Transfer, 76, 69–83) is implemented into a three-dimensional unstructured computational fluid dynamics (CFD) code for nongrey radiation modelling. The mixture full-spectrum k-distributions including nongrey absorbing soot particles are constructed from a narrow-band k-distribution database created for individual gas-phase species, and an efficient scheme is employed for their construction in CFD simulations. A detailed reaction mechanism including NO _x and soot kinetics is used to predict flame structure, and a detailed soot model using a method of moments is employed to determine soot particle size distributions. A spherical harmonic P₁ approximation is invoked to solve the radiative transfer equation. An oxygen-enriched, turbulent, nonpremixed jet flame is simulated, which features large concentrations of gas-phase radiating species and soot particles. Nongrey soot modelling is shown to be of greater importance than nongrey gas modelling in sooty flame simulations, with grey soot models producing large errors. The nongrey treatment of soot strongly influences flame temperatures in the upstream and the flame-tip region and is essential for accurate predictions of NO. The nongrey treatment of gases, however, weakly influences upstream flame temperatures and, therefore, has only a small effect on NO predictions. The effect of nongrey soot radiation on flame temperature is also substantial in downstream regions where the soot concentration is small. Limitations of the P₁ approximation are discussed for the jet flame configuration; the P₁ approximation yields large errors in the spatial distribution of the computed radiative heat flux for highly anisotropic radiation fields such as those in flames with localized, near-opaque soot regions.     

Strain rate and flame orientation statistics in the near-wall region for turbulent flame-wall interaction

Flame–wall interaction (FWI) in premixed turbulent combustion has been analysed based on a counter-flow like configuration at the statistically stationary state. For the present configuration, the two FWI sub-zones, i.e the influence zone and the quenching zone, can be identified from the DNS results. Detailed analysis of the important quantities, such as the flame temperature, flame–wall distance, wall heat flux, flame curvature and dilatation (including the flame normal and tangential strain rates), and some orientation relations between the flame normal and the principal strain rate directions, have been reported, together with the physical explanations. All these statistical results are determined by the relative strengths of the wall heat flux, thermal expansion and the flame–turbulence interaction.