相关论文: An Extension of the Athena++ Framework for Fully C…
We propose a new method to conserve the total energy to round-off error in grid-based codes for hydrodynamic simulations with self-gravity. A formula for the energy flux due to the work done by the the self-gravitational force is given, so…
We present a co-scaling grid formalism and its implementation in the magnetohydrodynamics code Athena++. The formalism relies on flow symmetries in astrophysical problems involving expansion, contraction, and center-of-mass motion. The grid…
We describe the implementation of sophisticated numerical techniques for general-relativistic magnetohydrodynamics simulations in the Athena++ code framework. Improvements over many existing codes include the use of advanced Riemann solvers…
High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $\Delta x$) with only a moderate increase to computational expense. Significant…
This paper proposes a hierarchy of numerical fluxes for the compressible flow equations which are kinetic-energy and pressure equilibrium preserving and asymptotically entropy conservative, i.e., they are able to arbitrarily reduce the…
We present a new general relativistic magnetohydrodynamics (GRMHD) code integrated into the Athena++ framework. Improving upon the techniques used in most GRMHD codes, ours allows the use of advanced, less diffusive Riemann solvers, in…
We present a self-contained overview of GR-Athena++, a general-relativistic magnetohydrodynamics (GRMHD) code, that incorporates treatment of dynamical space-time, based on the recent work of (Daszuta+, 2021)[49] and (Cook+, 2023)[45].…
We develop a flux-conservative formalism for a Newtonian multi-fluid system, including dissipation and entrainment (i.e. allowing the momentum of one fluid to be a linear combination of the velocities of all fluids). Maximum use is made of…
We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle…
We develop a second-order well-balanced central-upwind scheme for the compressible Euler equations with gravitational source term. Here, we advocate a new paradigm based on a purely conservative reformulation of the equations using global…
We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of…
A conservative finite-volume framework, based on a collocated variable arrangement, for the simulation of flows at all speeds, applicable to incompressible, ideal-gas and real-gas fluids is proposed in conjunction with a fully-coupled…
We describe the implementation of multigrid solvers in the Athena++ adaptive mesh refinement (AMR) framework and their application to the solution of the Poisson equation for self-gravity. The new solvers are built on top of the AMR…
(Abridged) We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics when the light-crossing time is only marginally shorter than dynamical timescales. The transfer…
The present work concerns the derivation of a numerical scheme to approximate weak solutions of the Euler equations with a gravitational source term. The designed scheme is proved to be fully well-balanced since it is able to exactly…
An algorithm for simulating self-gravitating cosmological astrophysical fluids is presented. The advantages include a large dynamic range, parallelizability, high resolution per grid element and fast execution speed. The code is based on a…
We extend the co-scaling formalism of Habegger & Heitsch (2021) implemented in Athena++ to magneto-hydrodynamics. The formalism relies on flow symmetries in astrophysical problems involving expansion, contraction, and center-of-mass motion.…
A relativistic self-gravitating equilibrium system with steady flow as well as spherical symmetry is discovered. The energy-momentum tensor contains the contribution of a current related to the flow and the metric tensor does an…
Numerical relativity is central to the investigation of astrophysical sources in the dynamical and strong-field gravity regime, such as binary black hole and neutron star coalescences. Current challenges set by gravitational-wave and…
Centered numerical fluxes can be constructed for compressible Euler equations which preserve kinetic energy in the semi-discrete finite volume scheme. The essential feature is that the momentum flux should be of the form $f^m_\jph =…