Related papers: Toward computing gravitational initial data withou…
We explore whether a new method to solve the constraints of Einstein's equations, which does not involve elliptic equations, can be applied to provide initial data for black holes. We show that this method can be successfully applied to a…
We present a novel implicit numerical implementation of the parabolic-hyperbolic formulation of the constraints of general relativity. The proposed method is unconditionally stable, has the advantage of not requiring the imposition of any…
We show how the basic idea of parabolic Jacobi relaxation can be modified to obtain a new class of hyperbolic relaxation schemes that are suitable for the solution of elliptic equations. Some of the analytic and numerical properties of…
In this paper we investigate the parabolic-hyperbolic formulation of the vacuum constraint equations introduced by R{\'a}cz with a view to constructing multiple black hole initial data sets without spin. In order to respect the natural…
An algebraic-hyperbolic method for solving the Hamiltonian and momentum constraints has recently been shown to be well posed for general nonlinear perturbations of the initial data for a Schwarzschild black hole. This is a new approach to…
Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudo-spectral code Elliptica, an…
Solving the 4-d Einstein equations as evolution in time requires solving equations of two types: the four elliptic initial data (constraint) equations, followed by the six second order evolution equations. Analytically the constraint…
Numerical solutions of Kepler's Equation are critical components of celestial mechanics software, and are often computation hot spots. This work uses symbolic regression and a genetic learning algorithm to find new initial guesses for…
The parabolic-hyperbolic form of the constraints is integrated numerically. The applied numerical stencil is $4^{th}$ order accurate (in the spatial directions) while 'time'-integration is made by using the method of lines with a $4^{th}$…
We solve the elliptic equations associated with the Hamiltonian and momentum constraints, corresponding to a system composed of two black holes with arbitrary linear and angular momentum. These new solutions are based on a Kerr-Schild…
It was shown recently that the constraints on the initial data for Einstein's equations may be posed as an evolutionary problem [9]. In one of the proposed two methods the constraints can be replaced by a first order symmetrizable…
A procedure for evolving hyperbolic systems of equations on compact computational domains with no boundary conditions was recently described in [arXiv:1905.08657]. In that proposal, the computational grid is expanded in spacelike directions…
We present approximate analytical solutions to the Hamiltonian and momentum constraint equations, corresponding to systems composed of two black holes with arbitrary linear and angular momentum. The analytical nature of these initial data…
We discuss the initial value problem of general relativity in its recently unified Lagrangian and Hamiltonian pictures and present a multi-domain pseudo-spectral collocation method to solve the resulting coupled nonlinear partial…
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…
Near-Kerr black hole initial datasets are constructed by applying either the parabolic-hyperbolic or the algebraic-hyperbolic form of the constraints. In both cases, strongly and weakly asymptotically flat initial datasets with desirable…
By applying a parabolic-hyperbolic formulation of the constraints and superposing Kerr-Schild black holes, a simple method is introduced to initialize time evolution of binary systems. As the input parameters are essentially the same as…
With the goal of taking a step toward the construction of astrophysically realistic initial data for numerical simulations of black holes, we for the first time derive a family of fully general relativistic initial data based on…
We introduce a new approach for finding high accuracy, free and closed-form expressions for the gravitational waves emitted by binary black hole collisions from ab initio models. More precisely, our expressions are built from numerical…
We construct initial data suitable for the Kerr stability conjecture, that is, solutions to the constraint equations on a spacelike hypersurface with boundary entering the black hole horizon that are arbitrarily decaying perturbations of a…