SEARCH · Search NASA

Results for “PLASMA FLUID EQUATIONS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

969 records · Page 54

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Towards QCD-assisted hydrodynamics for heavy-ion collision phenomenology

Heavy-ion collisions are well described by a dynamical evolution with a long hydrodynamical phase. In this phase the properties of the strongly coupled quark-gluon plasma are reflected in the equation of state (EoS) and the transport coefficients, most prominently by the shear and bulk viscosity over entropy density ratios η/s( T ) and ς/s( T ), respectively. While the EoS is by now known to a high accuracy, the transport coefficients and in particular their temperature and density dependence are not well known from first-principle computations yet, as well as the possible influence they can have once used in hydrodynamical simulations. In this work, the most recent QCD-based parameters are provided as input to the MUSIC framework. A ratio η/s( T ) computed with a QCD based approach is used for the first time [1, 2]. The IP-Glasma model is used to describe the initial energy density distribution, and UrQMD for the dilute hadronic phase. Simulations are performed for Pb–Pb collisions at √ s NN = 2.76 TeV, for different centrality intervals. The resulting kinematic distributions of the particles produced in the collisions are compared to data from the LHC, for several experimental observables. Furthermore, the high precision of the experimental results and the broad variety of observables considered allow to critically verify the quality of the description based on first-principle input to the hydrodynamic evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

On fast solid-body rotation of the solar core and differential (liquid-like) rotation of the solar surface

On the basis of a two-component (two-fluid) hydrodynamic model, it is shown that the probable phenomenon of solar core rotation with a velocity higher than the average velocity of global rotation of the Sun, discovered by the SOHO mission, can be related to fast solid-body rotation of the light hydrogen component of the solar plasma, which is caused by thermonuclear fusion of hydrogen into helium inside the hot dense solar core. Thermonuclear fusion of four protons into a helium nucleus (α-particle) creates a large free specific volume per unit particle due to the large difference between the densities of the solar plasma and nuclear matter. As a result, an efficient volumetric sink of one of the components of the solar substance—hydrogen—forms inside the solar core. Therefore, a steady-state radial proton flux converging to the center should exist inside the Sun, which maintains a constant concentration of hydrogen as it burns out in the solar core. It is demonstrated that such a converging flux of hydrogen plasma with the radial velocity v{sub r}(r) = −βr creates a convective, v{sub r}∂v{sub φ}/∂r, and a local Coriolis, v{sub r}v{sub φ}/r,φ nonlinear hydrodynamic forces in the solar plasma, rotating with the azimuthal velocity v{sub φ}. In the absence of dissipation, these forces should cause an exponential growth of the solid-body rotation velocity of the hydrogen component inside the solar core. However, friction between the hydrogen and helium components of the solar plasma due to Coulomb collisions of protons with α-particles results in a steady-state regime of rotation of the hydrogen component in the solar core with an angular velocity substantially exceeding the global rotational velocity of the Sun. It is suggested that the observed differential (liquid-like) rotation of the visible surface of the Sun (photosphere) with the maximum angular velocity at the equator is caused by sold-body rotation of the solar plasma in the radiation zone and strong turbulence in the tachocline layer, where the turbulent viscosity reaches its maximum value at the equator. There, the tachocline layer exerts the most efficient drag on the less dense outer layers of the solar plasma, which are slowed down due to the interaction with the ambient space plasma (solar wind).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Nonlinear instability of an Oldroyd elastico–viscous magnetic nanofluid saturated in a porous medium

Through viscoelastic potential theory, a Kelvin-Helmholtz instability of two semi-infinite fluid layers, of Oldroydian viscoelastic magnetic nanofluids (MNF), is investigated. The system is saturated by porous medium through two semi-infinite fluid layers. The Oldroyd B model is utilized to describe the rheological behavior of viscoelastic MNF. The system is influenced by uniform oblique magnetic field that acts at the surface of separation. The model is used for the MNF incorporated the effects of uniform basic streaming and viscoelasticity. Therefore, a mathematical simplification must be considered. A linear stability analysis, based upon the normal modes analysis, is utilized to find out the solutions of the equations of motion. The onset criterion of stability is derived; analytically and graphs have been plotted by giving numerical values to the various parameters. These graphs depict the stability characteristics. Regions of stability and instability are identified and discussed in some depth. Some previous studies are recovered upon appropriate data choices. The stability criterion in case of ignoring the relaxation stress times is also derived. To relax the mathematical manipulation of the nonlinear approach, the linearity of the equations of motion is taken into account in correspondence with the nonlinear boundary conditions. Taylor's theory is adopted to expand the governing nonlinear characteristic equation according to of the multiple time scales technique. This analysis leads to the well-known Ginzburg–Landau equation, which governs the stability criteria. The stability criteria are achieved theoretically. To simplify the mathematical manipulation, a special case is considered to achieve the numerical estimations. The influence of orientation of the magnetic fields on the stability configuration, in linear as well as nonlinear approaches, makes a dual role for the magnetic field strength in the stability graphs. Stability diagram is plotted for different sets of physical parameters. A new stability as well as instability region, in the parameter space, appears due to the nonlinear effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2011

IONIC COMPOSITION STRUCTURE OF CORONAL MASS EJECTIONS IN AXISYMMETRIC MAGNETOHYDRODYNAMIC MODELS

We present the ionic charge state composition structure derived from axisymmetric MHD simulations of coronal mass ejections (CMEs), initiated via the flux-cancellation and magnetic breakout mechanisms. The flux-cancellation CME simulation is run on the Magnetohydrodynamics-on-A-Sphere code developed at Predictive Sciences, Inc., and the magnetic breakout CME simulation is run on ARC7 developed at NASA GSFC. Both MHD codes include field-aligned thermal conduction, radiative losses, and coronal heating terms which make the energy equations sufficient to calculate reasonable temperatures associated with the steady-state solar wind and model the eruptive flare heating during CME formation and eruption. We systematically track a grid of Lagrangian plasma parcels through the simulation data and calculate the coronal density and temperature history of the plasma in and around the CME magnetic flux ropes. The simulation data are then used to integrate the continuity equations for the ionic charge states of several heavy ion species under the assumption that they act as passive tracers in the MHD flow. We construct two-dimensional spatial distributions of commonly measured ionic charge state ratios in carbon, oxygen, silicon, and iron that are typically elevated in interplanetary coronal mass ejection (ICME) plasma. We find that the slower CME eruption has relatively enhanced ionic charge states and the faster CME eruption shows basically no enhancement in charge states-which is the opposite trend to what is seen in the in situ ICME observations. The primary cause of the difference in the ionic charge states in the two simulations is not due to the different CME initiation mechanisms per se. Rather, the difference lies in their respective implementation of the coronal heating which governs the steady-state solar wind, density and temperature profiles, the duration of the connectivity of the CME to the eruptive flare current sheet, and the contribution of the flare-heated plasma associated with the reconnection jet outflow into the ejecta. Despite the limitations inherent in the first attempt at this novel procedure, the simulation results provide strong evidence in support of the conclusion that enhanced heavy ion charge states within CMEs are a direct consequence of flare heating in the low corona. We also discuss future improvements through combining numerical CME modeling with quantitative ionic charge state calculations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Sparsified time-dependent Fourier neural operators for fusion simulations

This paper presents a sparsified Fourier neural operator for coupled time-dependent partial differential equations (ST-FNO) as an efficient machine learning surrogate for fluid and particle-based fusion codes such as NIMROD (Non-Ideal Magnetohydrodynamics with Rotation - Open Discussion) and GTC (Gyrokinetic Toroidal Code). ST-FNO leverages the structures in the governing equations and utilizes neural operators to represent Green's function-like numerical operators in the corresponding numerical solvers. Once trained, ST-FNO can rapidly and accurately predict dynamics in fusion devices compared with first-principle numerical algorithms. In general, ST-FNO represents an efficient and accurate machine learning surrogate for numerical simulators for multi-variable nonlinear time-dependent partial differential equations, with the proposed architectures and loss functions. The efficacy of ST-FNO has been demonstrated using quiescent H-mode simulation data from NIMROD and kink-mode simulation data from GTC. The ST-FNO H-mode results show orders of magnitude reduction in memory and central processing unit usage in comparison with the numerical solvers in NIMROD when computing fields over a selected poloidal plane. The ST-FNO kink-mode results achieve a factor of 2 reduction in the number of parameters compared to baseline FNO models without accuracy loss.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Effect of causality constraints on Bayesian analyses of heavy-ion collisions

There have long been questions about the limits to the validity of relativistic fluid dynamics and whether it is being used outside its regime of validity in modern simulations of relativistic heavy-ion collisions. An important new tool for answering this question is a causality analysis in the nonlinear regime—if the solutions of the evolution equations do not respect relativistic causality, then they are not a faithful representation of the underlying relativistic theory (in this case, quantum chromodynamics). Using this nonlinear criterion, it has recently been shown that hydrodynamics is indeed being used outside its regime of validity in simulations, at least sometimes. Here we explore the phenomenological implications, particularly the quantitative effects of demanding limits on acausality in modern Bayesian parameter estimation. We find that, while typically only a small fraction of the system's energy is initially in an acausal regime, placing strict limits on the allowed energy fraction significantly changes the preferred properties of the initial condition, which in turn alters the extracted medium properties such as bulk viscosity, where large values are no longer favored. Furthermore, these findings highlight the importance of developing better theoretical descriptions of the early-time, out-of-equilibrium dynamics of relativistic heavy-ion collisions.

Bayesian methods↗

DOE OSTI2022

Performance Analysis of Speculative Parallel Adaptive Local Timestepping for Conservation Laws

Stable simulation of conservation laws, such as those used to model fluid dynamics and plasma physics applications, requires the satisfaction of the so-called Courant-Friedrichs-Lewy condition. By allowing regions of the mesh to advance with different timesteps that locally satisfy this stability constraint, significant work reduction can be attained when compared to a time integration scheme using a single timestep size. However, parallelizing this algorithm presents considerable difficulty. Since the stability condition depends on the state of the system, dependencies become dynamic and potentially non-local. In this article, we present an adaptive local timestepping algorithm using an optimistic (Timewarp-based) parallel discrete event simulation. We introduce waiting heuristics to limit misspeculation and a semi-static load balancing scheme to eliminate load imbalance as parts of the mesh require finer or coarser timesteps. Last, we outline an interface for separating the physics of the specific conservation law from the temporal integration allowing for productive adoption of our proposed algorithm. We present a misspeculation study for three conservation laws, demonstrating both the productivity of the local timestepping API, for which 74% of the lines of code are reused across different conservation laws, and the robustness of the waiting heuristics—at most 1.5% of element updates are rolled back. Our performance studies demonstrate up to a 2.8× speedup versus a baseline unoptimized local timestepping approach, a 4x improvement in per-node throughput compared to an MPI parallelization of synchronous timestepping, and scalability up to 3,072 cores on NERSC’s Cori Haswell partition.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

Weakly nonlinear Bell-Plesset effects for a uniformly converging cylinder

In this research, a weakly nonlinear (WN) model has been developed considering the growth of a small perturbation on a cylindrical interface between two incompressible fluids which is subject to arbitrary radial motion. We derive evolution equations for the perturbation amplitude up to third order, which can depict the linear growth of the fundamental mode, the generation of the second and third harmonics, and the third-order (second-order) feedback to the fundamental mode (zero-order). WN solutions are obtained for a special uniformly convergent case. WN analyses are performed to address the dependence of interface profiles, amplitudes of inward-going and outward-going parts, and saturation amplitudes of linear growth of the fundamental mode on the Atwood number, the mode number (m), and the initial perturbation. The difference of WN evolution in cylindrical geometry from that in planar geometry is discussed in some detail. It is shown that interface profiles are determined mainly by the inward and outward motions rather than bubbles and spikes. The amplitudes of inward-going and outward-going parts are strongly dependent on the Atwood number and the initial perturbation. For low-mode perturbations, the linear growth of fundamental mode cannot be saturated by the third-order feedback. For fixed Atwood numbers and initial perturbations, the linear growth of fundamental mode can be saturated with increasing m. The saturation amplitude of linear growth of the fundamental mode is typically 0.2λ–0.6λ for m < 100, with λ being the perturbation wavelength. Thus, it should be included in applications where Bell-Plesset [G. I. Bell, Los Alamos Scientific Laboratory Report No. LA-1321, 1951; M. S. Plesset, J. Appl. Phys. 25, 96 (1954)] converging geometry effects play a pivotal role, such as inertial confinement fusion implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Wave topology in Hall magnetohydrodynamics

Hall magnetohydrodynamics (HMHD) extends ideal MHD by incorporating the Hall effect via the induction equation, making it more accurate for describing plasma behavior at length scales below the ion skin depth. Despite its importance, a comprehensive description of the eigenmodes in HMHD has been lacking. In this work, we derive the complete spectrum and eigenvectors of HMHD waves and identify their underlying topological structure. We prove that the HMHD wave spectrum is homotopic to that of ideal MHD, consisting of three distinct branches: the slow magnetosonic-Hall waves, the shear Alfvén-Hall waves, and the fast magnetosonic-Hall waves, which continuously reduce to their ideal MHD counterparts in the limit of vanishing Hall parameter. Contrary to a recent claim [Mahajan, Sharma, and Lingam, Phys. Plasmas 31, 090701 (2024)], we find that HMHD does not admit any additional wave branches beyond those in ideal MHD. In conclusion, the key qualitative difference lies in the topological nature of the HMHD wave structure: it exhibits nontrivial topology characterized by a Weyl point—an isolated eigenmode degeneracy point—and associated nonzero Chern numbers of the eigenmode bundles over a 2-sphere in 𝐤-space surrounding the Weyl point.

Alfvén waves↗

DOE OSTI2011

Nonlinear ELM simulations based on a nonideal peeling–ballooning model using the BOUT++ code

A minimum set of equations based on the peeling–ballooning (P–B) model with nonideal physics effects (diamagnetic drift, E × B drift, resistivity and anomalous electron viscosity) is found to simulate pedestal collapse when using the BOUT++ simulation code, developed in part from the original fluid edge code BOUT. Linear simulations of P–B modes find good agreement in growth rate and mode structure with ELITE calculations. The influence of the E × B drift, diamagnetic drift, resistivity, anomalous electron viscosity, ion viscosity and parallel thermal diffusivity on P–B modes is being studied; we find that (1) the diamagnetic drift and E × B drift stabilize the P–B mode in a manner consistent with theoretical expectations; (2) resistivity destabilizes the P–B mode, leading to resistive P–B mode; (3) anomalous electron and parallel ion viscosities destabilize the P–B mode, leading to a viscous P–B mode; (4) perpendicular ion viscosity and parallel thermal diffusivity stabilize the P–B mode. With addition of the anomalous electron viscosity under the assumption that the anomalous kinematic electron viscosity is comparable to the anomalous electron perpendicular thermal diffusivity, or the Prandtl number is close to unity, it is found from nonlinear simulations using a realistic high Lundquist number that the pedestal collapse is limited to the edge region and the ELM size is about 5–10% of the pedestal stored energy. Furthermore, this is consistent with many observations of large ELMs. The estimated island size is consistent with the size of fast pedestal pressure collapse. In the stable α-zones of ideal P–B modes, nonlinear simulations of viscous ballooning modes or current-diffusive ballooning mode (CDBM) for ITER H-mode scenarios are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Inviscid evolution of large amplitude filaments in a uniform gravity field

The inviscid evolution of localized density stratifications under the influence of a uniform gravity field in a homogeneous, ambient background is studied. The fluid is assumed to be incompressible, and the stratification, or filament, is assumed to be initially isotropic and at rest. It is shown that the center of mass energy can be related to the center of mass position in a form analogous to that of a solid object in a gravity field g by introducing an effective gravity field g{sub eff}, which is less than g due to energy that goes into the background and into non-center of mass motion of the filament. During the early stages of the evolution, g{sub eff} is constant in time and can be determined from the solution of a 1D differential equation that depends on the initial, radially varying density profile of the filament. For small amplitude filaments such that ρ{sub 0} ≪ 1, where ρ{sub 0} is the relative amplitude of the filament to the background, the early stage g{sub eff} scales linearly with ρ{sub 0}, but as ρ{sub 0}→∞, g{sub eff}→g and is thus independent of ρ{sub 0}. Fully nonlinear simulations are performed for the evolution of Gaussian filaments, and it is found that the time t{sub max}, which is defined as the time for the center of mass velocity to reach its maximum value U{sub max}, occurs very soon after the constant acceleration phase and so U{sub max}≈g{sub eff}(t=0)t{sub max}. The simulation results show that U{sub max}∼1/t{sub max}∼√(ρ{sub 0}) for ρ{sub 0} ≪ 1, in agreement with theory and results from previous authors, but that U{sub max} and t{sub max} both scale approximately with √(ρ{sub 0}) for ρ{sub 0} ≫ 1. The fact that U{sub max} and t{sub max} have the same scaling with ρ{sub 0} for large amplitude filaments is in agreement with the theory presented in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Nonlinear fluid simulation of particle and heat fluxes during burst of ELMs on DIII-D with BOUT++ code [Fluid Simulation of Particle and Heat Fluxes during Burst of ELMs on DIID with BOUT++ code]

In order to study the distribution and evolution of the transient particle and heat fluxes during edge-localized mode (ELM) bursts, a BOUT++ six-field two-fluid model based on the Braginskii equations with non-ideal physics effects is used to simulate pedestal collapse in divertor geometry. We used the profiles from the DIII-D H-mode discharge #144382 with fast target heat flux measurements as the initial conditions for the simulations. Moreover, a flux-limited parallel thermal conduction is used with three values of the flux-limiting coefficient ${{\alpha}_{j}}$ , free streaming model with ${{\alpha}_{j}}=1$ , sheath-limit with ${{\alpha}_{j}}=0.05$ , and one value in between. The studies show that a 20 times increase in ${{\alpha}_{j}}$ leads to ~6 times increase in the heat flux amplitude to both the inner and outer targets, and the widths of the fluxes are also expanded. The sheath-limit model of flux-limiting coefficient is found to be the most appropriate one, which shows ELM sizes close to the measurements. The evolution of the density profile during the burst of ELMs of DIII-D discharge #144382 is simulated, and the collapse in width and depth of ${{n}_{\text{e}}}$ are reproduced at different time steps. The growing process of the profiles for the heat flux at divertor targets during the burst of ELMs measured by IRTV (infrared television) is also reproduced by this model. The widths of heat fluxes towards targets are a little narrower, and the peak amplitudes are twice the measurements possibly due to the lack of a model of divertor radiation which can effectively reduce the heat fluxes. The magnetic flutter combined with parallel thermal conduction is found to be able to increase the total heat loss by around 33% since the magnetic flutter terms provide the additional conductive heat transport in the radial direction. Finally, the heat flux profile at both the inner and outer targets is obviously broadened by magnetic flutter. The lobe structures near the X-point at LFS are both broadened and elongated due to the magnetic flutter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗

NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Shebalin, John V↗

DOE OSTI2021

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗