Related papers: Binary black hole simulation with an adaptive fini…
We study slowly rotating black hole solutions in the Einstein-Bel-Robinson gravity (EBR) in four dimensions. At the leading order in the rotation parameter, the only modification with respect to the static case is the appearance of a…
We present techniques for long-term, stable, and accurate evolutions of multiple-black-hole spacetimes using the `moving puncture' approach with fourth- and eighth-order finite difference stencils. We use these techniques to explore…
This work proposes a unified $hp$-adaptivity framework for hybridized discontinuous Galerkin (HDG) method for a large class of partial differential equations (PDEs) of Friedrichs' type. In particular, we present unified $hp$-HDG…
We analyze the excision strategy for simulating black holes. The problem is modeled by the propagation of quasi-linear waves in a 1-dimensional spatial region with timelike outer boundary, spacelike inner boundary and a horizon in between.…
The subject of this work is an adaptive stochastic Galerkin finite element method for parametric or random elliptic partial differential equations, which generates sparse product polynomial expansions with respect to the parametric…
In this work, we introduce a new space-time variational formulation of the second-order wave equation, where integration by parts is also applied with respect to the time variable, and a modified Hilbert transformation is used. For this…
We solve for single distorted black hole initial data using the puncture method, where the Hamiltonian constraint is written as an elliptic equation in R^3 for the nonsingular part of the metric conformal factor. With this approach we can…
Novel applications of Numerical Relativity demand for more flexible algorithms and tools. In this paper, I develop and test a multigrid solver, based on the infrastructure provided by the Einstein Toolkit, for elliptic partial differential…
In this thesis, we use numerical relativity to investigate gravitational waves from binary black holes in extensions of GR. We first study spherically symmetric gravitational collapse in cubic Horndeski theories of gravity. By varying the…
In this paper we study the finite element approximation of systems of second-order nonlinear hyperbolic equations. The proposed numerical method combines a $hp$-version discontinuous Galerkin finite element approximation in the time…
On the basis of a recently proposed strategy of finite element integration in time domain for partial differential equations with a singular source term, we present a fourth order algorithm for non-rotating black hole perturbations in the…
In this work, we introduce a spectral-infinite element method for solving Einstein's constraint equations in hyperbolic form. As an application of this, we use this method for computing asymptotically flat perturbations of a Kerr black hole…
In response to LIGO's observation of GW170104, we performed a series of full numerical simulations of binary black holes, each designed to replicate likely realizations of its dynamics and radiation. These simulations have been performed at…
We study gravitational collapse of a spherically symmetric scalar field in Einstein-aether theory (general relativity coupled to a dynamical unit timelike vector field). The initial value formulation is developed, and numerical simulations…
Several improvements in numerical methods and gauge choice are presented that make it possible now to perform simulations of the merger and ringdown phases of "generic" binary black-hole evolutions using the pseudo-spectral evolution code…
One of the most promising avenues to perform numerical evolutions in theories beyond General Relativity is the fixing-the-equations approach, a proposal in which new ``driver'' equations are added to the evolution equations in a way that…
An approximate solution to Einstein's equations representing two widely-separated non-rotating black holes in a circular orbit is constructed by matching a post-Newtonian metric to two perturbed Schwarzschild metrics. The spacetime metric…
A homogenization approach is one of effective strategies to solve multiscale elliptic problems approximately. The finite element heterogeneous multiscale method (FEHMM) which is based on the finite element makes possible to simulate such…
Singular solutions of the harmonic Einstein evolution equation are constructed which are related to spatially global and time-local solutions for a certain class of quasilinear hyperbolic systems of second order. The constructed…
This paper presents the first family of conforming finite element divdiv complexes on tetrahedral grids in three dimensions. In these complexes, finite element spaces of $H(\text{divdiv},\Omega;\mathbb{S})$ are from a current preprint [Chen…