An operator-splitting numerical scheme for relativistic magnetohydrodynamics

David Phillips,Serguei Komissarov
2023-11-07
Abstract:We describe a novel operator-splitting approach to numerical relativistic magnetohydrodynamics designed to expand its applicability to the domain of ultra-high magnetisation. In this approach, the electromagnetic field is split into the force-free component, governed by the equations of force-free degenerate electrodynamics (FFDE), and the perturbation component, governed by the perturbation equations derived from the full system of relativistic magnetohydrodynamics (RMHD). The combined system of the FFDE and perturbation equations is integrated simultaneously, for which various numerical techniques developed for hyperbolic conservation laws can be used. At the end of every time-step of numerical integration, the force-free and the perturbation components of the electromagnetic field are recombined and the result is regarded as the initial value of the force-free component for the next time-step, whereas the initial value of the perturbation component is set to zero. To explore the potential of this approach, we build a 3rd-order WENO code, which was used to carry out 1D and 2D test simulations. Their results show that this operator-splitting approach allows us to bypass the stiffness of RMHD in the ultra-high-magnetisation regime where the perturbation component becomes very small. At the same time, the cod
Instrumentation and Methods for Astrophysics,High Energy Astrophysical Phenomena,Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the difficulties in numerical simulations of relativistic magnetohydrodynamics (RMHD) under extremely high magnetization conditions. Specifically: 1. **Numerical Instability under High Magnetization Conditions**: Under ultra - high magnetization parameters \( \sigma \) (for example, \( \sigma \sim 10^3 - 10^6 \)), modern conservative numerical schemes start to fail in multi - dimensional problems. This is because the conserved variables cannot be transformed into physically meaningful fundamental variables. Moreover, in problems involving shocks and magnetic reconnection, numerical resistivity leads to excessive artificial plasma heating. 2. **Limitations of Existing Methods**: For highly magnetized fluids, existing numerical methods are difficult to handle. For example, in highly magnetized astrophysical environments, such as the magnetospheres around black holes and neutron stars, traditional methods cannot effectively simulate the physical processes under these extreme conditions. To solve these problems, the author proposes a new operator - splitting method, which divides the electromagnetic field into two parts: - **Force - free Field Component**: Controlled by the force - free degenerate electrodynamics (FFDE) equations. - **Perturbation Component**: Controlled by the perturbation equations derived from the complete RMHD system. Through this method, the author hopes to be able to effectively perform numerical simulations under ultra - high magnetization conditions and at the same time maintain good performance for medium - and low - magnetization problems. ### Formula Summary - The definition of the magnetization parameter \( \sigma \) is: \[ \sigma=\frac{\hat{B}^2}{4\pi(e + p)} \] where \( \hat{B} \) is the magnetic field strength in the plasma rest frame, \( e=\rho c^2+e_t \) is the total energy density of the plasma, \( \rho \) is the rest mass density, \( e_t \) is the thermal energy density, and \( p \) is the thermal pressure. - Decomposition of the electromagnetic field: \[ F_{\mu\nu}=f^{(0)}_{\mu\nu}+f^{(1)}_{\mu\nu} \] - Decomposition of the stress - energy - momentum tensor: \[ T_{\mu\nu}^{(e)}=t^{(0)}_{\mu\nu}+t^{(1)}_{\mu\nu} \] Through this decomposition method, the author can better handle the numerical simulation problems under high magnetization conditions and has demonstrated the effectiveness and wide applicability of this method in tests.