Related papers: SpECTRE: A Task-based Discontinuous Galerkin Code …
Elliptic partial differential equations must be solved numerically for many problems in numerical relativity, such as initial data for every simulation of merging black holes and neutron stars. Existing elliptic solvers can take multiple…
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,…
We have developed the formalism necessary to employ the discontinuous-Galerkin approach in general-relativistic hydrodynamics. The formalism is firstly presented in a general 4-dimensional setting and then specialized to the case of…
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…
Discontinuous Galerkin (DG) methods for solving elliptic equations are gaining popularity in the computational physics community for their high-order spectral convergence and their potential for parallelization on computing clusters.…
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.…
The goal of this paper is to create a fruitful bridge between the numerical methods for approximating partial differential equations (PDEs) in fluid dynamics and the (iterative) numerical methods for dealing with the resulting large linear…
We present a discontinuous Galerkin-finite difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics in dynamical spacetimes. We present several…
A considerable amount of attention has been given to discontinuous Galerkin methods for hyperbolic problems in numerical relativity, showing potential advantages of the methods in dealing with hydrodynamical shocks and other…
Hydrodynamical numerical methods that converge with high-order hold particular promise for astrophysical studies, as they can in principle reach prescribed accuracy goals with higher computational efficiency than standard second- or…
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…
We propose a new parallel Discontinuous Galerkin method for the approximation of hyperbolic systems of conservation laws. The method remains stable with large time steps, while keeping the complexity of an explicit scheme: it does not…
We present and analyze a discontinuous Galerkin method for the numerical modeling of the non-linear fully-coupled thermo-hydro-mechanic problem. We propose a high-order symmetric weighted interior penalty scheme that supports general…
The stable and accurate approximation of discontinuities such as shocks on a finite computational mesh is a challenging task. Detection of shocks or strong discontinuities in the flow solution is typically achieved through a priori troubled…
A novel parallel technique for Fourier-Galerkin pseudo-spectral methods with applications to two-dimensional Navier-Stokes equations and inviscid Boussinesq approximation equations is presented. It takes the advantage of the programming…
In this paper we discuss a new and very efficient implementation of high order accurate ADER discontinuous Galerkin (ADER-DG) finite element schemes on modern massively parallel supercomputers. The numerical methods apply to a very broad…
The recently developed semi-Lagrangian discontinuous Galerkin approach is used to discretize hyperbolic partial differential equations (usually first order equations). Since these methods are conservative, local in space, and able to limit…
We extend the entropy stable high order nodal discontinuous Galerkin spectral element approximation for the non-linear two dimensional shallow water equations presented by Wintermeyer et al. [N. Wintermeyer, A. R. Winters, G. J. Gassner,…
The study of uncertainty propagation poses a great challenge to design numerical solvers with high fidelity. Based on the stochastic Galerkin formulation, this paper addresses the idea and implementation of the first flux reconstruction…
We present a new class of high-order accurate numerical algorithms for solving the equations of general-relativistic ideal magnetohydrodynamics in curved spacetimes. In this paper we assume the background spacetime to be given and static,…