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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Magnetohydrodynamic Effects in Propagating Relativistic Ejecta: Reverse Shock and Magnetic Acceleration

We solve the Riemann problem for the deceleration of arbitrarily magnetized relativistic ejecta injected into a static unmagnetized medium. We find that for the same initial Lorentz factor, the reverse shock becomes progressively weaker with increasing magnetization s (the Poynting-to-kinetic energy flux ratio), and the shock becomes a rarefaction wave when s exceeds a critical value, sc, defined by the balance between the magnetic pressure in the ejecta and the thermal pressure in the forward shock. In the rarefaction wave regime, we find that the rarefied region is accelerated to a Lorentz factor that is significantly larger than the initial value. This acceleration mechanism is due to the strong magnetic pressure in the ejecta.

Mizuno, Y.↗
NASA NTRS2012

Failure Mode Analysis of V-Shaped Pyrotechnically Actuated Valves

Current V-shaped stainless steel pyrovalve initiators have rectified many of the deficiencies of the heritage Y-shaped aluminum design. However, a credible failure mode still exists for dual simultaneous initiator (NSI) firings in which low temperatures were detected at the booster cap and less consistent ignition was observed than when a single initiator was fired. In order to asses this issue, a numerical framework has been developed for predicting the flow through pyrotechnically actuated valves. This framework includes a fully coupled solution of the gas-phase equation with a non-equilibrium dispersed phase for solid particles as well as the capability to model conjugate gradient heat transfer to the booster cap. Through a hierarchy of increasingly complex simulations, a hypothesis for the failure mode of the nearly simultaneous dual NSI firings has been proven. The simulations indicate that the failure mode for simultaneous dual NSI firings may be caused by flow interactions between the flame channels. The shock waves from each initiator interact in the booster cavity resulting in a high pressure that prevents the gas and particulate velocity from rising in the booster cap region. This impedes the bulk of the particulate phase from impacting the booster cap and reduces the heat transfer to the booster cap since the particles do not impact it. Heat transfer calculations to the solid metal indicate that gas-phase convective heat transfer may not be adequate by itself and that energy transfer from the particulate phase may be crucial for the booster cap burn through.

Sachdev, Jai S.↗
DOE OSTI2026

Image Distinguishability Analysis Testing Through Principal Components and Its Application to Hot Spot Scale Invariance

Hot spots are spatial regions of intense energy localization that govern initiation of secondary high explosives. Studies that characterize or compare simulated hot spots are frequently either qualitatively descriptive or resort to quantitative distribution functions that neglect stochastic variations and spatial correlations—effects that are also neglected in common comparison tests like the Kolmogorov–Smirnov test. To this end, we develop an image distinguishability analysis (IDA) test based on principal component (PC) analysis that makes pixel-by-pixel comparisons between small, for example, O(<10), image data sets. The IDA test makes comparisons through a generalized distance metric in the PC space and a test statistic that is derived to calculate mathematical equation-values. Here, we derive a statistical distribution and criticality criterion to determine whether images are distinguishable from established baselines. We apply the IDA test on images generated from molecular dynamics simulations of hot spots from pore collapse in TATB to assess scale invariance in the complex patterns of hot spots that form in a representative high explosive crystal. The IDA test shows that TATB hot spot spatial temperature fields and their derived temperature histograms exhibit scale-invariant features over specific intervals of shock orientation, strength, and initial pore diameter. However, the IDA test also shows that qualitatively different conclusions regarding invariance can be reached depending on whether the hot spot is treated as a spatially correlated field as opposed to a distribution function that lacks spatial information.

organic↗
NASA NTRS2002

On the Exit Boundary Condition for One-Dimensional Calculations of Pulsed Detonation Engine Performance

In one-dimensional calculations of pulsed detonation engine (PDE) performance, the exit boundary condition is frequently taken to be a constant static pressure. In reality, for an isolated detonation tube, after the detonation wave arrives at the exit plane, there will be a region of high pressure, which will gradually return to ambient pressure as an almost spherical shock wave expands away from the exit, and weakens. Initially, the flow is supersonic, unaffected by external pressure, but later becomes subsonic. Previous authors have accounted for this situation either by assuming the subsonic pressure decay to be a relaxation phenomenon, or by running a two-dimensional calculation first, including a domain external to the detonation tube, and using the resulting exit pressure temporal distribution as the boundary condition for one-dimensional calculations. These calculations show that the increased pressure does affect the PDE performance. In the present work, a simple model of the exit process is used to estimate the pressure decay time. The planar shock wave emerging from the tube is assumed to transform into a spherical shock wave. The initial strength of the spherical shock wave is determined from comparison with experimental results. Its subsequent propagation, and resulting pressure at the tube exit, is given by a numerical blast wave calculation. The model agrees reasonably well with other, limited, results. Finally, the model was used as the exit boundary condition for a one-dimensional calculation of PDE performance to obtain the thrust wall pressure for a hydrogen-air detonation in tubes of length to diameter ratio (L/D) of 4, and 10, as well as for the original, constant pressure boundary condition. The modified boundary condition had no performance impact for values of L/D > 10, and moderate impact for L/D = 4.

Wilson, Jack↗
DOE OSTI2026

Reduced‐Order Modeling of Energetic Materials Using Physics‐Aware Recurrent Convolutional Neural Networks in a Latent Space (LatentPARC)

Physics-aware deep learning (PADL) has gained popularity for use in spatiotemporal dynamics simulations, such as those in computational modeling of energetic materials (EM). We show that the challenge PADL methods face while learning complex field evolution problems can be simplified and accelerated by decoupling it into two tasks: learning complex geometric features in evolving fields and modeling dynamics over these features in a lower-dimensional feature space. We build upon our previous work on physics-aware recurrent convolutional neural networks (PARC). PARC embeds knowledge of underlying physics into its neural network architecture for more robust and accurate prediction of evolving physical fields. PARC was shown to effectively learn complex nonlinear features such as the formation of hotspots and coupled shock fronts in various initiation scenarios of EMs, as a function of microstructures, serving effectively as a microstructure-aware burn model. Here, we further accelerate PARC and reduce its computational cost by projecting the original dynamics onto a lower-dimensional invariant manifold, or “latent space.” The projected latent representation encodes the complex geometry of evolving fields (e.g., temperature and pressure) in a set of data-driven features. The reduced dimension of this latent space allows us to learn the dynamics during the initiation of EM with a lighter and more efficient model. We observe a significant decrease in training and inference time while maintaining results comparable to PARC at inference. This work takes steps towards enabling rapid prediction of EM thermomechanics at larger scales and characterization of EM structure–property–performance linkages at a full application scale.

Mathematics and Computing↗
DOE OSTI2023

VIPIR: A High-Throughput Drop-Weight Impact Instrument for Imaging the Initiation and Propagation of Reactions in Energetic Materials

Characterizing the handling safety and sensitivity of explosives has been a challenging area of study for over 60 years. Historically one of the most accessible and widely utilized experiments has been the drop-weight impact test, which involves dropping a weight on a small sample sandwiched between two anvils. Because this experiment generally only utilizes sound thresholds to determine whether or not a sample reacted, the physical and chemical properties governing sensitivity remain convolved. Better understanding of chemical and material characteristics is needed to give the chemistry and engineering communities a predictive tool to determine the handling sensitivity of explosives prior to pursuing expensive and potentially hazardous synthesis and formulation operations. Here, we are developing a high throughput drop tower instrument capable of imaging the deformation and flow of energetic materials during impact and the resulting thermal ignition and propagation events. This instrument is based on previous designs but has been modified for higher throughput and tailorable modifications in the future. Herein, we present key design features that enable high-speed visible and thermal imaging of explosive initiation by sub-shock impacts, as well as preliminary results in which ignition sites were observed in an erythritol tetranitrate sample.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Invited article: X-ray phase contrast imaging in inertial confinement fusion and high energy density research

X-ray phase contrast imaging (XPCI) provides enhanced image contrast beyond absorption-based x-ray imaging alone due to refraction and diffraction from gradients in the object material density. It is sensitive to small variations in density, such as internal voids, cracks, grains, defects, and material flow, as well as to stronger density variations such as from a shock wave. Beyond its initial use in biology and materials science, XPCI is now routinely used in inertial confinement fusion (ICF) and high energy density (HED) research, first to characterize ICF capsules and targets, and later applied in dynamic experiments, where coherent x-ray sources, ultrafast x-ray pulses, and high temporal and spatial resolution are required. In this Review article, XPCI image formation theory is presented, its diverse use in ICF and HED research is discussed, the unique requirements for ultrafast XPCI imaging are given, as well as current challenges and issues in its use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

AT2025ulz and S250818k: Zooming in with the Hubble Space Telescope

AT2025ulz is an optical/near-infrared transient discovered during follow-up of the candidate gravitational wave (GW) event S250818k. Its young age (≲1 day), rapid decline, and strong color evolution over the first 48 hr classify it as a potential kilonova candidate. In this work, we present the results of our observing campaign, carried out with the Gran Telescopio Canarias and the Hubble Space Telescope (HST). Although the early time evolution of AT2025ulz resembles some aspects of a kilonova, its rapid onset (∼3 hr after the GW trigger) and luminosity (a factor of ∼5 brighter than AT2017gfo in g band) are difficult to reproduce. Only a small subset of our kilonova models matches its multicolor light curve, and the inferred ejecta mass is uncomfortably large given the low chirp mass (≲0.87 M ⊙ ) of the GW candidate. HST observations place the transient within a nearby (z = 0.08489) spiral galaxy with ongoing star formation and measure a color (F336W − F160W ≈ 1.4 mag) that is too blue to match with a kilonova. Our data support the classification of AT2025ulz as a supernova, initially undergoing a shock-cooling phase and later entering its photospheric phase and spectroscopically identified via its broad absorption features.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

A gas-dynamic calculation of type 2 shock propagation through the corona

An approximate analytic theory of acoustic shock propagation in nonuniform media is used to determine the motion of a flare-generated shock wave in the corona. The shock is followed from the time it strikes the chromosphere-corona transition region (density interface) out to 5 solar radii under the assumption that the corona in this region is approximately in hydrostatic equilibrium. The strength of the shock incident on the transition region from below determines the ejection velocity of eruptive prominence material, as well as the initial velocity of the coronal shock. The calculation is applied to one well-documented case of a related flare spray, moving type 4 isolated source, and type 2 burst. It is shown that a chromospheric shock of the appropriate strength to produce the observed prominence and type 4 velocities strengthens as it moves out in the corona by an amount sufficient to account for the observed high velocity of the type 2 burst.

Kopp, R. A.↗
NASA NTRS1981

The Cygnus 'superbubble' - A supernova explosion in a tenuous intercloud medium

Cash et al. (1980) observed an incomplete X-ray shell in the Cygnus region with a temperature of 2 x 10 to the 6th K and a density of 0.02 per cu cm. This feature closely resembles Loop I, a nearby supernova remnant expanding in a tenuous (less than 0.01 per cu cm) medium. A similar origin - remnant expansion in a tenuous intercloud medium with the bulk of the interstellar gas concentrated in clouds - is suggested for the Cygnus X-ray feature. This model greatly reduces the required initial thermal energy. Remnant shock interaction with ambient interstellar clouds produces the observed H-alpha filaments.

Higdon, J. C.↗
NASA NTRS1985

Active Pulse Suppression of Structural Vibrations

Corrective opposing vibrations reduce Earthquake damage. Algorithm for controlling reactions of buildings and other structures to Earthquake initiates counterpulses when incoming shock pulses exceed threshold level. By using "go/no-go" criterion rather than continuously adjusting response to every input excitation no matter how small, control system needs less computation effort than comparable closed-loop systems.

Udwadia, F.↗
NASA NTRS1985

Evidence for Non-radiative Activity in Hot Stars

Indications of non-radiative activity in hot stars were investigated. Data used as evidence of non-radiative activity includes: (1) evidence for high temperature regions, especially regions that are too hot to be explained by shocks that can be initiated and driven by radiative forces; (2) evidence for multi-plume flow. Non-spherically symmetric flow can be diagnosed by the effect it can have on the polarization of the light from a star; and (3) direct evidence for stellar magnetic fields.

Cassinelli, J. P.↗
NASA NTRS1984

Magnetic field amplification and generation in hypervelocity meteoroid impacts with application to lunar paleomagnetism

A one-dimensional numerical model for the expansion of impact-produced vapor clouds is used to investigate magnetic field generation mechanisms in events such as meteor collisions with the moon. The resulting cloud properties, such as ionization fraction, electrical conductivity, radial expansion velocity, mass density, and energy density are estimated. The model is initiated with the peak shock states and pressure thresholds for incipient and complete vaporization of anorthosite lunar surface materials by iron and GA composition meteorites. The expansion of the spherical gas cloud into a vacuum was traced with a one-dimensional explicit lagrangian hydrodynamic code. The hypervelocity impact plasmas produced are found to be significant in the amplitudes and orientations of the magnetic fields generated. An ambient magnetic field could have been provided by the core dynamo, which would have interacted with the expanding plasmas and formed induced paleomagnetic fields. Several other field-contribution mechanisms are discussed and discarded as potential remanent magnetism contributors.

Hood, L. L.↗
NASA NTRS1985

Electrical safety requirements: Implications for the module designer

Commercial photovoltaic array installations, which include residential and intermediate applications, are subject to building and electrical codes and to product safety standards. The National Electrical Code (NEC) Article 690, titled Solar Photovoltaic Systems, contains provisions defining acceptable levels of system safety and emphasizes the system design and its installation. The Underwriters Laboratories, Inc. (UL), document titled: Proposed First Edition of the Standard for Flat Plate Photovoltaic Modules and Panels, UL-1703, identifies module and panel construction requirements that ensure product safety. Together these documents describe requirements intended to minimize hazards such as shock and fire. Although initial focus of these requirements is on single crystal silicon modules, they are generic in nature, and are equally applicable to high voltage ( 30 Vdc), multikilowatt, thin film systems. A major safety concern is insulation breakdown within the module or array wiring system, or discontinuities within the electrical conductors. These failures can result in ground faults, in circuit arcs, or exposure to hazardous electrical parts. Safeguards are discussed.

Sugimura, R. S.↗
NASA NTRS1988

Supernovae and supernova remnants

The basic features and astrophysics of supernovae of types Ia, Ib, and II and of SNRs are discussed, summarizing the results of recent ground-based and space observations obtained at different wavelengths. Typical data are presented in graphs and sample images, and consideration is given to their implications for shock heating by the initial explosion, radioactive energy input, circumstellar interaction, and pulsar energy input. SNRs described include Tycho's remnant, Kepler's remnant, the Cyg loop, Pup A, the Crab Nebula, the Vela and Crab pulsars, and 0540-69.3 in the LMC. The need for high-spatial-resolution IR and X-ray spectra of these objects is indicated.

Chevalier, Roger A.↗
NASA NTRS1994

An approximate Riemann solver for magnetohydrodynamics (that works in more than one dimension)

An approximate Riemann solver is developed for the governing equations of ideal magnetohydrodynamics (MHD). The Riemann solver has an eight-wave structure, where seven of the waves are those used in previous work on upwind schemes for MHD, and the eighth wave is related to the divergence of the magnetic field. The structure of the eighth wave is not immediately obvious from the governing equations as they are usually written, but arises from a modification of the equations that is presented in this paper. The addition of the eighth wave allows multidimensional MHD problems to be solved without the use of staggered grids or a projection scheme, one or the other of which was necessary in previous work on upwind schemes for MHD. A test problem made up of a shock tube with rotated initial conditions is solved to show that the two-dimensional code yields answers consistent with the one-dimensional methods developed previously.

Powell, Kenneth G.↗
NASA NTRS1995

Transition of slow shocks to fast shocks

As a CME pushes its way through preceding slower solar wind, large disturbances in the interaction region may evolve to form transient MHD shocks. The shocks can be slow shocks in the coronal space but they appear as fast shocks near 1 AU. We use a polytropic MHD model to study the transition of slow shocks to fast shocks in an interaction region sandwiched between a faster solar wind and a slower solar wind near the equatorial plane. The polytropic index controls the radial increase of Beta outside the interaction region. The transition process is attributed chiefly to the increase of Beta and theta to a lesser degree to the radial increase of the shock angle 0. Under the initial condition of small Beta and theta near 0.1 AU, the interaction region evolves to form a pair of slow shocks inside 0.15 AU. As the interaction region convects outward. the increases of B and theta cause a transition of the shock system between 0.15 and 0.3 AU from a slow shock pair to a double shock pair consisting of both slow and fast shocks. As the system moves outward, and Beta continue to increase, the fast shock grows stronger, and the slow shock becomes weaker. Eventually. the slow shocks fade away, and the shock system finally evolves to a pair of fast shocks. Parametric studies of the transition process are carried out for shocks formed over a wide range of disturbances.

Whang, Y. C.↗