Related papers: A double-domain spectral method for black hole exc…
In calculations of the inspiral of binary black holes an intermediate approximation is needed that can bridge the post-Newtonian methods of the early inspiral and the numerical relativity computations of the final plunge. We describe here…
The focus of this work is on the construction of initial data including a neutron star on a hyperboloidal slice. As simplest scenario for this first step, spherical symmetry is considered together with a polytropic-like equation of state…
A new gridding technique for the solution of partial differential equations in cubical geometry is presented. The method is based on volume penalization, allowing for the imposition of a cubical geometry inside of its circumscribing sphere.…
At the beginning of binary black hole simulations, there is a pulse of spurious radiation (or junk radiation) resulting from the initial data not matching astrophysical quasi-equilibrium inspiral exactly. One traditionally waits for the…
Many evolution problems in physics are described by partial differential equations on an infinite domain; therefore, one is interested in the solutions to such problems for a given initial dataset. A prominent example is the binary black…
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…
We consider the 2D quasi-periodic scattering problem in optics, which has been modelled by a boundary value problem governed by Helmholtz equation with transparent boundary conditions. A spectral collocation method and a tensor product…
The standard approach to initial data for both analytic and numerical computations of black hole collisions has been to use conformally-flat initial geometry. Among other advantages, this choice allows the simple superposition of holes with…
The quasi-stationary method for black hole binary inspiral is an approximation for studying strong field effects while suppressing radiation reaction. In this paper we use a nonlinear scalar field toy model (i) to explain the underlying…
In this paper, we propose a numerical method to approximate the solution of partial differential equations in irregular domains with no-flux boundary conditions by means of spectral methods. The main features of this method are its…
We model the electromagnetic signatures of massive black hole binaries (MBHBs) with an associated gas component. The method comprises numerical simulations of relativistic binaries and gas coupled with calculations of the physical…
Given the absence of observations of black hole binaries, it is critical that the full range of accessible parameter space be explored in anticipation of future observation with gravitational wave detectors. To this end, we compile the…
Using the Reduced Basis approach, we efficiently compress and accurately represent the space of waveforms for non-precessing binary black hole inspirals, which constitutes a four dimensional parameter space (two masses, two spin…
We present details of a new numerical code designed to study the formation and evaporation of 2-dimensional black holes within the CGHS model. We explain several elements of the scheme that are crucial to resolve the late-time behavior of…
We describe a grid generation procedure designed to construct new classes of orthogonal coordinate systems for binary black hole spacetimes. The computed coordinates offer an alternative approach to current methods, in addition to providing…
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 construct approximate initial data for non-spinning black hole binary systems by asymptotically matching the 4-metrics of two tidally perturbed Schwarzschild solutions in isotropic coordinates to a resummed post-Newtonian 4-metric in…
We present an alternative approach to setting initial data in general relativity. We do not use a conformal decomposition, but instead express the 3-metric in terms of a given unit vector field and one unknown scalar field. In the case of…
When using black hole excision to numerically evolve a black hole spacetime with no continuous symmetries, most 3+1 finite differencing codes use a Cartesian grid. It's difficult to do excision on such a grid, because the natural $r =…
This paper presents a novel shooting method for solving two-point boundary value problems for second order ordinary differential equations. The method works as follows: first, a guess for the initial condition is made and an integration of…