Related papers: WhiskyMHD: a new numerical code for general relati…
We present a novel high-order nodal artificial viscosity approach designed for solving Magnetohydrodynamics (MHD) equations. Unlike conventional methods, our approach eliminates the need for ad hoc parameters. The viscosity is…
This paper describes the design and implementation of our new multi-group, multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a suite of code tests to validate its application in a wide range of physical regimes.…
Compressible ideal magnetohydrodynamics (MHD) is formulated in terms of the time evolution of potential vorticity and magnetic flux per unit mass using a compact Lie bracket notation. It is demonstrated that this simplifies analytic…
Realistic numerical simulations, i.e., those that make minimal use of ad hoc modeling, are essential for understanding the complex turbulent dynamics of the interiors and atmospheres of the Sun and other stars and the basic mechanisms of…
We consider the Cauchy problem to the three-dimensional isentropic compressible Magnetohydrodynamics (MHD) system with density-dependent viscosities. When the initial density is linearly equivalent to a large constant state, we prove that…
Binary neutron star mergers are among the most energetic events in our Universe, with magnetic fields significantly impacting their dynamics, particularly after the merger. While numerical-relativity simulations that correctly describe the…
Although barotropic matter does not constitute a realistic model for magnetic stars, it would be interesting to confirm a recent conjecture that states that magnetized stars with a barotropic equation of state would be dynamically unstable…
We present a new approach for stably evolving general relativistic magnetohydrodynamic (GRMHD) simulations in regions where the magnetization $\sigma=b^2/\rho c^2$ becomes large. GRMHD codes typically struggle to evolve plasma above…
In an attempt to investigate the structures of ultra-relativistic jets injected into the intracluster medium (ICM) and the associated flow dynamics, such as shocks, velocity shear, and turbulence, we have developed a new special…
Motivated by the desire for highly accurate numerical computations of compact binary spacetimes in the era of gravitational wave astronomy, we reexamine hyperbolicity and well-posedness of the initial value problem for popular models of…
We have recently constructed a numerical code that evolves a spherically symmetric spacetime using a hyperbolic formulation of Einstein's equations. For the case of a Schwarzschild black hole, this code works well at early times, but…
The 2D time dependent solution of thin accretion disk in a close binary system have been presented on the equatorial plane around the Schwarzschild black hole. To do that, the special part of the General Relativistic Hydrodynamical(GRH)…
We present the first results on the construction of an exascale hyperbolic PDE engine (ExaHyPE), a code for the next generation of supercomputers with the objective to evolve dynamical spacetimes of black holes, neutron stars and binaries.…
In this paper we propose a novel thermodynamically compatible finite volume scheme for the numerical solution of the equations of magnetohydrodynamics (MHD) in one and two space dimensions. As shown by Godunov in 1972, the MHD system can be…
We report on a new 3D numerical code designed to solve the Einstein equations for general vacuum spacetimes. This code is based on the standard 3+1 approach using cartesian coordinates. We discuss the numerical techniques used in developing…
Binary systems of compact objects with electromagnetic field are modeled by helically symmetric Einstein-Maxwell spacetimes with charged and magnetized perfect fluids. Previously derived thermodynamic laws for helically-symmetric…
Iharm3D is an open-source C code for simulating black hole accretion systems in arbitrary stationary spacetimes using ideal general-relativistic magnetohydrodynamics (GRMHD). It is an implementation of the HARM ("High Accuracy Relativistic…
A new multidimensional simulation code for relativistic two-fluid electrodynamics (RTFED) is described. The basic equations consist of the full set of Maxwell's equations coupled with relativistic hydrodynamic equations for separate two…
We present a new numerical tool for solving the special relativistic ideal MHD equations that is based on the combination of the following three key features: (i) a one-step ADER discontinuous Galerkin (DG) scheme that allows for an…
Astrophysical relativistic flow problems require high resolution three-dimensional numerical simulations. In this paper, we describe a new parallel three-dimensional code for simulations of special relativistic hydrodynamics (SRHD) using…