Species Separation and Other Kinetic Effects at Shock Wave Fronts in Plasma Mixtures Assessed using Molecular Dynamics
The propagation of a shock through a binary plasma mixture of hydrogen isotopes is simulated with molecular dynamics (MD) and a multispecies continuum hydrodynamics model. Separation of the two species due to their mass difference is observed at the shock front in both simulations. We analyze the shock structure in detail to determine which aspects of the wave propagation and shock profile obey single-fluid hydrodynamics and which require multi-component hydrodynamics with Fickian and barodiffusion. In order to reduce possible sources of uncertainty due to transport models used in the continuum simulations, we have performed separate, equilibrium MD simulations which verify the accuracy of the recently proposed Kinetic Molecular Dynamics (KMD) viscosity model, even at the relatively cold pre-shock temperature of 1 eV. We find the post-shock temperature and density are generally consistent between continuum and atomic-scale calculations, but the MD simulations predict a significantly larger thickness of the leading light-ion shock front, with some associated differences in the species separation. We analyze the shape and anisotropy of the MD velocity distributions in detail and compute the electric field at the shock interface from both the MD and continuum simulations. We find that the continuum hydrodynamics calculation captures the major features of the electric fi eld through the electron pressure term in the barodiffusion model and conclude that the differences observed in the shock structures from the two simulations arise largely from non-Maxwellian effects in the longitudinal velocities (the velocities in the direction that the shock propagates) that are not captured in the hydrodynamic simulations.