Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.
We present the development of a pulsed power experimental technique to infer the electrical conductivity of metals from ambient to high energy density conditions. The method is implemented on Thor, a moderate scale (1-2 MA) pulsed power driver. The electrical conductivity of copper at elevated temperature (>4000 K) and pressure (>10 GPa) is determined, and a new tabular material model is developed, guided by density functional theory, which preserves agreement with existing experimental data. Minor modifications (<10%) are found to be necessary to the previous Lee-More-Desjarlais model isotherms in the vicinity of the melt transition in order to account for observed discrepancies with the new experimental data. An analytical model for magnetic direct drive flyer acceleration and Joule heating induced vaporization based on the Tsiolkovsky "rocket equation"is presented to assess sensitivity of the method to minor changes in electrical conductivity.
The outer core of the Earth is composed primarily of liquid iron, and the inner core boundary is governed by the intersection of the melt line and the geotherm. While there are many studies on the thermodynamic equation of state for solid iron, the equation of state of liquid iron is relatively unexplored. We use dynamic compression to diagnose the high-pressure liquid equation of state of iron by utilizing the shock-ramp capability at Sandia National Laboratories’ Z-Machine. This technique enables measurements of material states off the Hugoniot by initially shocking samples and subsequently driving a further, shockless compression. Planetary studies benefit greatly from isentropic, off-Hugoniot experiments since they can cover pressure-temperature (P-T) conditions that are close to adiabatic profiles found in planetary interiors. We used this method to drive iron to P-T conditions similar to those of the Earth’s outer-inner core boundary, along an elevated-temperature isentrope in the liquid from 275 GPa to 400 GPa. We derive the equation of state using a hybrid backward integration – forward Lagrangian technique on particle velocity traces to determine the pressure-density history of the sample. Our results are in excellent agreement with SESAME 92141, a previously published equation of state table. With our data and previous experimental data on liquid iron we provide new information on the iron melting line and derive new parameters for a Vinet-based equation of state. The table and our parameterized equation of state are applied to provide an updated means of modeling the pressure, mass, and density of liquid iron cores in exoplanetary interiors.
Elucidating the mechanisms responsible for sub-microsecond desorption of water and other impurities from electrode surfaces at high heating rates is crucial for understanding pulsed power behavior. Ionization of desorbed impurities in the vacuum regions causes power or current loss; devising methods to limit such desorption during the short time scale of pulsed power is needed to improve corresponding applications. Previous molecular modeling studies have strongly suggested that, under high vacuum conditions, the amount of water impurity desorbing from oxide surfaces on metal electrodes is at a sub-monolayer level at room temperature, which appears insufficient to explain observed pulsed power energy losses at high current densities. In this work, we apply Density Functional Theory (DFT) techniques to show that hydrogen trapped inside iron metal can diffuse into hematite (α-Fe2O3) on the metal surface, ultimately reacting with the oxide to form Fe(II) and H2O. The latter desorbs at elevated temperature and may explain the anomalous amount of desorbed impurity inferred from pulsed-power experiments. We also apply a suite of characterization techniques to demonstrate that when iron metal is heated to 650 °C, the dominant surface oxide component becomes α-Fe2O3. The oxide facets exposed are found to be a mixture of (0001), (10-10), and others, in agreement with the DFT models used.
This report describes the high-level accomplishments from the Plasma Science and Engineering Grand Challenge LDRD at Sandia National Laboratories. The Laboratory has a need to demonstrate predictive capabilities to model plasma phenomena in order to rapidly accelerate engineering development in several mission areas. The purpose of this Grand Challenge LDRD was to advance the fundamental models, methods, and algorithms along with supporting electrode science foundation to enable a revolutionary shift towards predictive plasma engineering design principles. This project integrated the SNL knowledge base in computer science, plasma physics, materials science, applied mathematics, and relevant application engineering to establish new cross-laboratory collaborations on these topics. As an initial exemplar, this project focused efforts on improving multi-scale modeling capabilities that are utilized to predict the electrical power delivery on large-scale pulsed power accelerators. Specifically, this LDRD was structured into three primary research thrusts that, when integrated, enable complex simulations of these devices: (1) the exploration of multi-scale models describing the desorption of contaminants from pulsed power electrodes, (2) the development of improved algorithms and code technologies to treat the multi-physics phenomena required to predict device performance, and (3) the creation of a rigorous verification and validation infrastructure to evaluate the codes and models across a range of challenge problems. These components were integrated into initial demonstrations of the largest simulations of multi-level vacuum power flow completed to-date, executed on the leading HPC computing machines available in the NNSA complex today. These preliminary studies indicate relevant pulsed power engineering design simulations can now be completed in (of order) several days, a significant improvement over pre-LDRD levels of performance.
Hypervelocity impact-driven vaporization is characteristic of late-stage planet formation. Yet the behavior and properties of liquid-vapor mixtures of planetary materials of interest are typically unknown. Multiphase equations of state used in hydrodynamic simulations of planet impacts therefore lack reliable data for this important phenomenon. Here, we present the first constraints on the liquid-vapor critical point and coexistence phase boundary of Mg2SiO4 computed from ab initio molecular dynamics simulations. We found that the vapor is depleted in magnesium and enriched in silica and oxygen, while the coexisting liquid is enriched in magnesium and depleted in oxygen, from which we infer vaporization is incongruent. The critical point was estimated from an equation of state fit to the data. The results are in line with recent calculations of MgSiO3 and together confirm that extant multiphase equation of state (EOS) models used in planetary accretion modeling significantly underestimate the amount of supercritical material postimpact.
Mixing of cold, higher-Z elements into the fuel region of an inertial confinement fusion target spoils the fusion burn efficiency. This mixing process is driven by both "turbulent" and "atomic" mixing processes, the latter being modeled through transport corrections to the basic hydrodynamic models. Recently, there has been a surge in the development of dense plasma transport modeling and the associated transport coefficients; however, experimental validation remains in its infancy. To address this gap in our knowledge of interfacial mixing, Sandia National Laboratories is developing a new experimental platform at the Z-facility to investigate plasma transport in dense plasmas that span the entire warm dense matter regime. Specifically, this platform is being developed to measure species transport across a V/CH interface, using an x-ray driven hohlraum to drive the sample to [?] 190eV over 5ns. The heated sample is diagnosed using radiography optimized to measure the distribution of Vanadium perpendicular the interface. In order to interpret measurements made using this experimental platform, modeling tools that incorporate transport effects in strongly coupled plasmas are required. To this end, we utilize new advances in multi-species kinetic theory, collision models applicable to strongly coupled plasmas and modeling of degenerate electron plasmas to develop such a capability. The resulting kinetic transport code has been applied, along with state-of-the-art radiation hydrodynamic codes, to model the experiments. Results from this modeling effort highlight the importance of strong electric fields, which are present in the kinetic transport code, but absent in the radiation hydrodynamics code, in driving interfacial mixing. Synthetic radiography generated from all of these models reveals the ability of experimental diagnostics to distinguish interfacial mixing driven by a range of transport effects. We demonstrate that the spatial and temporal resolution of radiography diagnostics currently available at the Z-facility can distinguish between these different transport effects when multiple (3 [?] 4) radiographs, separated in time ( [?] 2 ns ) with accurate timing are captured per experiment.