The fast multi-frame X-ray diffraction detector at the Dynamic Compression Sector
Explore the source record for details and available documents.
Sinclair, N.W.↗
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.
Three-dimensional blade vortex interactions
A three-dimensional time dependent Navier-Stokes analysis was applied to the rotor blade vortex interaction problem. The numerical procedure is an iterative implicit procedure using three point central differences to represent spatial derivatives. A series of calculations were made to determine the time steps, pseudo-time steps, iterations, artificial dissipation level, etc. required to maintain a nondissipative vortex. Results show the chosen method to have excellent non-dissipative properties provided the correct parameters are chosen. This study was used to set parameters for both two- and three-dimensional blade vortex interaction studies. The case considered was the interaction between a vortex and a NACA0012 airfoil. The results showed the detailed physics during the interaction including the pressure pulse propagating from the blade. The simulated flow physics was qualitatively similar to that experimentally observed. The BVI phenomena is the result of the buildup and violent collapse of the shock waves and local supersonic pockets on the blade surfaces. The resulting pressure pulse build-up appears to be centered at the blade leading edge.
Davoudzadeh, Farhad↗
Beyond Basic Drag in Interplanetary CME Modeling: Effects of Solar Wind Pileup and High-Speed Streams
Coronal mass ejections (CMEs) cause severe space weather effects throughout our solar system. As a fast CME propagates through interplanetary space, it accumulates solar wind materials at its front. This pileup of materials, or CME-driven sheath, can be important in determining the geoeffectiveness of a CME. We take an existing arrival time model that includes expansion and deformation of the CME flux rope (ANTEATR; Kay & Gopalswamy, 2018, https://doi.org/10.1029/2018JA025780; Kay & Nieves-Chinchilla, 2021a, https://doi.org/10.1029/2020JA028911) and add a pileup procedure (PUP) as a physics-based approach to modeling the CME-driven sheath. ANTEATR-PUP solves the Rankine-Hugoniot equations for an oblique shock to determine the shock speed and sheath density, magnetic field, and temperature. The extra sheath mass affects the background drag calculation. Additionally, ANTEATR can now use any 1D profile for the background solar wind as opposed to the simple empirical models it previously relied upon. We present initial results from ANTEATR-PUP and compare with previous ANTEATR findings. Using results from an MHD simulation, we explore the effects of interactions with a static high-speed stream (HSS) on the CME's and sheath's interplanetary evolution. The drag forces essentially disappear while a CME remains within the HSS, but reappear stronger once the CME exits. The HSS-CME interaction produces the largest changes in the CME and sheath properties at 1 au when it occurs either close to the Sun near the inner simulation boundary at 0.1 au or right before the CME reaches 1 au. We estimate that these changes could significantly affect the geoeffectiveness.
C Kay↗
Use of computer vision for analysis of image datasets from high temperature plasma experiments
Great strides have been made in improving the quality of x-ray radiographs in high energy density plasma experiments, enabled in part by innovations in engineering and manufacturing of integrated circuits and materials. As a consequence, the radiographs of today are filled with a great deal of detail, but few of these features are extracted in a systematic way. Analysis techniques familiar to plasma physicists tend toward brittle 1D lineout or Fourier transform type analyses. The techniques applied to process our data have not kept pace with improvements in the quality of our data. Fortunately, the field of computer vision has a wealth of tools to offer, which have been widely used in industrial imaging and, more recently, adopted in biological imaging. We demonstrate the application of computer vision techniques to the analysis of x-ray radiographs from high energy density plasma experiments, as well as give a brief tutorial on the computer vision techniques themselves. These tools robustly extract 2D contours of shocks, boundaries of inhomogeneities, and secondary flows, thereby allowing for increased automation of analysis, as well as direct and quantitative comparisons with simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Calculation of the detonation state of HN 3 with quantum accuracy
HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Z Line VISAR Analysis with the LineVISAR SMASH Class
The Z line VISAR system (ZLV) is a spatially-resolved velocimeter that measures surface velocities in high-energy density experiments on the Z Pulsed Power Facility to facilitate the investigation of fusion, power flow, and dynamic material physics. The data measured in these experiments are analyzed with the LineVISAR SMASH class, a MATLAB software suite that provides tools for data importation, streak image correction, spatiotemporal registration, wrapped phase computation, phase unwrapping, shock handling, and velocity calculation. This report overviews the LineVISAR class, discusses its use, provides example implementations, and supplies analysis specifics not typically recorded in journal publications.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Planar, longitudinal, compressive waves in solids: Thermodynamics and uniaxial strain restrictions
A new tri-axial pressure-based constitutive expression has been found using Cauchy's stress tensor. This stress state emphasizes pressure and shear stress. The description is a pressure plus an effective shear stress allowing for a constitutive law based on atomic solid-state phase changes in crystalline cells due to pressure plus shear-based dislocation motion commonly associated with plasticity. Pressure has a new role in the material's constitutive response as it is separated from plasticity. The thermo-mechanical system describes third-order Gibbs’ expressions without specific volume restrictions placed upon the material. Isothermally, the ratio of heat to shear work in elastic copper is shown to approach zero at a very low temperature and become larger than one as temperature approaches melting. Wave compression models investigated are elastic and plastic: in fully elastic materials, the planar wave is restricted by Poisson's effect although plastic shear changes this constraint. Plastic deformation, dominated by dissipative shear stresses in uniaxial strain, heats the material while excluding phase changes from hydrostatic pressures. The material properties per se across Hugoniot shocks are described with entropy concepts. Shock waves are exceedingly complex since the constitutive laws are linked at extreme temperatures, pressures, and shear stresses. Isothermal, isentropic, isochoric, and iso-shear conditions are used throughout with Jacobian algebra.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Pre-main-sequence Brackett Emitters in the APOGEE DR17 Catalog: Line Strengths and Physical Properties of Accretion Columns
Very young (t ≲ 10 Myr) stars possess strong magnetic fields that channel ionized gas from the interiors of their circumstellar disks to the surface of the star. Upon impacting the stellar surface, the shocked gas recombines and emits hydrogen spectral lines. To characterize the density and temperature of the gas within these accretion streams, we measure equivalent widths of Brackett (Br) 11–20 emission lines detected in 1101 APOGEE spectra of 326 likely pre-main-sequence accretors. For sources with multiple observations, we measure median epoch-to-epoch line strength variations of 10% in Br11 and 20% in Br20. We also fit the measured line ratios to predictions of radiative transfer models by Kwan & Fischer. We find characteristic best-fit electron densities of n e = 10 11 –10 12 cm -3 , and excitation temperatures that are inversely correlated with electron density (from T ~ 5000 K for n e ~ 10 12 cm -3 to T ~ 12,500 K at n e ~ 10 11 cm -3 ). These physical parameters are in good agreement with predictions from modeling of accretion streams that account for the hydrodynamics and radiative transfer within the accretion stream. We also present a supplementary catalog of line measurements from 9733 spectra of 4255 Brackett emission-line sources in the APOGEE Data Release 17 data set.
79 ASTRONOMY AND ASTROPHYSICS↗
A temperature profile diagnostic for radiation waves on OMEGA-60
Predicting and matching radiation wave propagation with computational models has proven difficult. Information provided by experiments studying radiation flow has been limited when only radiation breakout is measured. We have developed the COAX (co-axial) diagnostic platform to provide spatial temperature profiles of a radiation wave through low density foams as a more detailed constraint for simulations. COAX uses a standard, laser-driven OMEGA-60 halfraum to drive radiation down a titanium-laden silicon oxide foam. Point-projection X-ray absorption spectroscopy perpendicular to the radiation flow measures the spatial profile of titanium ionization. The spectroscopic measurement utilizes a broadband capsule backlighter. Imaging and streak spectroscopy are used to characterize the size and spectrum of this source. Radiography provides an additional constraint by capturing the developing shock as the radiation flow becomes subsonic. The DANTE diagnostic is used to measure the halfraum temperature. Here, we provide a spectroscopic analysis of COAX data to determine temperature, and we describe experimental sources of uncertainty. The temperature is obtained by comparison to multi-temperature synthetic spectra post-processed from radiation-hydrodynamics simulations. Quantitative comparison between data and synthetic spectra generated from temperature profiles at relevant simulation times enable determination of a peak temperature of 114 ± 8 eV at 265 ± 22.4 μm from the halfraum. This represents an improvement over the temperature uncertainties of previous radiation flow experiments. Further refinements to the spectroscopic analysis could achieve ±4 eV. The combination between space-resolved spectroscopy and radiography enables us to determine the distance from the halfraum of both the radiation front and the shock front at the time of measurement. For the example shown in this paper the radiation front position is 600–630 μm at 3.43 ± 0.16 ns and the shock front position is 633 μm at 3.3 ± 0.24 ns.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Experimental shock lithification of water-bearing rock powders
The geology and morphology of the terrestrial planets and their moons have been profoundly affected by impact cratering. Some of these bodies contain substantial quantities of water or ice in their regoliths. It is pointed out that the effects of impacts into water-bearing target rocks are not well understood, and may be significantly different from those produced by identical impacts on a desiccated surface. The present investigation has the objective to determine the effects of water on targets of powdered rock and to seek evidence of impact-induced hydration or clay formation. Samples of andesitic basalt were crushed and sieved, and experiments were conducted on the material smaller than 150 micrometers. These experiments show that the water content of a powdered rock target can strongly affect its physical condition following an impact. A relatively small component (5-15 wt %) of water or ice in planetary surface material could inhibit or prevent lithification throughout much of the shocked volume at an impact site.
Allen, C. C.↗
Calculation of shocked flows by mathematical programming
A framework for using mathematical programming to solve Burgers' equation is presented. The steady-state inviscid Burgers' equation with given boundary conditions is considered, and the physically relevant solution is discretized using a four-point difference scheme for the viscous term and a two-point difference scheme for the inviscid term. The framework is then used to solve the Euler equations for quasi-one-dimensional flows on grids with variable spacing.
Lavery, John E.↗
Physical lumping methods for developing linear reduced models for high speed propulsion systems
In gasdynamic systems, information travels in one direction for supersonic flow and in both directions for subsonic flow. A shock occurs at the transition from supersonic to subsonic flow. Thus, to simulate these systems, any simulation method implemented for the quasi-one-dimensional Euler equations must have the ability to capture the shock. In this paper, a technique combining both backward and central differencing is presented. The equations are subsequently linearized about an operating point and formulated into a linear state space model. After proper implementation of the boundary conditions, the model order is reduced from 123 to less than 10 using the Schur method of balancing. Simulations comparing frequency and step response of the reduced order model and the original system models are presented.
Immel, S. M.↗
Refractive index of lithium fluoride at high dynamic stresses
Alkali halides are prototypical ionic solids whose refractive index is a fundamental property related to their lattice structure and ionic polarizability. In particular, lithium fluoride (LiF) has the largest band gap of any known transparent material and maintains transparency at high pressures, making it well suited for use as an optical window in dynamic compression experiments. While empirical fits to the measured density dependence of the refractive index have been provided, a model that is valid over the entire pressure range of experiments and based on a polarization-based theoretical description is lacking. We present a refractive index model for dynamically compressed LiF based on the Lorentz-Lorenz equation, where the molecular polarizability is determined using a single oscillator model with a strain polarizability parameter Λ=0.73. Here, we show that our modeling approach provides an excellent match to the LiF refractive index data for both shock and ramp compression experiments to 900 GPa (density compression greater than threefold). Additionally, updated fits are provided to determine the refractive index correction for LiF windows used in dynamic compression experiments.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Shock state distributions in porous tantalum and characterization with multipoint velocimetry
Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Numerical simulations of the Los Alamos gapstick experiment
The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime
With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.
3-dimensional systems↗
New Sesame Equations of State for Polypropylene, Polyvinylchloride (PVC), and Nylon 66
I built new Sesame equations of state (EOS) for polypropylene (PP), polyvinylchloride (PVC), and nylon 66. PP has the same stoichiometry as polyethylene and TPX but with a density intermediate between the two. Nylon is a polyamide manufactured chiefly in the form of fibers. PP and PVC are the second- and third-most widely produced polymers in the world. All three are thermoplastics. Their monomer units are depicted in Figure 1. I discuss the theory on which the EOS are based in the following section, then compare results to thermal and shock compression data in Section III.
36 MATERIALS SCIENCE↗
Using neutrons and x rays to measure plasma conditions in a solid sphere of deuterated polyethylene compressed to densities of 35 g/cc at temperatures of 2 keV and pressures of 40 Gbar
This paper describes an experiment that shock compresses the center of a solid deuterated polyethylene sphere, CD 2 , to densities of 35 g/cc and temperatures of 2 keV with corresponding pressure of 40 Gbar. The design employs a strong spherically converging shock launched through a solid ball of material using a Hohlraum radiation drive. As the shock coalesces at the center it produces a hot spot that we characterize by measuring the x-ray self-emission and 2.45 MeV neutrons emitted. Two-dimensional images and time-resolved measurements of the x rays emitted determine the size and time duration of the hot spot, leading to an estimated 2k eV electron temperature. The neutron time of flight spectrometer measures an average ion temperature of 1.06 ± 0.15 keV and neutron yield of 7.0 (±0.5) × 10 9 DD neutrons. Our new distribution function tool enables us to create a forward model of the experimental data based on 1D radiation-hydrodynamic simulations, leading to a better understanding of the plasma conditions that produce the measured neutrons and x rays. Furthermore, our simulations indicate that the x rays are produced in a short-lived hot-dense core over tens of picoseconds, whereas the neutron emission continues for about 200 ps, as the hot core starts to expand, thereby leading to a lower mean temperature of the plasma during neutron production. This finding is in agreement with the experimental data, and we therefore conclude that the forward-modeling is a useful tool for inferring the conditions of the hot spot in a laser-driven implosion during burn.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗