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 =0 for the magnetic field vector , is proposed. The suggested procedure is based on 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 (r+1)th order accurate =0 relation for any numerical field having a rth 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.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.