Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “detonation wave”

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 379 records · Page 21

Detonation of insensitive high explosives by a Q-switched ruby laser.

Immediate longitudinal detonations have been observed in confined small-diameter columns of PETN, RDX, and tetryl by using a focused Q-switched ruby laser. The energy ranged from 0.8 to 4.0 J in a pulse width of 25 nsec. A 1000-A-thick aluminum film deposited on a glass window was used to generate a shock wave at the window-explosive interface when irradiated by the laser. In some cases, steady-state detonations were reached in less than .5 microsec with less than 10% variation in the detonation velocity.

Yang, L. C.↗

Modeling for pRad Explosive Experiments

The hydrodynamics code CTH was used to simulate high-explosive detonations in various geometries to design a planar, high-pressure drive with Taylor wave release for experiments in the proton radiography containment vessel. The most successful geometry is based on a P25 lens with some of its outer diameter and output end removed. Shock planarity was significantly improved relative to an unmodified P25 lens for this case, where its full diameter was not needed, and the calculated timing was adequate to trigger the experiments successfully.

36 MATERIALS SCIENCE↗

Artificial stimulation of auroral electron acceleration by intense field aligned currents

A cesium-doped high explosion was detonated at 165 km altitude in the auroral ionosphere during quiet conditions. An Alfven wave pulse with a 200-mV/m electric field was observed, with the peak occurring 135 ms after the explosion at a distance of about 1 km. The count rate of fixed energy 2-keV electron detectors abruptly increased at 140 ms, peaked at 415 ms, and indicated a downward field-aligned beam of accelerated electrons. An anomalously high-field aligned beam of backscattered electrons was also detected. The acceleration is interpreted as due to production of an electrostatic shock or double layer between 300 and 800 km altitude. The structure was probably formed by an instability of the intense field-aligned currents in the Alfven wave launched by the charge-separation electric field due to the explosion.

Holmgren, G.↗

Noninvasive acoustic time-of-flight measurements in heated, hermetically-sealed high explosives using a convolutional neural network

In this work, we present a data-driven technique for measuring the time-of-flight through material sealed within a container. Time-of-flight measurement provides a noninvasive means of quantifying the sound speed profile within a material by transmitting an acoustic burst and then measuring the time required for the burst to arrive at an opposing receiver. In a hermetically-sealed cylindrical container, a portion of the acoustic energy propagates through the material as a bulk wave, while the remainder of the acoustic energy propagates around the container walls as guided waves. As a result, interference from the guided waves obscures the bulk arrival, inhibiting measurement of the sound speed. The technique uses a Convolutional Neural Network (CNN) to identify critical features in the measured waveforms and identify bulk wave arrivals. We demonstrate this time-of-flight measurement technique on high explosive-filled containers as they are heated from room temperature to detonation. This is a particularly challenging application for acoustic time-of-flight measurements as the high explosives have significant sound speed gradients as they undergo heating, and they lead to significant attenuation of the bulk wave, as opposed to the guided waves, which do not suffer significant attenuation. We characterize the performance of the CNN as a function of the high explosive temperature and as a function of the CNN hyperparameters. We then provide physical insight into the error trends.

47 OTHER INSTRUMENTATION↗

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↗

A numerical study of the NOL large-scale test

The large-scale gap test was modeled numerically using the reactive hydrodynamics code 2DE with forest fire burn rates. The model gives good agreement between calculated and experimental critical gap values. The calibration of peak pressure in the gap versus gap length is obtained from these calculations, and is in good agreement with the published experiment curve. The two dimensional nature of the gap test is evident from the curvature of the shock wave in the gap and test sample, and from the observed distance of run to detonation in the test samples.

Bowan, A. L.↗

A hierarchy for modeling high speed propulsion systems

General research efforts on reduced order propulsion models for control systems design are overviewed. Methods for modeling high speed propulsion systems are discussed including internal flow propulsion systems that do not contain rotating machinery, such as inlets, ramjets, and scramjets. The discussion is separated into four areas: (1) computational fluid dynamics models for the entire nonlinear system or high order nonlinear models; (2) high order linearized models derived from fundamental physics; (3) low order linear models obtained from the other high order models; and (4) low order nonlinear models (order here refers to the number of dynamic states). Included in the discussion are any special considerations based on the relevant control system designs. The methods discussed are for the quasi-one-dimensional Euler equations of gasdynamic flow. The essential nonlinear features represented are large amplitude nonlinear waves, including moving normal shocks, hammershocks, simple subsonic combustion via heat addition, temperature dependent gases, detonations, and thermal choking. The report also contains a comprehensive list of papers and theses generated by this grant.

Hartley, Tom T.↗

Trinity

At 5:29:45 am Mountain War Time, on July 16, 1945, the world’s first atomic bomb, codenamed Trinity, exploded in a blinding flash 100 feet over a portion of the southern New Mexican desert known as the Jornada del Muerto – the Journey of Death – temporarily blinding a future Nobel Laureate, Richard Feynman. Forty seconds after the detonation, Noble laureate Enrico Fermi dropped small pieces of paper before, during, and after the blast wave passed him. From the lateral dispersion of the paper, he calculated a yield of ten kilotons. Another laureate, I.I. Rabi, won the betting pool for Trinity’s yield. Arriving late, he bought the last available number, eighteen kilotons. The yield was later calculated to be twenty kilotons. On August 9th, a copy of the Trinity device, dubbed Fat Man, was dropped on Nagasaki.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗

Diffraction Measurement of Reaction Products in Shock Compressed TATB on the NIF

The shock-induced reaction of high explosives to gaseous and solid reaction products is a rapid, complex process. Understanding solid reaction product structure and formation kinetics is essential in determining the high-pressure equation of state of these multicomponent systems. We use the National Ignition Facility (NIF) to shock compress ~500-µm thick pressed powder high-explosive TATB samples to 70-135 GPa and collect in situ structural data on detonation byproducts using the TARDIS X-ray diffraction diagnostic. Velocimetry is used to record the transmitted compression wave profile, which is well described by Cheetah hydrocode simulations coupled with a reactive flow model, providing strong evidence of reaction in the TATB sample. While an unambiguous determination of the product phases was not possible owing to the low signal-to-noise quality of the diffraction signal, X-ray diffraction data of the product phases formed within the first 50 ns of the reaction is most consistent with a mixture of amorphous products and crystalline hexagonal diamond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CFD Applications in Support of the Space Shuttle Risk Assessment

The paper describes a numerical study of a potential accident scenario of the space shuttle, operating at the same flight conditions as flight 51L, the Challenger accident. The interest in performing this simulation is derived by evidence that indicates that the event itself did not exert large enough blast loading on the shuttle to break it apart. Rather, the quasi-steady aerodynamic loading on the damaged, unbalance vehicle caused the break-up. Despite the enormous explosive potential of the shuttle total fuel load (both liquid and solid), the post accident explosives working group estimated the maximum energy involvement to be equivalent to about five hundreds of pounds of TNT. This understanding motivated the simulation described here. To err on the conservative side, we modeled the event as an explosion, and used the maximum energy estimate. We modeled the transient detonation of a 500 lbs spherical charge of TNT, placed at the main engine, and the resulting blast wave propagation about the complete stack. Tracking of peak pressures and impulses at hundreds of locations on the vehicle surface indicate that the blast load was insufficient to break the vehicle, hence demonstrating likely crew survivability through such an event.

Baum, Joseph D.↗

Detection of the Large Surface Explosion Coupling Experiment by a Sparse Network of Balloon-Borne Infrasound Sensors

In recent years, high-altitude infrasound sensing has become more prolific, demonstrating an enormous value especially when utilized over regions inaccessible to traditional ground-based sensing. Similar to ground-based infrasound detectors, airborne sensors take advantage of the fact that impulsive atmospheric events such as explosions can generate low frequency acoustic waves, also known as infrasound. Due to negligible attenuation, infrasonic waves can travel over long distances, and provide important clues about their source. Here, we report infrasound detections of the Apollo detonation that was carried on 29 October 2020 as part of the Large Surface Explosion Coupling Experiment in Nevada, USA. Infrasound sensors attached to solar hot air balloons floating in the stratosphere detected the signals generated by the explosion at distances 170–210 km. Three distinct arrival phases seen in the signals are indicative of multipathing caused by the small-scale perturbations in the atmosphere. We also found that the local acoustic environment at these altitudes is more complex than previously thought.

47 OTHER INSTRUMENTATION↗

Parameterized Study of Heat Load Trends in a Subscale Rotating Detonation Rocket Engine

NASA has developed a subscale rotating detonation rocket engine platform, enabling rapid parameterization of components and their associated performances. A major performance metric of interest for RDREs is the total heat absorbed at a given operating condition and wave mode operation. To investigate this, several sets of hardware were produced to have variation of the contraction ratio, injector geometry, and length. A wide range of total mass flow rates and mixture ratios were also explored using gaseous methane/oxygen. All hardware was produced using laser powder bed fusion GRCop-42, GRX-810, or C-103 alloys depending on the component cooling requirements. Outer and inner body hardware were made from GRCop-42 where heat fluxes were expected to be high and were water cooled through integrated coolant channels. Nozzles were made from C-103 which are radiatively cooled. Finally, the injectors were produced using GRCop-42 or GRX-810 as passive cooling via propellant injection has previously found to be sufficient given the small gap width of exposed injection surface. Total heat load and bulk heat flux measurements to the hot walls are reported. Trends in heat load and heat flux are shown with relation to parameters such as chamber pressure, area ratio, and wave mode. Heat loads between the inner and outer bodies and the chamber and shroud section are compared.

Heat Flux↗

Parameterized Study of Heat Load Trends in a Subscale Rotating Detonation Rocket Engine

NASA has developed a subscale rotating detonation rocket engine platform, enabling rapid parameterization of components and their associated performances. A major performance metric of interest for RDREs is the total heat absorbed at a given operating condition and wave mode operation. To investigate this, several sets of hardware were produced to have variation of the contraction ratio, injector geometry, and length. A wide range of total mass flow rates and mixture ratios were also explored using gaseous methane/oxygen. All hardware was produced using laser powder bed fusion GRCop-42, GRX-810, or C-103 alloys depending on the component cooling requirements. Outer and inner body hardware were made from GRCop-42 where heat fluxes were expected to be high and were water cooled through integrated coolant channels. Nozzles were made from C-103 which are radiatively cooled. Finally, the injectors were produced using GRCop-42 or GRX-810 as passive cooling via propellant injection has previously found to be sufficient given the small gap width of exposed injection surface. Total heat load and bulk heat flux measurements to the hot walls are reported. Trends in heat load and heat flux are shown with relation to parameters such as chamber pressure, area ratio, and wave mode. Heat loads between the inner and outer bodies and the chamber and shroud section are compared.

Heat Loads↗

Propagation of wide bandwidth signals through strongly turbulent ionized media

Analytic and numerical techniques are presented which directly address the problem of propagation of wide bandwidth signals through random ionized media. This work is applicable to the problems of satellite communication and space based radar observation through a disturbed ionospheric propagation channel that would result from a high altitude chemical release or nuclear detonation. An analytic solution is presented for the two-position, two-frequency mutual coherence function for spherical wave propagation in the strong scattering limit. This solution is used to derive simple expressions for the power impulse response function and to determine its relationship to the important parameters of decorrelation distance and coherence bandwidth which describe the disturbed propagation channel. Results for mean time delay and time delay jitter are presented and compared to direct simulation results and to other theoretical calculations. A numerical/analytical solution to the parabolic wave equation is presented in the form of a multiple phase-screen (MPS) propagation simulation. In this very general solution technique, the random medium is divided into a finite number of layers. The field fluctuation through each layer is obtained by replacing the layer by a centrally located thin phase-screen, whose statistical properties are determined from the statistics of the electron-density irregularities. The waveform then propagates from phase-screen to phase-screen via an exact solution to the Fresnel-Kirchnoff equation. For wide bandwidth waveforms, numerical solutions are obtained at a number of discrete frequencies centered about the carrier and then time-domain solutions are obtained by Fourier transform techniques. Detailed results are presented for a number of interesting cases including propagation of a 20 MHz bandwidth signal through a finite barium cloud at a carrier frequency of 100 MHz. One of the major uses of the MPS propagation simulation has been to provide realizations of the received signal after propagation through a disturbed channel. The MPS simulation obtains a general solution of the parabolic wave equation under both weak and strong scattering conditions. A second technique to directly obtain realizations of wide bandwidth waveforms is presented here. This technique is limited to the case of strong scattering but requires only a fraction of the computer resources needed for MPS signal generation. Detailed comparisons of the two signal generation techniques are presented.

Random Media↗

Shallow Soil Response to a Buried Chemical Explosion With Geophones and Distributed Acoustic Sensing

Abstract Shallow sediments can respond non‐linearly to large dynamic strains and undergo a subsequent healing phase as the material gradually recovers following the passing of seismic waves. This study focuses on the physical changes in the subsurface caused by the shaking from a buried chemical explosion detonated in a borehole in Nevada, USA, as a part of the Source Physics Experiment Phase II. The explosion damaged the shallow subsurface and modified the frequency content recorded by 491 geophones and 2240 Distributed Acoustic Sensing (DAS) channels within 2.5 km from surface ground zero. We observe a gradual shift of resonance frequencies in the 10–25 Hz frequency band in the hours following the explosion and develop a method to characterize the related logarithm‐type healing process of the shallow (i.e., upper ∼25 m) subsurface. We find that stronger levels of ground motion increase the relative degree of damage and duration of the subsurface healing; with the spall region exhibiting the largest degree of damage and longest healing recovery time. We observe coherent spatial patterns of damage with the region located to the southeast of the explosion exhibiting more damage than the southwest region. This study demonstrates that both DAS and co‐located geophones capture similar temporal changes associated with the physical processes occurring in the subsurface, with the high‐density sampling of DAS measurements enabling a new capability to monitor the fine‐scale changes of the Earth's shallow subsurface following the detonation of a buried explosion.

Viens, Loïc↗

Shock initiation of the HMX-based explosive PBX 9012: Experiments, uncertainty analysis, and unreacted equation-of-state

In this study, shock initiation experiments have been carried out on the polymer-bonded explosive PBX 9012 [nominally 90% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazine (HMX), 10% vinylidene–hexafluoropropylene copolymer (Viton A) by weight] in order to provide calibration data for the explosive’s unreacted equation-of-state (EOS) and shock-to-detonation transition for reactive burn rate calibration. The input pressures covered the range of 1.86–4.43 GPa. This provided run-to-detonation depths ranging from > 22.3 mm for the lowest pressure shot to 4.91 mm at the highest pressure. The relative shock sensitivity of PBX 9012 is compared to other HMX-explosives in terms of the Pop-plot, showing that the studied explosive is more sensitive than other similar HMX-based counterparts (with notable exceptions). The uncertainty in the shock velocity determinations from the shock tracker measurements are also investigated, yielding new uncertainty measures in the generated Hugoniot data and run-to-detonation coordinates. Finally, the unreacted equation-of-state is determined using a linear U s – u p Mie–Grüneisen relation and the Davis reactants EOS analytical form, the latter being more suited for reactive burn modeling.

36 MATERIALS SCIENCE↗