Shijia Chen / National University of Defense Technology
A new algorithm to model resistive magnetohydrodynamics (MHD) in Magnetised Liner Inertial Fusion (MagLIF) has been developed for one-dimensional axisymmetric geometry. The corresponding differential equations are discretized in the cylindrical Lagrangian coordinates. The algorithm is fully conservative for mass, momentum, and energy. Matter energy and magnetic energy including Nernst effects are managed separately. The evolution of magnetic field and energy are solved using a derivative of semi-implicit scheme, where unconditional stability is obtained without to solve large sparse systems of equations. This MHD package has been integrated into the radiation-hydrodynamics code MULTI-IFE (R. Ramis, et al., Comp. Phys. Comm. 203 (2016) 226–237) that includes hydrodynamics of two-component plasmas, laser light ray-tracing, thermal diffusion, multigroup radiation transport, deuterium-tritium burning, and α particle diffusion. This setup allows to simulate MagLIF configurations relevant for Inertial Confinement Fusion.