General relativistic force-free electrodynamics with a discontinuous Galerkin-finite difference hybrid method

Yoonsoo Kim,Elias R. Most,William Throwe,Saul A. Teukolsky,Nils Deppe
DOI: https://doi.org/10.1103/PhysRevD.109.123019
2024-06-13
Abstract:Relativistic plasmas around compact objects can sometimes be approximated as being force-free. In this limit, the plasma inertia is negligible and the overall dynamics is governed by global electric currents. We present a novel numerical approach for simulating such force-free plasmas, which allows for high accuracy in smooth regions as well as capturing dissipation in current sheets. Using a high-order accurate discontinuous Galerkin method augmented with a conservative finite-difference method, we demonstrate efficient global simulations of black hole and neutron star magnetospheres. In addition to a series of challenging test problems, we show that our approach can-depending on the physical properties of the system and the numerical implementation-be up to 10x more efficient than conventional simulations, with a speedup of 2-3x for most problems we consider in practice.
Instrumentation and Methods for Astrophysics,High Energy Astrophysical Phenomena,General Relativity and Quantum Cosmology,Plasma Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in a strong gravitational field environment (such as around black holes and neutron stars), how to efficiently and accurately simulate the force - free electrodynamics (FFE) behavior of relativistic plasmas. Specifically, the author proposes a new numerical method - the discontinuous Galerkin - finite difference hybrid method to overcome the limitations of existing methods in dealing with smooth regions and discontinuous features (such as current sheets and magnetic reconnection points). ### Background and Problem Description of the Paper 1. **Research Objects** - Relativistic plasmas around black holes and neutron stars. - These celestial bodies have the most powerful magnetic fields in the universe, and the plasma environment around them is highly conductive. 2. **Physical Models** - Under these extreme conditions, the plasma can be approximated as a force - free state, that is, the plasma inertia can be ignored, and the evolution of the system is mainly driven by the overall current. - Force - free electrodynamics (FFE) is an effective theoretical framework for describing the evolution of electromagnetic fields under such extreme conditions. 3. **Limitations of Existing Methods** - Although traditional finite - difference or finite - volume methods are robust, a large number of grid points are required to capture wave solutions (such as Alfvén waves) during long - time integration, resulting in high computational costs. - Although spectral methods have exponential convergence in smooth regions, non - physical oscillations (Gibbs phenomenon) will occur in places with discontinuities (such as current sheets), and it is difficult to parallelize. ### Solutions The author proposes a new numerical method - the discontinuous Galerkin - finite difference hybrid method (DG - FD hybrid method), aiming to combine the advantages of the two methods: - **Discontinuous Galerkin Method (DG)**: Provides high precision in smooth regions. - **Finite - Difference Method (FD)**: Captures dissipation effects in discontinuous regions (such as current sheets). Through this method, the author hopes to achieve the following goals: 1. **Improve Computational Efficiency**: Compared with traditional methods, the new method can achieve a 2 - 3 times speed increase in most practical problems, and in some cases, it can even achieve a 10 - fold efficiency increase. 2. **Accurately Capture Discontinuous Features**: Such as magnetic reconnection points and current sheets to ensure the accuracy of simulation results. 3. **Large - Scale Parallel Computing**: Utilize the high parallelism of the DG method and be suitable for large - scale parallel computing environments. ### Numerical Implementation To achieve the above goals, the author describes in detail the specific implementation steps of the numerical method, including: - **Spatial Discretization**: Use the discontinuous Galerkin method for high - order polynomial expansion. - **Time Integration**: Use the implicit - explicit (IMEX) time integration method to handle stiff source terms. - **Hybrid Solver**: Automatically switch to the finite - difference method when discontinuities are detected, and switch back to the discontinuous Galerkin method after the smoothness is restored. Through these improvements, the author demonstrates the superior performance of this method in a series of challenging test problems and proves its effectiveness in simulating the magnetospheres of black holes and neutron stars.