Related papers: Primitive Variable Solvers for Conservative Genera…
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, we propose a conservative nonconforming virtual element method for the full stationary incompressible magnetohydrodynamics model. We leverage the virtual element satisfactory divergence-free property to ensure mass…
We present a multi-dimensional numerical code to solve isothermal magnetohydrodynamic (IMHD) equations for use in modeling astrophysical flows. First, we have built a one-dimensional code which is based on an explicit finite-difference…
A new implementation for magnetohydrodynamics (MHD) simulations in full general relativity (involving dynamical spacetimes) is presented. In our implementation, Einstein's evolution equations are evolved by a BSSN formalism, MHD equations…
We present a new numerical code, X-ECHO, for general relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes. This is aimed at studying astrophysical situations where strong gravity and magnetic fields are both supposed to play an…
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with…
The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO…
We present a class of non-standard numerical schemes which are modifications of the discrete gradient method. They preserve the energy integral exactly (up to the round-off error). The considered class contains locally exact discrete…
Relativistic magnetic reconnection is thought to power various multi-wavelength emission signatures from neutron stars and black holes. Relativistic resistive magnetohydrodynamics (RRMHD) offers the simplest model of reconnection. However,…
In this paper we discuss energy conservation issues related to the numerical solution of the nonlinear wave equation. As is well known, this problem can be cast as a Hamiltonian system that may be autonomous or not, depending on the…
In this work, we design an entropy stable, finite volume approximation for the shallow water magnetohydrodynamics (SWMHD) equations. The method is novel as we design an affordable analytical expression of the numerical interface flux…
We use the general framework of summation-by-parts operators to construct conservative, energy-stable, and well-balanced semidiscretizations of two different nonlinear systems of dispersive shallow water equations with varying bathymetry:…
This paper describes a numerical scheme for multi-fluid hydrodynamics in the limit of small mass densities of the charged particles. The inertia of the charged particles can then be neglected, which makes it possible to write an evolution…
Many systems of current interest in relativistic astrophysics require a knowledge of radiative transfer in a magnetized gas flowing in a strongly-curved, dynamical spacetime. Such systems include coalescing compact binaries containing…
Numerically solving magnetohydrodynamic (MHD) equations faces many challenges: avoiding divergence error, maintaining positivity, and satisfying entropy conditions. Among discontinuous Galerkin (DG) schemes, there has been a modal version…
We assemble the equations of general relativistic magnetohydrodynamics (MHD) in 3+1 form. These consist of the complete coupled set of Maxwell equations for the electromagnetic field, Einstein's equations for the gravitational field, and…
Existing model reduction techniques for high-dimensional models of conservative partial differential equations (PDEs) encounter computational bottlenecks when dealing with systems featuring non-polynomial nonlinearities. This work presents…
We present an improved version of the ECHO-QGP numerical code, which self-consistently includes for the first time the effects of electromagnetic fields within the framework of relativistic magnetohydrodynamics (RMHD). We discuss results of…
A conservative formulation of the drift-reduced fluid plasma model is constructed by analytically inverting the implicit relation defining the polarisation velocity as a function of the time-derivative of the electric field. The obtained…
Active galactic nuclei, x-ray binaries, pulsars, and gamma-ray bursts are all believed to be powered by compact objects surrounded by relativistic plasma flows driving phenomena such as accretion, winds, and jets. These flows are often…