A Divergence-free Upwind Code for Multidimensional Magnetohydrodynamic Flows

A description is given for preserving 𝛁 ∙ = 0 in a magnetohydrodynamic (MHD) code that employs the upwind, total variation diminishing (TVD) scheme and Strang type operator splitting for multidimensionality. The method is based on the staggered mesh technique to constrain the transport of magnetic field: the magnetic field components are defined at grid interfaces with their advective fluxes on grid edges, while other quantities are defined at grid centers. The magnetic field at grid centers for the upwind step is calculated by interpolating the values from grid interfaces. The advective fluxes on grid edges for the magnetic field evolution are calculated from the upwind fluxes at grid interfaces. Then the magnetic field can be maintained with 𝛁 ∙ = 0 exactly, if this is so initially, while the upwind scheme is used for the update of fluid quantities. The correctness of the code is demonstrated through tests comparing numerical solutions either with analytic solutions or with numerical solutions from a code using an explicit divergence-cleaning method. Also, the robustness is shown through tests involving realistic astrophysical problems.

[1]  Dinshaw S. Balsara,et al.  Total Variation Diminishing Scheme for Adiabatic and Isothermal Magnetohydrodynamics , 1998 .

[2]  Paul R. Woodward,et al.  On the Divergence-free Condition and Conservation Laws in Numerical Simulations for Supersonic Magnetohydrodynamical Flows , 1998 .

[3]  A. Frank,et al.  Effects of Cooling on the Propagation of Magnetized Jets , 1997, astro-ph/9711250.

[4]  A. Frank,et al.  The MHD Kelvin-Helmholtz Instability. II. The Roles of Weak and Oblique Fields in Planar Flows , 1997, astro-ph/9701044.

[5]  Department of Materials Science,et al.  The Magnetohydrodynamics of Supersonic Gas Clouds: MHD Cosmic Bullets and Wind-swept Clumps , 1996, astro-ph/9606106.

[6]  Dinshaw S. Balsara,et al.  Notes on the Eigensystem of Magnetohydrodynamics , 1996, SIAM J. Appl. Math..

[7]  A. Frank,et al.  The Magnetohydrodynamic Kelvin-Helmholtz Instability: A Two-dimensional Numerical Study , 1995, astro-ph/9510115.

[8]  D. Ryu,et al.  A MULTI-DIMENSIONAL MAGNETOHYDRODYNAMIC CODE IN CYLINDRICAL GEOMETRY , 1995 .

[9]  D. Ryu,et al.  Numerical magnetohydrodynamics in astrophysics: Algorithm and tests for multidimensional flow , 1994, astro-ph/9505073.

[10]  Paul R. Woodward,et al.  An approximate Riemann solver for ideal magnetohydrodynamics , 1994 .

[11]  Phillip Colella,et al.  A Higher-Order Godunov Method for Multidimensional Ideal Magnetohydrodynamics , 1994, SIAM J. Sci. Comput..

[12]  Phillip Colella,et al.  A higher-order Godunov method for the equations of ideal magnetohydrodynamics , 1992 .

[13]  R. Blandford,et al.  Numerical simulations of magnetized jets , 1989 .

[14]  Phillip Colella,et al.  Higher order Godunov methods for general systems of hyperbolic conservation laws , 1989 .

[15]  C. Richard DeVore,et al.  Flux-corrected transport techniques for multidimensional compressible magnetohydrodynamics , 1989 .

[16]  J. Hawley,et al.  Simulation of magnetohydrodynamic flows: A Constrained transport method , 1988 .

[17]  M. Brio,et al.  An upwind differencing scheme for the equations of ideal magnetohydrodynamics , 1988 .

[18]  Michael L. Norman,et al.  Numerical Simulations of a Magnetically Confined Jet , 1986 .

[19]  J. Brackbill,et al.  The Effect of Nonzero ∇ · B on the numerical solution of the magnetohydrodynamic equations☆ , 1980 .

[20]  S. Orszag,et al.  Small-scale structure of two-dimensional magnetohydrodynamic turbulence , 1979, Journal of Fluid Mechanics.

[21]  Paul R. Woodward,et al.  High-resolution simulations of compressible convection using the piecewise-parabolic method , 1994 .

[22]  M. Norman,et al.  ZEUS-2D : a radiation magnetohydrodynamics code for astrophysical flows in two space dimensions. II : The magnetohydrodynamic algorithms and tests , 1992 .