We present a new analysis methodology that allows for the self-consistent integration of multiple diagnostics including nuclear measurements, x-ray imaging, and x-ray power detectors to determine the primary stagnation parameters, such as temperature, pressure, stagnation volume, and mix fraction in magnetized liner inertial fusion (MagLIF) experiments. The analysis uses a simplified model of the stagnation plasma in conjunction with a Bayesian inference framework to determine the most probable configuration that describes the experimental observations while simultaneously revealing the principal uncertainties in the analysis. We validate the approach by using a range of tests including analytic and three-dimensional MHD models. An ensemble of MagLIF experiments is analyzed, and the generalized Lawson criterion χ is estimated for all experiments.
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.
Understanding the effects of contaminant plasmas generated within the Z machine at Sandia is critical to understanding current loss mechanisms. The plasmas are generated at the accelerator electrode surfaces and include desorbed species found in the surface and substrate of the walls. These desorbed species can become ionized. The timing and location of contaminant species desorbed from the wall surface depend non-linearly on the local surface temperature. For accurate modeling, it is necessary to utilize wall heating models to estimate the amount and timing of material desorption. One of these heating mechanisms is Joule heating. We propose several extended semi-analytic magnetic diffusion heating models for computing surface Joule heating and demonstrate their effects for several representative current histories. We quantitatively assess under what circumstances these extensions to classical formulas may provide a validatable improvement to the understanding of contaminant desorption timing.
Tin has been shock compressed to ∼69 GPa on the Hugoniot using Sandia's Z Accelerator. A shockless compression wave closely followed the shock wave to ramp compress the shocked tin and probe a high temperature quasi-isentrope near the melt line. A new hybrid backwards integration - Lagrangian analysis routine was applied to the velocity waveforms to obtain the Lagrangian sound velocity of the tin as a function of particle velocity. Surprisingly, an elastic wave was observed on initial compression from the shock state. The presence of the elastic wave indicates tin possess a small but finite strength at this shock pressure, strongly indicating a (mostly) solid state. High fidelity shock Hugoniot measurements on tin sound velocities in this stress range may be required to refine the shock melting stress for pure tin.
For over a decade, velocimetry based techniques have been used to infer the electrical current delivered to dynamic materials properties experiments on pulsed power drivers such as the Z Machine. Though originally developed for planar load geometries, in recent years, inferring the current delivered to cylindrical coaxial loads has become a valuable diagnostic tool for numerous platforms. Presented is a summary of uncertainties that can propagate through the current inference technique when applied to expanding cylindrical anodes. An equation representing quantitative uncertainty is developed which shows the unfold method to be accurate to a few percent above 10 MA of load current.
By deploying a photon Doppler velocimetry based plasma diagnostic, we have directly observed low density plasma in the load anode/cathode gap of cylindrically converging pulsed power targets. The arrival of this plasma is temporally correlated with gross current loss and subtle power flow differences between the anode and the cathode. The density is in the range where Hall terms in the electromagnetic equations are relevant, but this physics is lacking in the magnetohydrodynamics codes commonly used to design, analyze, and optimize pulsed power experiments. The present work presents evidence of the importance of physics beyond traditional resistive magnetohydrodynamics for the design of pulsed power targets and drivers.
Photonic Doppler velocimetry tracks motion during high-speed, single-event experiments using telecommunication fiber components. The same technology can be applied in situations where there is no actual motion, but rather a change in the optical path length. Migration of plasma into vacuum alters the refractive index near a fiber probe, while intense radiation modifies the refractive index of the fiber itself. These changes can diagnose extreme environments in a flexible, time-resolved manner.
Dynamic materials properties experiments on Sandia National Laboratories Z Machine require increasingly precise electrical current pulse shaping. In the experiment described here, a copper flyer plate is accelerated from rest to 1,600 m/s over a 40 micron flight gap in 50 ns. This flyer then impacts a cerium sample, shock melting the cerium, before subsequent quasi-isentropic ramping to mega-bar pressures. Through predictive simulations, postdicted analysis, and a new computational tool for characterizing inherent Z Machine timing accuracy, qualitative estimates of pulse controllability and experimental design robustness are arrived upon.