6533b857fe1ef96bd12b4647

RESEARCH PRODUCT

Numerically solving the relativistic Grad–Shafranov equation in Kerr spacetimes: numerical techniques

J. F. MahlmannM. A. AloyPablo Cerdá-durán

subject

High Energy Astrophysical Phenomena (astro-ph.HE)Physics010308 nuclear & particles physicsGeneralizationRotational symmetryFOS: Physical sciencesAstronomy and AstrophysicsGeneral Relativity and Quantum Cosmology (gr-qc)01 natural sciencesGeneral Relativity and Quantum CosmologyMagnetic fieldGrad–Shafranov equationQuality (physics)Space and Planetary Science0103 physical sciencesConvergence (routing)Applied mathematicsPoint (geometry)Astrophysics - High Energy Astrophysical Phenomena010303 astronomy & astrophysicsNumerical stability

description

The study of the electrodynamics of static, axisymmetric and force-free Kerr magnetospheres relies vastly on solutions of the so called relativistic Grad-Shafranov equation (GSE). Different numerical approaches to the solution of the GSE have been introduced in the literature, but none of them has been fully assessed from the numerical point of view in terms of efficiency and quality of the solutions found. We present a generalization of these algorithms and give detailed background on the algorithmic implementation. We assess the numerical stability of the implemented algorithms and quantify the convergence of the presented methodology for the most established setups (split-monopole, paraboloidal, BH-disk, uniform).

https://doi.org/10.1093/mnras/sty858