In a single process, the merger of binary neutron star systems combines extreme gravity, the
copious emission of gravitational waves, complex microphysics and electromagnetic …

Binary neutron star mergers observations are a unique way to constrain fundamental
physics and astrophysics at the extreme. The interpretation of gravitational-wave events and …

In 1917, the British scientist LF Richardson made the first reported attempt to predict the
weather by solving partial differential equations numerically, by hand! It is generally …

This paper describes the open-source code Enzo, which uses block-structured adaptive
mesh refinement to provide high spatial and temporal resolution for modeling astrophysical …

We present a new numerical code, PLUTO, for the solution of hypersonic flows in 1, 2, and 3
spatial dimensions and different systems of coordinates. The code provides a multiphysics …

In simulations of magnetohydrodynamic (MHD) processes the violation of the divergence
constraint causes severe stability problems. In this paper we develop and test a new …

Context. Computational astrophysics nowadays routinely combines grid-adaptive
capabilities with modern shock-capturing, high resolution spatio-temporal integration …

Seven schemes to maintain the∇· B= 0 constraint numerically are compared. All these
algorithms can be combined with shock-capturing Godunov type base schemes. They fall …

In this article, a conservative least-squares polynomial reconstruction operator is applied to
the discontinuous Galerkin method. In a first instance, piecewise polynomials of degree N …

We describe a single step, second-order accurate Godunov scheme for ideal MHD based on
combining the piecewise parabolic method (PPM) for performing spatial reconstruction, the …