Publications

Abstract not provided.

The shock Hugoniot for full-density and porous CeO2 was investigated in the liquid regime using ab initio molecular dynamics (AIMD) simulations with Erpenbeck's approach based on the Rankine-Hugoniot jump conditions. The phase space was sampled by carrying out NVT simulations for isotherms between 6000 and 100 000 K and densities ranging from ρ=2.5 to 20g/cm3. The impact of on-site Coulomb interaction corrections +U on the equation of state (EOS) obtained from AIMD simulations was assessed by direct comparison with results from standard density functional theory simulations. Classical molecular dynamics (CMD) simulations were also performed to model atomic-scale shock compression of larger porous CeO2 models. Results from AIMD and CMD compression simulations compare favorably with Z-machine shock data to 525 GPa and gas-gun data to 109 GPa for porous CeO2 samples. Using results from AIMD simulations, an accurate liquid-regime Mie-Grüneisen EOS was built for CeO2. In addition, a revised multiphase SESAME-Type EOS was constrained using AIMD results and experimental data generated in this work. This study demonstrates the necessity of acquiring data in the porous regime to increase the reliability of existing analytical EOS models.

Absolute Hugoniot measurements for empty plastic foams at ∼10% of solid polystyrene density and supporting rad-hydro simulation results are reported. Planar foam slabs, ∼400 μm thick and ∼500 μm wide, some of which were covered with a 10 μm solid plastic ablator, were directly driven by 4 ns long Nike krypton-fluoride 248 nm wavelength laser pulses that produced strong shock waves in the foam. The shock and mass velocities in our experiments were up to 104 km/s and 84 km/s, respectively, and the shock pressures up to ∼9 Mbar. The motion of the shock and ablation fronts was recorded using side-on monochromatic x-ray imaging radiography. The steadiness of the observed shock and ablation fronts within ∼1% has been verified. The Hugoniot data inferred from our velocity measurements agree with the predictions of the SESAME and CALEOS equation-of-state models near the highest pressure ∼9 Mbar and density compression ratio ∼5. In the lower pressure range 2-5 Mbar, a lower shock density compression is observed than that predicted by the models. Possible causes for this discrepancy are discussed.

Presented are the results from the liner ablation experiments conducted at 550 kA on the Michigan Accelerator for Inductive Z-Pinch Experiments. These experiments were performed to evaluate a hypothesis that the electrothermal instability (ETI) is responsible for the seeding of magnetohydrodynamic instabilities and that the cumulative growth of ETI is primarily dependent on the material-specific ratio of critical temperature to melting temperature. This ratio is lower in refractory metals (e.g., tantalum) than in non-refractory metals (e.g., aluminum or titanium). The experimental observations presented herein reveal that the plasma-vacuum interface is remarkably stable in tantalum liner ablations. This stability is particularly evident when contrasted with the observations from aluminum and titanium experiments. These results are important to various programs in pulsed-power-driven plasma physics that depend on liner implosion stability. Examples include the magnetized liner inertial fusion (MagLIF) program and the cylindrical dynamic material properties program at Sandia National Laboratories, where liner experiments are conducted on the 27-MA Z facility.

Abstract not provided.

In recent years, α-quartz has been used prolifically as an impedance matching standard in shock wave experiments in the multi-Mbar regime (1 Mbar = 100 GPa = 0.1 TPa). This is due to the fact that above ∼90-100 GPa along the principal Hugoniot α-quartz becomes reflective, and thus, shock velocities can be measured to high precision using velocity interferometry. The Hugoniot and release of α-quartz have been studied extensively, enabling the development of an analytical release model for use in impedance matching. However, this analytical release model has only been validated over a range of 300-1200 GPa (0.3-1.2 TPa). Here, we extend this analytical model to 200-3000 GPa (0.2-3 TPa) through additional α-quartz Hugoniot and release measurements, as well as first-principles molecular dynamics calculations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Included in this memo are the final values for the BSE and SiO₂ critical points. VASP 5.3.3 was used to calculate the values in each table. Standard methodology is to continue equilibrating the simulation until the block averaged standard deviation is less than 1%. However, due to the simulation sizes in these simulations, the average standard deviation of the BSE pressure is 5% (ranging between 2% and 7%, which is much higher than the usual under 1% considered acceptable) and the SiO₂ pressure is 6% with the lower density error bars being significantly larger than the highest density. The critical point is chosen where the dP/dρ >= 0 for an entire isotherm but is often difficult to distinguish because of the noise along each isotherm. As such, the actual values of the critical points are found by “best guess”. Further computer simulations creating a more pressure values along each isotherm would help but time constraints preclude this. Similarly, the ability to calculate pressures at lower densities that are currently feasible with VASP would also make determination of the critical point easier.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The electrical conductivity of materials under extremes of temperature and pressure is of crucial importance for a wide variety of phenomena, including planetary modeling, inertial confinement fusion, and pulsed power based dynamic materials experiments. There is a dearth of experimental techniques and data for highly compressed materials, even at known states such as along the principal isentrope and Hugoniot, where many pulsed power experiments occur. We present a method for developing, calibrating, and validating material conductivity models as used in magnetohydrodynamic (MHD) simulations. The difficulty in calibrating a conductivity model is in knowing where the model should be modified. Our method isolates those regions that will have an impact. It also quantitatively prioritizes which regions will have the most beneficial impact. Finally, it tracks the quantitative improvements to the conductivity model during each incremental adjustment. In this paper, we use an experiment on Sandia National Laboratories Z-machine to isentropically launch multiple flyer plates and, with the MHD code ALEGRA and the optimization code DAKOTA, calibrated the conductivity such that we matched an experimental figure of merit to +/-1%.

Abstract not provided.

Abstract not provided.

We report on a new technique for obtaining off-Hugoniot pressure vs. density data for solid metals compressed to extreme pressure by a magnetically driven liner implosion on the Z-machine (Z) at Sandia National Laboratories. In our experiments, the liner comprises inner and outer metal tubes. The inner tube is composed of a sample material (e.g., Ta and Cu) whose compressed state is to be inferred. The outer tube is composed of Al and serves as the current carrying cathode. Another aluminum liner at much larger radius serves as the anode. A shaped current pulse quasi-isentropically compresses the sample as it implodes. The iterative method used to infer pressure vs. density requires two velocity measurements. Photonic Doppler velocimetry probes measure the implosion velocity of the free (inner) surface of the sample material and the explosion velocity of the anode free (outer) surface. These two velocities are used in conjunction with magnetohydrodynamic simulation and mathematical optimization to obtain the current driving the liner implosion, and to infer pressure and density in the sample through maximum compression. This new equation of state calibration technique is illustrated using a simulated experiment with a Cu sample. Monte Carlo uncertainty quantification of synthetic data establishes convergence criteria for experiments. Results are presented from experiments with Al/Ta, Al/Cu, and Al liners. Symmetric liner implosion with quasi-isentropic compression to peak pressure ∼1000 GPa is achieved in all cases. These experiments exhibit unexpectedly softer behavior above 200 GPa, which we conjecture is related to differences in the actual and modeled properties of aluminum.

In this paper, we explore magnetized liner inertial fusion (MagLIF) [S. A. Slutz et al., Phys. Plasmas 17, 056303 (2010)] using a semi-analytic model [R. D. McBride and S. A. Slutz, Phys. Plasmas 22, 052708 (2015)]. Specifically, we present simulation results from this model that: (a) illustrate the parameter space, energetics, and overall system efficiencies of MagLIF; (b) demonstrate the dependence of radiative loss rates on the radial fraction of the fuel that is preheated; (c) explore some of the recent experimental results of the MagLIF program at Sandia National Laboratories [M. R. Gomez et al., Phys. Rev. Lett. 113, 155003 (2014)]; (d) highlight the experimental challenges presently facing the MagLIF program; and (e) demonstrate how increases to the preheat energy, fuel density, axial magnetic field, and drive current could affect future MagLIF performance.

Abstract not provided.

Abstract not provided.

Abstract not provided.