Numerical simulation of unsteady MHD flows and applications

We present a robust numerical method for solving the compressible Ideal Magneto-Hydrodynamic equations. It is based on the Residual Distribution (RD) algorithms already successfully tested in many problems. We adapted the scheme to the multi-dimensional unsteady MHD model. The constraint ∇ · B = 0 is enforced by the use a Generalized Lagrange Multiplier (GLM) technique. First, we present this complete system and the keys to get its eigensystem, as we may need it in the algorithm. Next, we introduce the numerical scheme built in order to get a compressible, unsteady and implicit solver which has good shock-capturing properties and is second-order accurate at the converged state. To show the efficiency of our method, we will then comment some 2D results. We will end by pointing out some issues and the extensions we plan for this solver