Abstract:We investigate how the accuracy and stability of numerical relativity simulations of 1D colliding plane waves depends on choices of equation formulations, gauge conditions, boundary conditions, and numerical methods, all in the context of a first-order 3+1 approach to the Einstein equations, with basic variables some combination of first derivatives of the spatial metric and components of the extrinsic curvature tensor. Hyperbolic schemes, specifically variations on schemes proposed by Bona and Masso and Anderson and York, are compared with variations of the Arnowitt-Deser-Misner formulation. Modifications of the three basic schemes include raising one index in the metric derivative and extrinsic curvature variables and adding a multiple of the energy constraint to the extrinsic curvature evolution equations. Redundant variables in the Bona-Masso formulation may be reset frequently or allowed to evolve freely. Gauge conditions which simplify the dynamical structure of the system are imposed during each time step, but the lapse and shift are reset periodically to control the evolution of the spacetime slicing and the longitudinal part of the metric. We show that physically correct boundary conditions, satisfying the energy and momentum constraint equations, generically require the presence of some ingoing eigenmodes of the characteristic matrix. Numerical methods are developed for the hyperbolic systems based on decomposing flux differences into linear combinations of eigenvectors of the characteristic matrix. These methods are shown to be second-order accurate, and in practice second-order convergent, for smooth solutions, even when the eigenvectors and eigenvalues of the characteristic matrix are spatially varying.

Chiral solitons from dimensional reduction of Chern-Simons gauged non-linear Schrödinger model of FQHE: classical and quantum aspects

Black Hole Excision for Dynamic Black Holes

The Evolution of Distorted Rotating Black Holes I: Methods and Tests

Foundations of multiple black hole evolutions

Initial Data and Eccentricity Reduction Toolkit for Binary Black Hole Numerical Relativity Waveforms

Advances in Simulations of Generic Black-Hole Binaries

Numerical-relativity simulations of the quasi-circular inspiral and merger of non-spinning, charged black holes: methods and comparison with approximate approaches

Simulating binary black hole mergers using discontinuous Galerkin methods

Nonlinear ringdown at the black hole horizon

Numerical modeling of black holes as sources of gravitational waves in a nutshell

Improved methods for simulating nearly extremal binary black holes

Non-linear evolutions of magnetised thick discs around black holes: dependence on the initial data

Inspiraling black-hole binary spacetimes: Challenges in transitioning from analytical to numerical techniques

Numerical Relativity Multimodal Waveforms using Absorbing Boundary Conditions

Numerical Tests of Evolution Systems, Gauge Conditions, and Boundary Conditions for 1D Colliding Gravitational Plane Waves

Recent Advances in the Numerical Simulations of Binary Black Holes

Late-time tails in nonlinear evolutions of merging black holes

Eccentric Binary Black Hole Simulations with Numerical Relativity

Towards adiabatic waveforms for inspiral into Kerr black holes: II. Dynamical sources and generic orbits

Second MAYA Catalog of Binary Black Hole Numerical Relativity Waveforms

An eccentric binary blackhole in post-Newtonian theory