Related papers: Black Hole Excision with Multiple Grid Patches
We apply the twistor construction for higher-dimensional black holes to known examples in five space-time dimensions. First the patching matrices are calculated from the explicit metric for these examples. Then an ansatz is proposed for…
Spherically symmetric (1D) black-hole spacetimes are considered as a test for numerical relativity. A finite difference code, based in the hyperbolic structure of Einstein's equations with the harmonic slicing condition is presented.…
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 present results of numerical simulations of the formation of black holes from the gravitational collapse of a massless, minimally-coupled scalar field in 2+1 dimensional, axially-symmetric, anti de-Sitter (AdS) spacetime. The geometry…
It is well known that multigrid methods are optimally efficient for solution of elliptic equations (O(N)), which means that effort is proportional to the number of points at which the solution is evaluated). Thus this is an ideal method to…
Binary black hole simulations become increasingly more computationally expensive with smaller mass ratios, partly because of the longer evolution time, and partly because the lengthscale disparity dictates smaller time steps. The program…
Detailed observations of phenomena involving black holes, be it via gravitational waves or more traditional electromagnetic means, can probe the strong field regime of the gravitational interaction. The prediction of features in such…
A paradigm deeply rooted in modern numerical relativity calculations prescribes the removal of those regions of the computational domain where a physical singularity may develop. We here challenge this paradigm by performing…
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…
We present numerical evolutions of three equal-mass black holes using the moving puncture approach. We calculate puncture initial data for three black holes solving the constraint equations by means of a high-order multigrid elliptic…
Spacetime singularities in numerical relativity can be avoided by excising a region of the computational domain from inside the apparent horizon. We report on results of such a scheme that is based on using ({\it i}) a horizon locking…
A numerical solution scheme for the Einstein field equations based on generalized harmonic coordinates is described, focusing on details not provided before in the literature and that are of particular relevance to the binary black hole…
We have developed a numerical code to study the evolution of self-gravitating matter in dynamic black hole axisymmetric spacetimes in general relativity. The matter fields are evolved with a high-resolution shock-capturing scheme that uses…
We experiment with modifications of the BSSN form of the Einstein field equations (a reformulation of the ADM equations) and demonstrate how these modifications affect the stability of numerical black hole evolution calculations. We use…
We consider the numerical evolution of dynamic black hole initial data sets with a full 3D, nonlinear evolution code. These data sets consist of single black holes distorted by strong gravitational waves, and mimic the late stages of…
We analyze the nonlinear evolution of spherically symmetric wormhole solutions coupled to a massless ghost scalar field using numerical methods. In a previous article we have shown that static wormholes with these properties are unstable…
In this paper a new double-domain spectral method to compute binary black hole excision initial data is presented. The method solves a system of elliptic partial differential equations in the exterior of two excised spheres. At the surface…
Many astrophysical simulations involve extreme dynamic range of timescales around 'special points' in the domain (e.g. black holes, stars, planets, disks, galaxies, shocks, mixing interfaces), where processes on small scales couple strongly…
We have developed a new numerical code to study the evolution of distorted, rotating black holes. We discuss the numerical methods and gauge conditions we developed to evolve such spacetimes. The code has been put through a series of tests,…
Discontinuous Galerkin (DG) methods for solving elliptic equations are gaining popularity in the computational physics community for their high-order spectral convergence and their potential for parallelization on computing clusters.…