Related papers: Performance-Portable Binary Neutron Star Mergers w…
We present the first neutron star merger simulations performed with the newly developed Numerical Relativity code SPHINCS_BSSN. This code evolves the spacetime on a mesh using the BSSN formulation, but matter is evolved via Lagrangian…
We present a code for solving the coupled Einstein-hydrodynamics equations to evolve relativistic, self-gravitating fluids. The Einstein field equations are solved in generalized harmonic coordinates on one grid using pseudospectral…
(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…
Numerical simulations are essential to understand the complex physics accompanying the merger of binary systems of neutron stars. However, these simulations become computationally challenging when they have to model the merger remnants on…
Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved…
We aim to study the possibility of a hadron-quark phase transition in the interior of neutron stars, taking into account different schematic evolutionary stages at finite temperature. We also discuss the strange quark matter stability in…
The force-free limit of magnetohydrodynamics (MHD) is often a reasonable approximation to model black hole and neutron star magnetospheres. We describe a general relativistic force-free (GRFFE) formulation that allows general relativistic…
Neutron stars and black holes in X-ray binaries are observed to host strong collimated jets in the hard spectral state. Numerical simulations can act as a valuable tool in understanding the mechanisms behind jet formation and its…
We present the first results on the construction of an exascale hyperbolic PDE engine (ExaHyPE), a code for the next generation of supercomputers with the objective to evolve dynamical spacetimes of black holes, neutron stars and binaries.…
We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of merging binary neutron stars incorporating neutrino transport and magnetic fields. Our new radiative transport module for neutrinos adopts a general relativistic,…
We present an extensive study of the effects of neutrino transport in 3-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino…
We describe a high-order numerical magnetohydrodynamics (MHD) solver built upon a novel non-linear entropy stable numerical flux function that supports eight travelling wave solutions. By construction the solver conserves mass, momentum,…
The numerical study of relativistic magnetohydrodynamics (MHD) plays a crucial role in high-energy astrophysics, but unfortunately is computationally demanding, given the complex physics involved (high Lorentz factor flows, extreme…
Large scale simulations are a key pillar of modern research and require ever-increasing computational resources. Different novel manycore architectures have emerged in recent years on the way towards the exascale era. Performance…
In this work we present a proof of concept of CUDA-capable, resistive, multi-fluid models of relativistic magnetohydrodynamics (RMHD). Resistive and multi-fluid codes for simulating models of RMHD suffer from stiff source terms, so it is…
The merger of binary neutron stars (BNSs) can lead to large amplifications of the magnetic field due to the development of turbulence and instabilities in the fluid, such as the Kelvin-Helmholtz shear instability, which drive small-scale…
We perform magnetohydrodynamic simulations in full general relativity (GRMHD) of quasi-circular, equal-mass, binary neutron stars that undergo merger. The initial stars are irrotational, $n=1$ polytropes and are magnetized. We explore two…
Black hole-neutron star (BHNS) binaries are expected to be among the leading sources of gravitational waves observable by ground-based detectors, and may be the progenitors of short-hard gamma ray bursts (SGRBs) as well. Here, we discuss…
We solve Einstein's field equations coupled to relativistic hydrodynamics in full 3+1 general relativity to evolve astrophysical systems characterized by strong gravitational fields. We model rotating, collapsing and binary stars by…
We present a discontinuous Galerkin-finite-difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple.…