English

High Order Upwind Schemes for Multidimensional Magnetohydrodynamics

Astrophysics 2009-10-31 v1

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\divb=0 for the magnetic field vector \bb\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\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 (r+1)th(r+1)^{th} order accurate \divb=0\divb=0 relation for any numerical \bb\bb field having a rthr^{th} 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\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.

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