Search NASA⌕ Search

SEARCH · Search NASA

Results for “ignition modeling”

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 19 records

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗

Application of a Comprehensive Lagrangian–Eulerian Spark-Ignition Model to Different Operating Conditions

Increasing engine efficiency is essential to reducing emissions, which is a priority for automakers. Unconventional modes such as boosted and highly dilute operation have the potential to increase engine efficiency but suffer from stability concerns and cyclic variability. To aid engineers in designing ignition systems that reduce cyclic variability in such engine operation modes, reliable and accurate spark-ignition models are necessary. Here, in this article, a Lagrangian–Eulerian spark-ignition (LESI) model is used to simulate electrical discharge, spark channel elongation, and ignition in inert or reacting crossflow within a combustion vessel, at different pressures, flow speeds, and dilution rates. First the model formulation is briefly revisited. Then, the experimental and simulations setups are presented. The results showcase the model’s ability to predict the secondary circuit voltage, current, and power signals, in addition to the spark channel elongation, for the inert cases, or flame front growth, for the reacting cases. The results also compare simulation spark channel and flame growth plots to experimental Schlieren images at different instants in time. This work serves to highlight LESI’s ability to predict the characteristics of discharge and ignition across a variety of operating conditions.

42 ENGINEERING↗

A novel ignition model for low velocity impact of heterogeneous explosives based on interacting hot spots

While numerous studies have focused on the ignition of explosives occurring in high velocity impact and the associated shock-to-detonation transition, there has been growing interest in developing computational models focused on low-velocity impact regimes. A predictive low-velocity impact ignition model will be important for analyzing high explosive safety and potential accident scenarios. This work introduces a novel ignition model based on the concept of thermally interacting hot spots to simulate low velocity impacted heterogeneous explosives where observed ignition times are on the order of milliseconds. The model asserts that relevant hot spots are micron-sized, the typical separation between neighboring hot spots is on the order of a hundred microns, and that neighbors interact thermally through heat conduction across the interstitial region between them. To achieve tractable numerical solutions, hot spots are assumed to form a periodic array as opposed to the highly irregular positioning in an actual explosive. This idealization allows a single two hotspot system to characterize the ignition process. Consequently, the model is referred to as the two hot spot Frank-Kamenetskii ignition model. In the present study, hot spots are modeled as constant heat sources terms, but this can be extended to include grain-scale phenomena like frictional heating of micron-sized growing cracks that are confined under high pressure. Because the micron-sized features are below the scale that can be efficiently resolved at a systems level, an efficient subscale scheme based on the Method of Weighted Residuals (MWR) is used to efficiently solve the equations. In conclusion, we carry out numerical examples and analytic predictions illustrating the accuracy and the functioning of the model.

97 MATHEMATICS AND COMPUTING↗

Commercialization of a Comprehensive Spark-Ignition Model for Automotive Engine Applications in CONVERGE CFD

The work completed under this CRADA commercializes an advanced hybrid Lagrangian-Eulerian spark-ignition (LESI) model for multi-dimensional engine simulations. In this report, LESI is coupled with combustion models (WSR, TFM, G-equation) to model spark channel elongation and kernel and flame development in a constant volume combustion vessel, and a four-valve direct-injection spark-ignition (4VDISI) engine in two operating conditions: stoichiometric E30 with no dilution, and stoichiometric PACE-20 with dilution using CONVERGE CFD solver.

42 ENGINEERING↗

Catalytic ignition model in a monolithic reactor with in-depth reaction

Two transient models have been developed to study the catalytic ignition in a monolithic catalytic reactor. The special feature in these models is the inclusion of thermal and species structures in the porous catalytic layer. There are many time scales involved in the catalytic ignition problem, and these two models are developed with different time scales. In the full transient model, the equations are non-dimensionalized by the shortest time scale (mass diffusion across the catalytic layer). It is therefore accurate but is computationally costly. In the energy-integral model, only the slowest process (solid heat-up) is taken as nonsteady. It is approximate but computationally efficient. In the computations performed, the catalyst is platinum and the reactants are rich mixtures of hydrogen and oxygen. One-step global chemical reaction rates are used for both gas-phase homogeneous reaction and catalytic heterogeneous reaction. The computed results reveal the transient ignition processes in detail, including the structure variation with time in the reactive catalytic layer. An ignition map using reactor length and catalyst loading is constructed. The comparison of computed results between the two transient models verifies the applicability of the energy-integral model when the time is greater than the second largest time scale of the system. It also suggests that a proper combined use of the two models can catch all the transient phenomena while minimizing the computational cost.

Tien, Ta-Ching↗

Catalytic ignition model in monolithic reactor with in-depth reaction

A transient model has been developed to study ignition in a monolithic catalytic reactor. Special features are the inclusion of thermal and species structures in the porous catalytic layer and the resolution of small-time-scale events. Catalyst loading is used as a parameter to study these structure variations before and after ignition. Regions where the catalytic reaction is kinetically controlled and regions where the reaction is diffusionally controlled are identified in both steady and transient states.

Tien, Ta-Ching↗

Modelling ignition characteristics of rich H2/O2 mixture in a monolithic catalytic reactor

Ignition of rich hydrogen-oxygen mixture in a monolithic catalytic reactor is studied using a transient combustion model. The model assumes a quasi-steady gas phase and a thermally-thin substrate. The ignition time lag is due to the thermal inertia of the substrate. One-step global chemical reaction is assumed both on the surface and in the gas phase. An effectiveness factor is introduced to account for the transition from a kinetically-limited catalytic surface reaction to a diffusion-limit one during ignition. Results presented include a catalytic ignition boundary, ignition delay time and the transient response in the catalytic bed.

Tien, Ta-Ching↗

LOX/Methane Main Engine Glow Plug Igniter Tests and Modeling

Ignition data for tests with a LOX/methane igniter that utilized a glow plug as the ignition source are presented. The tests were conducted in a vacuum can with thermally conditioned (cold) hardware. Data showing the effects of glow plug geometry, type, and igniter operating conditions are discussed. Comparisons between experimental results and multidimensional, transient computer models are also made.

Breisacher, Kevin↗

Coupling a Lagrangian–Eulerian Spark-Ignition (LESI) model with LES combustion models for engine simulations

In the United States transportation sector, Light-Duty Vehicles (LDVs) are the largest energy consumers and CO 2 emitters. Electrification of LDVs is posed as a potential solution, but SI engines can still contribute to decarbonization. Car manufacturers have turned to unconventional engine operation to increase the efficiency of Spark-Ignition (SI) engines and reduce the carbon emissions of their fleets. Dilute, lean, and stratified-charge engine operation has the potential for engine efficiency improvements at the expense of increased cyclic variability and combustion instability. At such demanding engine conditions, the spark ignition event is key for flame initiation and propagation and for enhanced combustion stability. Reliable and accurate spark ignition models can help design ignition systems that reduce cyclic variability. Multiple computational spark-ignition models exist that perform well under conventional conditions, but the underlying physics needs to be expanded, for unconventional engine operation. In this paper, a hybrid Lagrangian–Eulerian Spark-Ignition (LESI) model is coupled with different turbulent flame propagation models for engine simulations. LESI relies on Lagrangian arc tracking and Eulerian energy deposition. The LESI model is coupled with the Well-Stirred Reactor (WSR), Thickened Flame Model (TFM), and g-equation model and used to simulate several cycles of a Direct-Injection Spark-Ignition (DISI) engine using a commercial Computational Fluid Dynamics (CFD) engine solver. The results showcase the successful coupling of LESI with the combustion models. Global engine metrics, such as pressure and Apparent Heat Release Rate (AHRR), for each simulation setup are compared to experimental engine results, for validation. In addition, results highlight the successful prediction of spark channel movement by comparing simulation images to experimental optical engine images. Finally, the successful coupling of LESI to combustion models, making it a usable model in the engine modeling community, is emphasized and future development details are discussed.

33 ADVANCED PROPULSION SYSTEMS↗

Heat transfer model for predicting squib ignition times

A squib ignition model based on transient heat condition from the hot bridgewire to the pyrotechnic is described. No Arrhenius-type chemical reaction is included. Instead, a thermal contact resistance is postulated to exist between the hot bridgewire and the pyrotechnic. Ignition is assumed to occur when a 2.5 micron layer of pyrotechnic next to the bridgewire reaches a characteristic ignition temperature for that pyrotechnic. This model was applied to the JPL squib, which uses a 50 micron (0.002-in.) diameter Tophet A bridgewire to ignite a boron, potassium perchlorate mix. A computer program was utilized that solves the transient heat condition problem with the boundary conditions stipulated by the model. The thermal contact conductance at the interface was determined by trial and error so that the experimentally determined ignition time for one firing condition would be properly predicted by the model. The agreement was quite good for tests run between -129 C and +93.3 C at current levels of 3.5 and 5 A. Axial heat conduction along the bridgewire is shown to be negligible.

Sernas, V.↗

Effects of spatiotemporal plasma power distribution on the modeling of ignition kernel evolution in quiescent and turbulent methane/air mixtures

Abstract The present work improves a phenomenological plasma-assisted combustion model by integrating the spatiotemporal distribution of plasma power density, thereby considering the evolution of plasma streamers in the modeling, and subsequently, better predicting the ignition kernel evolution. The improved phenomenological model is validated against experiments representing the plasma discharge and post-discharge ignition kernel evolution. Specifically, the new model demonstrates a more accurate prediction of ultrafast gas heating and O 2 dissociation during the plasma discharge, compared to the original model. In addition, the new model is found to closely match the experimental pressure wave and heated channel profiles post-discharge without the need for tuning the energy deposition (unlike the original model), highlighting its accuracy of post-discharge ignition kernel dynamics. The improved phenomenological model is then employed to investigate ignition kernel evolution for a stoichiometric methane-air discharge across various discharge gap configurations. Simulations reveal a non-uniform temperature and streamer distribution progressing from the electrode tips toward the center, contrasting uniform cylindrical discharges previously described in the original model. Streamer propagation is observed to be faster for larger gaps when maintained at the same average electric field for different discharge gaps. The tendency of smaller gaps to produce detached toroidal ignition kernels is observed, while larger gaps promote cylindrical and attached ignition kernels. Interactions between successive ignition kernels from consecutive discharges varied significantly, with the smallest gap (1 mm) promoting the quenching of the preceding ignition kernel due to the initial kernel–kernel separation. The intermediate gap (2 mm) promotes detached kernel growth. In contrast, in the largest gap (4 mm), kernels consistently combine and expand attached to electrodes. The impact of homogeneous isotropic turbulence is also explored, showing the persistence of ignition kernels early on but eventually quenching due to enhanced radical and heat losses with pronounced turbulence intensity.

Johnson, Praise Noah↗

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗

A comprehensive model to capture electrical discharge and spark channel evolution during spark-ignition processes

Emissions reduction through engine efficiency improvements is a priority for automakers who have turned to unconventional engine operation such as highly dilute, boosted, and stratified charge. Given its importance to flame initiation and sustained turbulent flame propagation, reliable and accurate spark ignition models are necessary to design ignition systems that reduce cyclic indicated mean effective pressure (IMEP) variability and increase engine efficiency in these operation modes. In this paper, secondary electric circuit, short-circuit, blowout, and re-strike sub-models are added to the Lagrangian-Eulerian spark ignition (LESI) model to simulate electrical discharge and spark channel elongation in an inert cross-flow combustion vessel. First, the physics of the underlying sub-models are described and the governing equations discussed with the spark channel voltage expression playing a critical role. Then, the experimental and simulations setups are presented. The results section begins with the derivation of spark channel voltage from experimental results. Then, electrical discharge and spark channel elongation simulations in inert flows are carried out using LESI and compared against experimental results. Finally, the results validate the model's ability to accurately predict spark channel elongation, as well as the occurrence of short-circuits, blowouts, re-strikes, and end of discharge.

42 ENGINEERING↗

Comparing the shock sensitivity of insensitive energetic materials

We present a tabletop method to study the shock sensitivity of plastic-bonded explosives that are considered shock insensitive using high dynamic range optical emission spectroscopy with laser-launched km/s flyer plate impacts (2, 3, and 4 km/s), which measures the spectral radiance (the emission spectrum vs a calibrated intensity standard) every 0.8 ns in the nanosecond and microsecond regimes. The explosives were TATB (1,3,5 trinitro, 2,4,6 triamino benzene), FOX-7 (1,1-diamino-2,2-dinitroethylene), LLM-105 (2,6-diamino-3,5-dinitropyrazine-1-oxide), and NTO (nitrotriazolone), all with 20% Sylgard 182 binder. Time resolving emission from shocked explosives produces a unique fingerprint at each shock pressure, and the emission temporal profile can be used to understand each material's shock response and the underlying mechanisms within the framework of the well-known shock ignition and growth models. Ignition was characterized by the emission intensity during nanosecond hot spot ignition. The growth process, which occurs after the shock has unloaded, was characterized by the time derivative of microsecond emission intensity. Singular-value decomposition was used to determine a characteristic spectral radiance over the entire time range, and this radiance was fit to a graybody model to determine a characteristic temperature Tc. Even though these insensitive explosives have similar shock sensitivities, our method can clearly distinguish their shock response. FOX-7 and NTO were the most sensitive, with FOX-7 having the greatest shock ignition sensitivity, growth rate, and Tc. LLM-105 was much less sensitive than FOX-7 and NTO, and TATB consistently showed the lowest sensitivity by any measure.

Akhtar, Meysam (ORCID:0000000243897837)↗

CFD Evaluation of Lean-Direct Injection Combustors for Commercial Supersonics Technology

An overview is given of an effort that focused on using CFD analysis to complement design and configuration definition of Lean-Direct Injection (LDI) combustion concepts for NASA's Commercial Supersonic Transport (CST) program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for second and third generation LDI configurations at CST cruise conditions. All computations were performed with a consistent approach of mesh-generation, spray modeling, ignition and kinetics modeling. Emissions (EINOx) characteristics were predicted for CST cruise conditions, and compared with emissions data from experimental measurements to evaluate the fidelity of the CFD modeling approach to predict emissions changes in response to changes in supersonic cycle conditions.

nitrogen oxides injectors↗

Pilot Injector Redesign to Reduce N+3 Cycle Emissions for a Gas-Turbine Combustor

An overview is given of an effort for the use of CFD analysis to complement design and configuration definition of third generation Lean-Direct Injection combustion concepts (LDI-3) for NASA’s N+3 program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for a three-cup, nineteen-element flame tube array with redesigned pilot injectors to improve spray and emissions characteristics when compared to a previous LDI-3 design. All computations were performed with a consistent approach to mesh-generation, spray modeling, ignition and kinetics modeling for a ‘medium-power’ cycle condition. Computational predictions of the aerodynamics of a new pre-filming pilot injector were used to arrive at an optimized aerothermal design that meets effective area and fuel-air mixing criteria. The newly designed pilot injectors were shown to provide considerable improvements in aerodynamic stability, flame-tube pattern factor and NOx emissions, when compared to the original design.

nitrogen oxides injectors↗

Pilot Injector Redesign to Reduce N+3 Cycle Emissions for a Gas-Turbine Combustor

An overview is given of an effort for the use of CFD analysis to complement design and configuration definition of third generation Lean-Direct Injection combustion concepts (LDI-3) for NASA's N plus 3 program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for a three-cup, nineteen-element flametube array with redesigned pilot injectors to improve spray and emissions characteristics when compared to a previous LDI-3 design. All computations were performed with a consistent approach to mesh-generation, spray modeling, ignition and kinetics modeling for a "medium-power" cycle condition. Computational predictions of the aerodynamics of a new pre-filming pilot injector were used to arrive at an optimized aerothermal design that meets effective area and fuel-air mixing criteria. The newly designed pilot injectors were shown to provide considerable improvements in aerodynamic stability, flame-tube pattern factor and NOx emissions, when compared to the original design.

Nitrogen Oxides↗

CFD Evaluation Of Lean-Direct Injection Combustors for Commercial Supersonics Technology

An overview is given of an effort that focused on using CFD analysis to complement design and configuration definition of Lean-Direct Injection (LDI) combustion concepts for NASA's Commercial Supersonic Transport (CST) program. The National Combustion Code (OpenNCC) was used to perform non-reacting and two-phase reacting flow computations for second and third generation LDI configurations at CST cruise conditions. All computations were performed with a consistent approach of mesh-generation, spray modeling, ignition and kinetics modeling. Emissions (EINOx) characteristics were predicted for CST cruise conditions, and compared with emissions data from experimental measurements to evaluate the fidelity of the CFD modeling approach to predict emissions changes in response to changes in supersonic cycle conditions.

Ajmani, Kumud↗