Related papers: An adaptive mesh, GPU-accelerated, and error minim…
There is great interest in numerical relativity simulations involving matter due to the likelihood that binary compact objects involving neutron stars will be detected by gravitational wave observatories in the coming years, as well as to…
In this paper, we describe a numerical algorithm for the self-consistent simulations of surface water and sediment dynamics. The method is based on the original Lagrangian-Eulerian CSPH-TVD approach for solving the Saint-Venant and Exner…
The characteristic decomposition for GRMHD is not known in a form useful for current numerical simulations. This prevents us from using the most accurate known computational methods, such as full-wave Riemann solvers. In this paper, we…
We present the extension of the differentiable hydrodynamics code, diffhydro, enabling scalable PDE-constrained inference and integrated hybrid physics-ML models for a wide range of astrophysical applications. New physics additions include…
We have developed a one-dimensional code to solve ultra-relativistic hydrodynamic problems, using the Glimm method for an accurate treatment of shocks and contact discontinuities. The implementation of the Glimm method is based on an exact…
Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel…
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.…
Computational magnetohydrodynamics (MHD) for space physics has become an essential area in understanding the multiscale dynamics of geophysical and astrophysical plasma processes, partially motivated by the lack of space data. Full MHD…
We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement…
We investigate the ideal and incompressible magnetohydrodynamic (MHD) equations in three space dimensions for the development of potentially singular structures. The methodology consists in implementing the four-fold symmetries of the…
The accurate modelling of astrophysical scenarios involving compact objects and magnetic fields, such as the collapse of rotating magnetized stars to black holes or the phenomenology of gamma-ray bursts, requires the solution of the…
In certain astrophysical systems the commonly employed ideal magnetohydrodynamics (MHD) approximation breaks down. Here, we introduce novel explicit and implicit numerical schemes of ohmic resistivity terms in the moving-mesh code AREPO. We…
This paper develops high-order accurate entropy stable (ES) adaptive moving mesh finite difference schemes for the two- and three-dimensional special relativistic hydrodynamic (RHD) and magnetohydrodynamic (RMHD) equations, which is the…
Numerical methods for radiative transfer play a key role in modern-day astrophysics and cosmology, including study of the inhomogeneous reionization process. In this context, ray tracing methods are well-regarded for accuracy but notorious…
In this paper we test a special-relativistic formulation of Smoothed Particle Hydrodynamics (SPH) that has been derived from the Lagrangian of an ideal fluid. Apart from its symmetry in the particle indices, the new formulation differs from…
We present a new computational method for smoothly matching general relativistic ideal magnetohydrodynamics (MHD) to its force-free limit. The method is based on a flux-conservative formalism for MHD and its force-free limit, and a vector…
We have proposed a new High Resolution Shock Capturing (HRSC) scheme for Special Relativistic Hydrodynamics (SRHD) based on the semidiscrete central Godunov-type schemes and a modified Weighted Essentially Non-oscillatory (WENO) data…
This paper highlights first steps towards enabling graphics processing unit (GPU) acceleration of the task-parallel smoothed particle hydrodynamics (SPH) solver SWIFT. Novel combinations of algorithms are presented, enabling SWIFT to…
We consider the relativistic hydrodynamics of non-perfect fluids with the goal of determining a formulation that is suited for numerical integration in special-relativistic and general-relativistic scenarios. To this end, we review the…
We propose a new Harten-Lax-van Leer discontinuities (HLLD) approximate Riemann solver to improve the stability of shocks and the accuracy of low-speed flows in multidimensional magnetohydrodynamic (MHD) simulations. Stringent benchmark…