Shock-Induced Mix Across an Ideal Interface - Supplemental Material

Published: 07-10-2016| Version 1 | DOI: 10.17632/w273xgdtyy.1
Claudio Bellei,
Peter Amendt


This Mendeley Data repository contains: 1. Supplemental material for Fig 1. of the paper "Shock-Induced Mix Across an Ideal Interface" by C. Bellei and P. A. Amendt. The movies were obtained from a particle-in-cell simulation using the code Lsp. The simulation models a shock wave transiting an interface that separates a high-density deuterium-tritium (DT) plasma from a low density 3He plasma. It shows the formation of a DT ion bunch ahead of the 3He transmitted shock. For more details, please refer to the above paper. 2. Movie of a Noh test, used to benchmark the multi-ion plasma fluid code described in the paper (this is a test in spherical geometry).


Steps to reproduce

PIC simulation ---------------- The test problem presented here refers to a shock with Mach number $M=3.8$ that is about to cross an interface separating a DT plasma with ionization state $Z=1$ from a $^3$He plasma with $Z=2$. * The DT region is initialized with a single ion species, representing an average DT-ion $\langle$DT$\rangle$, with the unshocked number density $n_\mathrm{DT}=10^{20}$ cm$^{-3}$ ($\rho_\mathrm{DT}=4.2$ mg/cc) and temperature $T=100$ eV * the upstream $^3$He region is initialized with a density $n_\mathrm{3He}=10^{19}$ cm$^{-3}$ ($\rho_\mathrm{3He}=0.5$ mg/cc) and $T=1$ eV. * The spatial resolution is $\Delta x = 20 \ \mu$m, with 2500 particles/cell for each ion species and a simulation box of length $L=2.0$ cm. * The $\langle$DT$\rangle$ ions occupy the region $x=[0,1.2]$ cm at $t=0$. * The ions are interacting through a binary collision operator (Takizuka-Abe). * The temporal resolution is $\Delta t=0.27$ ps, several hundreds of times smaller than the slowing-down time of a shocked $\langle$DT$\rangle$-ion in the background $^3$He plasma. * The electrons are treated as a fluid species, with use of an equation of state of an ideal gas (with $\gamma=5/3$) and a flux limiter of 0.01. * An implicit electromagnetic field solver is used. Noh Test ---------- * one 3He ion species * no coupling to plasma electrons (pure hydro simulation with no electric field, no thermalization with electrons) * at t=0: particle density = 10^23 cm-3; temperature = 0.1 eV; velocity = -10^7 cm/s * ideal gas eos