Publications

Abstract not provided.

While magnetized turbulence is ubiquitous in many astrophysical and terrestrial systems, our understanding of even the simplest physical description of this phenomena, ideal magnetohydrodynamic (MHD) turbulence, remains substantially incomplete. In this work, we highlight the shortcomings of existing theoretical and phenomenological descriptions of MHD turbulence that focus on the joint (kinetic and magnetic) energy fluxes and spectra by demonstrating that treating these quantities separately enables fundamental insights into the dynamics of MHD turbulence. This is accomplished through the analysis of the scale-wise energy transfer over time within an implicit large eddy simulation of subsonic, super-Alfvénic MHD turbulence. Our key finding is that the kinetic energy spectrum develops a scaling of approximately k -4/3 in the stationary regime as magnetic tension mediates large-scale kinetic to magnetic energy conversion and significantly suppresses the kinetic energy cascade. This motivates a reevaluation of existing MHD turbulence theories with respect to a more differentiated modeling of the energy fluxes.

Abstract not provided.

Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn-Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼ 20 %.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This article considers two algorithms of a finite-volume solver for the MHD equations with a real-gas equation of state (EOS). Both algorithms use a multistate form of the Harten-Lax-Van Leer approximate Riemann solver as formulated for MHD discontinuities. This solver is modified to use the generalized sound speed from the real-gas EOS. Two methods are tested: EOS evaluation at cell centers and flux interfaces where the former is more computationally efficient. A battery of 1-D and 2-D tests is employed: convergence of 1-D and 2-D linearized waves, shock tube Riemann problems, a 2-D nonlinear circularly polarized Alfvén wave, and a 2-D magneto-Rayleigh-Taylor instability test. The cell-centered-EOS-evaluation algorithm produces unresolvable thermodynamic inconsistencies in the intermediate states leading to spurious solutions while the flux-interface EOS evaluation algorithm robustly produces the correct solution. The linearized wave tests show that this inconsistency is associated with the magnetosonic waves and the magneto-Rayleigh-Taylor instability test demonstrates simulation results, where the spurious solution leads to an unphysical simulation.

Abstract not provided.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.