A High-Order Relativistic Two-Fluid Electrodynamic Scheme with Consistent Reconstruction of Electromagnetic Fields and a Multidimensional Riemann Solver for Electromagnetism
arXiv:1603.06975 · doi:10.1016/j.jcp.2016.05.006
Abstract
In various astrophysics settings it is common to have a two-fluid relativistic plasma that interacts with the electromagnetic field. While it is common to ignore the displacement current in the ideal, classical magnetohydrodynamic limit, when the flows become relativistic this approximation is less than absolutely well-justified. In such a situation, it is more natural to consider a positively charged fluid made up of positrons or protons interacting with a negatively charged fluid made up of electrons. The two fluids interact collectively with the full set of Maxwell's equations. As a result, a solution strategy for that coupled system of equations is sought and found here. Our strategy extends to higher orders, providing increasing accuracy. Three important innovations are reported here. In our first innovation, the magnetic field within each zone is reconstructed in a divergence-free fashion while the electric field within each zone is reconstructed in a form that is consistent with Gauss' law. In our second innovation, a multidimensionally upwinded strategy is presented which ensures that the magnetic field can be updated via a discrete interpretation of Faraday's law and the electric field can be updated via a discrete interpretation of the generalized Ampere's law. Our third innovation consists of an efficient design for several popular IMEX schemes so that they provide strong coupling between the finite-volume-based fluid solver and the electromagnetic fields at high order. Several accuracy analyses are presented showing that our method meets its design accuracy in the MHD limit as well as in the limit of electromagnetic wave propagation. Several stringent test problems are also presented. We also present a relativistic version of the GEM problem, which shows that our algorithm can successfully adapt to challenging problems in high energy astrophysics.
arXiv commentResubmitted to Journal of Computational Physics after major revision
References in corpus (31)
- HARM: A Numerical Scheme for General Relativistic Magnetohydrodynamics
- Particle Acceleration in Relativistic Magnetized Collisionless Electron-Ion Shocks
- ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
- Second Order Accurate Schemes for Magnetohydrodynamics With Divergence-Free Reconstruction
- Divergence-Free Adaptive Mesh Refinement for Magnetohydrodynamics
- An efficient shock-capturing central-type scheme for multidimensional relativistic flows. II. Magnetohydrodynamics
- Efficient, High Accuracy ADER-WENO Schemes for Hydrodynamics and Divergence-Free Magnetohydrodynamics
- Multidimensional HLLE Riemann solver; Application to Euler and Magnetohydrodynamic Flows
- Very High Order $\PNM$ Schemes on Unstructured Meshes for the Resistive Relativistic MHD Equations
- A Two-dimensional HLLC Riemann Solver for Conservation Laws : Application to Euler and MHD Flows
- Divergence-Free Reconstruction of Magnetic Fields and WENO Schemes for Magnetohydrodynamics
- An HLLC Solver for Relativistic Flows -- II. Magnetohydrodynamics
- An Intercomparison Between Divergence-Cleaning and Staggered Mesh Formulations for Numerical Magnetohydrodynamics
- Multi-dimensional Numerical Scheme for Resistive Relativistic MHD
- A five-wave HLL Riemann solver for relativistic MHD
- A Second Order Godunov Method for Multidimensional Relativistic Magnetohydrodynamics
- Tearing instability in relativistic magnetically dominated plasmas
- Equation of State in Numerical Relativistic Hydrodynamics
- Relativistic Magnetohydrodynamics: Renormalized eigenvectors and full wave decomposition Riemann solver
- Two-dimensional MHD simulations of relativistic magnetic reconnection
- A discontinuous Galerkin method for solving the fluid and MHD equations in astrophysical simulations
- A Discontinuous Galerkin Method for Ideal Two-Fluid Plasma Equations
- A Subluminal Relativistic Magnetohydrodynamics Scheme with ADER-WENO Predictor and Multidimensional Riemann Solver-Based Corrector
- Entropy Stable Numerical Schemes for Two-Fluid Plasma Equations
- Total-Variation-Diminishing Implicit-Explicit Runge-Kutta Methods for the Simulation of Double-Diffusive Convection in Astrophysics
- The Role of Superluminal Electromagnetic Waves in Pulsar Wind Termination Shocks
- Relativistic Two-fluid Simulations of Guide Field Magnetic Reconnection
- A multi-dimensional numerical scheme for two-fluid Relativistic MHD
- Numerical Construction of Magnetosphere with Relativistic Two-fluid Plasma Flows
- A Second-order Divergence-constrained Multidimensional Numerical Scheme for Relativistic Two-Fluid Electrodynamics
- Simulations of pulsar wind formation
Cited by in corpus (19)
- Chiral magnetohydrodynamic turbulence in core-collapse supernovae
- A Subluminal Relativistic Magnetohydrodynamics Scheme with ADER-WENO Predictor and Multidimensional Riemann Solver-Based Corrector
- A Constrained Transport Method for the Solution of the Resistive Relativistic MHD Equations
- Exact and Locally Implicit Source Term Solvers for Multifluid-Maxwell Systems
- von Neumann Stability Analysis of Globally Constraint-Preserving DGTD and PNPM Schemes for the Maxwell Equations using Multidimensional Riemann Solvers
- Globally constraint-preserving FR/DG scheme for Maxwell's equations at all orders
- A Six-moment Multi-fluid Plasma Model
- A numerical approach to the non-uniqueness problem of cosmic ray two-fluid equations at shocks
- A Second-order Divergence-constrained Multidimensional Numerical Scheme for Relativistic Two-Fluid Electrodynamics
- An Energy Conserving Vlasov Solver That Tolerates Coarse Velocity Space Resolutions: Simulation of MMS Reconnection Events
- Resistive and multi-fluid RMHD on graphics processing units
- High-order finite-difference entropy stable schemes for two-fluid relativistic plasma flow equations
- High-order finite-volume integration schemes for subsonic magnetohydrodynamics
- Magnetically driven coupling in relativistic radiation-mediated shocks
- Relativistic resistive magnetohydrodynamics for a two-component plasma
- Multidimensional Hall magnetohydrodynamics with isotropic or anisotropic thermal pressure: numerical scheme and its validation using solitary waves
- An Alternative Finite Difference WENO-like Scheme with Physical Constraint Preservation for Divergence-Preserving Hyperbolic Systems
- Techniques, Tricks and Algorithms for Efficient GPU-Based Processing of Higher Order Hyperbolic PDEs
- MIR: a general-relativistic resistive-magneto-hydrodynamic code to study the effect of resistivity in Neutron Star dynamics