Harmonic initial-boundary evolution in general relativity

Maria C. Babiuc, Béla Szilágyi, and Jeffrey Winicour
Phys. Rev. D 73, 064017 – Published 15 March 2006

Abstract

Computational techniques which establish the stability of an evolution-boundary algorithm for a model wave equation with shift are incorporated into a well-posed version of the initial-boundary value problem for gravitational theory in harmonic coordinates. The resulting algorithm is implemented as a 3-dimensional numerical code which we demonstrate to provide stable, convergent Cauchy evolution in gauge wave and shifted gauge wave testbeds. Code performance is compared for Dirichlet, Neumann, and Sommerfeld boundary conditions and for boundary conditions which explicitly incorporate constraint preservation. The results are used to assess strategies for obtaining physically realistic boundary data by means of Cauchy-characteristic matching.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
7 More
  • Received 10 January 2006

DOI:https://doi.org/10.1103/PhysRevD.73.064017

©2006 American Physical Society

Authors & Affiliations

Maria C. Babiuc1, Béla Szilágyi2, and Jeffrey Winicour1,2

  • 1Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
  • 2Max-Planck-Institut für Gravitationsphysik, Albert-Einstein-Institut, 14476 Golm, Germany

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 73, Iss. 6 — 15 March 2006

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×
×

Images

×

Sign up to receive regular email alerts from Physical Review D

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×