Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published August 15, 2007 | public
Journal Article Open

Constraint damping in first-order evolution systems for numerical relativity

Owen, Robert

Abstract

A new constraint suppressing formulation of the Einstein evolution equations is presented, generalizing the five-parameter first-order system due to Kidder, Scheel, and Teukolsky (KST). The auxiliary fields, introduced to make the KST system first order, are given modified evolution equations designed to drive constraint violations toward zero. The algebraic structure of the new system is investigated, showing that the modifications preserve the hyperbolicity of the fundamental and constraint evolution equations. The evolution of the constraints for perturbations of flat spacetime is completely analyzed, and all finite-wavelength constraint modes are shown to decay exponentially when certain adjustable parameters satisfy appropriate inequalities. Numerical simulations of a single Schwarzschild black hole are presented, demonstrating the effectiveness of the new constraint-damping modifications.

Additional Information

©2007 The American Physical Society (Received 29 March 2007; published 23 August 2007) I thank Lee Lindblom for extensive discussions and advice, and also Harald Pfeiffer and Mark Scheel, for help on the manuscript and for assistance on the use of their code. Michael Boyle provided access to some input files and numerical results of Ref. [37], which were of great help in debugging and understanding the numerical issues. The numerical calculations presented in this paper were performed with the Caltech/Cornell Spectral Einstein Code (SpEC) written primarily by Lawrence E. Kidder, Harald P. Pfeiffer, and Mark A. Scheel. This work was supported in part by NSF grants No. PHY-0099568, No. PHY-0244906, No. PHY-0601459, and No. DMS-0553302, NASA grants No. NAG5-12834 and No. NNG05GG52G, and a grant from the Sherman Fairchild Foundation. Some of the numerical calculations leading to this paper were performed with the Tungsten cluster at NCSA.

Files

OWEprd07.pdf
Files (366.9 kB)
Name Size Download all
md5:6903ea617406c5b674e486095ebd2f26
366.9 kB Preview Download

Additional details

Created:
August 22, 2023
Modified:
October 16, 2023