Volume 9, Issue 2
Approximate Riemann Solvers and Robust High-Order Finite Volume Schemes for Multi-Dimensional Ideal MHD Equations

Franz Georg Fuchs, Andrew D. McMurry, Siddhartha Mishra, Nils Henrik Risebro & Knut Waagan

Commun. Comput. Phys., 9 (2011), pp. 324-362.

Published online: 2011-09

Preview Full PDF 125 1005
Export citation
  • Abstract

We design stable and high-order accurate finite volume schemes for the ideal MHD equations in multi-dimensions. We obtain excellent numerical stability due to some new elements in the algorithm. The schemes are based on three- and five-wave approximate Riemann solvers of the HLL-type, with the novelty that we allow a varying normal magnetic field. This is achieved by considering the semi-conservative Godunov-Powell form of the MHD equations. We show that it is important to discretize the Godunov-Powell source term in the right way, and that the HLL-type solvers naturally provide a stable upwind discretization. Second-order versions of the ENO- and WENO-type reconstructions are proposed, together with precise modifications necessary to preserve positive pressure and density. Extending the discrete source term to second order while maintaining stability requires non-standard techniques, which we present. The first- and second-order schemes are tested on a suite of numerical experiments demonstrating impressive numerical resolution as well as stability, even on very fine meshes.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-9-324, author = {}, title = {Approximate Riemann Solvers and Robust High-Order Finite Volume Schemes for Multi-Dimensional Ideal MHD Equations}, journal = {Communications in Computational Physics}, year = {2011}, volume = {9}, number = {2}, pages = {324--362}, abstract = {

We design stable and high-order accurate finite volume schemes for the ideal MHD equations in multi-dimensions. We obtain excellent numerical stability due to some new elements in the algorithm. The schemes are based on three- and five-wave approximate Riemann solvers of the HLL-type, with the novelty that we allow a varying normal magnetic field. This is achieved by considering the semi-conservative Godunov-Powell form of the MHD equations. We show that it is important to discretize the Godunov-Powell source term in the right way, and that the HLL-type solvers naturally provide a stable upwind discretization. Second-order versions of the ENO- and WENO-type reconstructions are proposed, together with precise modifications necessary to preserve positive pressure and density. Extending the discrete source term to second order while maintaining stability requires non-standard techniques, which we present. The first- and second-order schemes are tested on a suite of numerical experiments demonstrating impressive numerical resolution as well as stability, even on very fine meshes.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.171109.070510a}, url = {http://global-sci.org/intro/article_detail/cicp/7502.html} }
TY - JOUR T1 - Approximate Riemann Solvers and Robust High-Order Finite Volume Schemes for Multi-Dimensional Ideal MHD Equations JO - Communications in Computational Physics VL - 2 SP - 324 EP - 362 PY - 2011 DA - 2011/09 SN - 9 DO - http://dor.org/10.4208/cicp.171109.070510a UR - https://global-sci.org/intro/article_detail/cicp/7502.html KW - AB -

We design stable and high-order accurate finite volume schemes for the ideal MHD equations in multi-dimensions. We obtain excellent numerical stability due to some new elements in the algorithm. The schemes are based on three- and five-wave approximate Riemann solvers of the HLL-type, with the novelty that we allow a varying normal magnetic field. This is achieved by considering the semi-conservative Godunov-Powell form of the MHD equations. We show that it is important to discretize the Godunov-Powell source term in the right way, and that the HLL-type solvers naturally provide a stable upwind discretization. Second-order versions of the ENO- and WENO-type reconstructions are proposed, together with precise modifications necessary to preserve positive pressure and density. Extending the discrete source term to second order while maintaining stability requires non-standard techniques, which we present. The first- and second-order schemes are tested on a suite of numerical experiments demonstrating impressive numerical resolution as well as stability, even on very fine meshes.

Franz Georg Fuchs, Andrew D. McMurry, Siddhartha Mishra, Nils Henrik Risebro & Knut Waagan. (2020). Approximate Riemann Solvers and Robust High-Order Finite Volume Schemes for Multi-Dimensional Ideal MHD Equations. Communications in Computational Physics. 9 (2). 324-362. doi:10.4208/cicp.171109.070510a
Copy to clipboard
The citation has been copied to your clipboard