High-order finite-difference entropy stable schemes for two-fluid relativistic plasma flow equations
arXiv:2210.08568 · doi:10.1016/j.jcp.2023.112207
Abstract
In this article, we propose high-order finite-difference entropy stable schemes for the two-fluid relativistic plasma flow equations. This is achieved by exploiting the structure of the equations, which consists of three independent flux components. The first two components describe the ion and electron flows, which are modeled using the relativistic hydrodynamics equation. The third component is Maxwell's equations, which are linear systems. The coupling of the ion and electron flows, and electromagnetic fields is via source terms only. Furthermore, we also show that the source terms do not affect the entropy evolution. To design semi-discrete entropy stable schemes, we extend the RHD entropy stable schemes in Bhoriya et al. to three dimensions. This is then coupled with entropy stable discretization of the Maxwell's equations. Finally, we use SSP-RK schemes to discretize in time. We also propose ARK-IMEX schemes to treat the stiff source terms; the resulting nonlinear set of algebraic equations is local (at each discretization point). These equations are solved using the Newton's Method, which results in an efficient method. The proposed schemes are then tested using various test problems to demonstrate their stability, accuracy and efficiency.
References in corpus (19)
- Kinetic energy preserving and entropy stable finite volume schemes for compressible Euler and Navier-Stokes equations
- An efficient shock-capturing central-type scheme for multidimensional relativistic flows. II. Magnetohydrodynamics
- Global ill-posedness of the isentropic system of gas dynamics
- An HLLC Solver for Relativistic Flows -- II. Magnetohydrodynamics
- Multi-dimensional Numerical Scheme for Resistive Relativistic MHD
- Two-Fluid Magnetohydrodynamic Simulations of Relativistic Magnetic Reconnection
- ENO reconstruction and ENO interpolation are stable
- Resistive Magnetohydrodynamic Simulations of Relativistic Magnetic Reconnection
- A High-Order Relativistic Two-Fluid Electrodynamic Scheme with Consistent Reconstruction of Electromagnetic Fields and a Multidimensional Riemann Solver for Electromagnetism
- A Subluminal Relativistic Magnetohydrodynamics Scheme with ADER-WENO Predictor and Multidimensional Riemann Solver-Based Corrector
- Entropy Stable Numerical Schemes for Two-Fluid Plasma Equations
- The Role of Superluminal Electromagnetic Waves in Pulsar Wind Termination Shocks
- High-order accurate entropy stable finite difference schemes for the shallow water magnetohydrodynamics
- Relativistic Tearing and Drift-kink Instabilities in Two-fluid Simulations
- High-order accurate entropy stable finite difference schemes for one- and two-dimensional special relativistic hydrodynamics
- Relativistic Two-fluid Simulations of Guide Field Magnetic Reconnection
- A multi-dimensional numerical scheme for two-fluid Relativistic MHD
- A Second-order Divergence-constrained Multidimensional Numerical Scheme for Relativistic Two-Fluid Electrodynamics
- Richtmyer-Meshkov Instability of an Imploding Flow with a Two-Fluid Plasma Model