Understanding Low-Temperature First-Stage Ignition Delay: Propane
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL 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.
This work describes the development of a transported Livengood-Wu (L-W) integral model for computational fluid dynamics (CFD) simulation to predict auto-ignition and engine knock tendency. The currently employed L-W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition combustion for three different fuels, showing good accuracy in prediction of auto-ignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L-W integral model is coupled with G-equation model to indicate end-gas auto-ignition and knock tendency in CFD simulations of a direct injection spark ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely Co-Optima Alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multi-variable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency indicated by the fuel energy contained in the auto-ignited region of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.
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.
Abstract An iterative dynamic chemical stiffness removal method (IDCSR) based on quasi-steady-state approximation (QSSA) is proposed. The IDCSR method is built on a previously developed non-iterative method which has proved to work well for small timestep sizes. A novel iterative procedure is designed in IDCSR to enable explicit time integration of stiff chemistry at relatively large timestep sizes relevant to practical reacting flow simulations. The effectiveness of the iterative procedure is first demonstrated with a toy problem and homogeneous auto-ignition with fixed integration step sizes, showing that larger timestep sizes can be allowed for explicit time integration using IDCSR compared with the previous non-iterative method. IDCSR is then compared with existing explicit chemistry solvers for simulations of homogeneous auto-ignition and shows similar or lower computational cost but significantly higher accuracy across a wide range of timestep sizes. IDCSR is further combined with an automatic adaptive time-stepping scheme for simulations of 0-D homogeneous auto-ignition and a 2-D laminar lifted n -dodecane jet flame. For the 0-D auto-ignition simulations, IDCSR is shown to reduce both the error (by 43%–90%) and computational cost (by 6–15 times) compared with existing explicit solvers, while achieving speed-up factors of up to 400 compared with VODE for a wide range of timestep sizes and reaction mechanisms. For the 2-D jet flame simulations, speed-up factors of 15 and 31 for chemistry integration, and 5 and 9 for overall simulation, are achieved by IDCSR compared with CVODE with and without analytic Jacobian, respectively.
In this work, spray flames under diesel engine conditions are numerically investigated to understand the ignition and stabilization characteristics, based on state-of-the-art large eddy simulations (LES) and chemical explosive mode analysis (CEMA). In particular, CEMA is employed to identify such critical flame features as auto-ignition, premixed flame fronts, cool and strongly burning flames, and local extinction. In particular, a recently developed CEMA-based criterion for demarcating local combustion modes, that is, auto-ignition, diffusion-assisted ignition, and extinction, is applied to gain in-depth insights into the interactions of local chemical and transport processes. A series of one-dimensional non-premixed flamelets with varying air temperatures and scalar dissipation rates are first analyzed using CEMA with a focus on the role of each local mode in the transient ignition process. CEMA-based diagnostics are then applied to diesel spray flames at different ambient temperatures. The structure of local combustion modes is identified, while the contribution from each local mode to overall flame burning is quantified at different stages of the flame developing process. Effects of the ambient temperature on local combustion mode distribution, and subsequently on the ignition and stabilization mechanisms, are further examined. Results show that auto-ignition is the main mechanism for flame initiation and stabilization, while the role of diffusion is substantially enhanced as ambient temperature decreases.
A space-based experiment is currently under development to study diffusion-controlled, gas-phase, low temperature oxidation reactions, cool flames and auto-ignition in an unstirred, static reactor. At Earth's gravity (1g), natural convection due to self-heating during the course of slow reaction dominates diffusive transport and produces spatio-temporal variations in the thermal and thus species concentration profiles via the Arrhenius temperature dependence of the reaction rates. Natural convection is important in all terrestrial cool flame and auto-ignition studies, except for select low pressure, highly dilute (small temperature excess) studies in small vessels (i.e., small Rayleigh number). On Earth, natural convection occurs when the Rayleigh number (Ra) exceeds a critical value of approximately 600. Typical values of the Ra, associated with cool flames and auto-ignitions, range from 104-105 (or larger), a regime where both natural convection and conduction heat transport are important. When natural convection occurs, it alters the temperature, hydrodynamic, and species concentration fields, thus generating a multi-dimensional field that is extremely difficult, if not impossible, to be modeled analytically. This point has been emphasized recently by Kagan and co-workers who have shown that explosion limits can shift depending on the characteristic length scale associated with the natural convection. Moreover, natural convection in unstirred reactors is never "sufficiently strong to generate a spatially uniform temperature distribution throughout the reacting gas." Thus, an unstirred, nonisothermal reaction on Earth does not reduce to that generated in a mechanically, well-stirred system. Interestingly, however, thermal ignition theories and thermokinetic models neglect natural convection and assume a heat transfer correlation of the form: q=h(S/V)(T(bar) - Tw) where q is the heat loss per unit volume, h is the heat transfer coefficient, S/V is the surface to volume ratio, and (T(bar) - Tw ) is the spatially averaged temperature excess. This Newtonian form has been validated in spatially-uniform, well-stirred reactors, provided the effective heat transfer coefficient associated with the unsteady process is properly evaluated. Unfortunately, it is not a valid assumption for spatially-nonuniform temperature distributions induced by natural convection in unstirred reactors. "This is why the analysis of such a system is so difficult." Historically, the complexities associated with natural convection were perhaps recognized as early as 1938 when thermal ignition theory was first developed. In the 1955 text "Diffusion and Heat Exchange in Chemical Kinetics", Frank-Kamenetskii recognized that "the purely conductive theory can be applied at sufficiently low pressure and small dimensions of the vessel when the influence of natural convection can be disregarded." This was reiterated by Tyler in 1966 and further emphasized by Barnard and Harwood in 1974. Specifically, they state: "It is generally assumed that heat losses are purely conductive. While this may be valid for certain low pressure slow combustion regimes, it is unlikely to be true for the cool flame and ignition regimes." While this statement is true for terrestrial experiments, the purely conductive heat transport assumption is valid at microgravity (mu-g). Specifically, buoyant complexities are suppressed at mu-g and the reaction-diffusion structure associated with low temperature oxidation reactions, cool flames and auto-ignitions can be studied. Without natural convection, the system is simpler, does not require determination of the effective heat transfer coefficient, and is a testbed for analytic and numerical models that assume pure diffusive transport. In addition, mu-g experiments will provide baseline data that will improve our understanding of the effects of natural convection on Earth.
In this paper, the first flamelet analysis is conducted of a highly resolved DNS of a multi-injection flame with both auto-ignition and ignition induced by flame-flame interaction. A novel method is proposed to identify the different combustion modes of ignition processes using generalized flamelet equations. The state-of-the-art DNS database generated by Rieth et al. (US National Combustion Meeting, 2019) for a multi-injection flame in a Diesel engine environment is investigated. Three-dimensional flamelets are extracted from the DNS at different time instants with a focus on auto-ignition and interaction-ignition processes. The influences of mixture field interactions and the scalar dissipation rate on the ignition process are investigated by varying the species composition boundary conditions of the transient flamelet equations. Budget analyses of the generalized flamelet equations show that the transport along the mixture fraction iso-surface is insignificant during the auto-ignition process, but becomes important when interaction-ignition occurs, which is further confirmed through a flamelet regime classification method.
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.
We propose an anomaly detection method for multi-variate scientific data based on analysis of high-order joint moments. Using kurtosis as a reliable measure of outliers, we suggest that principal kurtosis vectors, by analogy to principal component analysis (PCA) vectors, signify the principal directions along which outliers appear. The inception of an anomaly, then, manifests as a change in the principal values and vectors of kurtosis. Obtaining the principal kurtosis vectors requires decomposing a fourth order joint cumulant tensor for which we use a simple, computationally less expensive approach that involves performing a singular value decomposition (SVD) over the matricized tensor. We demonstrate the efficacy of this approach on synthetic data, and develop an algorithm to identify the occurrence of a spatial and/or temporal anomalous event in scientific phenomena. The algorithm decomposes the data into several spatial sub-domains and time steps to identify regions with such events. Feature moment metrics, based on the alignments of the principal kurtosis vectors, are computed at each sub-domain and time step for all features to quantify their relative importance towards the overall kurtosis in the data. Accordingly, spatial and temporal anomaly metrics for each sub-domain are proposed using the Hellinger distance of the feature moment metric distribution from a suitable nominal distribution. Finally, we apply the algorithm to two turbulent auto-ignition combustion cases and demonstrate that the anomaly metrics reliably capture the occurrence of auto-ignition in relevant spatial sub-domains at the right time steps.
Here, this work introduces a new predictive multi-zone model for the description of combustion in Homogeneous Charge Compression Ignition (HCCI) engines. The model exploits the existing OpenSMOKE++ computational suite to handle detailed kinetic mechanisms, providing reliable predictions of the in-cylinder auto-ignition processes. All the elements with a significant impact on the combustion performances and emissions, like turbulence, heat and mass exchanges, crevices, residual burned gases, thermal and feed stratification are taken into account. Compared to other computational approaches, this model improves the description of mixture stratification phenomena by coupling a wall heat transfer model derived from CFD application with a proper turbulence model. Furthermore, the calibration of this multi-zone model requires only three parameters, which can be derived from a non-reactive CFD simulation: these adaptive variables depend only on the engine geometry and remain fixed across a wide range of operating conditions, allowing the prediction of auto-ignition, pressure traces and pollutants. This computational framework enables the use of detail kinetic mechanisms, as well as Rate of Production Analysis (RoPA) and Sensitivity Analysis (SA) to investigate the complex chemistry involved in the auto-ignition and the pollutants formation processes. In the final sections of the paper, these capabilities are demonstrated through the comparison with experimental data.
To meet current emission regulations and increasingly demanding global fleet CO2 standards on fuel economy and future trends towards life cycle GHG emissions, advanced combustion engines remain significant in the passenger vehicle sector to achieve high efficiency and low emissions over the full operating range. Lean burn gasoline compression ignition (GCI) technology has shown to have the most potential in reaching these goals, although it faces challenges in the operating range. In this study, an advanced GCI engine is considered with the capability to operate under two combustion modes, namely low temperature combustion (LTC) and GCI. They are enabled with the use of two Hyundai in-house developed technologies; an advanced valve control mechanism known as continuously variable valve duration (CVVD) and a high-pressure gasoline injection system. At low load, the engine utilizes dual CVVD and dual CVVT (continuously variable valve timing) mechanisms for both intake and exhaust valvetrains to enable NVO (negative valve overlap) to trap hot residuals. The hot residuals enable low load auto-ignition, and the phasing control of that auto-ignition is achieved by varying the degrees of NVO, and in turn, the amount of hot residuals. This is done in real time with respect to the changes in engine loads or operating conditions, significantly improving cylinder-phasing control in LTC mode, and subsequently is able to realize a fuel economy gain and reduced emission benefits. Early direct injection during NVO or intake stroke is used to form homogenous or quasi-homogenous air-fuel mixture for LTC. At mid to high loads, late direct injection of gasoline or GCI mode is employed to take advantage of the high volatility and high resistance to auto-ignition of gasoline fuel, so that diesel-engine like high combustion efficiency can be achieved with lower soot emissions for the same NOx level as diesel combustion. In this paper, the key control and fueling technologies that enable this multiple mode combustion are introduced, and the resultant analysis on engine test and CFD simulation at LTC and GCI operating conditions are presented in detail.
To this date, extensive research has been conducted to understand the low-temperature auto-ignition chemistry of gasoline. The detection of low-temperature chemical reactions under Spark Ignition (SI) combustion cannot be detected, as they are hidden by the flame propagation. Alternatively, Homogeneous Charge Compression Ignition (HCCI) combustion has a two-stage combustion involving low and high-temperature heat release (LTHR and HTHR respectively). Both Knocking SI and HCCI combustion involve auto-ignition and are governed by fuel characteristics and the pressure-temperature (P-T) history. Therefore, HCCI combustion might provide an alternative to understand the knocking behavior and LTHR in modern SI engines. A standard Cooperative Fuel Research (CFR) engine was operated at lean HCCI conditions (lambda 3), as well as SI conditions at stoichiometry. For SI combustion, the CFR engine was operated with RON-like conditions, but at late spark timing to induce LTHR prior to flame propagation. Three RON 90 binary fuel blends were investigated, being composed of n-heptane with isooctane, toluene, or ethanol. This work demonstrated that the CFR engine under stoichiometric SI with late spark timing and HCCI combustion mode can help to detect LTHR which is not possible in the standard RON test. The intake pressure and temperature sweeps showed similar effects on LTHR for both combustion modes. The linking of auto-ignition behavior of SI and HCCI was dependent primarily on intake valve closing (IVC) conditions. In conclusion, the high exhaust temperature in SI lead to high IVC temperatures. In order to match the IVC temperatures and to overlap the P-T trajectories, the intake temperature for HCCI was increased.
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.
Decades of work on low-carbon fuels have established their potential for substantial emissions reductions; however, their adoption is still limited by infrastructure concerns and engine efficiency deficits. As infrastructures have begun to evolve, research on strategies that maximize engine efficiency through interactions with fuel properties must also now take center place. This paper compares the performance, emissions, and combustion characteristics of two forefront low-carbon fuels: compressed natural gas (CNG) and liquefied petroleum gas (LPG) in a cooperative fuel research (CFR) engine over a range of compression ratios and engine loads. The effects of exhaust gas recirculation (EGR), end-gas auto-ignition, and a novel combustion control tool, the combustion intensity metric (CIM), were also evaluated at different stoichiometric engine operating conditions. In comparison to LPG, CNG operation demonstrated an extended knock-free regime, allowing engine operation at higher engine loads and compression ratios, but LPG operation exhibited enhanced combustion characteristics with higher peak pressures and faster apparent heat release rates (AHRR). LPG operation achieved higher brake thermal efficiencies and lower equivalent CO 2 emissions compared to CNG operation at the tested engine loads and compression ratios. LPG demonstrated significantly higher EGR tolerance limits compared to CNG, with a maximum of 28% EGR rate, compared to 23% for CNG. This improved EGR dilution tolerance was responsible for a 90% reduction in NOx emissions for LPG compared to a maximum of 70% with CNG. EGR dilution also exhibited more effective knock mitigation potential with LPG, suppressing knock intensity values by up to 98% and transitioning the engine operation towards normal combustion from heavy knocking conditions. As a result, the CIM was found to decrease burn durations and improve the quality of combustion by controlling the desired fraction of end-gas auto-ignition.
An n-dodecane spray flame (Spray A from Engine Combustion Network) was simulated using a detailed combustion model along with a dynamic structure LES model to evaluate its performance at engine-relevant conditions and understand the transient behavior of this turbulent flame. The liquid spray was treated with a traditional Lagrangian method and the gas-phase reaction was modeled using a detailed combustion model. A 103-species skeletal mechanism was used for the n-dodecane chemical kinetic model. Significantly different flame structures and ignition processes are observed for the LES compared to those of RANS predictions. The LES data suggests that the first ignition initiates in lean mixture and propagates to rich mixture, and the main ignition happens in rich mixture, preferable less than 0.14 in mixture fraction space. LES was observed to have multiple ignition spots in the mixing layer simultaneously while the main ignition initiates in a clearly asymmetric fashion. The temporal flame development also indicates the flame stabilization mechanism is auto-ignition controlled and modulated by flame propagation. Soot predictions by LES present much better agreement with experiments compared to RANS both qualitatively and quantitatively. Multiple realizations for LES were performed to understand the realization to realization variation and to establish best practices for ensemble-averaging diesel spray flames. The relevance index analysis suggests that an average of 2 and 5 realizations can reach 99\% of similarity to the target average of 16 realizations on the temperature and mixture fraction fields, respectively. However, more realizations are necessary for OH and soot mass fraction due to their high fluctuations.
The transportation sector accounts for almost a third of the United States energy consumption. We are developing predictive simulations of the complex in-cylinder processes of internal combustion engines to improve performance and decrease pollutant emissions. Using volume rendering, we provide a first look at the jet impingement and ignition inside a piston-cylinder chamber. Through the visualization, we were able to construct a complete 3D picture of the jet impingement on the piston-bowl walls leading to the formation of pockets of auto-ignition of fuel-air mixture. These auto-ignition pockets subsequently lead to flame stabilization releasing the chemical energy for conversion to useful work. A detailed physical picture of flame stabilization as provided by the visualization techniques presented in this study help in the design next-generation internal combustion engines.