Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLAME PROPAGATION”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Numerical Investigation of Fuel Property Effects on Mixed-Mode Combustion in a Spark-Ignition Engine

In this research, lean mixed-mode combustion is numerically investigated using computational fluid dynamics (CFD) in a spark-ignition engine. A new E30 fuel surrogate is developed using a neural network model with matched octane numbers. A skeletal mechanism is also developed by automated mechanism reduction and by incorporating a NO x submechanism. A hybrid approach that couples the G-equation model and the well-stirred reactor model is employed for turbulent combustion modeling. The developed CFD model is shown to well predict pressure and apparent heat release rate (AHRR) traces compared with experiment. Two types of combustion cycles (deflagration-only and mixed-mode cycles) are observed. The mixed-mode cycles feature early flame propagation and subsequent end-gas auto-ignition, leading to two distinctive AHRR peaks. The validated CFD model is then employed to investigate the effects of NO x chemistry. The NO x chemistry is found to promote auto-ignition through the residual gas, while the deflagration phase remains largely unaffected. Sensitivity analysis is finally performed to understand effects of fuel properties, including heat of vaporization (HoV) and laminar flame speed (S L ). An increased HoV tends to suppress auto-ignition through charge cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased S L is found to significantly promote both flame propagation and end-gas auto-ignition. The promoting effect of S L on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end-gas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗
DOE OSTI2024

Flame Flashback Investigations in Hydrogen-Enriched Low Swirl Burner using High-Speed Hydroxyl Planar Laser-Induced Fluorescence

Hydrogen/ Hydrogen-rich flame is a potential strategy to reduce carbon emissions for the development of next-generation lean premixed combustors. However, the higher flashback propensity with hydrogen-enrichment owing to faster chemical kinetics is one of the key challenges in existing natural gas turbine combustors. In this study, the fundamentals of stabilized flame dynamics and flashback events in a premixed low swirl burner (LSB) configuration are visualized experimentally using nanosecond (ns)-based high-speed hydroxyl radical planar laser-induced fluorescence (OH-PLIF). The LSB has an optically accessible pre-mixing tube allowing laser diagnostic investigation of flashback events. The inlet conditions are varied systematically with respect to pre-mixer velocities, flame equivalence ratio (ϕ), and hydrogen content and OH-PLIF characterization is performed in stable flame at 20 kHz repetition rate. An increase in flashback propensity is observed with increasing ϕ, hydrogen content, and decreasing V via detailed investigation of flame lift-off length (L) above the burner rim. Flashback ϕ investigations show an expected linearly increasing trend with decreasing hydrogen content and increasing V for each swirler studied and the results agree well with detailed L investigations in stable flame configuration. The propagating flame anchored likely on the back side of the perforated plate within 30-40 ms of the entrance forming the distorted conical flame front.

Parajuli, Pradeep↗
DOE OSTI2021

Large-Eddy Simulation of Laser-Ignited Direct Injection Gasoline Spray for Emission Control

Large-Eddy Simulations (LES) of a gasoline spray, where the mixture was ignited rapidly during or after injection, were performed in comparison to a previous experimental study with quantitative flame motion and soot formation data [SAE 2020-01-0291] and an accompanying Reynolds-Averaged Navier–Stokes (RANS) simulation at the same conditions. The present study reveals major shortcomings in common RANS combustion modeling practices that are significantly improved using LES at the conditions of the study, specifically for the phenomenon of rapid ignition in the highly turbulent, stratified mixture. At different ignition timings, benchmarks for the study include spray mixing and evaporation, flame propagation after ignition, and soot formation in rich mixtures. A comparison of the simulations and the experiments showed that the LES with Dynamic Structure turbulence were able to capture correctly the liquid penetration length, and to some extent, spray collapse demonstrated in the experiments. For early and intermediate ignition timings, the LES showed excellent agreement to the measurements in terms of flame structure, extent of flame penetration, and heat-release rate. However, RANS simulations (employing the common G-equation or well-stirred reactor) showed much too rapid flame spread and heat release, with connections to the predicted turbulent kinetic energy. With confidence in the LES for predicted mixture and flame motion, the predicted soot formation/oxidation was also compared to the experiments. The soot location was well captured in the LES, but the soot mass was largely underestimated using the empirical Hiroyasu model. An analysis of the predicted fuel–air mixture was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing.

42 ENGINEERING↗
DOE OSTI2020

Spherically expanding flame in silane–hydrogen–nitrous oxide–argon mixtures

The effect of silane addition on the laminar flame speed ($S^0_u$) of flames propagating in hydrogen–nitrous oxide–argon mixtures has been investigated experimentally for the first time using the spherically expanding flame technique in a constant volume combustion chamber. Replacing hydrogen by silane and maintaining the equivalence ratio constant, much higher flame speeds, explosion peak pressures, and pressure rise coefficients were measured. A previously developed detailed reaction model has been updated based on ab initio thermodynamic properties calculations and collision limit violation analysis. The improved reaction model demonstrates encouraging performance in predicting the flame speed, with a mean absolute error below 11%. To explain the effect of silane addition on the flame dynamics, a number of parameters have been calculated including OH and H rate of production, heat release rate per reaction, and sensitivity coefficient on . The dynamics of freely propagating flames in SiH 4 –H 2 –N 2 O–Ar mixtures is essentially controlled by reactions of the H–O–N chemical system: N 2 O+H=N 2 +OH, OH+H 2 =H 2 O+H, and N 2 O(+M)=N 2 +O(+M). Whereas silane addition does not influence much the rate of production of OH, it significantly modifies that of H with a number of pyrolytic chemical pathways of silicon hydrides, such as SiH+H 2 =SiH 2 +H and Si+H 2 =SiH+H, which act as sink of H atom as they proceed in the backward direction. The reactions forming SiO(s) and SiO 2 (s), such as SiO+OH=SiO 2 (s)+H and 2SiO=2SiO(s), are exothermic and significantly contribute to the temperature increase. The adiabatic, constant pressure flame temperature for mixture containing silane is significantly higher, up to several 100’s K. The increase of induced by silane addition seems to be mostly related to the large increase of the flame temperature which leads to higher energy release rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Effects of dilution and pressure on combustion characteristics within externally heated micro-tubes

Here this study investigates the impact of pressure and dilution on combustion behavior of of C 2H6/air/N2 mixtures within an externally heated micro-tube using both experimental and analytical techniques. Experiments are conducted for pressures up to 10 bar at multiple dilution levels and document three flame regimes, i.e. weak flames, FREI and normal flames. Characteristic temperatures for flame stabilization, as well as and ignition and extinction are obtained by applying a novel thin filament pyrometry method (TFP), which uses a thin silicon carbide (SiC) filament that enables temperature measurements for both reacting and non-reacting cases. Results show a decrease of characteristic temperatures with increasing pressure, which can be which can be partially compensated by a decrease of O2 content. Further, regime transitions occur at lower velocity as the pressure increases, whereas weak flame positions are insensitive to dilution at a given pressure. Predictions from the analytical model compare favorably with experimental results. The most prominent finding is that as dilution increases, the transition from weak flame to FREI becomes less abrupt and is completely lost for high dilution levels (i.e. low O2 content). Both experimental and analytical results indicate that the unstable FREI combustion regime does not exist at high dilution levels, i.e. weak flame and normal flame regimes are no longer clearly separated. As FREI are interpreted as a cyclic transition from an ignition zone to a propagating flame, results for highly diluted conditions are interpreted as a combustion regime where classical flame propagation ceases to exist.

42 ENGINEERING↗
DOE OSTI2025

Plasma assisted NH 3 /H 2 /air ignition in nanosecond discharges with non-equilibrium energy transfer

Ammonia (NH 3 ), with its high energy density and easiness to store and transport as a hydrogen carrier, has become a promising alternative green fuel. However, its adoption in power generation is hindered by challenges such as low burning velocity, slow low-temperature oxidation, high NO x emissions, and ignition difficulty. Here, this work computationally investigates the effects of non-equilibrium energy transfer by nanosecond discharges on NH 3 ignition and flame propagation in an NH 3 /H 2 /air flow at 700 K and 1 atm. The simulation results demonstrate that NH 3 /air mixtures require a large ignition energy due to their large critical ignition radius. It is shown that adding 30 % hydrogen significantly reduces the critical ignition radius and minimum ignition energy. Two-dimensional modeling further shows a non-monotonic dependence of ignition kernel volume on the applied voltage and reduced electric field. The optimum ignition enhancement occurs at 200 Td where the generation of electronically excited species and radicals including N 2 (B), O( 1 D) and OH becomes most efficient. Higher voltages divert electron energy toward ionization, which makes it less effective for NH 3 ignition. The study also identifies an optimal electrode gap size for a given pulse energy. Smaller gap sizes increase deposited energy density, raising temperature and radical concentrations. However, excessive reduction of the gap distance reduces flame propagation speed due to the flame stretch effect in rich mixtures with the effective Lewis number greater than unity. A nonlinear relationship between pulse repetition frequency and ignition kernel volume is observed in a nanosecond pulsed high frequency discharge (NPHFD). An optimal frequency range of 200 kHz to 2 MHz is found when two pulses are used. In addition, an optimal number of pulses exists for each pulse repetition frequency, with higher frequencies requiring more pulses to maximize the overlap region. These findings provide critical insights on developing controlled plasma discharge techniques for efficient NH 3 ignition in reactive flows within internal combustion engines and gas turbines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

THE EFFECT OF FUEL REACTIVITY AND EXHAUST GAS RECIRCULATION ON KNOCK PROPENSITY OF NATURAL GAS

The development of high efficiency, spark ignited natural gas engines is currently limited by engine knock at high compression ratio/elevated boost pressures and misfire at lean conditions/high exhaust gas recirculation (EGR) levels. The knock and misfire limits are further confounded by the wide variety in fuel reactivity observed in “pipeline quality” natural gas. In this study, a rapid compression machine was used to characterize the effects of EGR and variation in natural gas fuel reactivity on the homogeneous ignition delay, flame propagation rate, and end- gas autoignition propensity for stoichiometric natural gas/oxidizer/EGR blends. A reduced chemical kinetic mechanism was also developed to accurately model the homogeneous ignition delays measured in the Colorado State University rapid compression machine (CSU RCM). Pipeline quality natural gas with a range of chemical reactivity (68 < Methane Number < 95) was simulated using mixtures of CH4, C2H6, and C3H8. Exhaust gas recirculation gases were simulated with mixtures of Ar, CO2, CO, and NO at substitution rates of 0 to 30 mass percent. Ignition delay period under homogeneous autoignition conditions was measured at compressed pressures of 30.2 to 34.0 bar and compressed temperatures of 667 to 980 K. End-gas autoignition fraction and flame propagation rate were measured by initiating a laser spark in the center of the combustion chamber, after compression, at pressures of 30.7 to 32.7 bar and temperatures of 751 to 795 K. The results indicate that both fuel reactivity and the presence of reactive species (NO and CO) in the exhaust gas recirculation have a strong impact on end-gas autoignition fraction. A chemical kinetic mechanism was developed to predict homogeneous ignition delays for pipeline quality natural gas in a pressure and temperature range of 1-100 bar and 500-1000 K respectively. This mechanism accurately predicted measured homogeneous ignition delay in the RCM with a total average relative error of 11.0%.

Mohr, Jeffrey↗
DOE OSTI2023

Sensitivity of He Flames in X-Ray Bursts to Nuclear Physics

Through the use of axisymmetric 2D hydrodynamic simulations, we further investigate laterally propagating flames in X-ray bursts (XRBs). Our aim is to understand the sensitivity of a propagating helium flame to different nuclear physics. Using the Castro simulation code, we confirm the phenomenon of enhanced energy generation shortly after a flame is established by adding 12 C(p, γ) 13 N(α, p) 16 O to the network, in agreement with the past literature. This sudden outburst of energy leads to a short accelerating phase, causing a drastic alteration in the overall dynamics of the flame in XRBs. Furthermore, we investigate the influence of different plasma screening routines on the propagation of the XRB flame. We finally examine the performance of simplified spectral deferred correction, a novel approach to hydrodynamics and reaction coupling incorporated in Castro, as an alternative to operator splitting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The role of collisional energy transfer in the thermal and prompt dissociation of 1-methyl allyl

Collisional energy transfer plays a pivotal role in a-priori calculations of the pressure-dependent rate constants obtained from a master equation. However, accurate determinations of collisional energy transfer parameters are rare and so most kinetic studies rely on best-fits of these parameters to experimental measurements (if available) or estimates guided from literature studies. In this work, we have quantified the effect of the uncertainty in energy transfer parameters on the thermal and prompt dissociation kinetics of a resonance stabilized radical, 1-methyl allyl (1MA), of relevance to the combustion of 1- and 2-butene isomers. Simulations using literature kinetics models were performed to assess the impact of these uncertainties on flame propagation and speciation data in laminar flames of 1- and 2-butene. Analyses of the uncertainty propagated by the energy transfer parameters for 1MA dissociation to the flame simulations, in particular the laminar flame speed, indicate an intricate coupling between the kinetics of 1MA dissociation (chain propagation) and its reaction with H-atoms (chain-termination). Ab-initio based theoretical calculations were also performed to obtain pressure-dependent kinetics for the reaction of 1MA with H-atoms. Lastly, theoretically calculated energy transfer parameters were used to best characterize the kinetics and branching between the chain propagating 1MA dissociation and the chain-terminating reaction 1MA + H.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Influence of gas expansion on the propagation of a premixed flame in a spatially periodic shear flow

It has been previously demonstrated that thermal gas expansion might have a role in boundary layer flashback of premixed turbulent flames, inducing local flow-reversal in the boundary layer's low-velocity streaks on the reactants’ side of the flame and facilitating its upstream propagation. Here, we perform a two-dimensional numerical investigation of the interaction between a periodic shear flow and a laminar premixed flame. The periodic shear is a simplified model for the oncoming prolonged streamwise velocity streaks with alternating regions of high and low velocities found in turbulent boundary layers in the vicinity of the walls. The parametric study focuses on the amplitude and wavelength of the periodic shear flow and on the gas expansion ratio (unburnt-to-burnt density ratio). With the increase of the amplitudes of the periodic shear flow and of the gas expansion, the curved flame velocity increases monotonically. The flame velocity dependence on the periodic shear wavelength is non-monotonic, which is consistent with previous theoretical studies of curved premixed flame velocity. The flame shape that is initially formed by the oncoming periodic shear appears to be metastable. At a later stage of the flame propagation, the flame shape transforms into the stationary one dominated by the Darrieus-Landau instability.

42 ENGINEERING↗
DOE OSTI2024

A numerical study of ignition and flame development characteristics in GCI combustion using large eddy simulations and chemical explosive mode analysis

Here, this work investigates the ignition and flame development processes of low reactivity fuel combustion under compression ignition conditions based on the large eddy simulation approach. The chemical explosive mode analysis (CEMA) is employed to characterize the local combustion features, including gas-liquid fuel zone, auto-ignition, diffusion-assisted, extinction, cool flame and post-ignition zone, among which auto-ignition and post-ignition are found to play a key role in the overall heat release process. The local flame propagation modes in gasoline compression ignition (GCI) are determined by quantifying the relative magnitude of diffusion/chemistry at a representative progress variable in the pre-ignition zone. The results show that autoignition fronts and deflagration waves exist simultaneously in the ignition and intense high temperature heat release (HTHR) stages, but autoignition fronts dominate. In addition, the chemical kinetic processes of four heat release periods are analyzed. The heat release during the ignition period is found to be dominated by the reactions CH 3 + H (+M) <=> CH 4 (+M) and CH 3 CHO + H <=> CH 2 CHO + H 2 . The reaction CH 2 OH + OH <=> CH 2 O + H 2 O always plays an important role in the heat releases during the other three combustion stages including intense HTHR, moderate HTHR and post-combustion.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Numerical and Experimental Investigations on the Ignition Behavior of OME

On the path towards climate-neutral future mobility, the usage of synthetic fuels derived from renewable power sources, so-called e-fuels, will be necessary. Oxygenated e-fuels, which contain oxygen in their chemical structure, not only have the potential to realize a climate-neutral powertrain, but also to burn more cleanly in terms of soot formation. Polyoxymethylene dimethyl ethers (PODE or OMEs) are a frequently discussed representative of such combustibles. However, to operate compression ignition engines with these fuels achieving maximum efficiency and minimum emissions, the physical-chemical behavior of OMEs needs to be understood and quantified. Especially the detailed characterization of physical and chemical properties of the spray is of utmost importance for the optimization of the injection and the mixture formation process. The presented work aimed to develop a comprehensive CFD model to specify the differences between OMEs and dodecane, which served as a reference diesel-like fuel, with regards to spray atomization, mixing and auto-ignition for single- and multi-injection patterns. The simulation results were validated against experimental data from a high-temperature and high-pressure combustion vessel. The sprays’ liquid and vapor phase penetration were measured with Mie-scattering and schlieren-imaging as well as diffuse back illumination and Rayleigh-scattering for both fuels. To characterize the ignition process and the flame propagation, measurements of the OH* chemiluminescence of the flame were carried out. Significant differences in the ignition behavior between OMEs and dodecane could be identified in both experiments and CFD simulations. Liquid penetration as well as flame lift-off length are shown to be consistently longer for OMEs. Zones of high reaction activity differ substantially for the two fuels: Along the spray center axis for OMEs and at the shear boundary layers of fuel and ambient air for dodecane. Additionally, the transient behavior of high temperature reactions for OME is predicted to be much faster.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Multi-regime mixing modeling for local extinction and re-ignition in turbulent non-premixed flame by using LES/FDF method

Local extinction and re-ignition occur in turbulent non-premixed combustion when the Damköhler number is not large enough and the combustion is not fully mixing controlled. The occurrence of local extinction introduces locally extinguished flame holes followed by a premixed flame propagation toward the hole center to potentially reignite the flame. The co-existence of the non-premixed and premixed combustion regimes complicates the modeling since traditional combustion models are mostly for a single regime. In this work, we examine the effect of multi-regime mixing modeling in the transported filtered density function (FDF) method on the predictions of local extinction and reignition. Predictions of local extinction and reignition remain a challenge for the FDF method despite the progress made in the past. To account for the multi-regime combustion, two different mixing timescale models for non-premixed and premixed combustion are combined. A flame index based on the gradients of fuel and oxidizer is used to define a weighting factor to blend the two mixing timescale models. A turbulent jet non-premixed flame with substantial local extinction, the Sydney piloted jet flame L, is adopted as a test case to examine the performance of the multi-regime model in large-eddy simulation/FDF modeling. It is found that the traditional non-premixed mixing timescale model when combined with the modified Curl mixing leads to global extinction for the Sydney flame L without the presence of the premixed combustion regime. After accounting for the multi-regime combustion with proper detection of the different combustion regimes, the predictions for the flame statistics and the amount of local extinction are significantly improved. Furthermore, it suggests the need of including multi-regime combustion for the predictions of local extinction and re-ignition in turbulent non-premixed combustion configurations.

42 ENGINEERING↗
DOE OSTI2021

Optical Investigation of a Partial Fuel Stratification Strategy to Stabilize Overall Lean Operation of a DISI Engine Fueled with Gasoline and E30

This paper offers new insights into a partial fuel stratification (PFS) combustion strategy that has proven to be effective at stabilizing overall lean combustion in direct injection spark ignition engines. To this aim, high spatial and temporal resolution optical diagnostics were applied in an optically accessible engine working in PFS mode for two fuels and two different durations of pilot injection at the time of spark: 210 µs and 330 µs for E30 (gasoline blended with ethanol by 30% volume fraction) and gasoline, respectively. In both conditions, early injections during the intake stroke were used to generate a well-mixed lean background. The results were compared to rich, stoichiometric and lean well-mixed combustion with different spark timings. In the PFS combustion process, it was possible to detect a non-spherical and highly wrinkled blue flame, coupled with yellow diffusive flames due to the combustion of rich zones near the spark plug. The initial flame spread for both PFS cases was faster compared to any of the well-mixed cases (lean, stoichiometric and rich), suggesting that the flame propagation for PFS is enhanced by both enrichment and enhanced local turbulence caused by the pilot injection. Different spray evolutions for the two pilot injection durations were found to strongly influence the flame kernel inception and propagation. PFS with pilot durations of 210 µs and 330 µs showed some differences in terms of shapes of the flame front and in terms of extension of diffusive flames. Yet, both cases were highly repeatable.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Direct numerical simulation of a high-pressure hydrogen micromix combustor: Flame structure and stabilisation mechanism

A high-pressure hydrogen micromix combustor has been investigated using direct numerical simulation with detailed chemistry to examine the flame structure and stabilisation mechanism. The configuration of the combustor was based on the design by Schefer et al., using numerical periodicity to mimic a large square array. A precursor simulation of an opposed jet-in-crossflow was first conducted to generate appropriate partially-premixed inflow boundary conditions for the subsequent reacting simulation. The resulting flame can be described as an predominantly-lean inhomogeneously-premixed lifted jet flame. Five main zones were identified: a jet mixing region, a core flame, a peripheral flame, a recirculation zone, and combustion products. The core flame, situated over the jet mixing region, was found to burn as a thin reaction front, responsible for over 85% of the total fuel consumption. The peripheral flame shrouded the core flame, had low mean flow with high turbulence, and burned at very lean conditions (in the distributed burning regime). It was shown that turbulent premixed flame propagation was an order-of-magnitude too slow to stabilise the flame at these conditions. Stabilisation was identified to be due to ignition events resulting from turbulent mixing of fuel from the jet into mean recirculation of very lean hot products. Ignition events were found to correlate with shear-driven Kelvin-Helmholtz vortices, and increased in likelihood with streamwise distance. At the flame base, isolated events were observed, which developed into rapidly burning flame kernels that were blown downstream. Further downstream, near-simultaneous spatially-distributed ignition events were observed, which appeared more like ignition sheets. The paper concludes with a broader discussion that considers generalising from the conditions considered here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Large-eddy simulation study on cycle-to-cycle variation of knocking combustion in a spark-ignition engine

The cycle-to-cycle variation in the knock intensity is commonly encountered under abnormal combustion conditions. The severity of these abnormal combustion events can vary significantly, and the efficiency of engines at high loads is limited in practice by heavy knocking phenomena. Since, a thorough analysis of such recurrent but non-cyclic phenomena via experiments alone becomes highly cumbersome, in the present work, a multi-cycle large-eddy simulation study was performed to quantitatively predict cyclic variability in the combustion process and cyclic knock intensity variability in a direct injection spark-ignition engine. To account for the turbulence-chemistry interaction effects on flame propagation, the G-equation combustion model was used. Detailed chemistry was solved outside the flame front with a toluene primary reference fuel skeletal kinetic mechanism. For both the mild knock and heavy knock conditions, the numerical results were validated against experimental measurements. Based on the simulation results, a correlation analysis was performed considering combustion phasing, peak cylinder pressure and maximum amplitude of pressure oscillation. Furthermore, a detailed three-dimensional spatial analysis illustrated the evolution of auto-ignition kernel development and propagation of pressure waves during knocking combustion for three typical cycles with different knock intensities. In this process, it was found that an early occurrence of auto-ignition in the end gas was prone to high knock intensity. Although multiple auto-ignition kernels were observed in different cycles, the degree of coupling between chemical heat release and pressure waves varied, thereby leading to different maximum amplitude of pressure oscillation values.

42 ENGINEERING↗