Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock initiation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Annular arc accelerator shock tube

An annular arc accelerator shock tube employs a cold gas driver to flow a stream of gas from an expansion section through a high voltage electrode section to a test section, thus driving a shock wave in front of it. A glow discharge detects the shock wave and actuates a trigger generator which in turn fires spark-gap switches to discharge a bank of capacitors across a centered cathode and an annular anode in tandem electrode sections. The initial shock wave passes through the anode section from the cathode section thereby depositing energy into the flow gas without the necessity of any diaphragm opening in the gas flow from the expansion section through the electrode sections.

Leibowitz, L. P.↗
NASA NTRS1986

Shock-induced volatile loss from a carbonaceous chondrite: Implications for planetary accretion

Solid recovery impact-induced volatile loss experiments on the Murchison C2M meteorite indicate that for an impact of a given velocity, H2O and total volatiles are driven from the sample in the same proportion as present initially. The primitive surface volatile budget of a planet growing by accretion would have the same bulk elemental composition as the volatiles in the incident planetesimals. Incipient devolatilization of Murchison occurs at an initial shock pressure of about 11 GPa and complete devolatilization occurs at a pressure of about 30 GPa. For the Earth, incipient and complete devolatilization of accreting planetesimals would occur when the planet reached approximately 12% and 27%, respectively, of its present-day radius. Impact-induced devolatilization would profoundly affect the volatile distribution within the accreting planet. Prior to metallic core formation and internal differentiation the growing planet would have a very small core with the same volatile content as the incident material, a volatile depleted mantle, and an extremely volatile rich surface. In the case of the Earth, 99.4 wt% of the total incident volatile material would end up on or near the planetary surface.

Tyburczy, J. A.↗
NASA NTRS1986

Shock-induced volatile loss from a carbonaceous chondrite Implications for planetary accretion

Solid recovery impact-induced volatile loss experiments on the Murchison C2M meteorite indicate that for an impact of a given velocity, H2O and total volatiles are driven from the sample in the same proportion as present initially. The primitive surface volatile budget of a planet growing by accretion would have the same bulk elemental composition as the volatiles in the incident planetesimals. Incipient devolatilization of Murchison occurs at an initial shock pressure of about 11 GPa and complete devolatilization occurs at a pressure of about 30 GPa. For the earth, incipient and complete devolatilization of accreting planetesimals would occur when the planet reached approximately 12 percent and 27 percent, respectively, of its present-day radius. Impact-induced devolatilization would profoundly affect the volatile distribution within the accreting planet. Prior to metallic core formation and internal differentiation the growing planet would have a very small core with the same volatile content as the incident material, a volatile depleted mantle, and an extremely volatile rich surface. In the case of the earth, 99.4 wt pct of the total incident volatile material would end up on or near the planetary surface.

Tyburczy, James A.↗
DOE OSTI2025

Sample glue layer investigation and mitigation for laser induced prompt impulse experiments

Understanding longer timescale material reactions under dynamic stress loading is critical for applications in materials engineering, shock physics, and planetary science. Prompt impulse experiments generate lower pressures since the ablator—the material first removed by the laser—is thicker and farther from the diagnostic plane, capturing aggregate material responses from the initial shock wave, rarefaction waves, and later time effects. This complexity demands thorough material characterization and simulation support. Since traditional sample construction is specific to supported shock experiments, designing prompt impulse experiments requires reconsideration around target design and sample engineering. Here, we present sample preparation techniques, experimental investigations, and theoretical simulations to investigate glue layer impacts, aiming to standardize samples for consistent data at lower laser fluences. We find that glue layers <30 μm have a minimal impact on peak velocity and pulse shape. The peak velocity scales linearly with glue layer thickness until a glue layer of 75 μm. For glue layers >75 μm, the peak velocity no longer scales with thickness; however, the pulse shape continues to degrade as described by simulations.

Lasers↗
NASA NTRS1987

The effect of neutrino transport on the collapse of iron stellar cores

A multigroup flux-limited diffusion approximation to neutrino transport, correct to first order in material velocities, is described. The role of neutrino-electron scattering in determining core deleptonization is studied using a modified Fokker-Planck approximation in which adjacent neutrino energy groups are coupled. The moderate deleptonizaton obtained during infall restricts homologous core masses to less than 0.8 solar mass. With 0.6-0.8 solar mass of outer core to traverse, nuclear dissociation alone stalls the shock waves produced at core bounce. In addition, neutrino energy losses remove at least another 4 x 10 to the 51st ergs. In no case is matter ejected promptly as a direct result of the initial shock wave.

Myra, E. S.↗
NASA NTRS1989

Cyclic behavior at quasi-parallel collisionless shocks

Large scale one-dimensional hybrid simulations with resistive electrons have been carried out of a quasi-parallel high-Mach-number collisionless shock. The shock initially appears stable, but then exhibits cyclic behavior. For the magnetic field, the cycle consists of a period when the transition from upstream to downstream is steep and well defined, followed by a period when the shock transition is extended and perturbed. This cyclic shock solution results from upstream perturbations caused by backstreaming gyrating ions convecting into the shock. The cyclic reformation of a sharp shock transition can allow ions, at one time upstream because of reflection or leakage, to contribute to the shock thermalization.

Burgess, D.↗
DOE OSTI2020

Generation of shock trains in free liquid jets with a nanosecond green laser

Shock wave trains in liquid jets were previously generated only by ablation with femtosecond x-ray lasers. Here we show that shock trains in water microjets can be also generated using nanosecond green laser pulses with 1- to 10-mJ energy. Furthermore, the ablation of 15-, 20-, 30-, and 70-μm water microjets opened a gap in the jets and launched an initial shock wave. Fully developed shock trains were observed in the 30- and 70-μm jets up to 250-ns delays, and these trains were also transmitted inside the nozzles. A few tens of nanoseconds after the pulse, the shock dynamics and its pressure became similar to the ones generated by x-ray lasers, with a more rapid pressure decay in thinner jets. At time delays exceeding 100 ns in the 30-μm jets, the leading shock pressure stabilized to an approximately constant pressure of 40 MPa. The energy density deposited in the jets was estimated at 30 MJ/cm 3 by comparing the jet gaps in the green and x-ray laser experiments, and matched previous estimates for optical ablation in water. Here, the pressure decay in the 30-μm jets was modeled based on the pressure decay observed in x-ray laser experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

High Explosive Ignition through Chemically Activated Nanoscale Shear Bands

Shock initiation and detonation of high explosives is considered to be controlled through hot spots, which are local regions of elevated temperature that accelerate chemical reactions. Using classical molecular dynamics, we predict the formation of nanoscale shear bands through plastic failure in shocked TATB high explosive crystal. By scale-bridging with quantum-based molecular dynamics, we show that shear bands exhibit lower reaction barriers. While shear bands quickly cool,they remain chemically activated and support increased reaction rates without the local heating typically evoked by the hot spot paradigm. We describe this phenomenon as chemical activation through shear banding

Chemistry - Inorganic↗
DOE OSTI2023

Compressible pairwise interaction extended point-particle model for force prediction of shock-particle bed interaction

We propose a pairwise influence framework for the complex unsteady compressible particle-laden flow problem by accounting for the scattered hydrodynamic waves emitting from neighboring particles in a Euler-Lagrange simulation. It has been observed from particle-resolved (PR) simulations of randomly dispersed particle beds under a loading shock that the compressible pseudoturbulence dominates the flow system even after the primary shock has passed, which causes fluctuations observed in the forces experienced by the particle. Moreover, the fact that each particle exists in the vicinity of a random arrangement of other particles modifies the time history of the drag force experienced by each particle during and after the passage of the shock. First, the scattering flow field due to an incoming shock interacting with a single sphere is constructed using an analysis of the flow in the acoustic limit. Then we examine the validity of the compressible Maxey-Riley-Gatignol force model by comparing the force prediction against a PR simulation of two interacting particles for various particle arrangements and incoming shock strength. Subsequently, the neighboring influences are stored as a library of maps that can be used readily in the calculation of the perturbation force. Lastly, the pairwise interaction assumption is evaluated by comparing the force predicted with the model with PR simulations of a randomly packed particle bed of 10% volume fraction for both water and air as the fluid medium for an incoming shock Mach number 1.22. With a considerably lower cost for the implementation of the model compared to PR simulations, it is verified that the model is reasonably accurate in pinpointing particles whose peak force is significantly larger or smaller than the mean drag but also to capture the prolonged fluctuations after the initial shock.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Femtosecond diffraction studies of the sodium chloride phase diagram under laser shock compression

The phase diagram of sodium chloride (NaCl) under laser shock compression has been studied at Linac Coherent Light Source (LCLS) at the x-ray free-electron laser facility. Both solid–solid (B1 → B2) and solid–liquid (B2 → liquid) transitions have been observed along the Hugoniot over nanosecond time scales. By combining structural measurements through in situ x-ray diffraction, pressure determination through velocimetry, and a thermal equation-of-state, the shock-compressed data are used to constrain the phase diagram of NaCl. Transformation into the B2 phase is found to occur at 28(2) GPa, and B2–liquid coexistence is observed between 54(4) and 66(6) GPa, with near full melt at 66(6) GPa. Late-time pressure release from an initial shocked B2-state results in a B2 → B1 back transformation. Our results show agreement with previous static compression data, suggesting that the time scale for melting is very rapid and that equilibrium states in NaCl are being accessed over nanosecond time scales. Finally, a multiphase equation-of-state description of NaCl incorporated into a one-dimensional hydrocode is used to interpret pressure and temperature evolution over these rapid time scales.

36 MATERIALS SCIENCE↗
NASA NTRS1994

On the hydrodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds

The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.

Klein, Richard I.↗
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 OSTI2022

Observations on the deformation of metal microspheres in shock-driven polymer flows

We report that solid particles can be fragmented by a fast-moving fluid if their velocity difference is great enough, such as during the atmospheric entry of meteoroids or the shock compression of engineered particulate composites. The extent of particle deformation and breakup in such systems is poorly understood because the necessary extreme conditions make observation difficult and data scarce. To meet this need, experiments combining ultrafast synchrotron-based radiography with plate impact loading were performed at the dynamic compression sector at the advanced photon source. Metal microspheres of several densities and strengths (Au, Ta, and W) were placed inside a polymer matrix. A planar shock wave was then produced in the polymer by the impact of a gun-launched flyer plate. X-ray images of the resulting flow were collected at ~150ns intervals. These images document the progression of particle deformation across a range of flow conditions and particle materials. They show that the extent of deformation is sensitive to the ratio of drag stress to particle strength. The deforming particle's shape is determined by the initial shock–particle interaction, fluid stagnation pressure, and vorticity, each acting on its own timescale. A set of scaling relationships is presented to capture these observations and enable comparison with prior hydrodynamic data. The result is a framework for predicting the conditions under which strong particles are severely deformed by a shock-driven flow.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Enhancements to viscous-shock-layer technique

A solution procedure is presented which considerably improves the computational efficiency of the viscous-shock-layer technique, especially for long slender bodies. The 'predictor-corrector' procedure suggested for obtaining the shock shape beyond the nose region requires only a single global pass. The accuracy of the present method is demonstrated by comparison with globally iterated results over the entire body and with ground- and flight-test data. A good comparison of the results computed with the two methods is shown for different flowfield chemistry models and axisymmetric body shapes. The new procedure results in computer run times 1/3 to 1/2 of the times required for the full-body global iteration procedure. Further, the algebraic expressions used to specify the initial shock shape eliminate the need for a shock shape generated by external means and permit immediate introduction of the full viscous-shock-layer equations. Finally, the present method of solution for the VSL equations provides to the aerothermal designer a very efficient and accurate tool for detailed flowfield as well as future technology studies.

Gupta, Roop N.↗
DOE OSTI2020

ZND Verification Tests for Reactive Burn Models in FLAG

The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Experimental Pathways for Detecting Double Superionicity in Planetary Ices

The ice giant planets Uranus and Neptune are assumed to contain large amounts of planetary ices such as water, methane, and ammonia. The properties of mixtures of such ices at the extreme pressures and temperatures of planetary interiors are not yet well understood. Ab initio computer simulations have predicted that a number of ices exhibit a hydrogen superionic state and a doubly superionic state. Since the latter state has not yet been generated with experiments, we outline here two possible pathways for reaching and detecting such a state with dynamic compression experiments. Here, we suggest X-ray diffraction as the principal tool for detecting when the material becomes doubly superionic and the sublattice of one of the heavy nuclei melts. That would require a temperature of ~3500 K and pressures greater than ~200 GPa for H 3 NO 4 , which we use as an example material here. Such conditions can be reached with experiments that employ an initial shock that is followed by a ramp compression wave. Alternatively, one may use triple-shock compression because a single shock does not yield sufficiently high densities.

Neptune↗
DOE OSTI2022

A Bayesian Deep Learning Approach to Near-Term Climate Prediction

Since model bias and associated initialization shock are serious shortcomings that reduce prediction skills in state-of-the-art decadal climate prediction efforts, we pursue a complementary machine-learning-based approach to climate prediction. The example problem setting we consider consists of predicting natural variability of the North Atlantic sea surface temperature on the interannual timescale in the pre-industrial control simulation of the Community Earth System Model. While previous works have considered the use of recurrent networks such as convolutional LSTMs and reservoir computing networks in this and other similar problem settings, we currently focus on the use of feedforward convolutional networks. In particular, we find that a feedforward convolutional network with a Densenet architecture is able to outperform a convolutional LSTM in terms of predictive skill. Next, we go on to consider a probabilistic formulation of the same network based on Stein variational gradient descent and find that in addition to providing useful measures of predictive uncertainty, the probabilistic (Bayesian) version improves on its deterministic counterpart in terms of predictive skill. Finally, we characterize the reliability of the ensemble of machine learning models obtained in the probabilistic setting by using analysis tools developed in the context of ensemble numerical weather prediction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Ultrafast spectroscopic studies of vibrational energy transfer in energetic materials

Shock-induced detonation is a key property of energetic materials (EM) that remains empirically understood. One proposed mechanism of shock-initiation in EM is “phonon up-pumping” to initiate chemical reactions, where excitation of lattice phonon modes rapidly transfers energy into intramolecular vibrations, ultimately resulting in the breaking of chemical bonds. We are developing novel ultrafast laser spectroscopy techniques to study vibrational energy transfer from phonon modes to intramolecular vibrations (phonon up-pumping), as well as competing energy transfer pathways from intramolecular vibrations to phonon modes (vibrational cooling). Through combinations of plasma- generated supercontinuum infrared, tunable near- and mid-infrared, and terahertz pulses in pump-probe spectroscopy, supplemented with ab inito simulations, we can explore the energy transfer processes on a sub-picosecond time scale to elucidate vibrational energy transfer pathways and lifetimes in EM. In this work, we highlight recent progress, including the spectral and temporal characteristics of the infrared and THz sources as well as preliminary results on select EM.

36 MATERIALS SCIENCE↗