Regular Article
Hyperbolic Divergence Cleaning for the MHD Equations

https://doi.org/10.1006/jcph.2001.6961Get rights and content

Abstract

In simulations of magnetohydrodynamic (MHD) processes the violation of the divergence constraint causes severe stability problems. In this paper we develop and test a new approach to the stabilization of numerical schemes. Our technique can be easily implemented in any existing code since there is no need to modify the solver for the MHD equations. It is based on a modified system in which the divergence constraint is coupled with the conservation laws by introducing a generalized Lagrange multiplier. We suggest a formulation in which the divergence errors are transported to the domain boundaries with the maximal admissible speed and are damped at the same time. This corrected system is hyperbolic and the density, momentum, magnetic induction, and total energy density are still conserved. In comparison to results obtained without correction or with the standard “divergence source terms,” our approach seems to yield more robust schemes with significantly smaller divergence errors.

References (49)

  • W. Dai et al.

    A simple finite difference scheme for multidimensional magnetohydrodynamical equations

    J. Comput. Phys.

    (1998)
  • P. Degond et al.

    Polynomial upwind schemes for hyperbolic systems

    C. R. Acad. Sci. Paris Ser. I, Math.

    (1999)
  • C.R. DeVore

    Flux-corrected transport techniques for multidimensional compressible magnetohydrodynamics

    J. Comput. Phys.

    (1991)
  • B.-N. Jiang et al.

    The origin of spurious solutions in computational electromagnetics

    J. Comput. Phys.

    (1996)
  • B. Marder

    A method for incorporating Gauss' law into electromagnetic PIC codes

    J. Comput. Phys.

    (1987)
  • C.-D. Munz et al.

    Divergence correction techniques for Maxwell solvers based on a hyperbolic model

    J. Comput. Phys.

    (2000)
  • C.-D. Munz et al.

    Maxwell's equations when the charge conservation is not satisfied

    C. R. Acad. Sci. Paris Ser. I, Math.

    (1999)
  • K.G. Powell et al.

    A solution-adaptive upwind scheme for ideal magnetohydrodynamics

    J. Comput. Phys.

    (1999)
  • G. Tóth

    The ▿·B=0 constraint in shock-capturing magnetohydrodynamics codes

    J. Comput. Phys.

    (2000)
  • G. Tóth et al.

    Comparison of some flux corrected transport and total variation diminishing numerical schemes for hydrodynamic and magnetohydrodynamic problems

    J. Comput. Phys.

    (1996)
  • A.L. Zachary et al.

    Note: A higher-order Godunov method for the equations of ideal magnetohydrodynamics

    J. Comput. Phys.

    (1992)
  • N. Aslan

    Computational Investigations of Ideal Magnetohydrodynamic Plasmas with Discontinuities

    (1993)
  • D.S. Balsara

    Linearized formulation of the Riemann problem for adiabatic and isothermal magnetohydrodynamics

    Astrophys. J. Suppl.

    (1998)
  • D.S. Balsara

    Total variation diminishing scheme for adiabatic and isothermal magnetohydrodynamics

    Astrophys. J. Suppl.

    (1998)
  • Cited by (1307)

    View all citing articles on Scopus
    View full text