A robust and efficient solver based on kinetic schemes for Magnetohydrodynamics (MHD) equations

This paper is devoted to the simulation of compressible magnetohydrodynamic (MHD) flows with the Lattice Boltzmann Method (LBM). The usual LBM is limited to low-Mach flows. We propose a robust and accurate numerical method based on the vectorial kinetic construction of [5,25], which allows us to extend the LBM to arbitrary Mach flows. We also explain how to adjust the numerical viscosity in order to obtain stable and accurate results in smooth or discontinuous parts of the flow and reduced divergence errors. The method can handle shock waves and can be made second order in smooth regions. It is also very well adapted to computing with Graphics Processing Unit (GPU). Our GPU im-plementation in 2D achieves state-of-the-art accuracy, with near-optimal performance. We finally present numerical computations of a tilt instability that demonstrate the capability of the method to handle physically relevant simulations.