Related papers: Multi-dimensional Numerical Scheme for Resistive R…
Motivated by a geometric method employed for the derivation of the Nambu bracket for ideal two-dimensional incompressible hydrodynamics, we reconstruct the reduced magnetohydrodynamic (RMHD) model by a priori imposition of its conservation…
The equations of general relativistic magnetohydrodynamics (GRMHD) have become the standard mathematical framework for modeling high-energy plasmas in curved spacetimes. However, the fragility of the primitive variable reconstruction…
Many interesting terrestrial and astrophysical scenarios involving magnetic fields can be approached in axial geometry. Even though the Lagrangian smoothed particle hydrodynamics (SPH) technique has been successfully extended to handle…
Conservative numerical schemes for general relativistic magnetohydrodynamics (GRMHD) require a method for transforming between ``conserved'' variables such as momentum and energy density and ``primitive'' variables such as rest-mass…
It is known that HLL-type schemes are more dissipative than schemes based on characteristic decompositions. However, HLL-type methods offer greater flexibility to large systems of hyperbolic conservation laws because the eigenstructure of…
Resistive magnetohydrodynamics is thought to play a key role in transient astrophysical phenomena such as black hole flares and neutron star magnetospheres. When performing numerical simulations of resistive magnetohydrodynamics, one is…
We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux conservative, 2nd order, total variation diminishing (TVD), upwind scheme of Jin and Xin.…
In this paper, an energy-consistent finite difference scheme for the compressible hydrodynamic and magnetohydrodynamic (MHD) equations is introduced. For the compressible magnetohydrodynamics, an energy-consistent finite difference…
The characteristic decomposition for GRMHD in the comoving frame of the fluid has been known for a long time. However, it has not been known in the coordinate frame of the simulation and in terms of the conserved variables evolved in…
A flux-splitting method is proposed for the hyperbolic-equation system (HES) of magnetized electron fluids in quasi-neutral plasmas. The numerical fluxes are split into four categories, which are computed by using an upwind method which…
In this paper, we prove the non-uniqueness of three-dimensional magneto-hydrodynamic (MHD) system in $C([0,T];L^2(\mathbb{T}^3))$ for any initial data in $H^{\bar{\beta}}(\mathbb{T}^3)$~($\bar{\beta}>0$), by exhibiting that the total energy…
In this paper we present a two-dimensional, time dependent, numerical simulation of a reconnection current layer in incompressible resistive magnetohydrodynamics with uniform resistivity in the limit of very large Lundquist numbers. We use…
Magnetic fields play an important role in many astrophysical systems and a detailed understanding of their impact on the gas dynamics requires robust numerical simulations. Here we present a new method to evolve the ideal…
High-energy astrophysical systems and compact objects are frequently modeled using ideal relativistic magnetohydrodynamic (MHD) or force-free electrodynamic (FFE) simulations, with the underlying assumption that the discretisation from the…
The advantages of high-order finite difference scheme for astrophysical MHD and turbulence simulations are highlighted. A number of one-dimensional test cases are presented ranging from various shock tests to Parker-type wind solutions.…
This work provides a comparative assessment of several low-dissipation numerical schemes for hyperbolic conservation laws, highlighting their performance relative to the classical Harten-Lax-van Leer (HLL) schemes. The schemes under…
We investigate the large scale evolution of a relativistic magnetic reconnection in an electron-positron pair plasma by a relativistic two-fluid magnetohydrodynamic (MHD) code. We introduce an inter-species friction force as an effective…
The magnetohydrodynamic dynamo equation is derived within general relativity, using the covariant 1+3 approach, for a plasma with finite electric conductivity. This formalism allows for a clear division and interpretation of plasma and…
We present an extension of the Piecewise Parabolic Method to special relativistic fluid dynamics in multidimensions. The scheme is conservative, dimensionally unsplit, and suitable for a general equation of state. Temporal evolution is…
Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new,…