SEARCH · Search NASA

Results for “detonation velocity”

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.

SourceResearch connectionFrom yearThrough year

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

106 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Seismic velocity modeling and earthquake relocations for the southern Nevada National Security Site (NNSS), with a focus on Rock Valley

The Rock Valley Direct Comparison (RV/DC) is the third phase of the Source Physics Experiment (SPE; Snelson et al., 2014), a project aimed at improving seismic discrimination between explosive sources and natural earthquakes. Phases I and II of SPE focused on studying the generation of shear waves in endmember geologies (granite and alluvium, respectively), while Phase III aims to detonate a chemical explosion co-located with the source region of the 1993 Rock Valley earthquake sequence, located in the Rock Valley Fault Zone (RVFZ), in the southern portion of the Nevada Test Site (now the Nevada National Security Site, or NNSS).

58 GEOSCIENCES↗

DOE OSTI2026

AWSD Reactive Burn Model for the HMX‐Based High Explosive LX‐04

An Arrhenius–Wescott–Stewart–Davis (AWSD) reactive burn model is applied to describe shock initiation and detonation properties of the HMX-based high explosive LX-04. The parameters in the model are calibrated to data from multiple sources. The thermodynamic equations of state used in the model are calibrated to a combination of thermochemical calculations for HMX and LX-04 as well as experimentally-measured cylinder expansion results for LX-04. The kinetic parameters are calibrated to velocity data from gas gun experiments performed using EDC-32—a high explosive with the same chemical composition as LX-04 but different structural properties, and scaled rate stick data for PBX 9012. The AWSD model is shown to accurately describe the shock initiation and propagation of LX-04. Very good agreement is observed between the available experimental data and the AWSD model output. The presented results constitute an accurate LX-04 reactive burn model for use in engineering-scale models and simulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2023

A study of explosive-induced fracture in polymethyl methacrylate (PMMA)

The fracture response of geologic materials is of interest for applications, including geothermal energy harnessing and containment of underground explosions. To better understand the explosively induced fracture response of geomaterials, polymethyl methacrylate (PMMA) was used as a transparent rock surrogate to allow imaging of internal shock propagation and fracture growth processes. Experiments were conducted using high-speed shadowgraphy and photon Doppler velocimetry (PDV), which were compared to numerical simulations. Experiments measured fractures produced in 304.8mm × 304.8mm × 304.8mm PMMA cubes with two simultaneously initiated detonators. The cubes were subjected to varying amounts and directions of externally applied uniaxial stresses, including no stress, 2 MPa stress, and 20 MPa stress. The fracture radius as a function of time was extracted from the high-speed videos. Post-test images of the PMMA cubes aided in the determination of three-dimensional effects not directly imaged by the cameras. The surface velocity history and the shock response captured in PDV and the high-speed videos were compared to the simulated explosive-induced shock response. The simulation results indicate that the shock drives the fracture for the first 20 μs corresponding to a fracture radius of approximately 15 mm in the experiments. The gas-driven fracture extent was estimated analytically using an equilibrium stress distribution calculated after the shock wave propagation through the sample. Reduction in the gas pressure due to the leakage of the explosive products through the crack as a function of time was accounted for. In conclusion, the estimated fracture lengths were in agreement with the experimentally observed fracture lengths.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Numerical modeling of plasma assisted deflagration to detonation transition in a microscale channel

Here, this work numerically studies the plasma assisted deflagration to detonation transition (DDT) of H 2 /O 2 mixtures in a microscale channel with detailed chemistry and transport. The results show that the DDT onset time is non-monotonically dependent on the discharge pulse number. The DDT is accelerated with small pulse numbers, whereas retarded with large ones. Two different DDT regimes, respectively at a small and large plasma discharge number, via acoustic choking of the burned gas and plasma-enhanced reactivity gradient without acoustic choking, are observed. Without plasma discharge, pronounced pressure and temperature gradients in front of the flame are generated by acoustic compression after the choking of the burned gas, triggering DDT via autoignition. With small plasma pulse numbers, the plasma-generated species enhance the ignition kinetics and lead to an increased reactivity in the boundary layer. After the choking of the burned gas, the plasma-enhanced reactivity advances the sequence of autoignition near the wall, strengthens ignition-shock wave coupling, and accelerates DDT. However, with a large discharge pulse number, a direct autoignition initiating DDT can occur without the acoustic choking of the burned gas due to the strongly accelerated reactivity and elevated temperature. In this case, DDT onset is retarded because the elevated temperature increases sonic velocity and the increased reactivity accelerates fuel oxidation in front of the flame, decelerating the formation of a leading shock and subsequent pressure buildup ahead of the flame. The present modeling reveals that no matter with or without plasma discharge, DDT is initiated by autoignition in thermal, pressure, and reactivity gradient fields via the Zel'dovich gradient mechanism. The acoustic choking of the burned gas may not be the necessary condition of DDT with strong plasma-enhanced reactivity gradient. This work provides an answer to the experimentally observed non-monotonic DDT onset time by plasma, which provides guidance to control DDT in advanced detonation engines and fire safety of hydrogen-fueled catalytic reactors in microchannels by non-equilibrium plasma discharge.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2020

Numerical Considerations in the Modeling of a High Explosive Cylinder Experiment Using an ALE Continuum Mechanics Code

A series of experiments involving the detonation of PBX 9501 encased in a copper cylinder are modeled with the objective of evaluating a proposed set of phenomenological parameters for the Wescott–Stewart–Davis reactive burn model. The numerical analysis is conducted using the Los Alamos continuum mechanics code FLAG. Numerical considerations pertaining to various aspects of modeling the experiments using FLAG are discussed. It is shown that use of the proposed set of phenomenological parameters results in predictions of free-surface velocity that match empirically measured velocities reasonably well.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Modification of Jet Velocities in an Explosively Loaded Copper Target with a Conical Cavity

Here, in this work, the design and execution of an experiment with the goal of demonstrating control over the evolution of a copper jet is described. Simulations show that when using simple multi-material buffers placed between a copper target with a conical cavity and a cylinder of high-explosive, a variety of jetting behaviors occur based on material placement, including both jet velocity augmentation and mitigation. A parameter sweep was performed to determine optimal buffer designs in two configurations. Experiments using the optimal buffer designs verified the effectiveness of the buffers at altering jet velocities. Similar trends were shown between the experimental results and the modeling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

A generalization of the shock invariant relationship

Shock invariant relationship, which was conceived for inert shock waves to derive the 4th power relationship between shock pressure and maximum strain rate, is generalized for reactive shock waves such as Chapman–Jouget detonation and shock-induced vaporization. The generalization, based on the first-order reaction models, is a power function relationship between overall dissipated energy (Δe dis ) and reaction time Δτ such that Δe dis Δτ 1/α = constant, where the power coefficient α is found to be in the range of 2/3–4. Experimental data, though scarce, are consistent with the generalization. Implication of the generalization for inert shocks is also considered and suggests a broad range of the 4th power coefficient including an inequality equation that constrains the shock and particle velocity relationship.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Violent mergers revisited: The origin of the fastest stars in the Galaxy

Binary systems of two carbon-oxygen white dwarfs are one of the most promising candidates for the progenitor systems of Type Ia supernovae. Violent mergers, where the primary white dwarf ignites when the secondary white dwarf smashes into it while being disrupted on its last orbit, were the first double degenerate merger scenario proposed that ignites dynamically. However, violent mergers likely contribute only a few percent to the total Type Ia supernova rate and do not yield normal Type Ia supernova light curves. Here we revisit the scenario, simulating a violent merger with better methods and, in particular, a more accurate treatment of the detonation. We find good agreement with previous simulations but with one critical difference: The secondary white dwarf being disrupted and accelerated towards the primary white dwarf and impacted by its explosion does not fully burn, and its core survives as a bound object. The explosion leaves behind a 0.16 M ⊙ star travelling 2800 km/s, making it an excellent (and so far the only) candidate to explain the origin of the fastest observed hypervelocity stars. We also show that before the explosion, 5 × 10 −3 M ⊙ of material predominantly consisting of helium, carbon, and oxygen had already been ejected at velocities above 1000 km/s. Finally, we argue that if a violent merger made the hypervelocity stars D6-1 and D6-3 and violent mergers require the most massive primary white dwarfs in binaries of two carbon-oxygen white dwarfs, there has to be a much larger population of white dwarf mergers with slightly lower mass primary white dwarfs. Because this population likely represents ≫10% of the Type Ia supernovae rate, it can essentially only give rise to normal Type Ia supernovae.

Astronomy and AstroPhysics↗

DOE OSTI2022

Date Release Report for Large Surface Explosion Coupling Experiment (LSECE) Nevada National Security Site

The DTRA sponsored 2020-2022 Large Surface Explosion Coupling Experiment (LSECE) consists of two large ground surface chemical explosions, data collection, analysis and modeling carried out in 2020-2022. The LSECE chemical explosions were carried out at the site of prior NNSA sponsored buried chemical explosions that were part of the Source Physics Experiment (SPE) Phase II in Dry Alluvium Geology (DAG). This report described the data collected under LSECE that is being made publicly available. The prior buried explosion DAG data are described in a separate report (https://www.osti.gov/biblio/1825534) and are also publicly available. The LSECE explosions, data analysis and numerical modeling work sponsored by DTRA were intended to address three objectives: 1) Generate seismo-acoustic data to test and improve numerical models of explosion energy coupling in dry alluvium geology; 2) To improve seismo-acoustic yield estimation techniques as a function of depth and medium properties; 3) To study and improve acoustic propagation modeling under two different atmospheric conditions. The LSECE data that were collected have expanded the prior DAG buried explosion dataset to include surface chemical explosions recorded on a common set of seismo-acoustic stations. They have allowed a more detailed study of above and below surface seismo-acoustic energy coupling with applications to explosion monitoring and assessment. The two approximately 1-ton TNT equivalent yield LSECE explosions were detonated at different times of the day to explore the effects of the different atmospheric conditions. The first chemical explosion “Artemis” was conducted before dawn when temperature inversions were present. The second chemical explosion “Apollo” was conducted on a sunny afternoon when the temperature gradient was more linear. The LSECE chemical explosion data were also collected across a variety of different sensor types to allow evaluation of their effectiveness in recovering useful information. The LSECE instrumentation included fiber optic or Distributed Acoustic Sensing (DAS), a dense array or Large-N array of seismometers, borehole accelerometers, DAS and a velocity meter, airborne acoustic instruments and a variety of visual and remote sensing data.

58 GEOSCIENCES↗

DOE OSTI2020

The Spectacular Ultraviolet Flash from the Peculiar Type Ia Supernova 2019yvq

Early observations of Type Ia supernovae (SNe Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. Here, we present exquisite observations of SN 2019yvq, the second observed SN Ia, after iPTF 14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN 2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for a SN Ia ($M_g \approx -18.5$ mag at peak) yet featured very high absorption velocities ($v \approx 15,0000$ km s -1 for Si II λ6355 at peak). We find that many of the observational features of SN 2019yvq, aside from the flash, can be explained if the explosive yield of radioactive 56 Ni is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN 2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of 56 Ni in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN 2019yvq. In closing, we predict that the nebular spectra of SN 2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Flow and Performance Characterization of Rotating Detonation Combustor Integrated with Various Convergent Nozzles

In this study, convergent nozzles of various area ratios (ARs) are used downstream of an annular rotating detonation combustor (RDC) to increase the operating pressure and approach sonic conditions at the nozzle throat. Reactant methane and oxygen-enriched air (67% [Formula: see text] and 33% [Formula: see text] by volume) are supplied in counterflow arrangement from two separate plenums located at the base of the RDC annulus. Based on experimentation, a total mass flow rate of [Formula: see text] was chosen to achieve stable, single-wave mode RDC operation for all test cases, allowing for one-to-one comparisons. The internal performance of the RDC was characterized by ion probes and pressure measurements (wall static and oscillating) in supply plenums and across different axial locations of the combustor. Particle image velocimetry (PIV) at 100 kHz was utilized to measure axial and circumferential velocity components within a two-dimensional region of interest located downstream of the converging nozzle exit. Results show higher internal performance of the RDC with increasing AR of the convergent nozzle. PIV measurement illustrated that the flow oscillation amplitudes decrease with an increasing AR of the converging nozzle. The exit flow contained significant nonuniformity and unsteadiness even with a converging nozzle of AR 2.0, indicating incomplete choking of the flow at the nozzle throat.

Engineering↗

DOE OSTI2025

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Automated nuclear cloud feature extraction from film

Chemical, biological, radiological, nuclear, and explosives incidents require rapid detection and characterization for appropriate response. For a nuclear detonation, visible-light cameras may be used to locate the cloud and characterize fallout deposition when coupled with numerical models. Films from the United States’ nuclear testing era compose the only sizeable collection of imagery depicting high-yield detonations. These films offer unique insights into characteristics of flows involving scales that are difficult to replicate experimentally, and they are a valuable source of data for the validation of models for nuclear fallout transport, either as part of emergency response or forensic activities. In this work, we implement modern computer vision and machine learning techniques to identify and track the cloud automatically and subsequently determine the time dependence of some of its features. We trained a ResNet-18 image classifier on hundreds of images to categorize nuclear cloud morphology. Each category or cloud regime is determined by early cloud evolution and is associated to constitutive properties of the flow, such as distribution of vorticity. Next, we identified keypoint features using the KAZE algorithm and tracked these keypoints in the images, allowing us to determine the dimensions and velocities of the cloud across film frames. These measurements converted to real-world units provide valuable experimental data that can be used in the development and validation of nuclear cloud models. We compared the results of this method against manual cloud rise measurements from two different films. In one, our automated method accelerated the feature extraction process without sacrificing measurement accuracy.

Khristy, Joel [ORNL] (ORCID:0000000209963060)↗

DOE OSTI2023

Performance Characteristics of a Rotating Detonation Combustor Exiting Into a Pressurized Plenum to Simulate Gas Turbine Inlet

Abstract The present study aims to experimentally characterize the performance of a rotating detonation combustion (RDC) system integrated with a pressurized downstream plenum to simulate the high-pressure inlet conditions of power-generating gas turbines. A thorough understanding of the operational behavior including wave mode behavior, static pressure profile along the combustor length, and dynamic features of pressure fluctuations is crucial for successful integration of RDC with the turbine. In this study, two RDC configurations are investigated, RDC with a constant area annulus and RDC with a converging nozzle. In both cases, the RDC flow exited into a plenum chamber kept at pressures varying from 155 kPa to 330 kPa. RDC was operated on methane and oxygen-enriched air to represent reactants used in land-based power generation. Experiments were conducted for the two RDCs configurations operated at three reactant mass flow rates (0.23, 0.32, and 0.46 kg/s). The RDC performance is characterized by time-averaged static pressure measurements, and wave velocity determined by ionization probes. In addition, dynamic pressure measurements were recorded both inside and near the exit of RDC channel to investigate wave interactions between RDC and downstream plenum. Results show that the RDC with the converging nozzle achieved superior performance while minimizing detrimental interactions with the reflected shock and/or acoustic waves from the downstream plenum.

Engineering↗

DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Shock initiation of low density polymer bonded explosive LX-14: A study of two morphologies

A series of six shock initiation experiments have been carried out on low density LX-14 powder in order to simulate a shock insult on the heavily damage polymer bonded explosive, LX-14. Two distinct morphologies were studied (tap density molding powder and machined swarf), both at the same density of 0.942 g/cm 3 , or 50.1% theoretical maximum density. The purpose of these experiments was to provide data to help make an assessment of the effects that damage has on the material sensitivity to a planar shock. This was achieved primarily by providing shock sensitivity data in the form of a Pop plot, and also reactants equation of state data, which aids in the determination of input conditions for both the experiments performed in this work and also future experiments of this material type. The experiments were of a cut-back format, consisting of four sample heights on each shot and diagnosed with optical velocimetry. The experiments were carried out at the Technical Area 40 Chamber 9 gas gun facility at the Los Alamos National Laboratory, where the LX-14 targets were subjected to Al 6061 and Oxygen Free High Conductivity copper impactors launched to velocities up to 2.14 km/s. Time corrected reactive growth wave profiles are presented in this paper along with the derivation of the following Hugoniot parameters for this explosive, where the molding powder and machine swarf are represented by the following linear equations, respectively, U s = 1.64 (±0.64)u p + 2.25 (±0.49) and U s = 0.87 (±.0.60)u p + 4.08 (±0.49). The results show that the steady increase in shock sensitivity with increasing void fraction reaches an inflection point beyond which the shock sensitivity begins to decease. This inflection point lies between 65% and 50% of the theoretical maximum density of the LX-14. In both cases, the damaged LX-14 was not as sensitive as expected, with a relative shock sensitivity of the molding powder being less than the pressed LX-14, and the machined swarf having a shock sensitivity that is comparable to pressed LX-14.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗