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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multidimensional Tests of a Finite-Volume Solver for MHD With a Real-Gas Equation of State

This work considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and at flux interfaces where the former is more computationally efficient. A battery of 1D and 2D tests are employed: convergence of 1D and 2D linearized waves, shock tube Riemann problems, a 2D nonlinear circularly polarized Alfvén wave, and a 2D magneto-Rayleigh-Taylor instability test. The cell-centered EOS evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation findings where the spurious solution leads to an unphysical simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Investigation of finite-volume methods to capture shocks and turbulence spectra in compressible flows

The aim of the present paper is to provide a comparison between several finite-volume methods of different numerical accuracy: the second-order Godunov method with PPM interpolation and the high-order finite-volume WENO method. In this work, the results show that while on a smooth problem the high-order method performs better than the second-order one, when the solution contains a shock all the methods collapse to first-order accuracy. In the context of the decay of compressible homogeneous isotropic turbulence with shocklets, the actual overall order of accuracy of the methods reduces to second-order, despite the use of fifth-order reconstruction schemes at cell interfaces. Most important, results in terms of turbulent spectra are similar regardless of the numerical methods employed, except that the PPM method fails to provide an accurate representation in the high-frequency range of the spectra. It is found that this specific issue comes from the slope-limiting procedure and a novel hybrid PPM/WENO method is developed that has the ability to capture the turbulent spectra with the accuracy of a high-order method, but at the cost of the second-order Godunov method. Overall, it is shown that virtually the same physical solution can be obtained much faster by refining a simulation with the second-order method and carefully chosen numerical procedures, rather than running a coarse high-order simulation. Our results demonstrate the importance of evaluating the accuracy of a numerical method in terms of its actual spectral dissipation and dispersion properties on mixed smooth/shock cases, rather than by the theoretical formal order of convergence rate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Non-contact velocity measurement instruments and systems, and related methods

Methods for calibration of non-contact velocity measurements and systems for implementing the same are described. Generally, the method comprises inducing a shock wave into a sample at a stress intensity that varies across the sample's elastic limit, which corresponds to the elastic-plastic state transition of the sample. That transition state may be at the sample's Hugoniot elastic limit. The velocity of the sample is measured using a non-contact velocity measurement instrument such as a velocimeter. The measurement may be compared to a predicted velocity or a velocity measurement made by another system to determine calibration parameters.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Dynamic Measurements of the Structural Evolution of Material Defects at the Mesoscale

This is a technical report to be published as part of the FY2021 annual overview package. We have designed experiments and developed analytical methods using images from x-ray and electron microscopy in order to quantify the behavior of crystal dislocation defects at their native nanometer-length and micro-length scales. The results of these analyses are used to directly inform theoretical physics models and guide interpretation of dark field x-ray microscopy (DFXM) images that offer spatial resolution surpassing the resolution of beamline optics. Additionally, we developed new analytical methods to automate the alignment of optics hardware used in both synchrotron and x-ray free electron laser beamline facilities. The tools and workflows developed for this project are now being deployed as proof-of-concept in reaction history analysis (RHA) for future programmatic integration.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗
DOE OSTI2020

Investigating off-Hugoniot states using multi-layer ring-up targets

Abstract Laser compression has long been used as a method to study solids at high pressure. This is commonly achieved by sandwiching a sample between two diamond anvils and using a ramped laser pulse to slowly compress the sample, while keeping it cool enough to stay below the melt curve. We demonstrate a different approach, using a multilayer ‘ring-up’ target whereby laser-ablation pressure compresses Pb up to 150 GPa while keeping it solid, over two times as high in pressure than where it would shock melt on the Hugoniot. We find that the efficiency of this approach compares favourably with the commonly used diamond sandwich technique and could be important for new facilities located at XFELs and synchrotrons which often have higher repetition rate, lower energy lasers which limits the achievable pressures that can be reached.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Tungsten doped diamond shells for record neutron yield inertial confinement fusion experiments at the National Ignition Facility

We report on fabrication and characterization of layered, tungsten doped, spherical about 2 mm diameter microcrystalline diamond ablator shells for inertial confinement fusion (ICF) experiments at the National Ignition Facility. As previously reported, diamond ICF ablator shells can be fabricated by chemical vapor deposition (CVD) on solid spherical silicon mandrels using an ellipsoidal microwave plasma reactor. In the present work, we further developed these ablator shells by embedding a W -doped diamond layer sandwiched between two undoped diamond regions. W incorporation in diamond was achieved by adding tungsten hexacarbonyl to the CH 4 /H 2 CVD feed gas. We observe that the W doping concentration decreases with increasing deposition rate which, in turn, is controlled by adjusting the total gas pressure. Cross sectional microstructural analysis reveals sharp interfaces between doped and undoped regions of the diamond shell and uniform W distribution with concentrations up to about 0.3 at.%. At higher W concentrations (>0.3 at.%) formation of tungsten carbide precipitates is observed. Using a 3-shock 1.6 MJ laser pulse, the targets described in this work produced the first laser driven implosion to break the 1 × 10 16 neutron yield barrier, followed by experiments (described in future publications) with similar targets and slightly more laser energy producing yields as high as 4 × 10 17 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Shock-compressed silicon: Hugoniot and sound speed up to 2100 GPa

In this work, the high-pressure equation of state and isentropic sound speed of fluid silicon were studied using laser-driven shock waves. Principal Hugoniot measurements to 2100 GPa were performed using impedance matching to an α-quartz reference. Sound speed was determined by time correlating imposed shock-velocity perturbations in both the sample (Si) and reference material (α-quartz). A change in shock velocity versus particle velocity (u s –u p ) slope on the fluid silicon principal Hugoniot was detected at 200 GPa. Density functional theory-based quantum molecular dynamics simulations suggest that an increase in ionic coordination and 50% increase in average ionization is coincident with the experimentally observed change in slope.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Neutrino oscillations in supernovae: Angular moments and fast instabilities

Recent theoretical work indicates that the neutrino radiation in core-collapse supernovae may be susceptible to flavor instabilities that set in far behind the shock, grow extremely rapidly, and have the potential to profoundly affect supernova dynamics and composition. Here we analyze the nonlinear collective oscillations that are prefigured by these instabilities. We demonstrate that a zero crossing in $n_{ν_e}$ – $n_{\overline ν_e}$ as a function of propagation angle is not sufficient to generate instability. Our analysis accounts for this fact and allows us to formulate complementary criteria. Using FORNAX simulation data, we show that fast collective oscillations qualitatively depend on how forward peaked the neutrino angular distributions are.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Melting of magnesium oxide up to two terapascals using double-shock compression

Constraining the melting behavior of magnesium oxide, a major constituent of gaseous and rocky planets, is key to benchmarking their evolutionary models. Using a double-shock technique, we extended the MgO melt curve measurements to 2 TPa; this is twice the pressure achieved by previous melting experiments on any material. A temperature plateau is observed between 1218 and 1950 GPa in the second shock states, which is attributed to latent heat of melting. At 1950 GPa, the measured melting temperature is 17,600 K, 17% lower than recent theoretical predictions. Furthermore, the melting curve is steeper than that of MgSiO 3 , indicating that MgO is likely solid in the interior of Saturn-sized gas giants and extra-solar super-Earth planets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2021

The Physics of the Coronal-line Region for Galaxies in Mapping Galaxies at Apache Point Observatory

The fundamental nature and extent of the coronal-line region (CLR), which may serve as a vital tracer for active galactic nucleus (AGN) activity, remain unresolved. Previous studies suggest that the CLR is produced by AGNdriven outflows and occupies a distinct region between the broad-line region and the narrow-line region, which places it tens to hundreds of parsecs from the galactic center. Here, we investigate 10 coronal line (CL; ionization potential ≳100 eV) emitting galaxies from the Sloan Digital Sky Survey IV Mapping Galaxies at Apache Point Observatory (MaNGA) catalog with emission from one or more CLs detected at ≥5$σ$ above the continuum in at least 10 spaxels—the largest such MaNGA catalog. We find that the CLR is far more extended, reaching out to 1.3–23 kpc from the galactic center. We crossmatch our sample of 10 CL galaxies with the largest existing MaNGA AGN catalog and identify seven in it; two of the remaining three are galaxy mergers and the final one is an AGN candidate. Further, we measure the average CLR electron temperatures as ranging between 12,331 and 22,530 K, slightly above the typical threshold for pure AGN photoionization (~20,000 K) and indicative of shocks (e.g., merger induced or from supernova remnants) in the CLR. In conclusion, we reason that ionizing photons emitted by the central continuum source (i.e., AGN photoionization) primarily generate the CLs, and that energetic shocks are an additional ionization mechanism that likely produce the most extended CLRs we measure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Modeling radiative-shocks created by laser–cluster interactions

Radiative-shocks induced by laser–cluster interactions are modeled using radiation-hydrodynamic simulations. A good agreement—in both shock velocity and density profiles—is obtained between experiment and simulations, indicating that non-local thermodynamic equilibrium (NLTE) radiative effects are important in the experimental regime examined, particularly at early times (≤30 ns) due to the elevated temperatures (≥35 eV). The enhanced NLTE radiative emission causes the shock to be reduced in amplitude, increased in width, and reduced in propagation velocity, while the amplitude of the radiative precursor is increased. As the density and temperature conditions are relatively modest, this potentially has important implications for the scalings that are used in laboratory–astrophysics to transform between laboratory and astrophysical scales, which do not hold for non-LTE systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Unraveling the role of temperature on the onset of ejecta formation at atomic scales

Shock wave interactions with perturbations on a free surface can lead to the inversion and growth of the perturbation and eventual ejection of a jet of material, referred to as an ejecta microjet, from shocked surfaces. Here, this study carries out large-scale molecular dynamics (MD) simulations using single-crystal Cu and Sn systems with a pre-existing groove to characterize the localized gradients in temperatures and pressures generated that render the microjet formation. MD simulations are carried out for three loading orientations (along the [001], [011], and [111] directions) and shock pressures ranging from 16 GPa to 100 GPa to understand the role of temperatures generated at the groove vertex on the formation of ejecta. The simulations suggest that the interaction of the shock wave with the groove results in a localized increase in temperatures, leading to localized softening at the groove vertex and the generation of an ejecta microjet. For Cu systems, the simulations suggest that the orientation effects on shock wave structures, velocities, and localized softening affect the ejecta formation at low pressures. Jetting is only observed when temperatures at the groove vertex are high enough to induce localized softening (close to melting temperature), and the jet velocity increases with shock pressure. In contrast, the loading orientation is rendered inconsequential for Sn systems due to the melting of the material during shock compression at the pressures chosen. The jet velocities are similar, regardless of crystal orientation in Sn systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗
DOE OSTI2021

A magnetic particle time-of-flight (MagPTOF) diagnostic for measurements of shock and compression-bang time at the NIF [National Ignition Facility]

A magnetic particle time-of-flight (MagPTOF) diagnostic has been designed to measure shock- and compression-bang time using DHe-fusion protons and DD-fusion neutrons, respectively, at the National Ignition Facility (NIF). This capability, in combination with shock-burn weighted areal density measurements, will significantly constrain the modeling of the implosion dynamics. This design is an upgrade to the existing particle time-of-flight (pTOF) diagnostic, which records bang times using DD or DT neutrons with an accuracy better than ±70 ps [H. G. Rinderknecht et al., Rev. Sci. Instrum. 83, 10D902 (2012)]. The inclusion of a deflecting magnet will increase DHe-proton signal-to-background by a factor of 1000, allowing for the first time simultaneous measurements of shock and compression-bang times in DHe-filled surrogate implosions at the NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Pressure-induced phase transformations in Nb 2 O 5 from ab initio molecular dynamics simulations

Here, the equations of state for the H-, L- and T-Nb 2 O 5 polymorphs were investigated with density functional theory. The shock Hugoniot curves were predicted from ab initio molecular dynamics simulations along 300–2500 K isotherms. Enthalpies derived from simulations at 300 K show that the H → L transition occurs at 6.2 GPa, in agreement with the measured value of ~5.9 GPa, and is more favorable than the H → T transition. Pressure-induced amorphization (PIA) of L-Nb 2 O 5 is in the range ~18-25 GPa, consistent with experiments reporting PIA above 21.4 GPa. The simulated PIA of T-Nb 2 O 5 between ~16 and 26 GPa reproduces PIA observed above 19.2 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗