Early Flame Kernel Growth Research Articles

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. The present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 µm can predict the LFL as ϕCFD=0.61 and match the displacement speed of the kernel's boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 µm) is needed to capture the consumption of formaldehyde (CH2O) in kernel's reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (ϕ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

Numerical investigation of a fueled pre-chamber spark-ignition natural gas engine

Pre-chamber spark-ignition (PCSI) is a leading advanced ignition concept for internal combustion engines with the potential to enable diesel-like efficiency in medium-duty/heavy-duty (MD/HD) natural gas (NG) engines. By leveraging distributed ignition sources from multiple turbulent jets, the PCSI technology can deliver extremely short combustion duration in ultra-lean mixtures and significantly improve the engine thermal efficiency. However, in the automotive industry there is a lack of adequate science base and predictive simulation tools required for commercial development of PCSI engines. In this study, Reynolds-Average Navier-Stokes simulations are carried out to describe the combustion process in lean-burn NG engines, focusing on the combustion modeling approach. Two combustion models, multi-zone well-stirred reactor (MZ-WSR) and G-equation, are used to simulate the combustion process in an MD NG engine equipped with a fueled-PCSI system for four operating conditions close to the lean operating limit. A skeletal chemical mechanism and a laminar flame speed tabulation are used to compute the combustion accurately. Simulation results are compared with experimental data regarding measured cylinder pressure, heat release rate, and combustion duration. By dividing the PCSI combustion process into four distinct phases, the difference between the two models’ results for each phase is analyzed in detail. The MZ-WSR model overestimates the combustion duration for early flame kernel growth in the pre-chamber due to the lack of a specific formulation to take turbulence-chemistry interaction into account. Despite the prolonged combustion duration and low pressure built-up inside the pre-chamber, the model matches the combustion rate in the main-chamber. In contrast, the G-equation model delivers good agreements for the pre-chamber combustion and turbulent jet-driven combustion processes. However, the model starts to underestimate the combustion rate in the main-chamber, especially under ultra-lean mixture conditions. Finally, improvements are needed for both models to simulate the later combustion stage that occurred in the near-wall regions.

G-equation based ignition model for direct injection spark ignition engines

To accurately model the Direct Injection Spark Ignition (DISI) combustion process, it is important to account for the effects of the spark energy discharge process. The proximity of the injected fuel spray and spark electrodes leads to steep gradients in local velocities and equivalence ratios, particularly under cold-start conditions when multiple injection strategies are employed. The variations in the local properties at the spark plug location play a significant role in the growth of the initial flame kernel established by the spark and its subsequent evolution into a turbulent flame. In the present work, an ignition model is presented that is compatible with the G-Equation combustion model, which responds to the effects of spark energy discharge and the associated plasma expansion effects. The model is referred to as the Plasma Velocity on G-surface (PVG) model, and it uses the G-surface to capture the early kernel growth. The model derives its theory from the Discrete Particle Ignition (DPIK) model, which accounts for the effects of electrode heat transfer, spark energy, and chemical heat release from the fuel on the early flame kernel growth. The local turbulent flame speed has been calculated based on the instantaneous location of the flame kernel on the Borghi-Peters regime diagram. The model has been validated against the experimental measurements given by Maly and Vogel,1 and the constant volume flame growth measurements provided by Nwagwe et al.2 Multi-cycle simulations were performed in CONVERGE3 using the PVG ignition model in combination with the G-Equation-based GLR4 model in a RANS framework to capture the combustion characteristics of a DISI engine. Good agreements with the experimental pressure trace and apparent heat-release rates were obtained. Additionally, the PVG ignition model was observed to substantially reduce the sensitivity of the default G-sourcing ignition method employed by CONVERGE.

LES investigation of cycle-to-cycle variation in a SI optical access engine using TFM-AMR combustion model

Multi-cycle large-eddy simulations (LES) are performed to investigate combustion cycle-to-cycle variability (CCV) in a gasoline spark ignited optical access engine operating under homogeneous stoichiometric conditions. Combustion is addressed with the Thickened Flame Model (TFM) and finite rate chemistry is accounted for through a reduced oxidation reaction mechanism. In view of the fact that computational costs of LES engine simulations are still very high today, this work investigates the use of adaptive mesh refinement (AMR) in the flame zone in conjunction with the artificial flame thickening applied by the TFM model. The paper discusses how the resulting coupled TFM-AMR combustion model allows good resolution of the flame, maintaining accuracy at acceptable costs. First, the details of the coupled model are presented and the effects of the parameters are explored, highlighting their impact on the combustion prediction. Then, computational fluid dynamics (CFD) simulation results are validated against experimental data collected in a low-speed low-load engine point, by comparing 20 LES cycles and 100 measured cycles, for mass fraction burned, combustion phasing, flame images and CCV indices. Lastly, a detailed investigation on the fastest and slowest numerical cycles is presented, analyzing instantaneous flame structures, ignition behaviors, propagation speeds, and probability density function (PDF) of the instantaneous velocity fluctuation around the spark region. The results show that combustion variability is highly correlated to the resolved velocity field and the resolved turbulence intensity, which is found to be the main cause of CCV and affects the early flame kernel growth. This work is an early attempt to use TFM-AMR combustion model for LES simulations of internal combustion engines.

Large-eddy simulation study of combustion cyclic variation in a lean-burn spark ignition engine

Multi-cycle large-eddy simulation (LES) was performed to investigate combustion cyclic variability (CCV) in a single cylinder spark ignition engine with a homogeneous lean (λ = 1.25) isooctane-air mixture. The aim was to obtain physical insights into the early stage of combustion and its influence on CCV. Propagation of the flame was modeled by a transport equation for the filtered flame surface density within the LES framework. The ignition process was represented by the imposed stretch spark ignition model (ISSIM-LES). Ten consecutive cold flow LES cycles followed by two initialization cycles (12 cycles in total) were used to perform the reactive simulations concurrently. The simulation results were compared with experimental data. Although the number of computed cycles was fairly low, the LES was able to reproduce the cyclic variability observed in experiments both quantitatively and qualitatively. Firstly, validation of the simulation was done by comparing measured pressure traces. Secondly, correlations between the timing of the 10% fuel burnt mass fraction with early flame kernel growth and initial-to-turbulent transition period (in which there was an asymmetric flame kernel that persisted through the early development periods) were determined. Thirdly, calculated results of the flame propagation were analyzed at two cross-sections (in swirl and tumble planes) of the combustion chamber, which highlighted differences in instantaneous flame structures and propagation characteristics between the fastest and slowest cycles. Good overall agreement was obtained between the measurements and simulation data. The results revealed that the instantaneous velocity and fluctuation of flows around the spark vicinity affect growth of the early flame kernel and cause combustion cyclic variability.

Control of Early Flame Kernel Growth by Multi-Wavelength Laser Pulses for Enhanced Ignition

The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

The role of spray-enhanced swirl flow for combustion stabilization in a stratified-charge DISI engine

Implementation of spray-guided stratified-charge direct-injection spark-ignited (DISI) engines is inhibited by the occurrence of misfire and partial burns. Engine-performance tests demonstrate that increasing engine speed induces combustion instability, but this deterioration can be prevented by generating swirling flow during the intake stroke. In-cylinder pressure-based heat-release analysis reveals that the appearance of poor-burn cycles is not solely dependent on the variability of early flame-kernel growth. Cycles can experience burning-rate regression during later combustion stages and may or may not recover before the end of the cycle. Thermodynamic analysis and optical diagnostics are used here to clarify why swirl improves the combustion repeatability from cycle to cycle.The fluid dynamics of swirl/spray interaction was previously demonstrated using high-speed PIV measurements of in-cylinder motored flow. It was found that the sprays of the multi-hole injector redistribute the intake-generated swirl flow momentum, thereby creating a better-centered higher angular-momentum vortex with reduced variability. The engine operation with high swirl was found to have significant improvement in cycle-to-cycle variations of both flow pattern and flow momentum.This paper is an extension of the previous work. Here, PIV measurements and flame imaging are applied to fired operation for studying how the swirl flow affects variability of ignition and subsequent combustion phases. PIV results for fired operation are consistent with the measurements made of motored flow. They demonstrate that the spark-plasma motion is highly correlated with the direction of the gas flow in the vicinity of the spark-plug gap. Without swirl, the plasma is randomly stretched towards either side of the spark plug, causing variability in the ignition of the two spray plumes that are straddling the spark plug. In contrast, swirl flow always convects the spark plasma towards one spray plume, causing a more repeatable ignition. The swirl decreases local RMS velocity, consistent with an observed reduction of early-burn variability. Broadband flame imaging demonstrates that with swirl, the flame consistently propagates in multiple directions to consume fuel–air mixtures within the piston bowl. In contrast, operation without swirl displays higher variability of flame-spread patterns, occasionally causing the appearance of partial-burn cycles.

Efficiency improvement of a spark-ignition engine at full load conditions using exhaust gas recirculation and variable geometry turbocharger – Numerical study

The numerical analysis of performance of a four cylinder highly boosted spark-ignition engine at full load is described in this paper, with the research focused on introducing high pressure exhaust gas recirculation for control of engine limiting factors such as knock, turbine inlet temperature and cyclic variability. For this analysis the cycle-simulation model which includes modeling of the entire engine flow path, early flame kernel growth, mixture stratification, turbulent combustion, in-cylinder turbulence, knock and cyclic variability was applied. The cylinder sub-models such as ignition, turbulence and combustion were validated by using the experimental results of a naturally aspirated multi cylinder spark-ignition engine. The high load operation, which served as a benchmark value, was obtained by a standard procedure used in calibration of engines, i.e. operation with fuel enrichment and without exhaust gas recirculation. By introducing exhaust gas recirculation and by optimizing other engine operating parameters, the influence of exhaust gas recirculation on engine performance is obtained. The optimum operating parameters, such as spark advance, intake pressure, air to fuel ratio, were found to meet the imposed requirements in terms of fuel consumption, knock occurrence, exhaust gas temperature and variation of indicated mean effective pressure. By comparing the results of the base point with the results that used exhaust gas recirculation the improvement in fuel consumption of 8.7%, 11.2% and 1.5% at engine speeds of 2000rpm, 3500rpm and 5000rpm is obtained. Additionally, by using the presented numerical methodology the influence of the specific operating parameter on the overall behavior of the complex charging system was shown.

Modelling of early flame kernel growth towards a better understanding of cyclic combustion variability in SI engines

The analysis of particular effects that influence cyclic combustion variability in spark-ignition engines is presented in this paper. The cycle-simulations are performed over 16 knock-free full load operating conditions of a single cylinder engine fuelled with gasoline. For the modelling of early flame kernel growth a newly developed quasi-dimensional ignition model is applied. The ignition model includes a detailed description of the electric circuit, the electric spark length, the spark plug geometry and the flame kernel growth. The main turbulent combustion is calculated using the extended quasi-dimensional fractal combustion model requiring the correct prediction of in-cylinder turbulence level. The cyclic combustion variability is simulated by a variation of the in-cylinder turbulence level, the flow angle at spark plug and the stratification of air equivalence ratio from cycle-to-cycle. The statistical analysis of the indicated mean effective pressure (IMEP) gives a more detailed quantitative picture of the causes of the cyclic combustion variability in SI engines over the different operating conditions. The variation of in-cylinder turbulence is found to be the dominant factor that affects the cyclic variability (up to 66%), while the variation in the early flame kernel growth rate is responsible for up to 18% of the overall cyclic combustion variability in the modelled SI engine at full load conditions.

Influence of flow and ignition fluctuations on cycle-to-cycle variations in early flame kernel growth

Engine Cycle-to-Cycle Variations (CCVs) constitute a restriction for numerous concepts aimed at improving fuel consumption in spark-ignited engines. In consequence, engine CCV modeling has become an active research topic. Nevertheless, few fundamental studies have thus far examined the individual influences of the key parameters. In the present study, we performed 2D Direct Numerical Simulations (DNSs) of early flame kernel growth with semi-detailed chemistry solving to identify the controlling parameters of heat release fluctuations that are relevant to engine CCVs. Our objective was to assess how certain turbulence and flame parameters – namely turbulent flow motion pattern, turbulence intensity, integral length scale, mixture strength and initial kernel size – influence combustion variability. Two extreme mixture conditions were selected (stable operation at stoichiometry and unstable operation at the lean flammability limit), and computation conditions were derived from actual low load engine operating points. Results provide a classification of the key parameters, indicating initial kernel size and turbulent flow motion pattern as the main contributors. Additionally, DNSs highlighted the crucial role of mixture strength in the occurrence of local flame quenching and in the related amplification of heat release fluctuations.

LES study of cycle-to-cycle variations in a spark ignition engine

Multi-cycle Large Eddy Simulation (LES) of flow and combustion in a single cylinder engine set-up are presented in order to validate the ability of this technique to reproduce cyclic combustion variability in a spark ignition engine. The simulated engine operating point is part of a larger experimental database acquired at IFP and specifically designed to validate LES techniques in terms of cycle-to-cycle prediction. To characterize individual cycles, the database combines acquisitions of pressure and temperature in different locations of the engine set-up with optical diagnostics in the intake pipe and the cylinder. This new experimental database combined with LES is a powerful tool to study cycle-to-cycle variations (CCV). The operating points include points with low and high CCV levels. The LES includes the whole engine set-up, and covers 25 consecutive four-stroke engine cycles for a reference operating point with low CCV. The results show that LES is able to reproduce the experimental findings, and illustrate how LES can help explain the sources of CCV. In the present case it is shown that CCV are essentially due to velocity fluctuations at the spark plug, which induce variations of the early flame kernel growth and of the overall combustion duration.

Investigation of pre-mixed flame-kernel/vortex interactions via high-speed imaging

To reduce cycle-to-cycle variations in SI-IC engines, knowledge of early flame kernel growth in a turbulent flow field is required. Understanding the interaction between a flame kernel and a vortex is an important fundamental step toward this goal. This paper presents high-speed movies of combustion luminosity during the interaction of a laminar vortex with a spark-generated pre-mixed flame kernel in a quiescent combustion chamber. The resulting time evolution of the perturbed flame kernel shows that laminar vortices of various sizes and vortex strengths can increase the kernel growth rate by at least a factor of 3 and significantly increase combustion reaction rates by involving additional highly curved and stretched flame fronts.

Time-resolved imaging of flame kernels: Laser spark ignition of H 2/O 2/Ar mixtures

The shape and structure of developing flame kernels in laser-induced spark ignited hydrogen/air mixtures is investigated as a function of gas composition and time. Using planar laser-induced fluorescence (PLIF) to measure the spatial distribution of OH radicals produced inside the reaction zone, we have recorded the evolution of the nascent flame kernel in a series of images following the laser-induced spark. This series provides the rate of flame growth, the evolution of the flame shape, and the intensity of the PLIF signal as a function of time for both igniting flames and nonignition events. The reaction zones grow quickly at early times, but slowly decrease in propagation rate as the energy density within the flame kernel decreases. A distinct anisotropy is observed in the expanding spark and flame kernel. At short times ( t < 100 μs), a toroidal shape is observed similar to that seen previously for electrode-spark ignitions and for laser ignitions in methane/air. There is also a tendency for the flame to grow back toward the ignition laser. Successful ignitions appear virtually identical to failed ignitions during the first 100 μs. Significant differences, notably in intensity, appear between 100 and 500 μs following the spark. These observations imply that early flame kernel growth is dominated by gas motion induced by the short-duration spark. The ultimate fate of an ignition lies with the chemistry of the reactions which determines whether the gas undergoes a transition from hot plasma to propagating flame.