Transport and the Equation of State for Asymmetric Plasma Mixtures

Robert Rudd (15-ERD-052)

Abstract

We plan to develop a new capability to calculate thermodynamic and transport properties of asymmetric plasma mixtures——that is, plasmas with two or more components differing significantly in charge and mass in the presence of strong gradients. Inertial Confinement Fusion with ablator materials mixing into the deuterium–—tritium plasma is such an example. The molecular dynamics plasma simulation code ddcMD developed in the Cimarron project— (which was dedicated to the prediction and measurement of properties of dense, highly collisional, and strongly coupled plasmas—) will form the basis of this capability. Specifically, we will deliver a high-fidelity description of plasmas, including cases where there is strong ion––ion coupling for at least one component. We will develop new molecular dynamics techniques for transport calculations including accurate equations of state of asymmetric plasma mixtures, even in the presence of strong gradients such as sharp interfaces. Our initial simulations will focus on an interface at rest, and then proceed to simulations of an accelerated interface including diffusion brought about by pressure gradients and the Rayleigh––Taylor hydrodynamic instability.

Upon successful completion of this project, we will develop the capability to conduct plasma simulations with strong gradients, including those in the strongly coupled regime. We expect to develop simulations of diffusion in asymmetric plasmas with strong gradients by properly treating the plasma equation of state (not making an ideal gas approximation) as well as the screening and electric fields associated with variations in electron density. We will have calculated the reflecting power of the interface relevant to questions about Knudsen molecular diffusion effects, and have demonstrated the performance of the code in massively parallel simulations on Livermore supercomputers.

Mission Relevance

The new capability will address LLNL's strategic need to predict transport processes in plasma and warm dense matter systems that are heterogeneous, such as at an interface. The new capability will be most relevant to the Laboratory's strategic focus area in stockpile stewardship science, in which new approaches are sought to model plasmas to better understand boost physics in weapons. In addition, this project supports core competencies in high-performance computing, simulation, and data science as well as high-energy-density science. The project also has applications to DOE Office of Science strategy, which specifically identifies nonequilibrium plasma properties as a priority research area.

FY15 Accomplishments and Results

In FY15 we (1) implemented and tested, in ddcMD, the capability for large-scale (million-ion) molecular dynamics simulation of screened Coulomb plasmas in which the electron density (and hence the screening) varies in space, and included the electron pressure; (2) conducted simulations of diffusion at a constant pressure interface between the deuterium––tritium plasma and a plastic ablator, as found in conditions for Inertial Confinement Fusion; (3) improved the equation-of-state fidelity in ddcMD; and (4) began to develop a unique fast multipole capability for systems with spatially varying screening.

Plasma viscosity can impact inertial confinement fusion hot-spot structure and shock-wave profiles. the viscosity of a deuterium–argon plasma mixture at 100, 200, and 500 ev and ion density of 10<sup>25</sup> ion/cm<sup>3</sup> is shown calculated in two ways. the symbols denote values of the shear viscosity computed using a green–kubo analysis of molecular dynamics simulations, and the solid lines denote the predictions of a hybrid kinetic molecular dynamics viscosity model, which combines the predictions
 
Plasma viscosity can impact Inertial Confinement Fusion hot-spot structure and shock-wave profiles. The viscosity of a deuterium–argon plasma mixture at 100, 200, and 500 eV and ion density of 1025 ion/cm3 is shown calculated in two ways. The symbols denote values of the shear viscosity computed using a Green–Kubo analysis of molecular dynamics simulations, and the solid lines denote the predictions of a hybrid kinetic molecular dynamics viscosity model, which combines the predictions of two models from the literature. The kinetic molecular dynamics model agrees well with the molecular dynamics values (within 16%) over the range of conditions studied here, providing an analytic model of the viscosity plasma mixtures suitable for hydrodynamic simulations.

Publications and Presentations

  • Haxhimali, T., et al., “Shear viscosity for dense plasmas by equilibrium molecular dynamics in asymmetric Yukawa ionic mixtures.” Phys. Rev. E 92, 053110 (2015). LLNL-JRNL-674655.
  • Haxhimali, T., et al., Study of shear viscosity for dense plasmas by equilibrium molecular dynamics in asymmetric Yukawa ionic mixtures. 57th Ann. Mtg. APS Division of Plasma Physics, Savannah, GA, Nov. 16–20, 2015. LLNL-ABS-675170.