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.
The on-road natural-gas (NG) fueled transportation relies on stoichiometric spark-ignition engines for the advantages of simple after-treatment system despite the efficiency penalty relative to lean-burn combustion strategies. Exhaust gas recirculation (EGR) has the potential to reduce this efficiency gap at low to moderate loads without the need for complex lean-exhaust aftertreatment systems. However, EGR dilution leads to reduced combustion stability and increased cycle-to-cycle variability. A promising technology that has the potential to achieve reliable operation under diluted conditions is the prechamber ignition (or turbulent jet ignition) which uses chemically active turbulent jets generated from combustion inside a prechamber to initiate, stabilize and accelerate combustion of the mixture inside the main chamber. The present work focusses on developing a RANS-based CFD approach to accurately reproduce in-cylinder phenomena in a stoichiometric NG prechamber-assisted heavy-duty engine without relying on complex combustion models that account for turbulence-chemistry interactions. This is necessary because reactive prechamber jets at high EGR dilution tend to extinguish while emerging into the main chamber, which is followed by a phase of re-ignition — a phenomenon that conventional G-equation or well-stirred reactor combustion models cannot reproduce. With addition of a damping multiplier to the well-stirred reactor model, the predictions are seen to show good agreement with experimental pressure evolution and combustion images acquired from a single cylinder Cummins N-14 optical diesel engine retrofitted with a prechamber ignition system. Model predictions of local heat release in the flame and temperature evolution inside the flame are used to investigate combustion dynamics in the prechamber and the main chamber. It is seen that the well-stirred reactor model with the inclusion of damping is able to reproduce the temporary reduction in heat release within the flame, which can be considered equivalent to quenching of jets, and the subsequent re-ignition of the flame inside the main chamber. The delay between quenching and re-ignition depends on the amount of dilution, as explained by an illustration of flame evolution in a Borghi diagram.
Equation of state (EOS) measurements in the 100 s Mbar range are needed to underwrite models employed in the simulation of high energy density plasmas. To this end, a platform has been developed for fielding on the National Ignition Facility, capable of producing high-quality impedance match EOS data, wherein a planar, high-pressure, steady shock is driven into a sample package, and sample and reference standard shock velocities are measured. This platform, dubbed planar high pressure, or PHP, was fielded with an initial proof-of-concept shot in January 2023. The first PHP shot, aiming to study gold, demonstrated a pressure close to 400 Mbar, two orders of magnitude higher than previously reported gold EOS data.
Access to reliable, clean energy sources is a major concern for national security. Much research is focused on the “grand challenge” of producing energy via controlled fusion reactions in a laboratory setting. For fusion experiments, specifically inertial confinement fusion (ICF), to produce sufficient energy, the fusion reactions in the ICF fuel need to become self-sustaining and burn deuterium-tritium (DT) fuel efficiently. The recent record-breaking NIF ignition shot was able to achieve this goal as well as produce more energy than used to drive the experiment. This achievement brings self-sustaining fusion-based power systems closer than ever before, capable of providing humans with access to secure, renewable energy. In order to further progress toward the actualization of such power systems, more ICF experiments need to be conducted at large laser facilities such as the United States's National Ignition Facility (NIF) or France's Laser Mega-Joule. The high cost per shot and limited number of shots that are possible per year make it prohibitive to perform large numbers of experiments. As such, experimental design relies heavily on complex predictive physics simulations for high-fidelity “preshot” analysis. These multidimensional, multi-physics, high-fidelity simulations have to account for a variety of input parameters as well as modeling the extreme conditions (pressures and densities) present at ignition. Such simulations (especially in 3D) can become computationally prohibitive to turn around for each ICF experiment. In this work, we explore using Bayesian optimization with Gaussian processes (GPs) to find optimal designs for ICF double shell targets, while keeping computational costs to manageable levels. These double shell targets have an inner shell that grades from beryllium on the outer surface to the higher Z material molybdenum, as opposed to the nominally used tungsten, on the inside in order to trade off between the high performance associated with high density inner shells and capsule stability. We describe our results for “capsule-only” xRAGE simulations to study the physics between different capsule designs, inner shell materials, and potential for future experiments.
Mountain ecosystems typically serve as carbon (C) sinks. However, studies also suggest that they could be C sources due to climate warming, drought and insect-related mortality, wildfires, and management actions. We applied the Lund-Potsdam-Jena General Ecosystem Simulator (LPJ-GUESS), a process-based dynamic vegetation-ecosystem model, to investigate the role of ecosystem management in C storage under Mediterranean climate over the 21 st century. We modified LPJ-GUESS to include implementing mechanical thinning by vegetation size classes, components, and types along with a new mechanistic fire-occurrence model that accounts for wind speed and lightning ignition. Simulations show that mechanical thinning or prescribed fire performed 5-20 years in advance of a high-severity wildfire reduced direct wildfire C emissions by 38-66 %. Our results also show that long-term management actions repeated every 5-20 years, including thinning relatively small trees (diameters up to 7 inches or ∼178 mm), can maintain stable C levels in the forest and lower dead-fuel amounts. We found that, although prescribed fire mitigated wildfire severity, ecosystem C storage from reduced wildfire emissions can be outweighed by the added emissions from the prescribed fire themselves. Thinning plus removing and sequestering the thinned biomass can ensure that forests act as net C sinks through the end of the 21 st century. However, addition of prescribed fire is needed to reduce understory and lower the projected extent of high-severity wildfire. Achieving the competing goals of reducing wildfire and making the Sierra Nevada long-term C sink can be advanced through carefully coordinated thinning, sequestration of thinned biomass, and prescribed fire.
Inertial confinement fusion implosions for power production need to be robust. Degradation sources need to be understood and considered for an economic power plant. In this work, the impact of low-mode asymmetries on yield is studied using 2D radiation hydrodynamics simulations for inertial fusion energy relevant implosions. First, a high-gain indirectly driven capsule design is developed. This design is subsequently degraded with low-mode drive asymmetries. The results show that yield is constant with increasing amounts of asymmetry until a cliff is reached where the yield quickly drops below ignition. Robustness arises because the degraded hotspot still has enough time to get to ignition before capsule decompression. This robustness occurs even though the asymmetry couples less energy to the hotspot and causes decompression to occur sooner. The robustness is linked to the generalized Lawson parameter, χ. The value of χ steadily decreases with increasing asymmetry. Once χ is close to one, the hotspot is no longer able to ignite before the capsule decompresses. A simple piston model is used to predict where the robustness cliff is for any design. The trade-off between yield and robustness is also demonstrated.
The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.
Context.Nuclear networks are widely used coupled with hydrodynamical simulations of explosive scenarios to account for the change of nuclear species and energy generation rate due to nuclear reactions. In this way, there is a feedback mechanism between the hydrodynamical state and the nuclear processes. Unfortunately, the timescale of nuclear reactions is orders of magnitude smaller than the dynamical timescale that drives hydrodynamical simulations. Therefore, these nuclear networks are usually very small, reduced in most cases to a dozen elements, especially when simulations are carried out in more than one dimension. Aims.We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. Methods.We perform several tests that are relevant to simulate explosive scenarios, such as Type Ia supernovae and core-collapse supernovae. We compare the results of the 90 nuclei network with a standardα-chain network with 14 elements to evaluate the differences in the energy generation rate. We also evaluate the relevance of including the electrons in the network in terms of generated yields and how it affects the pressure of a degenerate fluid such as that of white dwarfs. The results obtained with the 90-nuclei network have been verified with a much larger 2000-nuclei network built from REACLIB (WinNet), in terms of nuclear energy generation rate, pressure, and produced yields. Results.The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed withWinNetin scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer toWinNetthan typical reduced alpha networks. Conclusions.A nuclear network withN= 90 species, including electrons, aimed at multidimensional calculations of supernova explosions is described and verified. The proposed network is suitable for the study of Type Ia supernova explosions because it provides better values of pressure and electron abundance than other existing networks with smaller or even a similar size but without including electron capture processes.
There is an imminent need to displace fossil diesel fuel with cleaner burning, domestically produced, renewable fuels for use in heavy-duty engines. Bioethanol is a prime candidate as it widely adopted in the U.S. as a gasoline additive ranging in volume percentage from 10% (E10) up to 85% (E85). Direct substitution of market available ethanol-gasoline blends for diesel fuel is not plausible as the stark reactivity differences would not constitute the same ignition quality nor achieve auto-ignition at all. This work focuses on the development of prechamber enabled mixing-controlled combustion (PC-MCC) as an advanced combustion strategy to facilitate reliable ignition and diffusion style combustion ethanol-gasoline fuel blends. PC-MCC involves integration of an actively fueled prechamber (PC) into a conventional compression ignition combustion system. When ignited, the PC ejects hot turbulent jets into the main combustion chamber that then interact with the direct injected fuel, prompting immediate ignition. The PC jet flames provide a robust thermal ignition source that allows the engine to operate agnostic of fuel composition, or flex-fuel. Computational fluid dynamics (CFD) modeling was used to assess critical design features of the PC while garnering insights into the ignition strategies that facilitate robust performance. A key finding was the ignition performance benefits of fuel-rich PC operation which yield exothermic jets. Based on the numerical findings, a prototype igniter was tested experimentally on both single and multi-cylinder engine platforms at a variety of operating conditions. The experimental results indicate flex-fuel PC-MCC is well capable of diesel-like combustion processes by demonstrating matched or improved gross thermal efficiencies and load variability within 2%. Fuel grade ethanol (E98) exhibited consistently lower NOx and immeasurable soot across the load space. E98 also demonstrated a significant improvement in thermal efficiency at light loads.
We conduct capsule surrogate experiments at the National Ignition Facility to calibrate radiation hydrodynamic simulations to infer hydrodynamic conditions that are not observable in indirect drive ignition implosions. We tune the simulations by applying laser power and cross beam energy transfer (CBET) saturation multipliers to match the observables from capsule surrogate experiments. Shock timing, velocity, and symmetry are measured in liquid D 2 filled Keyhole capsule surrogate experiments and implosion trajectory, stagnation time, and shape time history are measured in in-flight 2D backlit x-ray radiography experiments (“2DConA”) of D 2 gas filled capsule implosions. Calibrated simulations suggest that the N210808 ignition implosion (fusion target gain = 0.7) had a shell mass remaining at stagnation of less than the nominal %5 (3.8%) and resulted in less confinement. For N221204, the shell was made 5.75 μm thicker to trade implosion velocity for increased confinement and resulted in a target gain = 1.5 with a shell mass remaining of 5.7%. Furthermore, a single adjusted model can reproduce all shock timing data as changes are made to shell thickness (79–85 μm) and laser wavelength separation (1.8–4.0 Å). However, for the 2DConA implosions, a 5% variation in the peak power laser multipliers and a 30% variation in late-time CBET between experiments are needed to match the observed stagnation times, in-flight $P_2$ shape, and hot-spot $P_2$ shape. While progress is being made to improve the models in simulations using focused experiments, capsule surrogate experiments will continue to be needed to optimize future ignition designs.
This study investigated the behavior of gasoline direct injection (GDI) sprays using computational fluid dynamics (CFD). The authors developed an approach to identify optimal spray breakup model parameters by evaluating an error function across numerous simulations, with the goal of minimizing discrepancies from experimental data. Using the optimal setup, the simulated spray matched well with projected liquid volume distributions, liquid penetration, and spray width measured in a constant-pressure continuous-flow chamber. To further validate the approach, the same setup was tested across various fuels, injectors, and operating conditions. Subsequently, the optimal setup, along with a recently developed spray-wall interaction model, were applied to a direct-injection spark-ignited engine under late-injection conditions to predict and evaluate fuel film formation and evolution at varying engine coolant temperatures. Here, with the centrally mounted injector directing the spray toward the piston, simulations indicated that the spray tends to impinge on the piston surface. The proposed simulation framework also accurately captured the aggregate film area on the piston surface, aligning with previously published experimental results. Moreover, simulations showed that increasing the coolant temperature from cold start conditions (333 K) to warm conditions (363 K) reduced the fuel mass deposited on the piston by roughly 50%. Furthermore, for the spray-guided engine configuration studied in this work, the CFD model predicted minimal film deposition on the spark plug electrodes regardless of the coolant temperatures due to a relatively weak in-cylinder flow during the compression phase.
The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.
Bayesian inference applied to x-ray spectroscopy data analysis enables uncertainty quantification necessary to rigorously test theoretical models. However, when comparing to data, detailed atomic physics and radiation transfer calculations of x-ray emission from non-uniform plasma conditions are typically too slow to be performed in line with statistical sampling methods, such as Markov Chain Monte Carlo sampling. Furthermore, differences in transition energies and x-ray opacities often make direct comparisons between simulated and measured spectra unreliable. Here, we present a spectral decomposition method that allows for corrections to line positions and bound–bound opacities to best fit experimental data, with the goal of providing quantitative feedback to improve the underlying theoretical models and guide future experiments. In this work, we use a neural network (NN) surrogate model to replace spectral calculations of isobaric hot-spots created in Kr-doped implosions at the National Ignition Facility. The NN was trained on calculations of x-ray spectra using an isobaric hot-spot model post-processed with Cretin, a multi-species atomic kinetics and radiation code. The speedup provided by the NN model to generate x-ray emission spectra enables statistical analysis of parameterized models with sufficient detail to accurately represent the physical system and extract the plasma parameters of interest.
There is interest in adapting renewable and low-carbon intensity fuels to heavy-duty engines to help displace criteria and greenhouse gas emissions associated with diesel combustion. Low-carbon fuels have inherently low cetane numbers and thus pose a significant challenge when considering direct substitution in diesel combustion systems. Use of an actively fueled prechamber as an ignition source to initiate mixing-controlled combustion (PC-MCC) of low-cetane fuels is becoming an attractive combustion mode to alleviate the identified reactivity deficit. This work focuses on the utilization of an active prechamber to facilitate diesel-like combustion of bioethanol-gasoline fuel blends in a heavy-duty engine. Recent results in this emerging technology suggest that ignition quality is uniquely coupled to prechamber equivalence ratio. In this numerical study, the fundamental implications of prechamber equivalence ratio on the ignition performance of direct injected fuels ranging from E10 to E100 are investigated using CONVERGE. Parametric studies of the prechamber operating strategy were assessed at diesel-like conditions to characterize the performance trends relative to a diesel baseline at the same boundary conditions. Simulation results indicate that PC-MCC is flex-fuel capable and achieves diesel-like ignition qualities and combustion processes for all fuels considered under stoichiometric and rich prechamber conditions. To characterize the equivalence ratio trends observed, a novel inflow boundary modeling technique was utilized to prescribe turbulent jets in place of the prechamber where selective speciation of the jet composition was conducted to isolate the ignition contributions of radical and reacting species in lean, stoichiometric, and rich jets. The inflow boundary modeling suggests that excess fuel and combustion intermediates present in jets produced from stratified and fuel-rich prechamber operation promotes higher jet temperatures and as such are superior ignition sources. Relative to a lean prechamber jet, the peak temperature of a reacting jet from a rich prechamber was up to 600 K hotter at fixed distances from the orifice exit. Furthermore, radicals also demonstrated an influence on the ignition process, but the combustion mode was identified to be thermally dominant.
Primordial black holes (PBHs) in the asteroid-mass window are compelling dark matter candidates, made plausible by the existence of black holes and by the variety of mechanisms of their production in the early Universe. If a PBH falls into a white dwarf (WD), the strong tidal forces can generate enough heat to trigger a thermonuclear runaway explosion, depending on the WD’s mass and the PBH’s orbital parameters. In this work, we investigate the WD explosion triggered by the passage of a PBH. We perform 2D simulations of the WD undergoing thermonuclear explosion in this scenario, with the predicted ignition site as a parameter assuming the deflagration–detonation transition model. We study the explosion dynamics, and predict the associated light curves and nucleosynthesis. We find that the model sequence predicts light curves which align with the Phillips relation (B max versus ΔM 15 ). Our models hint at a unifying approach in triggering Type Ia supernovae without involving two distinctive evolutionary tracks.
Wildfire behavior is shaped by complex fire dynamics, with firebrands playing a critical role in spot fire ignition and fire spread. While previous studies have explored firebrand generation and transport, the specific role of towers and troughs from wildland fires in the lofting of firebrands remains unquantified. This study addresses that gap by using physics-based coupled fire-atmosphere model simulations to examine how wildfire towers (updrafts) and troughs (downdrafts) influence firebrand lofting. Our results show that the majority of firebrands (78.85%) are lofted from towers, where strong updrafts drive long-range transport. In contrast, only 21.15% of firebrands are lofted within troughs, where downdrafts cause most firebrands to fall near the fireline. We also find that firebrand size significantly influences lofting behavior, with smaller particles (1 mm radius) exhibiting the strongest correlation with updraft intensity. These findings highlight the dominant role of wildfire towers in promoting long-distance firebrand dispersal—an essential factor in rapid wildfire growth and wildland-urban interface (WUI) fire risks. By quantifying the relationship between firebrand lofting and fire-induced atmospheric features, this study provides critical insights to improve spot fire modeling, support mitigation planning, and enhance firefighter and WUI community safety in spot fire-prone regions.
Production of low carbon gasoline-like fuels such as methanol-to-gasoline (MTG) is a promising approach to achieve rapid greenhouse gas emission reduction of the transportation sector. Despite the fact that gasoline that meets the ASTM D4814 standard for automotive spark-ignition engine fuel can be readily produced from these processes, it is unclear how the composition of MTG may affect engine performance and emissions. Here, in this paper, a surrogate for an MTG is used to numerically study the effects of gasoline composition on knock propensity and on the sensitivity of knock to thermal and fuel stratification, to oxygen dilution and to nitric oxide from exhaust gas recirculation of residual gases. Simulations were performed in ANSYS CHEMKIN-PRO using a comprehensive chemical kinetic mechanism for gasoline surrogates, and results of the MTG surrogate were compared against those of a petroleum-based regular E10 gasoline, termed PACE-20. A premium-grade MTG fuel was also formulated by adding ethanol to the MTG surrogate, and results were compared against those of four premium-grade, gasoline-like fuels representative of future alternative gasoline formulations. Surrogates and mechanism were evaluated by comparison against experimental engine data, and the model showed high accuracy at stoichiometric conditions (mean absolute error of ignition timing equal to 1.46 crank angle degrees) but larger deviations at lean conditions (mean absolute error of ignition timing equal to 5.52 crank angle degrees). Despite the fact that the MTG surrogate has a RON 1.1 units higher than that of PACE-20, it may show higher knock propensity at medium temperature conditions due to a less intense NTC behavior. MTG autoignition was more temperature- and equivalence ratio-sensitive than that of PACE20, suggesting that MTG can benefit more from naturally-occurring thermal stratification or from induced fuel stratification of the end gas to mitigate knock intensity. The sensitivity of autoignition reactivity to oxygen dilution and to NO concentration was higher for MTG than for regular gasoline at medium loads, but the opposite trend was observed at high loads due to the effect of pressure on the low-temperature chemistry of regular gasoline. Approximately 14 % vol ethanol content was required to upgrade the octane rating of MTG from regular grade to premium grade. Adding 13.6 % vol ethanol made the fuel autoignition less sensitive to both oxygen dilution and NO content (ignition time varies approx. 17 % and 50 % less with oxygen dilution and NO addition, respectively, when adding ethanol at high engine loads).
The effect of transient plasma modes on ignition kernel development are discussed here for a quiescent CH 4 -air combustion model system. A 10 ns high-voltage pulse was applied to a pin-to-pin electrode in lean fuel-air mixtures at room temperature and atmospheric pressure. High-impedance streamer, transient spark and low-impedance spark discharges were identified based on pulse waveforms of voltage and current. A sustained ignition kernel expansion was observed when the plasma discharge transitioned into a transient spark or spark discharge. The minimum ignition energy was obtained at the transient spark mode, which has less than a third of the energy or Coulomb transfer compared to the low-impedance spark. Employing repetitive 10-pulse sequence at 10 kHz, the lean-fuel limit was extended from an equivalence ratio of 0.6 for the single pulse ignition to 0.5. The use of repetitive pulses also allowed streamer breakdown or spark initiation to occur at a lower voltage.