Search NASA⌕ Search

SEARCH · Search NASA

Results for “shock sensitivity”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

DME-Propane Blends Ignition Experiments and Modeling for Heavy-Duty Mixing Controlled Compression Engines

A blend of dimethyl ether (DME) and propane (C3H8) is being studied in a shock tube at heavy-duty engine conditions at 110 bar. Due to its intrinsic combustion properties, DME/propane blend can potentially replace diesel in mixing controlled compression ignition engines. A blend of DME/propane can reduce emissions in mixing controlled compression ignition in heavy-duty engines through modifications, which require simulations using a high-fidelity chemical kinetics model that can accurately predict the chemistry of the blend. An essential aspect of testing the chemical kinetics model is doing baseline fundamental chemistry studies on neat DME and propane, which include ignition delay time measurements. In this work, using a high-pressure shock tube, ignition delay times were gathered for DME/Propane blends at 110 bar diluted with AR to test chemical kinetic models published in the literature. These models include Aramco 3.0, NUIG V1.1, C3mech V3.3, and Dames et al. Comparisons with the experimental IDTs and models were conducted, and general agreement was observed. A sensitivity analysis was conducted, and important reactions were outlined.

Mohammed, Zuhayr Pasha [University of Central Flor↗
DOE OSTI2020

A measurement of the equation of state of carbon envelopes of white dwarfs

White dwarfs represent the final state of evolution for most stars. Certain classes of white dwarfs pulsate, leading to observable brightness variations, and analysis of these variations with theoretical stellar models probes their internal structure. Modelling of these pulsating stars provides stringent tests of white dwarf models and a detailed picture of the outcome of the late stages of stellar evolution. However, the high-energy-density states that exist in white dwarfs are extremely difficult to reach and to measure in the laboratory, so theoretical predictions are largely untested at these conditions. In this paper we report measurements of the relationship between pressure and density along the principal shock Hugoniot (equations describing the state of the sample material before and after the passage of the shock derived from conservation laws) of hydrocarbon to within five per cent. The observed maximum compressibility is consistent with theoretical models that include detailed electronic structure. This is relevant for the equation of state of matter at pressures ranging from 100 million to 450 million atmospheres, where the understanding of white dwarf physics is sensitive to the equation of state and where models differ considerably. The measurements test these equation-of-state relations that are used in the modelling of white dwarfs and inertial confinement fusion experiments, and we predict an increase in compressibility due to ionization of the inner-core orbitals of carbon. We also find that a detailed treatment of the electronic structure and the electron degeneracy pressure is required to capture the measured shape of the pressure–density evolution for hydrocarbon before peak compression. Our results illuminate the equation of state of the white dwarf envelope (the region surrounding the stellar core that contains partially ionized and partially degenerate non-ideal plasmas), which is a weak link in the constitutive physics informing the structure and evolution of white dwarf stars.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Influence of the double bond position in combustion chemistry of methyl butene isomers: A shock tube and laser absorption study

Shock tube experiments have been carried out on 2–methyl–1–butene (2M1B), 2–methyl–2–butene (2M2B), and 3–methyl–1–butene (3M1B)—the three isomers of methyl butene compound. Here, carbon monoxide (CO) time–histories and ignition delay times are obtained behind reflected shockwaves over the temperature range of 1350–1630 K and pressures of 8.3–10.5 atm with stoichiometric mixtures of 0.075% fuel in O 2 /Ar. Comparative ignition study reveals that 3M1B ignites significantly faster than the other two isomers, while 2M1B dissociates earlier but ignites later than 2M2B. Possible mechanisms for this behavior are discussed with ignition delay time sensitivity and reaction path analysis. In addition, time–resolved CO measurements are compared with three different reaction mechanisms from the literature. Sensitivity analyses have been carried out to identify important reactions that need attention to accurately predict the chemistry of these isomers. Further investigation into the rates of unimolecular fuel decomposition reactions and C 3 H 3 + O 2 = CH 2 CO + HCO reaction are suggested based on the current investigation.

09 BIOMASS FUELS↗
DOE OSTI2020

Strength of porous α-SiO 2 in a shock loaded environment: Calibration via Richtmyer–Meshkov instability and validation via Mach lens

The strength of brittle porous media is of concern in numerous applications, for example, earth penetration, crater formation, and blast loading. Thus, it is of importance to possess techniques that allow for constitutive model calibration within the laboratory setting. The goal of the current work is to demonstrate an experimental technique allowing for strength assessment of porous media subjected to shock loading, which can be implemented into pressure-dependent yield surfaces within numerical simulation schemes. As a case study, the deviatoric response of distended α-SiO 2 has been captured in a tamped Richtmyer–Meshkov instability (RMI) environment at a pressure regime of 4–10 GPa. Hydrocode simulations were used to interpret RMI experimental data, and a resulting pressure-dependent yield surface akin to the often employed modified Drucker–Prager model was calibrated. Simulations indicate that the resulting jet length generated by the RMI is sensitive to the porous media strength, thereby providing a feasible experimental platform capable of capturing the pressurized granular deviatoric response. Furthermore, in efforts to validate the RMI-calibrated strength model, a set of Mach-lens experiments was performed and simulated with the calibrated pressure-dependent yield surface. Excellent agreement between the resulting Mach-lens length in experiment and simulation provides additional confidence to the RMI yield-surface calibration scheme.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Search for sub-TeV Neutrino Emission from Novae with IceCube-DeepCore

The understanding of novae, the thermonuclear eruptions on the surfaces of white dwarf stars in binaries, has recently undergone a major paradigm shift. Though the bolometric luminosity of novae was long thought to arise directly from photons supplied by the thermonuclear runaway, recent gigaelectronvolt (GeV) gamma-ray observations have supported the notion that a significant portion of the luminosity could come from radiative shocks. More recently, observations of novae have lent evidence that these shocks are acceleration sites for hadrons for at least some types of novae. In this scenario, a flux of neutrinos may accompany the observed gamma rays. As the gamma rays from most novae have only been observed up to a few GeV, novae have previously not been considered as targets for neutrino telescopes, which are most sensitive at and above teraelectronvolt (TeV) energies. Here, we present the first search for neutrinos from novae with energies between a few GeV and 10 TeV using IceCube-DeepCore, a densely instrumented region of the IceCube Neutrino Observatory with a reduced energy threshold. We search both for a correlation between gamma-ray and neutrino emission as well as between optical and neutrino emission from novae. We find no evidence for neutrino emission from the novae considered in this analysis and set upper limits for all gamma-ray detected novae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Infrared absorption spectroscopy of dynamically compressed water

Streaked visible spectroscopy is well established in dynamic compression research. Infrared measurements remain problematic, however, due to the diminishing sensitivity of streak camera photocathodes beyond 800 nm. Time-stretch techniques offer an alternative method for probing infrared features during single-event experiments. This paper discusses the development of a time-stretch spectroscopy diagnostic using dispersed supercontinuum laser pulses. The technique is applied to near-infrared measurements of liquid water during multiple shock compression.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

An Eulerian multimaterial framework for simulating high-explosive aquarium tests

Aquarium tests of cylindrical high-explosive charges provide optical data of the detonation front velocity and shape, propagation of the shock wave in the surrounding water, and expansion rates of the detonation products behind the front. Data from aquarium experiments is often used for calibration of reactive burn models based on phenomenological equations of state (EOS) and reaction rate laws. This paper presents a multimaterial numerical modeling framework to solve the 2D axisymmetric reactive Euler equations for high-explosive aquarium tests, in particular for ammonium nitrate - fuel oil (ANFO) explosives. An extension of the Ghost Fluid Method (GFM) is used to handle the dynamic material interfaces for the ANFO explosion products, the charge-confining material (polymethyl methacrylate PMMA), and the surrounding water. This study analyzes the sensitivity of calculations (both computational efficiency and numerical accuracy) to different algorithms for the material interface models including the original GFM versus Riemann solver-based strategies. A novel method for defining the left and right states in the interfacial Riemann problem eliminates the need for sorting or nodal interpolation during the projection along the material interface. Numerical tests indicate that populating the interface node values using the Riemann solution mitigate the overheating error observed in steady-state calculations. Solution convergence and computational efficiency are explored as a function of the spatial and temporal order of the schemes. Results from the computational model with analytical equations of state and fitted reaction rate parameters show very good quantitative agreement with experimentally observed detonation front velocity, reaction products expansion, and shock wave propagation in the surrounding water for a cylindrical ANFO charge encased in PMMA. Finally, the proposed modeling framework, in conjunction with experimental tests, provides a reliable tool to assess equations of state and reaction rate expressions for reactive burn models of confined high explosives.

42 ENGINEERING↗
DOE OSTI2022

Radiation hydrodynamics modelling of kilonovae with SNEC

ABSTRACT We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A&L Annual Report: GRM-MAPS

Exploding Bridge Wire (EBW) is an important design for initiating detonators relevant to several LLNL systems. In such designs, a metal bridge wire is in direct contact with a porous secondary high explosive (HE), typically pressed to around 50% of its crystalline density. Initiation is believed to occur through shock waves generated when the bridge wire is made to explode via the rapid passing of a strong current pulse. Permeametry is an important characterization tool for the porous powdered HE (PETN, HMX, etc.) component in such devices. Such powders are well known to coarsen with age, with an accompanying deterioration in performance, which includes increased function time and sensitivity. Among various material characteristics, one that has been shown to correlate the most with such performance parameters is the flow-permeable surface area (FSSA), measured using air-permeametry. Commercially available permeametry apparatuses, such as that from Fisher Scientific or a more modern version from Micromeritics suffer from a few drawbacks, most notably as related to monitoring age-related changes in FSSA. This is because it requires re-compaction of the aged powder into the permeametry tube, which in most cases leads to an unwanted increase in FSSA. To address this and a few other shortcomings, we have been developing a modern version of the permeameter, Gee-Reinstein-Maiti Modern Air-Permeametry System (GRM-MAPS). In this report, we summarize progress made in FY25 toward perfecting and calibrating the instrument, which we aim at deploying within surveillance in FY26.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulation and analysis of the 3D double shock-induced recompression of spall using a hybrid PAGOSA with FLIP+MPM

In this work, an innovative hybrid algorithm PAGOSA with FLIP + MPM is first presented by coupling PAGOSA with the particle FLIP + MPM, and applied to systematically exploring the 3D recompression of spall induced by double-shock waves in ductile materials as the few existing studies are currently limited to the pure longitudinal recompression of spall, and ignore the 3D effects. This algorithm solves the difficulties the grid-based methods encounter when applied to capturing the fracture in material. We first validate the capabilities of PAGOSA with FLIP + MPM to predict different fracture/fragmentation cases by solving two benchmark problems and comparing with the analytical solutions or the experiment results. The convergences and the infinity-norm errors are also investigated. Subsequently, the 3D effects are demonstrated by simulating the spallation in material using PAGOSA with FLIP + MPM and comparing with those in the longitudinal cases. Additionally, some factors that affect the spallation with 3D effects are also discussed. Finally, the recompressions of spall in material with/without 3D effects are systematically simulated and analyzed. The recompression of spall driven by a high explosive is also simulated to show the ability of PAGOSA with FLIP + MPM to handle this kind of problem. The conditions that can lead to a recompression of spall are explored. Also, numerical results show that PAGOSA with FLIP + MPM can accurately predict different fracture cases, that the 3D effects play an important role in fracture/fragmentation, and the recompression of spall is very sensitive to the occurrence conditions, shows great differences when considering 3D effects in real applications. Moreover, the present hybrid method can be easily extended to other grid-based techniques employed for fracture in materials.

3D effects↗
DOE OSTI2022

Experimental Ignition Delay Time Measurements and Chemical Kinetics Modeling of Hydrogen/Ammonia/Natural Gas Fuels

In recent years, hydrogen-carrying compounds have accrued interest as an alternative to traditional fossil fuels due to their function as zero-emission fuels. As such, there is interest in investigating hydrogen-carrying compounds to improve understanding of the fuels' characteristics for use in high pressure systems. In the current study, the oxidation of ammonia/natural gas/hydrogen mixtures was carried out to study carbon monoxide (CO) formation profiles as well as the ignition delay times (IDTs) behind reflected shock waves in order to refine chemical kinetic models. Experiments were carried out in the University of Central Florida's shock tube facility by utilizing chemiluminescence to obtain OH* emission and laser absorption spectroscopy to obtain CO profiles over a temperature range between 1200 K and 1800 K with an average pressure of 2.2 atm. Experimental mixtures included both neat and combination natural gas/hydrogen with ammonia addition, with all mixtures except one having an equivalence ratio of 1. Results were then compared with the GRI 3.0 mechanism, as well as the newly developed UCF 2022 mechanism utilizing chemkin-pro software. In general, both models were able to capture the trend in auto-ignition delay times and CO time histories for natural gas and ammonia mixtures. However, for ammonia–hydrogen mixtures, GRI 3.0 failed to predict ignition delay times, whereas the UCF 2022 mechanism was able to capture the IDTs within the uncertainty limits of the experiments. A sensitivity analysis was conducted for different mixtures to understand the important reactions at the experimental conditions. Lastly, a reaction pathway analysis was carried out to understand important ammonia decomposition pathways in the presence of hydrogen and natural gas.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Comments on the SURFplus reactive burn model

The SURFplus reactive burn model is intended for heterogeneous solid HE that produce an excess amount of carbon. It utilizes a fast hotspot reaction (SURF model) plus a second slow reaction for the energy release from carbon clustering. The fast reaction dominates shock initiation, though the fast rate parameters are slightly affected by the late energy release from the slow reaction. The slow reaction has a significant affect on propagating detonation waves. Propagating detonation waves are characterized by the curvature effect; i.e., detonation speed $D_n$ as function of the local front curvature κ. The $D_n(κ)$ curves have a qualitative different shape. This results from the sensitivity of the slope of the $D_n(κ)$ curve to the reaction-zone width up to the sonic point, and the sonic point moving as κ increases from the end of the slow reaction to the neighborhood of the end of the fast reaction. For PBX 9502, the steep slope for small curvature is determined by the slow reaction, while the lower slope for larger curvature is determined by the fast reaction rate in the high detonation pressure regime. The remainder of this paper is organized as follows. In section 2 we review the assumptions of the Shaw and Johnson carbon cluster model used by the SURFplus model for the slow reaction. Section 3 describes in detail the mixture model for a partly burned HE. The form of the slow rate used by the SURFplus model is described in section 4. We conclude with a brief discussion of applying a burn model with 2 rates to either a metalized or a composite PBX.

42 ENGINEERING↗
DOE OSTI2022

Path Integral Monte Carlo Simulations of Iron Plasmas (Final Technical Report)

This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

HED Hydrodynamics in the Common Modeling Framework

Within the Common Modeling Framework (CMF), the HED Hydrodynamics Project has implemented a common methodology, or "Authority", for modeling high energy density (HED) experiments by extending development of the Eulerian Applications Project code Authority, "EAPA," to include relevant physics models, customization tools and templates, and in situ and post-processing capabilities. Additionally, several models for small-scale HED experiments were successfully added to the CMF and simulations reproduce the results of preexisting, benchmarked input decks. We expect use of the CMF to improve our effectiveness in developing predictive capabilities for experiments ranging from small-scale planar single-interface, single-shock to multi-interface, multi-shock configurations and up through multi-shell inertial confinement fusion (ICF) implosions. Additionally, this will strengthen our ability to design future targeted experiments. Using the CMF infrastructure and Authorities for HED modeling has three key impacts. First, it enables a version-controlled, pedigreed and archived base model for each experiment with tools that allow easy setup, execution and analysis of simulations to assess and/or expose sensitivities to various parameters. Second, having a common base set of inputs and physics definitions helps identify systematic differences between models in the codes and experimental data to support code validation, and build confidence in our predictive capabilities. Finally, designed to work in concert with multiple programs and projects including PEM, IC, OES, V&V and DSW, the CMF supports improved collaboration and integration through shared model definitions, data, simulations and results, enhancing the workflow required to underwrite the knowledge base and capabilities ultimately required for stockpile stewardship

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A coupled hydrologic-agroeconomic modeling framework to evaluate adaptive irrigation strategies under groundwater withdrawal restrictions

Growing groundwater scarcity requires integrated tools to capture interactions among hydrology, agricultural production, markets, and land use. This study presents an iterative modeling framework that couples hydrologic, crop-yield, and economic models to capture two-way feedback among water availability, agricultural production, and market responses under groundwater constraints. The primary goal of this paper is to describe the methodological development of the coupled framework and demonstrate the significance of iterative model interaction. Applied to the western United States, we evaluated adaptive responses to restricting groundwater use beyond recharge levels, represented through changes in irrigation management and expansion or shrinkage of crop markets through land reallocation. Results demonstrate that the iterative coupling converges to stable equilibrium responses within 10 iterations. At equilibrium, deficit irrigation emerges as the dominant adaptation strategy in California, with irrigation levels stabilizing at approximately 70% of full irrigation demand, while Arizona and New Mexico experience stronger yield sensitivities. Early iterations produce commodity price increases of up to 10% for fruit and vegetable crops; however, these responses moderate as land allocation and production patterns adjust across regions. Deficit irrigation and spatial reallocation of irrigated land partially offset production losses, with variability observed across different states: California maintains yields primarily via deficit irrigation, whereas Arizona and New Mexico will rely mainly on reducing irrigated area to absorb the shock. By capturing feedback between biophysical and economic processes, this approach highlights how irrigation strategies and land-use decisions evolve under water stress and provides a transferable platform for evaluating water management policies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Development and modeling for a small-scale, rapidly heated high explosives initiation time (HEIT) experiment

Most small-scale assessments of explosive sensitivity, including the popular drop-weight impact test, convolute thermal and mechanical phenomena to the extent that it has been extremely challenging to decipher how an explosive ignites and propagates reactions. For instance, an impact generates heat through a range of dissipation mechanisms, which can in turn, depending on the reaction rates of the explosive, lead to chemical decomposition. To deconvolute the various contributions to the sub-shock initiation and propagation of explosive reactions, we describe the development and modeling of the High Explosives Initiation Time (HEIT) test - a new, small-scale, high-throughput experiment designed to rapidly heat milligram quantities of energetic materials confined within small diameter steel needles. Specifically, we have modeled and designed a 250 joule pulsed power system capable of rapidly delivering electrical current to the needles, resulting in rapid heat delivery to the sample. In conclusion, the energy deposition rate into the sample is controlled by different transmission line topologies. Modeling in COMSOL is performed to understand the energy required to heat up the explosive sample.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗