Related papers: A resistive MHD module in the GPU-accelerated GRMH…
In this paper, the general procedure to solve the General Relativistic Hydrodynamical(GRH) equations with Adaptive-Mesh Refinement (AMR) is presented. In order to achieve, the GRH equations are written in the conservation form to exploit…
In this era of precision cosmology, a detailed physical understanding on the evolution of cosmic baryons is required. Cosmic magnetic fields, though still poorly understood, may represent an important component in the global cosmic energy…
The advent of robust, reliable and accurate higher order Godunov schemes for many of the systems of equations of interest in computational astrophysics has made it important to understand how to solve them in multi-scale fashion. This is so…
We present the implementation of general-relativistic resistive magnetohydrodynamics solvers and three divergence-free handling approaches adopted in the General-relativistic multigrid numerical (Gmunu) code. In particular,…
We find two families of analytic solutions to the ideal magnetohydrodynamics (iMHD) equations, in a class of 4-dimensional (4D) curved spacetimes. The plasma current is null, and as a result, the stress-energy tensor of the plasma itself…
Several MHD works and, in particular, the recent one by Medina-Torrejon et al. (2021) based on three-dimensional MHD simulations of relativistic jets, have evidenced that particle acceleration by magnetic reconnection driven by the…
This paper first studies the admissible state set $\mathcal G$ of relativistic magnetohydrodynamics (RMHD). It paves a way for developing physical-constraints-preserving (PCP) schemes for RMHD equations with the solutions in $\mathcal G$.…
Short gamma-ray bursts originate when relativistic jets emerge from the remnants of binary neutron star mergers. Both the jet and the remnant are believed to be strongly magnetized, and the presence of magnetic fields is known to influence…
This paper describes the main features of a pioneering unsteady solver for simulating ideal two-fluid plasmas on unstructured grids, taking profit of GPGPU (General-purpose computing on graphics processing units). The code, which has been…
The constraint of incompressibility is often used to simplify the magnetohydrodynamic (MHD) description of linearized plasma dynamics because it does not affect the ideal MHD marginal stability point. In this paper two methods for…
Fast reconnection in magnetically dominated plasmas is widely invoked in models of dissipation in pulsar winds, gamma-ray flares in the Crab nebula, and to explain the radio nanoshots of pulsars. When current sheets evolve reaching a…
We present Horizon, a new graphics processing unit (GPU)-accelerated code to solve the equations of general relativistic magnetohydrodynamics in a given spacetime. We evaluate the code in several test cases, including magnetized Riemann…
We introduce an extension to the AthenaK code for general-relativistic magnetohydrodynamics (GRMHD) in dynamical spacetimes using a 3+1 conservative Eulerian formulation. Like the fixed-spacetime GRMHD solver, we use standard finite-volume…
A new computational model is presented suitable for exploring the self-consistent production of energetic electrons during magnetic reconnection in macroscale systems. The equations are based on the recent discovery that parallel electric…
Incorporating the radiation reaction force into two-fluid plasma in curved space, we get a set of one-fluid general relativistic magnetohydrodynamics (GRMHD) equations with the Landau-Lifshitz radiation reaction force. We analyze the…
Magnetised plasma turbulence pervades the universe and is likely to play an important role in a variety of astrophysical settings. Magnetohydrodynamics (MHD) provides the simplest theoretical framework in which phenomenological models for…
We present an extension of the multi-region relaxed magnetohydrodynamics (MRxMHD) equilibrium model that includes plasma flow. This new model is a generalization of Woltjer's model of relaxed magnetohydrodynamics equilibria with flow. We…
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…
GX is a code for solving the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for…
We present nonlinear magnetohydrodynamic (MHD) simulations of high-$\beta$ Wendelstein 7-X plasmas using the stellarator extension of the M3D-$C^1$ code, building on the recent work that shows benign saturation of ideal ballooning modes…