Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effect of ionic disorder on the principal shock Hugoniot

Here, the effect of ionic disorder on the principal Hugoniot is investigated using multiple scattering theory to very high pressure (Gbar). Calculations using molecular dynamics to simulate ionic disorder are compared to those with a fixed crystal lattice, for both carbon and aluminum. For the range of conditions considered here we find that ionic disorder has a relatively minor influence. It is most important at the onset of shell ionization and we find that, at higher pressures, the subtle effect of the ionic environment is overwhelmed by the larger number of ionized electrons with higher thermal energies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Neutrino halo effect on collective neutrino oscillation in iron core-collapse supernova model of a 9.6 M ⊙ star

We extend the multi-angle computational framework and investigate the time evolution of the neutrino halo on collective neutrino oscillation in the core collapse of an iron core progenitor. We find that in the case of the 9.6 M progenitor adopted in this work, there are windows of time when the effects of neutrino halo and collective neutrino oscillation are not simultaneously large. Inside the shock, the impact of the inward-scattered halo neutrino cannot in general be neglected compared to the outward-propagating neutrino flux. However, during early epochs, collective neutrino oscillation is effectively shut down by multi-angle matter suppression. During the intermediate epoch, collective neutrino oscillation is not suppressed, but its onset radius is beyond the still relatively small explosion shock front where the halo is prominent. We also find in the case of the 9.6 M progenitor the halo neutrinos induce a delay in the onset of collective neutrino oscillations. This causes novel flavor conversions which sharpen collective neutrino oscillation spectral features. We predict that the inclusion of neutrino halo effects makes neutrino signals that are more clearly distinct from thermal emission that when halo neutrinos are omitted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Dynamics of millicharged dark matter in supernova remnants

Dark matter could have a small electromagnetic charge, provided the charge-to-mass ratio is much less than that of electrons or protons. This candidate, commonly known as millicharged dark matter (mDM), would form a plasma and interact with the interstellar medium and electromagnetic fields within galaxies. In general, understanding the dynamics of mDM requires consideration of collective plasma effects. It has been proposed that mDM can be accelerated in supernova remnants, forming a dark cosmic ray population that would leave distinct experimental signatures. In this work, we study a microphysical model where mDM is shocked by a supernova remnant and isotropized in the frame of the expanding fluid. We find that for | q χ / m χ | ≳ 10 – 13 e / MeV , the isotropization length for electromagnetic plasma instabilities is much shorter than the size of the supernova remnant. This is a necessary, though not sufficient, first step for formation of a Fermi-accelerated mDM component, and determining the size of this component requires further study. Finally, we discuss additional implications of mDM interactions in supernova remnants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Equation of State of CO 2 Shock Compressed to 1 TPa

Equation-of-state (pressure, density, temperature, internal energy) and reflectivity measurements on shock-compressed CO 2 at and above the insulating-to-conducting transition reveal new insight into the chemistry of simple molecular systems in the warm-dense-matter regime. CO 2 samples were precompressed in diamond-anvil cells to tune the initial densities from 1.35 g/cm 3 (liquid) to 1.74 g/cm 3 (solid) at room temperature and were then shock compressed up to 1 TPa and 93 000 K. Lastly, variation in initial density was leveraged to infer thermodynamic derivatives including specific heat and Gruneisen coefficient, exposing a complex bonded and moderately ionized state at the most extreme conditions studied.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Physical Properties of the Star-forming Regions in the Interacting Galaxies NGC 3395/NGC 3396

We report long-slit spectroscopy of the early major merger galaxies NGC 3395/NGC 3396. The spectra are consistent with those for star-forming galaxies, but there is some indication of LINER-like active galactic nucleus activity in the center of NGC 3396. The total star formation rate in the regions observed is 2.83 M {sub ⊙} yr{sup −1}, consistent with estimates for the entire galaxies. The highest abundances are in the centers of the galaxies, with the abundances decreasing with distance. There is a correlation between high abundance and high ionization parameter, both of which can be attributed to the presence of massive stars. Modeling with SB 99 indicates the star-forming regions are younger than 10 Myr. There are 1000–2000 WNL stars in the system, along with several thousand O stars, consistent with the ages of the star-forming regions. The highest electron densities are found in young regions with high star formation rates. The electron temperatures are higher than results for non-interacting galaxies, which is probably due to shock waves produced by the galaxy–galaxy interaction, the outflow of gas from massive stars, and/or collisions between gas clouds in the galaxies. There is star formation in the bridge of material between the galaxies. These regions are among the youngest in the system and have low abundances, suggesting the gas was pulled from the outer parts of the galaxies. X-ray point sources, probably high-mass X-ray binaries, are associated with several star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Thermal decoupling of deuterium and tritium during the ICF shock-convergence phase

We report that a series of thin glass-shell shock-driven DT gas-filled capsule implosions was conducted at the OMEGA laser facility. These experiments generate conditions relevant to the central plasma during the shock-convergence phase of ablatively driven inertial confinement fusion (ICF) implosions. The spectral temperatures inferred from the DTn and DDn spectra are most consistent with a two-ion-temperature plasma, where the initial apparent temperature ratio, $\frac{T_T}{T_P}$, is 1.5. This is an experimental confirmation of the long-standing conjecture that plasma shocks couple energy directly proportional to the species mass in multi-ion plasmas. The apparent temperature ratio trend with equilibration time matches expected thermal equilibration described by hydrodynamic theory. This indicates that deuterium and tritium ions have different energy distributions for the time period surrounding shock convergence in ignition-relevant ICF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Species Separation and Hydrogen Streaming upon Shock Release from Polystyrene under Inertial Confinement Fusion Conditions

Shock release from inertial confinement fusion (ICF) shells poses a great challenge to single-fluid hydrodynamic equations, especially for describing materials composed of different ion species. This has been evidenced by a recent experiment [Haberberger et al., Phys. Rev. Lett. 123, 235001 (2019)], in which low-density plasmas ( 10 19 to 10 20 cm – 3 ) are measured to move far ahead of what radiation-hydrodynamic simulations predict. To understand such experimental observations, in this work we have performed large-scale nonequilibrium molecular-dynamics simulations of shock release in polystyrene (CH) at experimental conditions. These simulations revealed that upon shock releasing from the back surface of a CH foil, hydrogen can stream out of the bulk of the foil due to its mass being lighter than carbon. This released hydrogen, exhibiting a much broader velocity distribution than carbon, forms low-density plasmas moving in nearly constant velocities ahead of the in-flight shell, which is in quantitative agreement with the experimental measurements. Such kinetic effect of species separation is currently missing in single-fluid radiation-hydrodynamics codes for ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Solar Energetic Particle Acceleration at a Spherical Shock with the Shock Normal Angle $θ$ $B_n$ Evolving in Space and Time

We present a 2D kinematic model to study the acceleration of solar energetic particles (SEPs) at a shock driven by a coronal mass ejection. The shock is assumed to be spherical about an origin that is offset from the center of the Sun. This leads to a spatial and temporal evolution of the angle between the magnetic field and the shock-normal direction (θ Bn ) as it propagates through the Parker spiral magnetic field from the lower corona to 1 au. We find that the high-energy SEP intensity varies significantly along the shock front due to the evolution of θ Bn . Generally, the west flank of the shock preferentially accelerates particles to high energies compared to the east flank and shock nose. This can be understood in terms of the rate of acceleration, which is higher at the west flank. Double power-law energy spectra are reproduced in our model as a consequence of the local acceleration and transport effects. These results will help us to better understand the evolution of SEP acceleration and provide new insights into large SEP events observed by multiple spacecraft, especially those close to the Sun, such as Parker Solar Probe and Solar Orbiter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Non-thermal neutrinos created by shock acceleration in successful and failed core-collapse supernova

ABSTRACT We present a comprehensive study of neutrino shock acceleration in core-collapse supernova (CCSN). The leading players are heavy leptonic neutrinos, νμ and ντ; the former and latter potentially gain the energy up to ∼100 and ∼200 MeV, respectively, through the shock acceleration. Demonstrating the neutrino shock acceleration by Monte Carlo neutrino transport, we make a statement that it commonly occurs in the early post-bounce phase (≲50 ms after bounce) for all massive stellar collapse experiencing nuclear bounce and would reoccur in the late phase (≳100 ms) for failed CCSNe. This opens up a new possibility to detect high energy neutrinos by terrestrial detectors from Galactic CCSNe; hence, we estimate the event counts for Hyper(Super)-Kamiokande, DUNE, and JUNO. We find that the event count with the energy of ≳80 MeV is a few orders of magnitude higher than that of the thermal neutrinos regardless of the detectors, and muon production may also happen in these detectors by νμ with the energy of ≳100 MeV. The neutrino signals provide a precious information on deciphering the inner dynamics of CCSN and placing a constraint on the physics of neutrino oscillation; indeed, the detection of the high energy neutrinos through charged current reaction channels will be a smoking gun evidence of neutrino flavour conversion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Impacts of Jets and winds from primordial black holes

ABSTRACT Primordial black holes (PBHs) formed in the early Universe constitute an attractive candidate for dark matter. Within the gaseous environment of the interstellar medium, PBHs with accretion discs naturally launch outflows such as winds and jets. We discuss for the first time how PBHs with significant spin can sustain powerful relativistic jets and generate associated cocoons. Jets and winds can efficiently deposit their kinetic energies and heat the surrounding gas through shocks. Focusing on the Leo T dwarf galaxy, we demonstrate that these effects form novel tests and set new limits on PBHs over a significant ∼10−2 –106 M⊙ mass range, including the parameter space associated with gravitational wave observations by the LIGO and VIRGO Collaborations. Observing the morphology of emission will allow to distinguish between jet and wind contributions, and hence establishes a new method for identifying spinning PBHs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Impact of late-time neutrino emission on the diffuse supernova neutrino background

In the absence of high-statistics supernova neutrino measurements, estimates of the diffuse supernova neutrino background (DSNB) hinge on the precision of simulations of core-collapse supernovae. Understanding the cooling phase of protoneutron star (PNS) evolution (≳ 1 s after core bounce) is crucial, since approximately 50% of the energy liberated by neutrinos is emitted during the cooling phase. We model the cooling phase with a hybrid method by combining the neutrino emission predicted by 3D hydrodynamic simulations with several cooling-phase estimates, including a novel two-parameter correlation depending on the final baryonic PNS mass and the time of shock revival. We find that the predicted DSNB event rate at Super-Kamiokande can vary by a factor of ~2–3 depending on the cooling-phase treatment. Here we also find that except for one cooling estimate, the range in predicted DSNB events is largely driven by the uncertainty in the neutrino mean energy. With a good understanding of the late-time neutrino emission, more precise DSNB estimates can be made for the next generation of DSNB searches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

$Z$-Pinch Interferometry Analysis With the Fourier-Based TNT Code

Here, we present the analysis of interferometry diagnostics with the user-friendly Talbot Numerical Tool (TNT), a Fourier-based postprocessing code that enables real-time assessment of plasma systems. TNT performance was explored with visible and infrared interferometry in pulsed-power-driven Z -pinch configurations to expand its capabilities beyond Talbot X-ray interferometry in the high-intensity laser environment. TNT enabled accurate electron density characterization of magnetically driven plasma flows and shocks through phase-retrieval methods that did not require data modification or masking. TNT demonstrated enhanced resolution, detecting below 4 % fringe shift, which corresponds to 8.7 × 10 15 cm –2 within 28 μ m, approaching the laser probing system limit. TNT was tested against a well-known interferometry analysis software, delivering an average resolving power nearly ten times better (~28 μ m versus ~ 210 μ m) when resolving plasma ablation features. TNT demonstrated higher sensitivity when probing sharp electron density gradients in supersonic shocks. A maximum electron areal density of 4.1 × 10 17 cm –2 was measured in the shocked plasma region, and a minimum electron density detection of ~ 1.0 × 10 15 cm –2 was achieved. When probing colliding plasma flows, the calculations of the effective adiabatic index and the associated errors were improved from γ* = 2.6 ± 1.6 – 1.4 ± 0.2 with TNT postprocessing, contributing valuable data for the interpretation of radiative transport. Additional applications of TNT in the characterization of pulsed-power plasmas and beyond are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dehalogenation reactions between halide salts and phosphate compounds

Reactions between phosphoric acid [H 3 PO 4 ] or ammonium hydrogen phosphates [i.e., NH 4 H 2 PO 4 , (NH 4 ) 2 HPO 4 ] and halide salts can be used to dehalogenate (remove halides from) salt-based waste streams, where the process of removing halides yields products that have more efficient disposal pathways for repository storage. In this context, the term efficiency is defined as higher waste loadings and simplified immobilization processes with potential for recycle of certain salt components (e.g., 37 Cl as H 37 Cl or NH 4 37 Cl). The main streams identified for these processes are nuclear wastes generated during electrochemical reprocessing of used nuclear fuel as well as used halide salts from molten salt reactor operation. The potential byproducts of these reactions are fairly consistent across the range of halide species (i.e., F, Cl, Br, I) where the most common are hydrogen halides [e.g., HCl (g) ] or ammonium halides (e.g., NH 4 Cl). However, trihalide compounds (e.g., NCl 3 ), nitrogen triiodide ammine adducts [NI 3 ·(NH 3 ) x ], and ammonium triiodide (NH 4 I 3 ) are also possible. Several of these byproducts (i.e., NCl 3 , NBr 3 , NI 3 , and NH 4 I 3 ) are shock-sensitive contact explosives so their production in these processes must be tracked and carefully controlled, which includes methods of immediate neutralization upon production such as direct transport to a caustic scrubber for dissolution. Several benefits arise from utilizing H 3 PO 4 as the phosphate additive during dehalogenation reactions for making iron phosphate waste forms including more oxidized iron (higher Fe 3+ :Fe 2+ ratios), higher chemical durabilities, and the avoidance of trihalides, but the byproducts are hydrogen halides, which are corrosive and require special handling.

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

Standing shock prevents propagation of sparks in supersonic explosive flows

Abstract Volcanic jet flows in explosive eruptions emit radio frequency signatures, indicative of their fluid dynamic and electrostatic conditions. The emissions originate from sparks supported by an electric field built up by the ejected charged volcanic particles. When shock-defined, low-pressure regions confine the sparks, the signatures may be limited to high-frequency content corresponding to the early components of the avalanche-streamer-leader hierarchy. Here, we image sparks and a standing shock together in a transient supersonic jet of micro-diamonds entrained in argon. Fluid dynamic and kinetic simulations of the experiment demonstrate that the observed sparks originate upstream of the standing shock. The sparks are initiated in the rarefaction region, and cut off at the shock, which would limit their radio frequency emissions to a tell-tale high-frequency regime. We show that sparks transmit an impression of the explosive flow, and open the way for novel instrumentation to diagnose currently inaccessible explosive phenomena.

58 GEOSCIENCES↗
DOE OSTI2020

Implementation of Talbot–Lau x-ray deflectometry in the pulsed power environment using a copper X-pinch backlighter

A Talbot–Lau x-ray deflectometer can map electron density gradients in high energy density plasmas, as well as provide information about plasma elemental composition through single-image x-ray refraction and attenuation measurements. A new adaptation to a pulsed power environment used backlighting from copper X-pinches, allowing for electron density mapping of a low-Z object. Even though the X-pinch backlighter is not properly optimized for emitting x-rays in terms of source size and photon fluence, Moire fringe patterns with contrast up to 14% and fringe shift due to refraction on a beryllium object are obtained. Due to the proximity of the deflectometer with the X-pinch (~6 cm), it is shown that a protective filter is required to avoid damage in the closest (i.e., source) grating due to both plasma debris and mechanical shock. Regarding grating survival, these did not show any damage due to the intense magnetic field or heating induced by plasma radiation. Electron density on beryllium was measured with a difference lower than 16%. The areal electron density mapping of the sample was limited by source size characteristics, in similarity to transmission radiography. These results show the potential of plasma electron density as well as material mapping through Talbot–Lau x-ray deflectometry in a pulsed power environment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Quantum Monte Carlo determination of the principal Hugoniot of deuterium

Here, we present coupled electron-ion Monte Carlo results for the principal Hugoniot of deuterium together with an accurate study of the initial reference state of shock-wave experiments. We discuss the influence of nuclear quantum effects, thermal electronic excitations, and the convergence of the potential energy surface by wave-function optimization within variational Monte Carlo and projection quantum Monte Carlo methods. Compared to a previous study, our calculations also include low pressure-temperature (P,T) conditions resulting in close agreement with experimental data, while our revised results at higher (P,T) conditions still predict a more compressible Hugoniot than experimentally observed.

36 MATERIALS SCIENCE↗