Related papers: Towards high-precision inspiral gravitational wave…
The Harten-Lax-van Leer with contact (HLLC) scheme is known to be plagued by various forms of numerical shock instabilities. In this paper, we propose a new framework for developing shock stable, contact and shear preserving approximate…
We present a new flexible, fast and accurate way to implement massive neutrinos, warm dark matter and any other non-cold dark matter relics in Boltzmann codes. For whatever analytical or numerical form of the phase-space distribution…
The joint observation of the gravitational-wave and electromagnetic signal from the binary neutron-star merger GW170817 allowed for a new independent measurement of the Hubble constant $H_0$, albeit with an uncertainty of about 15\% at…
Accurate and reliable gravitational waveform models are crucial in determining the properties of compact binary mergers. In particular, next-generation gravitational-wave detectors will require more accurate waveforms to avoid biases in the…
Accurate numerical solutions of the equations of hydrodynamics play an ever more important role in many fields of astrophysics. In this work, we reinvestigate the accuracy of the moving-mesh code \textsc{Arepo} and show how its convergence…
Godunov-type methods, which obtain numerical fluxes through local Riemann problems at cell interfaces, are among the most fundamental and widely used numerical methods in computational fluid dynamics. Exact Riemann solvers faithfully solve…
We compare a particular selection of approximate solutions of the Riemann problem in the context of ideal relativistic magnetohydrodynamics. In particular, we focus on Riemann solvers not requiring a full eigenvector structure. Such solvers…
Gravitational waves are produced by orbiting massive binary objects, such as black holes and neutron stars, and propagate as ripples in the very fabric of spacetime. As the waves carry off orbital energy, the two bodies spiral into each…
We compare the results for a set of hydrodynamical tests performed with the AMR finite volume code, MG and the SPH code, SEREN. The test suite includes shock tube tests, with and without cooling, the non-linear thin-shell instability and…
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…
We present initial results of our study of numerical methods for modeling neutron star mergers (NSMs) with simulations that perform the full hydrodynamic evolution required to capture tidal effects, particularly in the last several orbits.…
Gravitational-wave observations of binary neutron star systems can provide information about the masses, spins, and structure of neutron stars. However, this requires accurate and computationally efficient waveform models that take <1s to…
The rich phenomenology of binary neutron star mergers offers a unique opportunity to test general relativity, investigate matter at supranuclear densities, and learn more about the origin of heavy elements. As multi-messenger sources, they…
We present a five-wave Riemann solver for the equations of ideal relativistic magnetohydrodynamics. Our solver can be regarded as a relativistic extension of the five-wave HLLD Riemann solver initially developed by Miyoshi and Kusano for…
Approximate Riemann solvers are widely used for solving hyperbolic conservation laws, including those of magnetohydrodynamics (MHD). However, due to the nonlinearity and complexity of MHD, obtaining accurate and robust numerical solutions…
Gravitational waves from binary neutron star mergers provide critical insights into dense matter physics and strong-field gravity, yet accurate waveform modeling remains computationally intensive. We present a deep generative model for…
We compare results from numerical simulations of spinning binaries in the "orbital hangup" case, where the binary completes at least nine orbits before merger, with post-Newtonian results using the approximants TaylorT1, T4 and Et. We find…
The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate…
Equilibria of binary neutron stars in close circular orbits are computed numerically in a waveless formulation: The full Einstein-relativistic-Euler system is solved on an initial hypersurface to obtain an asymptotically flat form of the…
For the analysis of gravitational-wave signals, fast and accurate gravitational-waveform models are required. These enable us to obtain information on the system properties from compact binary mergers. In this article, we introduce the…