Related papers: Multi-Frequency General Relativistic Radiation-Hyd…
We present a new approach for achieving high-order convergence in fully general-relativistic hydrodynamic simulations. The approach is implemented in WhiskyTHC, a new code that makes use of state-of-the-art numerical schemes and was key in…
We implement general relativistic hydrodynamics in the moving-mesh code AREPO. We also couple a solver for the Einstein field equations employing the conformal flatness approximation. The implementation is validated by evolving isolated…
This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the…
We derive a consistent set of moments equations for CR-magnetohydrodynamics, assuming a gyrotropic distribution function (DF). Unlike previous efforts we derive a closure, akin to the M1 closure in radiation hydrodynamics (RHD), that is…
This paper describes the first steps of development of a new multidimensional time implicit code devoted to the study of hydrodynamical processes in stellar interiors. The code solves the hydrodynamical equations in spherical geometry and…
We present a new code, CASTRO, that solves the multicomponent compressible hydrodynamic equations for astrophysical flows including self-gravity, nuclear reactions and radiation. CASTRO uses an Eulerian grid and incorporates adaptive mesh…
We present a numerical implementation of radiative transfer based on an explicitly photon-conserving advection scheme, where radiative fluxes over the cell interfaces of a structured or unstructured mesh are calculated with a second-order…
We describe the development of a flux-limited gray radiation solver for the compressible astrophysics code, CASTRO. CASTRO uses an Eulerian grid with block-structured adaptive mesh refinement based on a nested hierarchy of…
We present a numerical method to solve the equations of general relativistic hydrodynamics in a given external gravitational field. The method is based on a generalization of Roe's approximate Riemann solver for the non relativistic Euler…
Mixed-frame formulations of radiation-hydrodynamics (RHD), where the radiation quantities are computed in an inertial frame but matter quantities are in a comoving frame, are advantageous because they admit algorithms that conserve energy…
We describe a photon-conserving radiative transfer algorithm, using a spatially-adaptive ray tracing scheme, and its parallel implementation into the adaptive mesh refinement (AMR) cosmological hydrodynamics code, Enzo. By coupling the…
We present general-relativistic radiation magnetohydrodynamics simulations of binary neutron star mergers performed with GR-Athena++. Neutrino transport is treated using a moment-based, energy-integrated scheme (M1), augmented by neutrino…
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 present a conservative numerical method for radiation magnetohydrodynamics with frequency-dependent full transport in stationary spacetimes. This method is stable and accurate for both large and small optical depths and radiation…
We present the new general-relativistic magnetohydrodynamics (GRMHD) capabilities of the Einstein Toolkit, an open-source community-driven numerical relativity and computational relativistic astrophysics code. The GRMHD extension of the…
We derive new equations using the mixed-frame approach for one- and two-dimensional (axisymmetric) time-dependent radiation transport and the associated couplings with matter. Our formulation is multi-group and multi-angle and includes…
We present numerical hydrodynamical evolutions of rapidly rotating relativistic stars, using an axisymmetric, nonlinear relativistic hydrodynamics code. We use four different high-resolution shock-capturing (HRSC) finite-difference schemes…
We modify an existing magnetohydrodynamics algorithm to make it more compatible with a dimensionally-split (DS) framework. It is based on the standard reconstruct-solve-average strategy (using a Riemann solver), and relies on constrained…
The relativistic magnetohydrodynamics (RMHD) set of equations has recently seen increased use in astrophysical computations. Even so, RMHD codes remain fragile. The reconstruction can sometimes yield superluminal velocities in certain parts…
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…