Related papers: Construction of high precision numerical single an…
This work consists of two distinct parts. In the first part we present a new method for solving the initial value problem of general relativity. Given any spatial metric with a surface orthogonal Killing field and two freely specified…
We present improvements to construction of binary black hole initial data used in SpEC (the Spectral Einstein Code). We introduce new boundary conditions for the extended conformal thin sandwich elliptic equations that enforce the excision…
First-order fully implicit as well as implicit--explicit schemes for coupled elliptic-parabolic systems are discussed in [Ern and Meunier, ESAIM: M2AN, 2009] and [Altmann et al., Math.\ Comp., 2021], respectively. The extension of the…
We consider the initial data problem for several black holes in vacuum with arbitrary momenta and spins on a three space with punctures. We compactify the internal asymptotically flat regions to obtain a computational domain without inner…
A new method is proposed to improve the numeri- cal simulation of time dependent problems when the initial and boundary data are not compatible. Unlike earlier methods limited to space dimension one, this method can be used for any space…
Systematic numerical investigations of the asymptotics of near Schwarzschild vacuum initial data sets is carried out by inspecting solutions to the parabolic-hyperbolic and to the algebraic-hyperbolic forms of the constraints, respectively.…
We present a new initial data formulation to solve the full set of Einstein equations for spacetimes that contain a black hole under general conditions. The method can be used to construct complete initial data for spacetimes (the full…
We review some of the recent dramatic developments in the fully nonlinear simulation of generic, highly-precessing, black-hole binaries, and introduce a new approach for generating hybrid post-Newtonian / Numerical waveforms for these…
We establish optimal order a priori error estimates for implicit-explicit BDF methods for abstract semilinear parabolic equations with time-dependent operators in a complex Banach space settings, under a sharp condition on the…
The 3+1 Hamiltonian formulation in the gauge $D_tN=-K$ on the lapse function fixes the direction of time associated with the trace $K$ of the extrinsic curvature tensor. The Hamiltonian equations hereby become hyperbolic. We study this new…
We construct initial data corresponding to a single perturbed Kerr black hole in vacuum. These data are defined on specific hyperboloidal ("ACMC-") slices on which the mean extrinsic curvature K asymptotically approaches a constant at…
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…
We present a formalism for constructing quasi-equilibrium binary black hole initial data suitable for numerical evolution. We construct quasi-equilibrium models by imposing an approximate helical Killing symmetry appropriate for…
We study the method for generating the initial data of black hole systems in Gauss-Bonnet (GB) gravity. The initial data are assumed to be momentarily static and conformally flat. Although the equation for the conformal factor is highly…
We present a family of extensions of spherically symmetric Einstein-Lanczos-Lovelock gravity. The field equations are second order and obey a generalized Birkhoff's theorem. The Hamiltonian constraint can be written in terms of a…
When using the black hole exclusion (horizon boundary condition) technique, $K$ is usually nonzero and spatially variable, so none of the special cases of York's conformal-decomposition algorithm apply, and the full 4-vector nonlinear York…
In numerical evolutions of binary black holes (BBH) it is desirable to easily control the orbital eccentricity of the BBH, and the number of orbits completed by the binary through merger. This paper presents fitting formulae that allow to…
We consider unconstrained evolution schemes for the hyperboloidal initial value problem in numerical relativity as a promising candidate for the optimally efficient numerical treatment of radiating compact objects. Here, spherical symmetry…
Simulations of the dynamics generated by partial differential equations (PDEs) provide approximate, numerical solutions to initial value problems. Such simulations are ubiquitous in scientific computing, but the correctness of the results…
We present a new approach for setting initial Cauchy data for multiple black hole spacetimes. The method is based upon adopting an initially Kerr-Schild form of the metric. In the case of non-spinning holes, the constraint equations take a…