Related papers: Binary neutron star mergers using a discontinuous …
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.…
Simulations of relativistic hydrodynamics often need both high accuracy and robust shock-handling properties. The discontinuous Galerkin method combines these features --- a high order of convergence in regions where the solution is smooth…
Binary black holes are the most abundant source of gravitational-wave observations. Gravitational-wave observatories in the next decade will require tremendous increases in the accuracy of numerical waveforms modeling binary black holes,…
Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to…
We introduce a new hybridized discontinuous Galerkin method for the incompressible magnetohydrodynamics equations. If particular velocity, pressure, magnetic field, and magnetic pressure spaces are employed for both element and trace…
We develop a hybrid spatial discretization for the wave equation in second order form, based on high-order accurate finite difference methods and discontinuous Galerkin methods. The hybridization combines computational efficiency of finite…
A hybridized discontinuous Galerkin method is proposed for solving 2D fractional convection-diffusion equations containing derivatives of fractional order in space on a finite domain. The Riemann-Liouville derivative is used for the spatial…
We present here a new hybrid scheme that combines a discontinuous Galerkin (DG) method with compact finite volume (FV) and finite difference (FD) methods. The computational mesh is divided into smaller elements that touch but do not…
We introduce and analyze a space-time hybridized discontinuous Galerkin method for the evolutionary Navier--Stokes equations. Key features of the numerical scheme include point-wise mass conservation, energy stability, and pressure…
We introduce a hybrid method to couple continuous Galerkin finite element methods and high-order finite difference methods in a nonconforming multiblock fashion. The aim is to optimize computational efficiency when complex geometries are…
Interpreting gravitational wave observations and understanding the physics of astrophysical compact objects such as black holes or neutron stars requires accurate theoretical models. Here, we present a new numerical relativity computer…
We report simulations of the inspiral and merger of binary neutron stars performed with \texttt{WhiskyTHC}, the first of a new generation of numerical relativity codes employing higher than second-order methods for both the spacetime and…
The detection of GW170817, together with its electromagnetic counterparts, has proven that binary neutron star mergers are of central importance to the field of nuclear astrophysics, e.g., through a better understanding of the formation of…
Relativistic plasmas around compact objects can sometimes be approximated as being force-free. In this limit, the plasma inertia is negligible and the overall dynamics is governed by global electric currents. We present a novel numerical…
The merger of two neutron stars has been proposed as a source of gamma-ray bursts, r-process elements, and detectable gravitational waves. Extracting information from observations of these phenomena requires fully relativistic simulations.…
We develop a high order accurate numerical method for solving the elastic wave equation in second-order form. We hybridize the computationally efficient Cartesian grid formulation of finite differences with geometrically flexible…
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…
The mixed form of the Cahn-Hilliard equations is discretized by the hybridizable discontinuous Galerkin method. For any chemical energy density, existence and uniqueness of the numerical solution is obtained. The scheme is proved to be…
We present a new numerical relativity code designed for simulations of compact binaries involving matter. The code is an upgrade of the BAM code to include general relativistic hydrodynamics and implements state-of-the-art…
We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics…