A positive conservative method for magnetohydrodynamics based on HLL and Roe methods.

*(English)*Zbl 0967.76061Summary: The exact Riemann problem solutions of the usual equation of ideal magnetohydrodynamics (MHD) can have negative pressure, if the initial data has \(\nabla\cdot C\neq 0\). This creates a problem for numerical solving because in a first-order finite-volume conservative Godunov-type method one cannot avoid jumps in the normal magnetic field component even if the magnetic field was divergence-free in the three-dimensional sense. We show that by allowing magnetic monopoles in MHD equations and properly taking into account the magnetostatic contribution to the Lorentz force, an additional source term appers in Faraday’s law only. Using the Harten-Lax-van Leer (HLL) Riemann solver and discretizing the source term in a specific manner, we obtain a method which is positive and conservative. We show positivity by extensive numerical experimentation. This MHD-HLL method is positive and conservative but rather diffusive; thus we show how to hybridize this method with the Roe method to obtain a much higher accuracy while still retaining positivity. The result is a fully robust positive conservative scheme for ideal MHD, whose accuracy and efficiency properties are similar to the first-order Roe method and which keeps \(\nabla\cdot{\mathbf B}\) small in the same sense as Powell’s method. As a special case, a method with similar characteristics for accuracy and robustness is obtained for the Euler equations as well.

##### MSC:

76M12 | Finite volume methods applied to problems in fluid mechanics |

76W05 | Magnetohydrodynamics and electrohydrodynamics |

##### Keywords:

Harten-Lax-van Leer Riemann solver; Riemann problem; ideal magnetohydrodynamics; magnetic monopoles; Roe method; positive conservative scheme; Euler equations##### Software:

HLLE
PDF
BibTeX
XML
Cite

\textit{P. Janhunen}, J. Comput. Phys. 160, No. 2, 649--661 (2000; Zbl 0967.76061)

Full Text:
DOI

##### References:

[1] | Brio, M.; Wu, C.C., An upwind differencing scheme for the equations of ideal magnetohydrodynamics, J. comput. phys., 75, 400, (1988) · Zbl 0637.76125 |

[2] | Croisille, J.; Khanfir, R.; Chanteur, G., Numerical simulation of the MHD equations by a kinetic-type method, J. sci. comput., 10, 81, (1995) · Zbl 0839.76062 |

[3] | Einfeldt, B., On Godunov-type methods for gas dynamics, SIAM J. numer. anal., 25, 294, (1988) · Zbl 0642.76088 |

[4] | Einfeldt, B.; Munz, C.D.; Roe, P.L.; Sjögreen, B., On Godunov-type methods near low densities, J. comput. phys., 92, 273, (1991) · Zbl 0709.76102 |

[5] | Godunov, S.K., On symmetric form of the equations of magnetohydrodynamics, (1972) |

[6] | Harten, A.; Lax, P.D.; van Leer, B., On upstream differencing and Godunov-type schemes for hyperbolic conservation laws, SIAM rev., 25, 35, (1983) · Zbl 0565.65051 |

[7] | Jackson, J.D., Classical electrodynamics, (1975) · Zbl 0114.42903 |

[8] | Janhunen, P., GUMICS-3-A global ionosphere-magnetosphere coupling simulation with high ionospheric resolution, Proc. environmental modelling for space-based applications, sept. 18-20, 1996, (1996) |

[9] | Janhunen, P., On the possibility of using an electromagnetic ionosphere in global MHD simulations, Ann. geophys., 16, 397, (1998) |

[10] | Janhunen, P.; Amm, O.; Greenwald, R., Adaptive mesh in global magnetosphere simulation: comparisons with global convection, Portugal Cambridge symposium workshop Proceedings, (1998) |

[11] | LeVeque, R.J., Numerical methods for conservation laws, (1992) · Zbl 0847.65053 |

[12] | Linde, T.J., A three-dimensional adaptive multifluid MHD model of the heliosphere, (1998) |

[13] | Moschetta, J.-M.; Pullin, D.I., A robust low diffusive kinetic scheme for the navier – stokes/Euler equations, J. comput. phys., 133, 193, (1997) · Zbl 0882.76066 |

[14] | Myong, R.S., Theoretical and computational investigations of nonlinear waves in magnetohydrodynamics, (1996) |

[15] | Powell, K.G., An approximate Riemann solver for magnetohydrodynamics (that works in more than one dimension), (1994) |

[16] | Powell, K.G., Solution of the Euler and magnetohydrodynamic equations on solution-adaptive Cartesian grids, Computational fluid dynamics, (March 1996) |

[17] | Raeder, J.; Walker, R.J.; Ashour-Abdalla, M., The structure of the distant geomagnetic tail during long periods of northward IMF, Geophys. res. lett., 22, 349, (1995) |

[18] | Roe, P.L.; Balsara, D.S., Notes on the eigensystem of magnetohydrodynamics, SIAM J. appl. math., 56, 57, (1996) · Zbl 0845.35092 |

[19] | Schmidt-Voigt, M., Time-dependent MHD simulations for cometary plasmas, Astron. astrophys., 210, 433, (1989) · Zbl 0661.76129 |

[20] | Tanaka, T., Finite volume TVD scheme on an unstructured grid system for three-dimensional MHD simulation of inhomogeneous systems including strong background potential fields, J. comput. phys., 111, 381, (1994) · Zbl 0797.76061 |

[21] | Tanaka, T., Generation mechanisms for magnetosphere-ionosphere current systems deduced from a three-dimensional MHD simulation of the solar wind-magnetosphere-ionosphere coupling processes, J. geophys. res., 100, 12,057, (1995) |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.