Flame stabilization with a tubular flame

A new technique to stabilize a flame in a high-velocity stream with use of a tubular flame has been proposed. To elucidate the validity of this technique, an experiment has been conducted by mounting a tubular flame burner on the nozzle. Flame stability limits and temperature distributions around the burner port have been determined, and experiments have been extended to the ducted combustion to measure pressure fluctuations and to analyze the burned gases. Results show that the tubular flame can successfully stabilize the main flame up to 130 m/s, which is the upper limit of the present supply facility. The main flame is well anchored at the exit of the nozzle, and the tubular flame efficiently supplies heat and radicals to the main flame. In the ducted combustion, the pressure fluctuations are reduced significantly. The exhaust gas analyses, however, indicate that an almost chemical equilibrium condition can be achieved at 50 m/s, but not at 90 and 130 m/s. Since the energy input relative to the main flame is just 6.1% at 130 m/s, the present tubular burner is not enough to burn all the unburned gas completely at high velocities, although the main flame can be anchored. The slit length and/or the slit width of the tubular flame burner should be larger to overcome this shortage. From the above results, it is concluded that the tubular flame has a potential for stabilizing a flame in a high-speed stream.     

Numerical study of transient evolution of lifted jet flames: partially premixed flame propagation and influence of physical dimensions

Three-dimensional (3D) unsteady Reynolds-averaged Navier–Stokes simulations of a spark-ignited turbulent methane/air jet flame evolving from ignition to stabilisation are conducted for different jet velocities. A partially premixed combustion model is used involving a correlated joint probability density function and both premixed and non-premixed combustion mode contributions. The 3D simulation results for the temporal evolution of the flame's leading edge are compared with previous two-dimensional (2D) results and experimental data. The comparison shows that the final stabilised flame lift-off height is well predicted by both 2D and 3D computations. However, the transient evolution of the flame's leading edge computed from 3D simulation agrees reasonably well with experiment, whereas evident discrepancies were found in the previous 2D study. This difference suggests that the third physical dimension plays an important role during the flame transient evolution process. The flame brush's leading edge displacement speed resulting from reaction, normal and tangential diffusion processes are studied at different typical stages after ignition in order to understand the effect of the third physical dimension further. Substantial differences are found for the reaction and normal diffusion components between 2D and 3D simulations especially in the initial propagation stage. The evolution of reaction progress variable scalar gradients and its interaction with the flow and mixing field in the 3D physical space have an important effect on the flame's leading edge propagation.     

Numerical simulation of flame dynamics associated with negative velocity induced by deformed flame shape

In this study, the influence of the negative velocity field formed ahead of an abruptly deformed flame tip on the propagation behaviour of a laminar premixed flame is numerically investigated. A strong deformation in the flame front is induced by imposing a very narrow, in-line pre-heating zone in the unburned region. The simulation is performed under low Mach number approximation by using a multi-scale multi-physics Computational Fluid Dynamics (CFD) solver FrontFlow/Red with one-step finite rate chemistry in order to track the time-dependent flame dynamics. The computed results unveil that the flame front is deformed significantly within a short time due to the narrow in-line pre-heating effect. The flame deformation induces a strong negative velocity field ahead of the deformed flame tip, acting in the direction of propagation, which gives rise to a strong pair vortex. This strong pair vortex interacts with the flame tip and then slides down along the flame surface in the upstream direction during propagation. This flame-vortex interaction causes further deformation in the flame surface in the upstream direction, and consequently, the flame exhibits a wave-like surface, which enhances the flame propagation speed. The auto-generation of a strong pair vortex ahead of the flame front due to the localised thermal input could be applied as one of the methods to control the combustion externally. It is also expected that the results obtained in this study could have a significant impact on the detailed understanding of the local thermo-fluid dynamical interaction process of turbulent combustion in practical combustors.     

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.     

LES/FGM investigation of ignition and flame structure in a gasoline partially premixed combustion engine

This paper presents a joint numerical and experimental study of the ignition process and flame structures in a gasoline partially premixed combustion (PPC) engine. The numerical simulation is based on a five-dimension Flamelet-Generated Manifold (5D-FGM) tabulation approach and large eddy simulation (LES). The spray and combustion process in an optical PPC engine fueled with a primary reference fuel (70% iso-octane, 30% n-heptane by volume) are investigated using the combustion model along with laser diagnostic experiments. Different combustion modes, as well as the dominant chemical species and elementary reactions involved in the PPC engines, are identified and visualized using Chemical Explosive Mode Analysis (CEMA). The results from the LES-FGM model agree well with the experiments regarding the onset of ignition, peak heat release rate and in-cylinder pressure. The LES-FGM model performs even better than a finite-rate chemistry model that integrates the full-set of chemical kinetic mechanism in the simulation, given that the FGM model is computationally more efficient. The results show that the ignition mode plays a dominant role in the entire combustion process. The diffusion flame mode is identified in a thin layer between the ultra fuel-lean unburned mixture and the hot burned gas region that contains combustion intermediates such as CO. The diffusion flame mode contributes to a maximum of 27% of the total heat release in the later stage of combustion, and it becomes vital for the oxidation of relatively fuel-lean mixtures.     

Experimental and LES investigation of premixed methane/air flame propagating in a tube with a thin obstacle

In this paper, an experimental and numerical investigation of premixed methane/air flame dynamics in a closed combustion vessel with a thin obstacle is described. In the experiment, high-speed video photography and a pressure transducer are used to study the flame shape changes and pressure dynamics. In the numerical simulation, four sub-grid scale viscosity models and three sub-grid scale combustion models are evaluated for their individual prediction compared with the experimental data. High-speed photographs show that the flame propagation process can be divided into five stages: spherical flame, finger-shaped flame, jet flame, mushroom-shaped flame and bidirectional propagation flame. Compared with the other sub-grid scale viscosity models and sub-grid scale combustion models, the dynamic Smagorinsky–Lilly model and the power-law flame wrinkling model are better able to predict the flame behaviour, respectively. Thus, coupling the dynamic Smagorinsky–Lilly model and the power-law flame wrinkling model, the numerical results demonstrate that flame shape change is a purely hydrodynamic phenomenon, and the mushroom-shaped flame and bidirectional propagation flame are the result of flame–vortex interaction. In addition, the transition from “corrugated flamelets” to “thin reaction zones” is observed in the simulation.     

DME/Oxygen wall-stabilized premixed cool flame

Wall-stabilized cool flames have been studied through numerical analysis and a series of experiments. One- and two-dimensional numerical simulations were performed to estimate the characteristics of the wall-stabilized cool flames, such as flammability and temperature/species distributions. Based on the computational results, the ignition condition of the cool flame at a fixed wall temperature has been identified with the strain rate between the van't Hoff point and the cool flame extinction point. In experiments with an impinging jet burner and a heated plate, the spontaneous ignition of the cool flame on the heated wall has been successfully established under the conditions predicted by the present numerical simulations. Spatial distributions of the HCHO concentration and flame temperature were measured through formaldehyde Planar Laser-Induced Fluorescence (HCHO-PLIF) and thermocouple measurements, respectively. It is found that the measurement data show a reasonable accordance with the simulation results with reduced low-temperature reactivity.     

Large-eddy simulations of a stratified-charge direct-injection spark-ignition engine: Comparison with experiment and analysis of cycle-to-cycle variations

Multiple-cycle large-eddy simulations (LES) have been performed for an optically accessible, single-cylinder, four-stroke-cycle, spray-guided direct-injection spark-ignition (SG-DISI) engine operating in a stratified globally fuel-lean mode. The simulations combine a standard Smagorinsky turbulence model, a stochastic Lagrangian parcel method for liquid fuel injection and fuel spray modeling, a simple energy-deposition spark-ignition model, and a modified thickened flame model for turbulent flame propagation through highly stratified reactant mixtures. Comparisons between simulations and experiments include individual-cycle and ensemble-average pressure and apparent-heat-release-rate traces, individual-cycle and ensemble-average indicated mean effective pressures (IMEP), and instantaneous two-dimensional vapor-equivalence-ratio contours. Although the number of LES cycles is small (35), the results show that the simulations are able to capture the global combustion behavior that is observed in the experiments, including cycle-to-cycle variations. The simulation results are then analyzed further to provide insight into the conditions that lead to misfire versus robust combustion. As has been reported in earlier experimental and LES studies for homogeneous-charge SI engines, local conditions in the vicinity of the spark gap at the time of ignition largely determine the subsequent flame development. However, in contrast to homogeneous-charge engines, no single local or global quantity correlates as strongly with the eventual peak pressure or IMEP for each cycle. Rather, it is the interplay among the early flame kernel, the velocity field that it experiences, and the fuel distribution that it encounters that ultimately determines the fate of each combustion event. Deeper analysis and quantitative statistical comparisons between experiments and simulations will require the simulation of larger numbers of engine cycles.     

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.     

Effects of confinement on premixed turbulent swirling flame using large Eddy simulation

This paper utilises large eddy simulation (LES) to study swirling reacting flows by comparison with experimental observations. The purpose is to provide further insights in engineering designs, as well as to improve modelling. A reduced-scale swirl burner has been developed for the experiments. Comparison of particle image velocimetry (PIV) measurements with LES results using finite rate chemistry shows that LES captures all the salient features of an unconfined flame including velocity and temperature distributions. However, when the flame is confined within a cylindrical combustor, the simulated flame shape is initially not consistent with experimental observation. Investigations show that the discrepancy is caused by the often practised assumption of adiabatic wall temperature. With the use of an assumed wall temperature distribution guided by laboratory observation, results of LES are consistent with experiments. Although the latter LES approach requires more computational resources, the improvement is found to be justified.     

Numerical simulation and validation of SI-CAI hybrid combustion in a CAI/HCCI gasoline engine

SI-CAI hybrid combustion, also known as spark-assisted compression ignition (SACI), is a promising concept to extend the operating range of CAI (Controlled Auto-Ignition) and achieve the smooth transition between spark ignition (SI) and CAI in the gasoline engine. In this study, a SI-CAI hybrid combustion model (HCM) has been constructed on the basis of the 3-Zones Extended Coherent Flame Model (ECFM3Z). An ignition model is included to initiate the ECFM3Z calculation and induce the flame propagation. In order to precisely depict the subsequent auto-ignition process of the unburned fuel and air mixture independently after the initiation of flame propagation, the tabulated chemistry concept is adopted to describe the auto-ignition chemistry. The methodology for extracting tabulated parameters from the chemical kinetics calculations is developed so that both cool flame reactions and main auto-ignition combustion can be well captured under a wider range of thermodynamic conditions. The SI-CAI hybrid combustion model (HCM) is then applied in the three-dimensional computational fluid dynamics (3-D CFD) engine simulation. The simulation results are compared with the experimental data obtained from a single cylinder VVA engine. The detailed analysis of the simulations demonstrates that the SI-CAI hybrid combustion process is characterised with the early flame propagation and subsequent multi-site auto-ignition around the main flame front, which is consistent with the optical results reported by other researchers. Besides, the systematic study of the in-cylinder condition reveals the influence mechanism of the early flame propagation on the subsequent auto-ignition.     

Triple flame: Inherent asymmetries and pentasectional character

A two-dimensional triple-flame numerical model of a laminar combustion process reflects flame asymmetric structural features that other analytical models do not generate. It reveals the pentasectional character of the triple flame, composed of the central pure diffusion-flame branch and the fuel-rich and fuel-lean branches, each of which is divided into two sections: a near-stoichiometric section and a previously unreported near-flammability-limits section with combined diffusion and premixed character. Results include propagation velocity, fuel and oxidiser mass fractions, temperature and reaction rates. Realistic stoichiometric ratios and reaction orders match experimental planar flame characteristics. Constant density, a one-step reaction, and a mixture fraction gradient at the inlet as the simulation parameter are imposed. The upstream equivalence ratio or the upstream reactant mass fractions are linear or hyperbolic functions of the transverse coordinate. The use here of experimental kinetics data differs from previous analytical works and results in flame asymmetry and different flammability limits. Upstream mixture composition gradient affects propagation velocity, flame curvature, diffusion flame reaction rate, and flammability limits. Flammability limits extend beyond those of a planar flame due to transverse heat and mass diffusion causing the pentasectional character.     

Three-dimensional direct simulations and structure of expanding turbulent methane flames

Direct numerical simulations (DNS) are ideally suited to investigate in detail turbulent reacting flows in simple geometries. For an increasing number of applications, detailed models must be employed to describe the chemical processes with sufficient accuracy. Despite the huge cost of such simulations, recent progress has allowed the direct numerical simulation of turbulent premixed flames while employing complete reaction schemes. We briefly describe our own developments in this field and use the resulting DNS code to investigate more extensively the structure of premixed methane flames expanding in a three-dimensional turbulent velocity field, initially homogeneous and isotropic. This situation typifies, for example, the initial flame development after spark ignition in a gas turbine or an internal combustion engine. First investigation steps have been carried out at low turbulence levels on this same configuration in the past Symposium, and we build on top of these former results. Here, a considerably higher Reynolds number is considered, the simulation has been repeated twice in to limit the possibility of spurious, very specific results, and several complementary post-processing steps are carried out. Characteristic features concerning the observed combustion regime are presented. We then investigate in a quantitative manner the evolution of flame surface area, global stretch-rate, flame front curvature, flame thickness, and correlation between thickness and curvature. The possibility of obtaining reliable information on flame front curvature from two-dimensional slices is checked by comparison with the exact procedure.     

Effects of pressure rise rate on laminar flame speed under normal and engine-relevant conditions

Laminar flame speed (LFS) is one of the most important physicochemical properties of a combustible mixture. At normal and elevated temperatures and pressures, LFS can be measured using propagating spherical flames in a closed chamber. LFS is also used in certain turbulent premixed flame modelling for combustion in spark ignition engines. Inside the closed chamber or engine, transient pressure rise occurs during the premixed flame propagation. The effects of pressure rise rate (PRR) on LFS are examined numerically in this study. One-dimensional simulations are conducted for spherical flame propagation in a closed chamber. Detailed chemistry and transport are considered. Different values of PRR at the same temperature and pressure are achieved through changing the spherical chamber size. It is found that the effect of PRR on LFS is negligible under the normal and engine-relevant conditions considered in this study. This observation is then explained through the comparison between the unsteady and convection terms in the energy equation for a premixed flame.     

一种测量发光火焰温度的数据处理方法研究

火焰温度是燃烧诊断的重要参数之一,它对研究各种燃烧过程具有重要价值。根据多光谱辐射测温所用的参考温度数学模型,提出了一种基于遗传算法的新的数据处理方法。该方法对火焰发射率与波长的关系依次进行了三点直线拟合,并通过遗传算法进行了优化,从而得到了发光火焰的温度和发射率。采用多波长高温计测量了某种固体推进剂燃烧火焰的自辐射光谱,并进行了数据处理。计算结果表明:火焰温度计算值与理论值之差在±100K以内;这种基于遗传算法的新的数据处理方法是测量发光火焰温度的一种可行性方法。     

Simple model of inhibition of chain-branching combustion processes

A simple kinetic model has been suggested to describe the inhibition and extinction of flame propagation in reaction systems with chain-branching reactions typical for hydrocarbon systems. The model is based on the generalised model of the combustion process with chain-branching reaction combined with the one-stage reaction describing the thermal mode of flame propagation with the addition of inhibition reaction steps. Inhibitor addition suppresses the radical overshoot in flame and leads to the change of reaction mode from the chain-branching reaction to a thermal mode of flame propagation. With the increase of inhibitor the transition of chain-branching mode of reaction to the reaction with straight-chains (non-branching chain reaction) is observed. The inhibition part of the model includes a block of three reactions to describe the influence of the inhibitor. The heat losses are incorporated into the model via Newton cooling. The flame extinction is the result of the decreased heat release of inhibited reaction processes and the suppression of radical overshoot with the further decrease of the reaction rate due to the temperature decrease and mixture dilution. A comparison of the results of modelling laminar premixed methane/air flames inhibited by potassium bicarbonate (gas phase model, detailed kinetic model) with the results obtained using the suggested simple model is presented. The calculations with the detailed kinetic model demonstrate the following modes of combustion process: (1) flame propagation with chain-branching reaction (with radical overshoot, inhibitor addition decreases the radical overshoot down to the equilibrium level); (2) saturation of chemical influence of inhibitor, and (3) transition to thermal mode of flame propagation (non-branching chain mode of reaction). The suggested simple kinetic model qualitatively reproduces the modes of flame propagation with the addition of the inhibitor observed using detailed kinetic models.     

Microgravity droplet combustion: effect of tethering fiber on burning rate and flame structure

Droplets tethering on fibers has become a well established technique for conducting droplet combustion experiments in microgravity conditions. The effects of these supporting fibers are frequently assumed to be negligible and are not considered in the experimental analysis or in numerical simulations. In this work, the effect of supporting fibers on the characteristics of microgravity droplet combustion has been investigated numerically; a priori predictions have then been compared with published experimental data. The simulations were conducted using a transient one-dimensional spherosymmetric droplet combustion model, where the effect of the supporting fiber was implicitly taken into account. The model applied staggered convective flux finite volume method combined with high-order implicit time integration. Thermal radiation was evaluated using a statistical narrow band radiation model. Chemical kinetics and thermophysical properties were represented in rigorous detail. Tether fiber diameter, droplet diameter, ambient pressure and oxygen concentration were varied over a range for n-decane droplets in the simulations. The results of the simulations were compared to previously published experiments conducted in the Japan Microgravity Center (JAMIC) 10 second drop tower and the NASA Glenn Research Center (GRC) 5.2 second drop tower. The model reproduces closely nearly all aspects of tethered n-decane droplet burning phenomena, which included droplet burning history, transient and average burning rate, and flame standoff ratio. The predictions show that the presence of the tethering fiber significantly influences the observed burning rate, standoff ratio, and extinction.     

On combustion in a closed rectangular channel with initial vorticity

This article examines the detailed combustion process in a theoretical model with applicability to combustion in a wave rotor or wave disc engine. The model comprises a single channel into which an initial loading of methane and air is admitted and ignited after all inlet and exit ports have been closed. Combustion takes place at constant volume. However, the initial gaseous mixture in the channel is not at rest. The initial opening and closing of the ports generates significant vorticity which influences the evolution of the combustion process. Numerical evaluations are provided for the detailed flame shape for simplified (one-step) chemistry and a simulation using the detailed 235-step San Diego scheme. Quantities examined are the evolution of vorticity, pressure fluctuations, mass consumption rate, flame surface area and the influences on combustion of adiabatic and non-adiabatic channel walls. Combustion regimes are identified and compared with simpler model studies (no initial flow). Pointwise quantities are examined to describe the various stages of burning in the channel. The focus of the study is directed towards quantities that influence overall burning rate and completeness of combustion.     

Bayesian analysis of syngas chemistry models

Syngas chemistry modelling is an integral step toward the development of safe and efficient syngas combustors. Although substantial effort has been undertaken to improve the modelling of syngas combustion, models nevertheless fail in regimes important to gas turbine combustors, such as low temperature and high pressure. In order to investigate the capabilities of syngas models, a Bayesian framework for the quantification of uncertainties has been used. This framework, given a set of experimental data, allows for the calibration of model parameters, determination of uncertainty in those parameters, propagation of that uncertainty into simulations, as well as determination of model evidence from a set of candidate syngas models. Here, three syngas combustion models have been calibrated using laminar flame speed measurements from high pressure experiments. After calibration the resulting uncertainty in the parameters is propagated forward into the simulation of laminar flame speeds. The model evidence is then used to compare candidate models for the given set of experimental conditions and results. Additionally, the technique MUM-PCE, an interesting uncertainty minimisation method for kinetics models, has been compared to the Bayesian method for this application to the prediction of syngas laminar flame speeds. This comparison shows the importance of model form error and experimental error representations in the uncertainty quantification context, for these choices significantly affect uncertainty quantification results.     

Turbulent flame simulation taking advantage of tabulated chemistry self-similar properties

Predicting the flame shape, its stabilization process, and pollutant emissions in practical combustion devices requires to incorporate complex chemistry features. As detailed chemical schemes are too voluminous for practical numerical simulations, tabulated chemistry techniques have been proposed to account for the complexity of kinetics in turbulent flame simulations. Unfortunately, the size of these databases may become a crucial issue for efficient implementation on massively parallel computers. A reduction strategy that takes advantage of self-similar properties of tabulated chemistry is proposed for turbulent combustion modeling. A reduction of the database size by a factor of 1000 is achieved. This procedure is successfully applied to a RANS simulation of a turbulent jet flame.