Related papers: Numerical simulations with a first order BSSN form…
In this second part of our two-part paper, we extend to multiple spatial dimensions the one-dimensional, fully conservative, positivity-preserving, and entropy-bounded discontinuous Galerkin scheme developed in the first part for the…
A new modified Galerkin / Finite Element Method is proposed for the numerical solution of the fully nonlinear shallow water wave equations. The new numerical method allows the use of low-order Lagrange finite element spaces, despite the…
An $hp$-discontinuous Galerkin (DG) method is applied to a class of second order linear hyperbolic integro-differential equations. Based on the analysis of an expanded mixed type Ritz-Volterra projection, {\it a priori} $hp$-error estimates…
We investigate stability issues for steady states of the spherically symmetric Einstein-Vlasov system numerically in Schwarzschild, maximal areal, and Eddington-Finkelstein coordinates. Across all coordinate systems we confirm the…
We present numerical solutions of several spacetimes of physical interest, including binary black hole mergers, in shift-symmetric Einstein-scalar-Gauss-Bonnet (ESGB) gravity, and describe our methods for solving the full equations 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…
We propose and rigorously analyse semi- and fully discrete discontinuous Galerkin methods for an initial and boundary value problem describing inertial viscoelasticity in terms of elastic and viscoelastic stress components, and with mixed…
Recent perturbative studies have shown the existence of long-lived, quasi-stationary configurations of scalar fields around black holes. In particular, such configurations have been found to survive for cosmological timescales, which is a…
BSSN-type evolution equations are discussed. The name refers to the Baumgarte, Shapiro, Shibata, and Nakamura version of the Einstein evolution equations, without introducing the conformal-traceless decomposition but keeping the three…
This thesis is concerned with formulations of the Einstein equations in axisymmetric spacetimes which are suitable for numerical evolutions. We develop two evolution systems based on the (2+1)+1 formalism. The first is a (partially)…
Current spectral simulations of Einstein's equations require writing the equations in first-order form, potentially introducing instabilities and inefficiencies. We present a new penalty method for pseudo-spectral evolutions of second order…
In this manuscript, we put forth a general scheme for defining initial value problems from Einstein's equations of General Relativity constrained by homogeneous and isotropic expansion. The cosmological models arising as solutions are…
The high-order numerical solution of the non-linear shallow water equations (and of hyperbolic systems in general) is susceptible to unphysical Gibbs oscillations that form in the proximity of strong gradients. The solution to this problem…
In this letter we present the first multiparticle solutions to Einstein's field equations in the presence of matter. These solutions are iteratively obtained via the perturbiner method, which can circumvent gravity's infinite number of…
This is the first in a series of articles on the numerical solution of Friedrich's conformal field equations for Einstein's theory of gravity. We will discuss in this paper why one should be interested in applying the conformal method to…
In this paper, a uniformly high-order discontinuous Galerkin gas kinetic scheme (DG-HGKS) is proposed to solve the Euler equations of compressible flows. The new scheme is an extension of the one-stage compact and efficient high-order GKS…
The radiative transfer equation is a fundamental equation in transport theory and applications, which is a 5-dimensional PDE in the stationary one-velocity case, leading to great difficulties in numerical simulation. To tackle this…
We investigate the well-posedness of the characteristic initial-boundary value problem for the Einstein equations in Bondi-like coordinates (including Bondi, double-null and affine). We propose a definition of strong hyperbolicity of a…
The long-time evolution of extreme mass-ratio inspiral systems requires minimal phase and dispersion errors to accurately compute far-field waveforms, while high accuracy is essential near the smaller black hole (modeled as a Dirac delta…
In this paper, we propose a unified and high order accurate fully-discrete one-step ADER Discontinuous Galerkin method for the simulation of linear seismic waves in the sea bottom that are generated by the propagation of free surface water…