High Order Upwind Schemes for Multidimensional Magnetohydrodynamics
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 for the magnetic field vector , 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 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 relation for any numerical 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 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.
Keywords
Cite
@article{arxiv.astro-ph/9910086,
title = {High Order Upwind Schemes for Multidimensional Magnetohydrodynamics},
author = {P. Londrillo and L. Del Zanna},
journal= {arXiv preprint arXiv:astro-ph/9910086},
year = {2009}
}
Comments
34 pages, including 14 figures