Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock chemistry”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Shock-induced chemistry and high strain-rate viscoelastic behavior of a phenolic polymer

We use impact experiments and a finite element model (up to 1.2 GPa), and molecular dynamics simulations (up to 60 GPa), to examine the behavior of a phenolic polymer under shock compression, spanning both nonreactive and reactive regimes. In the nonreactive regime, relaxation following compression at strain rates of ∼105 s−1 can be explained by viscoelasticity observed at ordinary laboratory rates (≲1 s−1) by accounting for the temperature dependence of the phenolic β-transition. Reasonable agreement is found between the measured shock Hugoniot up to 1.2 GPa and molecular dynamics simulation for cross-linked structures of comparable density. We also observed a first-order mechanical transition near 0.36 GPa shock stress and estimated a spall strength of 0.102 GPa and Hugoniot elastic limit of 1–2 GPa. The shock stress is found to vary up to 24% among phenolics made with different resin and/or cure processes. Finally, molecular dynamics simulations are used to identify a reactive regime at shock pressures ≳20 GPa that is characterized by chemically driven, rate-dependent relaxation processes, including dehydrogenation and dehydration reactions that promote the formation of a dense, highly cross-linked carbonaceous solid and the release of light volatiles.

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

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗
DOE OSTI2023

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. Finally, this suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

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

Stability of the fcc phase in shocked nickel up to 332 GPa

Despite making up 5-20 wt.% of Earth’s predominantly iron core, the melting properties of elemental nickel at core conditions remain poorly understood, due largely to a dearth of experimental data. We present here an in situ X-ray diffraction study performed on laser shock-compressed samples of bulk nickel, reaching pressures up to ~ 500 GPa. Hugoniot states of nickel were targeted using a flat-top laser drive, with in situ X-ray diffraction data collected using the Linac Coherent Light Source. Rietveld methods were used to determine the densities of the shocked states from the measured diffraction data, while peak pressures were determined using a combination of measured particle velocities, shock transit times, hydrodynamic simulations, and laser intensity calibrations. We observed solid compressed face-centered cubic (fcc) Ni up to at least 332 ± 30 GPa along the Hugoniot—significantly higher than expected from the majority of melt lines that have been proposed for nickel. We also bracket the partial melting onset to between 377 ± 38 GPa and 486 ± 35 GPa.

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

Chemical templates of the Central Molecular Zone

Context . The Central Molecular Zone (CMZ) of the Milky Way exhibits extreme conditions, including high gas densities, elevated temperatures, enhanced cosmic-ray ionization rates, and large-scale dynamics. This makes it a perfect laboratory for astrochemical studies. With large-scale molecular surveys revealing increasing chemical and physical complexity in the CMZ, it is essential to develop robust methods to decode the chemical information embedded in this extreme region. Aims . A key step to interpreting the molecular richness found in the CMZ is building chemical templates tailored to its diverse conditions. In particular, understanding how CMZ environments affect shock and protostellar chemistry is crucial. The combined impact of high ionization, elevated temperatures, and dense gas remains insufficiently explored for observable tracers. Methods . For this study, we utilized UCLCHEM , a gas-grain time-dependent chemical model, to link physical conditions with their corresponding molecular signatures and identify key tracers of temperature, density, ionization, and shock activity. To achieve this, we ran a grid of models of shocks and protostellar objects representative of typical CMZ conditions, focusing on 24 species, including complex organic molecules. Results . Shocked and protostellar environments show distinct evolutionary timescales (≲10 4 vs. ≳10 4 years); 300 K emerges as a key temperature threshold for chemical differentiation. We find that cosmic-ray ionization and temperature are the main drivers of chemical trends. HCO + , H 2 CO, and CH 3 SH trace ionization, while HCO, HCO + , CH 3 SH, CH 3 NCO, and HCOOCH 3 show consistent abundance contrasts between shocks and protostellar regions over similar temperature ranges. Conclusions . We characterized the behavior of 24 species in protostellar and shock-related environments. While our models underpredict some complex organics in shocks, they reproduce observed trends for most species, supporting scenarios involving a need for recurring shocks in Galactic Center clouds and enhanced ionization toward Sgr B2(N2). Future work should assess the role of shock recurrence and metallicity in shaping chemistry.

Galaxy: center↗
DOE OSTI2023

Polyimide dynamically compressed to decomposition pressures: Two-wave structures captured by velocimetry and modeling

We performed a series of six plate impact experiments on polyimide and modeled them using new reactant and products equations of state combined with an Arrhenius rate model. The first experiment was diagnosed with embedded electromagnetic velocity gauges through which we directly observed attenuation of the lead shock to an approximately constant state over a propagation distance of roughly 4 mm. Simulated gauge profiles were in excellent qualitative agreement with experiment and suggested a sluggish chemical reaction that did not proceed to completion. The remaining five experiments were conducted in a transmission geometry and diagnosed velocimetrically at the sample/window interface. All five of these yielded profiles with a sharp shock followed by a more gradual approach to maximum interface velocity that was “rounded” to varying degree. These profiles proved difficult to interpret unambiguously due to the convolution of the reactive wave upon first shock with reflection of the lead wave and reshock or release by the window. Comparison with thermochemical calculations strongly suggests that the point of maximum interface velocity corresponds to the equilibrium reshock or release locus. We discuss the implications of this point for the practice of impedance matching based on the reflected Hugoniot of reactive materials such as polymers. The reactant and thermochemical products equations of state are developmental SESAME tables 97710 and 97720, respectively.

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

Rapid activation of non-oriented mechanophores via shock loading and spallation

Mechanophores (MPs), stimuli-responsive molecules that respond chromatically to mechanochemical reactions, are important for understanding the coupling between mechanics and chemistry as well as in engineering applications. However, the atomic-level understanding of their activation originates from gas phase studies or under simple linear elongation forces directly on molecules or polymer chains containing MPs. The effect of many-body distortions, pervasive in condensed-phase applications, is not understood. Therefore, we performed large-scale molecular dynamics (MD) simulations of a poly(methyl methacrylate)-spiropyran copolymer under dynamic mechanical loading and studied the activation of the MP under various conditions from dynamical compression to tension during unloading. Detailed analysis of the all-atom MD trajectories shows that the MP blocks experience significant many-body intramolecular distortion that can significantly decrease the activation barrier as compared with when deformation rates are slow relative to molecular relaxation time scales. Here we find that the reactivity of MPs under material compression states is governed by many-body effects of intramolecular torsions, whereas under tension, the reactions are governed by tensile stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Initiation and Carbene Induced Radical Chain Reactions in CH 2 F 2 Pyrolysis

High temperature dissociations of organic molecules typically involve a competition between radical and molecular processes. In this work, we use a modeling, experiment, theory (MET) framework to characterize the high temperature thermal dissociation of CH 2 F 2 , a flammable hydrofluorocarbon (HFC) that finds widespread use as a refrigerant. Initiation in CH 2 F 2 proceeds via a molecular elimination channel; CH 2 F 2 →CHF+HF. Here we show that the subsequent self-reactions of the singlet carbene, CHF, are fast multichannel processes and a facile source of radicals that initiate rapid chain propagation reactions. These have a marked influence on the decomposition kinetics of CH 2 F 2 . The inclusion of these reactions brings the simulations into better agreement with the present and literature experiments. Additionally, flame simulations indicate that inclusion of the CHF+CHF multichannel reaction leads to a noticeable enhancement in predictions of laminar flame speeds, a key parameter that is used to determine the flammability of a refrigerant.

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

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

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

Chemical decomposition in shock-compressed 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals: Time-resolved Raman spectroscopy

Understanding the evolution of chemical decomposition—at the molecular level—in shock-compressed insensitive high explosive (IHE) single crystals is an important need. Toward this goal, time-resolved Raman spectra were measured in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals shock compressed to stresses above 20 GPa. At 22 GPa and higher stresses, the Raman peak intensities were significantly reduced, showing the onset and extent of chemical changes. At 33 GPa, no Raman peaks were observed, suggesting complete decomposition of FOX-7. The 22 GPa onset from Raman data is in marked contrast to the much higher (32 GPa) onset determined previously from wave profile measurements—showing the significantly greater sensitivity of Raman measurements for gaining insight into shock-induced molecular-level changes in IHE single crystals. The close match between the complete loss of Raman peaks at 33 GPa and the 32 GPa reaction threshold determined from continuum data makes a good case that the two results are related and they reflect significant energy release due to decomposition.

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

Direct observation of diamond formation in a shock-compressed high explosive

Understanding the formation timescale and structure of carbonaceous reaction products is critical for modeling the high-pressure equation-of-state of organic materials. We use the National Ignition Facility to shock-compress polycrystalline TATB (C 6 H 6 N 6 O 6 ) samples to ~70–130 GPa and ~4000–5500 K, employing in situ nanosecond X-ray diffraction to probe reaction products and velocimetry to measure transmitted compression wave profiles. Our diffraction data is consistent with the formation of diamond over timescales less than ~60 ns. This represents carbon condensation from a molecular explosive on timescales three times faster than previously reported and the earliest observation of diamond produced from reacting TATB. Reactive flow simulations with explicit chemistry reproduce the observed temporal structure within wave profiles to inform the distribution of P-T states. These findings provide direct evidence of ultrafast diamond formation in a reactive system at extreme conditions and provide new constraints for models of shock and detonation chemistry.

Clarke, Samantha M. [Lawrence Livermore National L↗
DOE OSTI2024

TCF - Sustainable Well Cement for Geothermal, Thermal Recovery, and Carbon Storage wells

Primary cementing is the most important operation performed on a subterranean well. Cement, placed in the annulus between the casing and the formations, serves as a hydraulic seal preventing fluids and gas migrations, protecting steel casing from corrosion, and supporting the well structure. Poor cementing jobs can be directly responsible for wells not reaching their full capacity, casing corrosion, compromised well integrity and, in the worst-case scenarios, well collapse. Geothermal, thermal recovery, and carbon storage wells offer environments that are especially difficult for cements to survive while the required lifespan of these wells can be years. In these wells currently used cements cannot provide durable well integrity, and new solutions are necessary. The necessity of stabilizing the electric grid, increasing its flexibility, and providing energy on demand will further expand the market of durable cement solutions for applications in high-temperature underground reservoir thermal energy storage wells (HT RTES). The design challenges of special cement systems for such wells originate from chemistry limitations of currently used well cements, aggressive environments, temperature and repeated shock conditions associated with them, and very weak formations that are not uncommon in these wells. During the last 70 years the most common systems for thermal-well construction have been Ordinary Portland Cement (OPC) (absolute majority of the wells), silica-lime system, and high aluminum cement-containing formulations. The major issues of calcium-silicate hydrates, that form during the hydration of OPC resulting in hardened material, are their poor chemical resistance to acids due to the calcium interactions with acid anions followed by its eventual dissolution, even in mild acids, such as carbonic acid from CO 2 dissolution, and inadequate bonding to the steel causing serious casing corrosion problems (Sugama & Pyatina, 2019). Performance of OPC-based cement may be significantly improved with organic additives; however, those are limited by their temperature stability. An alternative to calcium-silicate cement was developed by BNL and commercialized by Halliburton as ThermaLock TM cement. Calcium-aluminate-phosphate-based chemistry of the material allowed overcoming acid-resistance problems of OPC, especially in CO 2 -rich geothermal environments. In 2012-2015, building upon this experience, BNL developed Thermal Shock Resistant Cement (TSRC) with the support of Geothermal Technology Office (GTO) of DOE. TSRC has high-temperature stable chemistry, superior properties of cement-metal casing bond and corrosion protection (4-times better than currently used OPC-based formulation), significantly outperforms common well-cements in thermal-shock tests, has self-healing properties, and like ThermaLock TM is CO 2 resistant incorporating carbonate ions into stable hydrates (Gill et al., 2012).

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Pyrolysis of Cyclohexane and 1-Hexene at High Temperatures and Pressures—A Photoionization Mass Spectrometry Study

Cycloalkanes are important components of a wide range of fuels. However, there are few experimental data at simultaneously high temperatures and pressures similar to those found in practical systems. Such data are necessary for developing and testing chemical kinetic models. In this study, data relevant to cycloalkane pyrolysis were obtained from high repetition rate shock tube experiments coupled with synchrotron-based photoionization mass spectrometry diagnostics. The pyrolysis of cyclohexane was studied over 1270–1550 K and ~9 bar, while the more reactive primary decomposition product, 1-hexene, was studied at 1160–1470 K and ~5 bar. Insights into the decomposition of the parent molecules, the formation of primary products and the production of aromatic species were gained. Simulations were performed with models for cyclohexane and 1-hexene that were based on literature models. The results indicate that over several hundred microseconds reaction time at high pressures and temperatures the pyrolysis of cyclohexane is largely dominated by reactions initiated by cyclohexyl radicals. Furthermore, good agreement between the simulations and the experiments were observed for cyclohexane and 1-hexene with a modified version of the cyclohexane model. Conversely, the 1-hexene model did not reproduce the experimental observations.

1-alkenes↗
DOE OSTI2022

Effect of lot microstructure variations on detonation performance of the triaminotrinitrobenzene (TATB)-Based insensitive high explosive PBX 9502

PBX 9502 is an important insensitive high explosive due to its combination of safety properties and detonation performance. It is a polymer-bonded formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the high explosive crystal, bound with Kel F-800 (FK-800), a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Two different types are used, one known as virgin that uses only pristine manufactured TATB, and a second known as recycled that has 50 wt.% of its TATB reclaimed from machining scraps of previously pressed virgin PBX 9502. Recycled lots have a higher percentage of fine particles compared to virgin lots, due to the fracturing and damage sustained by TATB crystals during pressing. We examine the influence that this TATB microstructure difference between virgin and recycled PBX 9502 has on detonation performance properties. New rate-stick geometry diameter effect, detonation front shape and cylinder expansion test data are obtained for two previously uncharacterized virgin lots and one recycled lot of PBX 9502. This is combined with previously published data for one virgin and one recycled lot for evaluation purposes. Detonation shock dynamics model calibrations are conducted on each of the five lots to provide an assessment of the detonation timing characteristics of virgin versus recycled PBX 9502 lots. For two of the virgin lots, detonations propagate slower in the rate-stick geometry than those in the two recycled lots. However, the other virgin lot tested has propagation rates comparable to that of the recycled lots for larger diameter rate-sticks. The latter result, though, is shown to be geometry dependent and depends on the range of detonation curvatures accessed in different geometries. New copper-confined, cylinder expansion tests are conducted on each of the three virgin and two recycled lots to obtain detonation product Jones-Wilkins-Lee equations of state, enabling an assessment of the metal push capabilities for each of the five lots. We find that the metal push capabilities, characterized by the evolution of the heat of detonation with volume, are similar between the virgin and recycled lots. Thus, changes in microstructure between different PBX 9502 lots seemingly affect the rate of reaction, but not the overall energy content.

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

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗