Very high order \(P_NP_M\) schemes on unstructured meshes for the resistive relativistic MHD equations.

*(English)*Zbl 1261.76028Summary: We propose the first better than second order accurate method in space and time for the numerical solution of the resistive relativistic magnetohydrodynamics (RRMHD) equations on unstructured meshes in multiple space dimensions. The nonlinear system under consideration is purely hyperbolic and contains a source term, the one for the evolution of the electric field, that becomes stiff for low values of the resistivity.For the spatial discretization we propose to use high order \(P_NP_M\) schemes as introduced in Dumbser et al. [M. Dumbser, D. Balsara, E. F. Toro, C. D. Munz, A unified framework for the construction of one-step finite volume and discontinuous Galerkin schemes, J. Comput. Phys. 227, No. 18, 8209–8253 (2008; Zbl 1147.65075)] for hyperbolic conservation laws and a high order accurate unsplit time-discretization is achieved using the element-local space-time discontinuous Galerkin approach proposed in Dumbser et al. [M. Dumbser, C. Enaux, E.F. Toro, Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws, J. Comput. Phys. 227, No. 8, 3971–4001 (2008; Zbl 1142.65070)] for one-dimensional balance laws with stiff source terms. The divergence-free character of the magnetic field is accounted for through the divergence cleaning procedure of Dedner et al. [A. Dedner, F. Kemm, D. Kröner, C.-D. Munz, T. Schnitzer, M. Wesenberg, Hyperbolic divergence cleaning for the MHD equations, J. Comput. Phys. 175, No. 2, 645–673 (2002; Zbl 1059.76040)]. To validate our high order method we first solve some numerical test cases for which exact analytical reference solutions are known and we also show numerical convergence studies in the stiff limit of the RRMHD equations using \(P_NP_M\) schemes from third to fifth order of accuracy in space and time. We also present some applications with shock waves such as a classical shock tube problem with different values for the conductivity as well as a relativistic MHD rotor problem and the relativistic equivalent of the Orszag-Tang vortex problem. We have verified that the proposed method can handle equally well the resistive regime and the stiff limit of ideal relativistic MHD. For these reasons it provides a powerful tool for relativistic astrophysical simulations involving the appearance of magnetic reconnection.

##### MSC:

76M25 | Other numerical methods (fluid mechanics) (MSC2010) |

65M60 | Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs |

76W05 | Magnetohydrodynamics and electrohydrodynamics |

85A30 | Hydrodynamic and hydromagnetic problems in astronomy and astrophysics |

##### Keywords:

resistive relativistic magnetohydrodynamics; unstructured meshes; stiff source terms; high order finite volume and discontinuous Galerkin methods; \(P_NP_M\) schemes##### Software:

ECHO
PDF
BibTeX
XML
Cite

\textit{M. Dumbser} and \textit{O. Zanotti}, J. Comput. Phys. 228, No. 18, 6991--7006 (2009; Zbl 1261.76028)

Full Text:
DOI

##### References:

[1] | Balsara, D., Total variation diminishing scheme for relativistic magneto-hydrodynamics, The astrophysical journal supplement series, 132, 83-101, (2001) |

[2] | Balsara, D.; Spicer, D., A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations, Journal of computational physics, 149, 270-292, (1999) · Zbl 0936.76051 |

[3] | Balsara, D.S.; Rumpf, T.; Dumbser, M.; Munz, C.D., Efficient high accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics, Journal of computational physics, 228, 2480-2516, (2009) · Zbl 1275.76169 |

[4] | T.J. Barth, P.O. Frederickson, Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction, in: AIAA Paper No. 90-0013, 28th Aerospace Sciences Meeting, January 1990. |

[5] | G. Cardano, Artis magnae sive de regulis algebraicis liber unus, Petreius, Nürnberg, Germany, 1545. · Zbl 1216.01001 |

[6] | Dahlburg, R.B.; Picone, J.M., Evolution of the orszagtang vortex system in a compressible medium. I. initial average subsonic flow, Physics of fluids B, 1, 2153-2171, (1989) |

[7] | Dedner, A.; Kemm, F.; Kröner, D.; Munz, C.-D.; Schnitzer, T.; Wesenberg, M., Hyperbolic divergence cleaning for the MHD equations, Journal of computational physics, 175, 645-673, (2002) · Zbl 1059.76040 |

[8] | Del Zanna, L.; Zanotti, O.; Bucciantini, N.; Londrillo, P., ECHO: a Eulerian conservative high-order scheme for general relativistic magnetohydrodynamics and magnetodynamics, Astronomy & astrophysics, 473, October, 11-30, (2007) |

[9] | di Matteo, T., Magnetic reconnection: flares and coronal heating in active galactic nuclei, Monthly notices of the royal astronomical society, 299, August, L15+, (1998) |

[10] | Dumbser, M.; Balsara, D.; Toro, E.F.; Munz, C.D., A unified framework for the construction of one-step finite-volume and discontinuous Galerkin schemes, Journal of computational physics, 227, 8209-8253, (2008) · Zbl 1147.65075 |

[11] | Dumbser, M.; Enaux, C.; Toro, E.F., Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws, Journal of computational physics, 227, 3971-4001, (2008) · Zbl 1142.65070 |

[12] | Dumbser, M.; Käser, M., Arbitrary high order non-oscillatory finite volume schemes on unstructured meshes for linear hyperbolic systems, Journal of computational physics, 221, 693-723, (2007) · Zbl 1110.65077 |

[13] | Dumbser, M.; Käser, M.; A Titarev, V.; Toro, E.F., Quadrature-free non-oscillatory finite volume schemes on unstructured meshes for nonlinear hyperbolic systems, Journal of computational physics, 226, 204-243, (2007) · Zbl 1124.65074 |

[14] | Dumbser, M.; Munz, C.D., Building blocks for arbitrary high order discontinuous Galerkin schemes, Journal of scientific computing, 27, 215-230, (2006) · Zbl 1115.65100 |

[15] | Giacomazzo, B.; Rezzolla, L., The exact solution of the Riemann problem in relativistic magnetohydrodynamics, Journal of fluid mechanics, 562, 223-259, (2006) · Zbl 1097.76073 |

[16] | Gruzinov, A., Power of an axisymmetric pulsar, Physical review letters, 94, 2, 21101, (2005) |

[17] | Harten, A.; Engquist, B.; Osher, S.; Chakravarthy, S., Uniformly high order essentially non-oscillatory schemes, III, Journal of computational physics, 71, 231-303, (1987) · Zbl 0652.65067 |

[18] | Jaroschek, C.H.; Lesch, H.; Treumann, R.A., Relativistic kinetic reconnection as the possible source mechanism for high variability and flat spectra in extragalactic radio sources, Astrophysical journal letters, 605, April, L9-L12, (2004) |

[19] | Jiang, G.S.; Wu, C.C., A high-order WENO finite difference scheme for the equations of ideal magnetohydrodynamics, Journal of computational physics, 150, 561-594, (1999) · Zbl 0937.76051 |

[20] | Kandus, A.; Tsagas, C.G., Generalized ohm’s law for relativistic plasmas, Monthly notices of the royal astronomical society, 385, April, 883-892, (2008) |

[21] | Klaij, C.; Van der Vegt, J.J.W.; Van der Ven, H., Space – time discontinuous Galerkin method for the compressible navier – stokes equations, Journal of computational physics, 217, 589-611, (2006) · Zbl 1099.76035 |

[22] | Komissarov, S.S., On the properties of Alfvén waves in relativistic magnetohydrodynamics, Physics letters A, 232, February, 435-442, (1997) · Zbl 1053.85500 |

[23] | Komissarov, S.S., Multidimensional numerical scheme for resistive relativistic magnetohydrodynamics, Monthly notices of the royal astronomical society, 382, December, 995-1004, (2007) |

[24] | Komissarov, S.S.; Barkov, M.; Lyutikov, M., Tearing instability in relativistic magnetically dominated plasmas, Monthly notices of the royal astronomical society, 374, 415-426, (2007) |

[25] | Lyubarsky, Y.E., On the relativistic magnetic reconnection, Monthly notices of the royal astronomical society, 358, March, 113-119, (2005) |

[26] | Lyutikov, M., Magnetar giant flares and afterglows as relativistic magnetized explosions, Monthly notices of the royal astronomical society, 367, April, 1594-1602, (2006) |

[27] | Lyutikov, M.; Uzdensky, D., Dynamics of relativistic reconnection, Astrophysical journal, 589, June, 893-901, (2003) |

[28] | Orszag, S.A.; Tang, C.M., Small-scale structure of two-dimensional magnetohydrodynamic turbulence, Journal of fluid mechanics, 90, 129, (1979) |

[29] | Palenzuela, C.; Lehner, L.; Reula, O.; Rezzolla, L., Beyond ideal MHD: towards a more realistic modeling of relativistic astrophysical plasmas, Monthly notices of the royal astronomical society, (2009) |

[30] | Picone, J.M.; Dahlburg, R.B., Evolution of the orszag – tang vortex system in a compressible medium. II. supersonic flow, Physics of fluids B, 3, 29-44, (1991) |

[31] | Schopper, R.; Lesch, H.; Birk, G.T., Magnetic reconnection and particle acceleration in active galactic nuclei, Astronomy & astrophysics, 335, July, 26-32, (1998) |

[32] | Stroud, A.H., Approximate calculation of multiple integrals, (1971), Prentice-Hall Inc. Englewood Cliffs, NJ · Zbl 0379.65013 |

[33] | Titarev, V.A.; Toro, E.F., ADER schemes for three-dimensional nonlinear hyperbolic systems, Journal of computational physics, 204, 715-736, (2005) · Zbl 1060.65641 |

[34] | Uzdensky, D.A., On the axisymmetric force-free pulsar magnetosphere, Astrophysical journal, 598, November, 446-457, (2003) |

[35] | van der Vegt, J.J.W.; van der Ven, H., Space – time discontinuous Galerkin finite element method with dynamic grid motion for inviscid compressible flows I. general formulation, Journal of computation al physics, 182, 546-585, (2002) · Zbl 1057.76553 |

[36] | Del Zanna, L.; Bucciantini, N.; Londrillo, P., An efficient shock-capturing central-type scheme for multidimensional relativistic flows II. magnetohydrodynamics, Astronomy and astrophysics, 400, 397-413, (2003) · Zbl 1222.76122 |

[37] | S. Zenitani, M. Hesse, A. Klimas. Two-Fluid MHD Simulations of Relativistic Magnetic Reconnection, ArXiv e-prints, February 2009. |

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.