paper

High Order Upwind Schemes for Multidimensional Magnetohydrodynamics

arXiv:astro-ph/9910086 · doi:10.1086/308344

Abstract

A general method for constructing high order upwind schemes for multidimensional magnetohydrodynamics (MHD), having as a main built-in condition the divergence-free constraint $\divb=0$ for the magnetic field vector $\bb$, is proposed. The suggested procedure is based on {\em consistency} arguments, by taking into account the specific operator structure of MHD equations with respect to the reference Euler equations of gas-dynamics. This approach leads in a natural way to a staggered representation of the $\bb$ field numerical data where the divergence-free condition in the cell-averaged form, corresponding to second order accurate numerical derivatives, is exactly fulfilled. To extend this property to higher order schemes, we then give general prescriptions to satisfy a order accurate $\divb=0$ relation for any numerical $\bb$ field having a order interpolation accuracy. Consistency arguments lead also to a proper formulation of the upwind procedures needed to integrate the induction equations, assuring the exact conservation in time of the divergence-free condition and the related continuity properties for the $\bb$ vector components. As an application, a third order code to simulate multidimensional MHD flows of astrophysical interest is developed using ENO-based reconstruction algorithms. Several test problems to illustrate and validate the proposed approach are finally presented.

34 pages, including 14 figures

References in corpus (1)

Cited by in corpus (97)