Search NASA⌕ Search

SEARCH · Search NASA

Results for “jet fuel”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DNS of ignition and flame stabilization in a simplified gas turbine premixer

With the increasing need for fuel flexibility, mitigation of auto-ignition (AI) inside gas turbine (GT) premixers becomes crucial. They must be designed to yield a sufficiently homogeneous fuel-air mixture to achieve low emissions while at the same time avoiding the occurrence of AI and subsequent flame stabilization. This challenge requires a detailed understanding of turbulent mixing and chemistry interactions. In the present work, a direct numerical simulation (DNS) of an array of jets in crossflow (JICF), representative of an industrial GT premixer, is reported to shed light on these complex phenomena. It is found that AI kernels form in the aft part of the premixer and coalesce into a flame front that then propagates upstream, mainly through the boundary layer, and successively engulfs the jets. This, therefore, suggests a significant role of the jet array pattern on the flame stabilization. It is noted that AI kernels continue to form independently during the whole time of the simulation. To clarify the contribution of AI and diffusion in the ignition kernels and the main flame, chemical explosive mode analysis (CEMA) is employed jointly with a kernel tracking algorithm. It is found that during the initial formation of the flame, many ignition kernels form in mixtures with low scalar dissipation rate and large contribution from AI mode. As they quickly grow, they merge into a single flame front that becomes increasingly more diffusion-assisted over time, balancing the AI mode. Turbulence is shown to have a significant enhancing effect in lean premixed flames, but further analysis is required to fully characterize it. These findings are relevant for the industrial premixer studied, and also for novel micromix concepts that may be used in the next generation of GT combustion systems.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Assessment of turbulent combustion models for simulating pre-chamber ignition in a natural gas engine

Lean combustion in an internal combustion engine is a promising strategy to increase thermal efficiency by leveraging a more favorable specific heat ratio of the fresh mixture and simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition events and slow flame propagation at fuel-lean condition lead to high cycle-to-cycle variability and hence limit the high-efficiency engine operating range. Pre-chamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling perspective, to date, still the science base and available simulation tools are inadequate for understanding and predicting the combustion processes in pre-chamber ignited engines. In this paper, conceptually different RANS combustion models widely adopted in the engine modeling community were used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a pre-chamber spark-ignition system. A flamelet-based turbulent combustion model, i.e., G-equation, and a multi-zone well-stirred reactor model were employed for the multi-dimensional study. Simulation results were compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two formulations respectively.

combustion model evaluation↗
DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 1

We report di-isobutylene has received significant attention as a promising fuel blendstock, as it can be synthesized via biological routes and is a short-listed molecule from the Co-Optima initiative. Di-isobutylene is also popularly used as an alkene representative in multi-component surrogate models for engine studies of gasoline fuels. However, there is limited experimental data available in the literature for neat di-isobutylene under engine-like conditions. Hence, most existing di-isobutylene models have not been extensively validated, particularly at lower temperatures (< 1000 K). Most gasoline surrogate models include the di-isobutylene sub-mechanism published by Metcalfe et al. with little or no modification. The current study is undertaken to develop a detailed kinetic model for di-isobutylene and validate the model using a wide range of relevant experimental data. Part 1 of this study exclusively focuses on the low- to intermediate temperature kinetics of di-isobutylene. An upcoming part 2 discusses the high-temperature model development and validation of the relevant experimental targets. Ignition delay time measurements for the di-isobutylene isomers were performed at pressures ranging from 15 - 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in air and in the temperature range 650 - 900 K using two independent rapid compression machine facilities. In addition, measurements of species identified during the oxidation of these isomers were performed in a jet-stirred reactor and in a rapid compression machine. A detailed kinetic model for the di-isobutylene isomers is developed to capture the wide range of new experimental targets. For the first time, a comprehensive low-temperature chemistry submodel is included. The differences in the important reaction pathways for the accurate prediction of the oxidation of the two DIB isomers are compared using reaction path analysis. In conclusion, the most sensitive reactions controlling the ignition delay times of the DIB isomers under the pressure and temperature conditions necessary for autoignition in engines are identified.

09 BIOMASS FUELS↗
DOE OSTI2020

Elucidating the differences in oxidation of high-performance α- and β- diisobutylene biofuels via Synchrotron photoionization mass spectrometry

Abstract Biofuels are a promising ecologically viable and renewable alternative to petroleum fuels, with the potential to reduce net greenhouse gas emissions. However, biomass sourced fuels are often produced as blends of hydrocarbons and their oxygenates. Such blending complicates the implementation of these fuels in combustion applications. Variations in a biofuel’s composition will dictate combustion properties such as auto ignition temperature, reaction delay time, and reaction pathways. A handful of novel drop-in replacement biofuels for conventional transportation fuels have recently been down selected from a list of over 10,000 potential candidates as part of the U.S. Department of Energy’s (DOE) Co-Optimization of Fuels and Engines (Co-Optima) initiative. Diisobutylene (DIB) is one such high-performing hydrocarbon which can readily be produced from the dehydration and dimerization of isobutanol, produced from the fermentation of biomass-derived sugars. The two most common isomers realized, from this process, are 2,4,4-trimethyl-1-pentene (α-DIB) and 2,4,4-trimethyl-2-pentene (β-DIB). Due to a difference in olefinic bond location, the α- and β- isomer exhibit dramatically different ignition temperatures at constant pressure and equivalence ratio. This may be attributed to different fragmentation pathways enabled by allylic versus vinylic carbons. For optimal implementation of these biofuel candidates, explicit identification of the intermediates formed during the combustion of each of the isomers is needed. To investigate the combustion pathways of these molecules, tunable vacuum ultraviolet (VUV) light (in the range 8.1–11.0 eV) available at the Lawrence Berkeley National Laboratory’s Advanced Light Source (ALS) has been used in conjunction with a jet stirred reactor (JSR) and time-of-flight mass spectrometry to probe intermediates formed. Relative intensity curves for intermediate mass fragments produced during this process were obtained. Several important unique intermediates were identified at the lowest observable combustion temperature with static pressure of 93,325 Pa and for 1.5 s residence time. As this relatively short residence time is just after ignition, this study is targeted at the fuels’ ignition events. Ignition characteristics for both isomers were found to be strongly dependent on the kinetics of C 4 and C 7 fragment production and decomposition, with the tert-butyl radical as a key intermediate species. However, the ignition of α-DIB exhibited larger concentrations of C 4 compounds over C 7 , while the reverse was true for β-DIB. These identified species will allow for enhanced engineering modeling of fuel blending and engine design.

09 BIOMASS FUELS↗
DOE OSTI2024

Next Gen High Efficiency Boosted Engine Development

This work represents an advanced engineering research project partially funded by the U.S. Department of Energy (DOE). Ford Motor Company, FEV North America, and Oak Ridge National Laboratory collaborated to develop a next generation boosted spark ignited engine concept. The project goals, specified by the DOE, were 23% improved fuel economy and 15% reduced weight relative to a 2015 or newer light-duty vehicle. The fuel economy goal was achieved by designing an engine incorporating high geometric compression ratio, high dilution tolerance, low pumping work, and low friction. The increased tendency for knock with high compression ratio was addressed using early intake valve closing (EIVC), cooled exhaust gas recirculation (EGR), an active pre-chamber ignition system, and careful management of the fresh charge temperature. Engine weight reduction measures were implemented throughout the engine system making use of composite materials, advanced manufacturing techniques, and architectural choices. This report outlines the analytical, design, fabrication, and test work conducted for the duration of the project. The combustion system stability, EGR tolerance, and knock resistance were validated on a single cylinder engine. An inline six-cylinder engine was then designed targeting application in the Ford F150. Multi-cylinder engines were produced and tested achieving the target vehicle fuel economy improvement of 23% assessed using measured engine fuel consumption combined with a vehicle drive cycle simulation. Actions were identified and designs were demonstrated to achieve the 15% weight reduction target. This project included items covering a range of technology readiness levels. Some of the technologies explored are production ready, while others were investigated to understand the limitations for what can be achieved in a stoichiometric, gasoline-fueled, spark-ignited internal combustion engine.

42 ENGINEERING↗
DOE OSTI2023

Representative Phenomena of Cyclic Turbulent Combustion in High-Pressure Fuel Sprays

Abstract Cyclic variations in conventional diesel combustion engines can lead to large differences in engine out emissions even at steady operation. This study uses an optically accessible constant-pressure flow chamber to acquire fuel injections in quick succession to analyze mixing, auto-ignition, and combustion of diesel-surrogate n-heptane using multiple high-speed optical diagnostics. Prior studies have utilized fewer injections and/or they rely on analysis of ensemble average behavior. These approaches do not yield information on injection-to-injection variation or provide confidence in utilizing individual injection measurements for high-fidelity computational fluid dynamics(CFD) model validation. In this study, a large set of 500 injections is used to obtain global parameters including liquid length, vapor penetration length, ignition delay time, and lift-off length. Results for multiple injections are presented to illustrate large injection to injection variations. Potential sources for these variations are analyzed to conclude localized, small scale turbulence and rate of injection variations as the likely sources. Then, a statistical method based on z-scores is proposed and implemented to identify instantaneous injections that best represent the bulk data-set of jet boundaries measured independently by three different diagnostics. This synthesis of statistics-guided screening of data set and ensemble-average analysis offers higher confidence for CFD model validation relying upon both a representative single and average injection results.

42 ENGINEERING↗
DOE OSTI2025

Quantification of regional net CO 2 flux errors in the Orbiting Carbon Observatory-2 (OCO-2) v10 model intercomparison project (MIP) ensemble using airborne measurements

Inverse model intercomparison projects (MIPs) provide a chance to assess the uncertainties in inversion estimates arising from various sources. However, accurately quantifying ensemble CO 2 flux errors remains challenging and often relies on the ensemble spread. This study proposes a method for quantifying the errors in regional net surface–atmosphere CO 2 flux estimates from models taken from the Orbiting Carbon Observatory-2 (OCO-2) v10 MIP by using independent airborne CO 2 measurements for the period 2015–2017. We first calculate the root mean square error (RMSE) between the ensemble mean of posterior CO 2 concentrations and airborne observations and then isolate the CO 2 concentration errors caused solely by the ensemble mean of posterior net fluxes by subtracting the observation, representation, and transport errors from seven regions. Our analysis reveals that the flux errors projected onto CO 2 space account for 55 %–85 % of the regional average RMSE over the 3 years, ranging from 0.88 to 1.91 ppm. In five regions, the error estimates based on observations exceed those computed from the ensemble spread of posterior fluxes by a factor of 1.33–1.93, implying an underestimation of the actual flux errors, while their magnitudes are comparable in two regions. The adjoint sensitivity analysis identifies that the underestimation of flux errors is prominent where the magnitudes of fossil fuel emissions exceed those of terrestrial-biosphere fluxes by a factor of 3–31 over the 3 years. This suggests the presence of systematic biases in the inversion estimates associated with errors in the prescribed fossil fuel emissions common to all models. Our study emphasizes the value of airborne measurements for quantifying regional errors in ensemble net CO 2 flux estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Quantitative Risk Assessment for Fuel Cell Electric Bus Hydrogen Storage and Refueling Facility

It is necessary to understand the safety implications and risk mitigation options for fuel cell electric bus fleet deployment, especially for related facilities responsible for operations such as production, storage, compression, and dispensing of hydrogen for use by the buses. In this report, we present a quantitative risk assessment for a potential fuel cell electric bus fleet that was motivated by efforts to improve resilience at the Portland International Airport but can be applicable to a range of hydrogen case studies and use cases. We estimated risk for a facility that produces, stores, compresses, and dispenses hydrogen for the fleet of buses, with a focus on individual risk to people in terms of annual frequency of fatality. We considered the frequency of hydrogen leaks that could result in harmful physical outcomes like jet fires or explosions, and the consequences of those outcomes for people. We created customized fault trees to calculate the frequencies of different sizes of leaks and event sequence diagrams to calculate ignition probabilities for the various leak sizes. We also leveraged the HyRAM+ toolkit to use these inputs to calculate overall risk for the facility, which we separated into one section responsible for producing, storing, and compressing hydrogen, and one section responsible for dispensing the hydrogen to the buses. We found that the dispensing area seemed to have a higher risk than the production/storage/compression area of the facility, largely because of the inclusion of a component with a high leak frequency (the heat exchanger used to cool the hydrogen before entering the vehicle, to prevent overheating and expansion of hydrogen in the onboard tank). For the example production and refueling facility we evaluated and the data we used for the analysis, the leak frequency had a larger impact on the risk differences between the two sections on the facility, compared to the physical outcome consequence, which was slightly different due to the varying fuel conditions, but not substantially different. Actions can be taken to prevent these hazards (e.g., lowering leak frequencies in system components) or to mitigate the consequences if they do occur (e.g., installing barriers to protect people if ignition events occur). The choice of which actions to take depends not only on safety considerations but also on space, time, staffing, feasibility, and financial constraints. Therefore, the quantitative risk assessment approach can help understand relative risk contributions from different components, leak sizes, consequences, and human actions, to prioritize risk reduction strategies and balance these parameters. The outcomes of this report may be useful for a variety of stakeholders working in the hydrogen, transportation, vehicle, and aviation sector, including those responsible for aspects like facility design, operations, and regulations. There is not a single value of risk that determines whether a hypothetical system is “safe” or not. The insights about risk mitigations may be leveraged, and the quantitative risk assessment approach can be applied to other case studies to understand risk priorities and contributions specific to different FCEB and hydrogen facility uses.

08 HYDROGEN↗
DOE OSTI2023

Parameter dependencies of the separatrix density in low triangularity L-mode and H-mode JET-ILW plasmas

Abstract The midplane electron separatrix density, n e,sep , in JET-ILW L-mode and H-mode low triangularity deuterium fuelled plasmas exhibits a strong explicit dependence on the averaged outer divertor target electron temperature, n e,sep ∼ T e,ot −1/2 . This dependence is reproduced by analytic reversed two point model (rev-2PM), and arises from parallel pressure balance, as well as the ratio of the power and momentum volumetric loss factors, (1 − f cooling )/(1 − f mom-loss ). Quantifying the influence of the (1 − f cooling ) and (1 − f mom-loss ) loss factors on n e,sep has been enabled by measurement estimates of these quantities from L-mode density (fueling) ramps in the outer horizontal, VH(C), and vertical target, VV, divertor configurations. Rev-2PM n e,sep estimates from the extended H-mode and more limited L-mode datasets are recovered to within ±25% of the measurements, with a scaling factor applied to account for use of T e,ot , an averaged quantity, rather than flux tube resolved target values. Both the (1 − f cooling ) and (1 − f mom-loss ) trends and recovery of n e,sep using the rev-2PM formatting are reproduced in EDGE2D-EIRENE L-mode-like and H-mode-like density scan simulations. The general lack of a divertor configuration effect in the JET-ILW n e,sep trends can be attributed to a significant influence of main chamber recycling, which has been shown in the EDGE2D-EIRENE results to moderate n e,sep with respect to changes in divertor neutral leakage imposed by changes in the divertor configuration. The unified n e,sep vs T e,ot trends can, however, be broken if large modifications to the divertor geometry (e.g. complete removal of the outer divertor baffle structure) are introduced in the model. The more pronounced high-field side high density region formation in the VH(C) configuration with reduced clearance to the separatrix does not appear to have a significant influence on the outer midplane separatrix and pedestal parameters when mapped to T e,ot , although conditions at the inner midplane could not be assessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

An arbitrary Lagrange-Eulerian investigation of HRAM shallow jet pre-spurt formation and time sensitivities to impact plate dynamics

With dry bay fires persisting as a significant contributor to aircraft vulnerability despite chronicled developments in survivability technologies, an accurate fire prediction capability remains paramount for credible vulnerability assessments. Physics based modeling of the hydrodynamic ram (HRAM) fluid deposition process is a key component of such capability, wherein capturing the first instance of fluid spurt, referred to herein as shallow jet spurts, is a core focus. Such pre-spurts as they have been formerly identified have only been witnessed sporadically in HRAM spurt experiments. In this work, ALE3D, a first-principles multi-physics code was employed to model the shallow jet spurt phenomenon with spherical projectiles impacting water-filled tanks faced with aluminum panels, such that the underlying physics and sensitivities could be explored. Development and verification of the 2D axisymmetric model is described relative to trends observed in a prior experimental campaign. Results from the verified model suggest shallow jet spurts have at least a quadratic sensitivity to the fundamental vibrational mode of the impact plate across impact velocities 610–1829 m/s (2000–6000 ft/s). It is further explained how shallow jet spurts are arrested for impacts at the two extremes of plate rigidity. This research constitutes the first fluid-structure modeling of shallow jet spurts which future three-dimensional analyses will expound upon.

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

Mixing-Controlled Combustion of Ethanol Enabled by Prechamber Ignition (PC-MCC): A Preliminary Experimental Demonstration

This experimental study presents preliminary investigations of prechamber-enabled mixing-controlled combustion (PC-MCC) at −2 bar brake mean effective pressure (BMEP) and 2200 rpm with fuel-grade ethanol (E98). Experimental results are conducted on a prechamber retrofitted single-cylinder Caterpillar C9.3B test engine. First, a series of prechamber-only experiments were conducted with a motored engine to evaluate the salient combustion trends in response to relevant prechamber operating parameters. Under firing conditions, the prechamber operating strategy was assessed with respect to the impact on ignition assistance of direct-injected E98 and overall engine performance. The preliminary results indicate the jet-induced ignition process is robust and prompts diffusion combustion of E98 at diesel-like boundary conditions. Here, the effect of external exhaust gas recirculation (EGR) on the residual tolerance of the prechamber combustion process was also investigated and showed stable combustion in both the main chamber and prechamber up to 30% EGR. Experiments were also conducted with the stock diesel engine for baseline comparison. At matched combustion phasing, mixing-controlled combustion of ethanol enabled by prechamber ignition was able to achieve heightened gross thermal efficiency while simultaneously reducing NOx and practically eliminating smoke emissions relative to diesel combustion. In addition, the covariance of load and standard deviation of combustion phasing was diesel-like and less than 2% and 1 CAD, respectively.

PC-MCC↗
DOE OSTI2021

Increasing large wildfires over the western United States linked to diminishing sea ice in the Arctic

Abstract The compound nature of large wildfires in combination with complex physical and biophysical processes affecting variations in hydroclimate and fuel conditions makes it difficult to directly connect wildfire changes over fire-prone regions like the western United States (U.S.) with anthropogenic climate change. Here we show that increasing large wildfires during autumn over the western U.S. are fueled by more fire-favorable weather associated with declines in Arctic sea ice during preceding months on both interannual and interdecadal time scales. Our analysis (based on observations, climate model sensitivity experiments, and a multi-model ensemble of climate simulations) demonstrates and explains the Arctic-driven teleconnection through regional circulation changes with the poleward-shifted polar jet stream and enhanced fire-favorable surface weather conditions. The fire weather changes driven by declining Arctic sea ice during the past four decades are of similar magnitude to other leading modes of climate variability such as the El Niño-Southern Oscillation that also influence fire weather in the western U.S.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗
DOE OSTI2020

High pressure/high temperature multiphase simulations of dodecane injection to nitrogen: Application on ECN Spray-A

The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.

42 ENGINEERING↗
DOE OSTI2020

Influence of ozone addition on the low-temperature oxidation of dimethyl ether in a jet-stirred reactor

The influence of ozone addition on the low-temperature oxidation of dimethyl ether (DME) was investigated experimentally in an atmospheric-pressure jet-stirred reactor, over the temperature range of 400–800 K. Detailed speciation information was obtained by employing synchrotron vacuum ultraviolet photoionization mass spectrometry. Experimental results revealed that the ozone addition had a positive influence on the production of the highly reactive intermediates. Moreover, the low-temperature reactivity of DME was significantly enhanced, which resulted in the broadening of the temperature window of fuel consumption and intermediates formation at lower temperatures. Therefore, novel experimental data of the low temperature regime (400–500 K) could be obtained. The data set of this special temperature regime yielded insights into the DME low-temperature kinetics, which were further supported with modeling analysis based on two existing DME models (Metcalfe et al., 2013; Wang et al., 2015) combined with an ozone sub-mechanism (Zhao et al., 2016). The analysis showed that temperature-sensitive reactions such as the second oxygen channel could be nearly “frozen” at this low temperature (T < 440 K). Furthermore, the production of some intermediates was found to be strongly governed by reaction pairs, such as CH 3 OCH 2 + O 2 = CH 3 OCH 2 O 2 and CH 3 OCH 2 + O 2 = 2 CH 2 O + OH for the CH 2 O formation. This finding could be useful for examining branching ratios in both models, and the analysis suggested the further modification of the branching ratios for the oxygen addition to CH 3 OCH 2 O 2 pathways and the CH 3 OCH 2 O 2 self-reactions were required. Finally, the influences of the O 3 addition in the sensitive reactions of the fuel initial low-temperature oxidation were investigated in this work. It was interesting to note that O 3 addition could change the dominating reactions in the initial low-temperature oxidation, by the addition of some O 3 -related pathways with relatively high sensitivity.

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

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↗
DOE OSTI2021

Effects of 18 O Exchange on Neutron Emission Rates of Aging 238 PuO 2 LWRHUs

The Lightweight Radioisotope Heater Unit (LWRHU) provides about 1 Watt of thermal power to critical electrical and mechanical systems in the extreme temperatures of space, maintaining the operating temperature of the spacecraft’s scientific and general operating equipment. Los Alamos National Laboratory (LANL) designed and tested the LWRHU in the 1970s through early 1980s. The first production campaign at LANL was in the early 1980s in support of the Galileo mission in which over 100 LWRHUs were manufactured. In the mid-1990s, almost 200 LWRHUs were manufactured in support of the Cassini mission. After the primary Cassini manufacturing campaign, the Department of Energy (DOE) approved the fabrication of more than 10 additional LWRHUs as spares in 1998. The LWRHU assembly contains the following components: a fuel pellet, a vented capsule, a pyrolytic graphite insulator, and a fine-weave pierced fabric (FWPF) graphite aeroshell. The fuel pellet is a hot-pressed 238 PuO 2 cylindrical pellet that is sintered to create a ceramic pellet. The pellet is encapsulated in a platinum-rhodium encapsulation. The welded fueled clad is placed inside a pyrolytic graphite thermal insulation to protect the fuel from the heat produced in reentry events. The exterior layer of the LWRHU assembly is the FWPF aeroshell. Before sintering, PuO 2 granules are heated in a furnace at temperatures greater than or equal to 700°C under an enriched 16 O atmosphere in order to reduce trace quantities of 18 O. Because 18 O undergoes an (α,n) reaction, 18 O needs to be removed as much as possible to meet radiation emission requirements set by the Jet Propulsion Laboratory (JPL) and the National Aeronautics and Space Administration (NASA), as well as other agencies. The relative abundances of the oxygen isotopes are indicated by neutron emission rate (NER) measurements. LWRHUs are vented with a sintered platinum frit vent (George 1986). During extraterrestrial operation, the vent in the cladding allows helium from the decay of plutonium to escape and relieves pressure from the system while preventing the solids from escaping (Rinehart 1996, Tate 1982 & 1985). This vent is protected during manufacture by a platinum-30 rhodium cover. The capsule vent is activated by milling a 0.025-inch diameter hole through the protective cover to a nominal depth of 0.015 inches, generally immediately before the LWRHU is loaded into the aeroshell (George 1986). Prior to long-term storage, LWRHUs are vented and thus potentially allow infiltration of 18 O from the ambient atmosphere. Although NER decreases with age due to the reduction in emitted alpha particles, it might be expected that NER reduction is mitigated by the exchange of 18 O. However, this is not necessarily the case and no trend was observed relating NER to age other than the natural decay of 238 Pu (Mulford 2021). The most probable and dominant factor for determining specific NER was manufacturing discrepancies between individual units. It is important to understand the lack of observed 16 O- 18 O exchange and to compare LWRHU pellets as much as possible.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗