Journal
SIAM JOURNAL ON SCIENTIFIC COMPUTING
Volume 22, Issue 3, Pages 917-950Publisher
SIAM PUBLICATIONS
DOI: 10.1137/S1064827599356171
Keywords
black hole; convolution spline; Einstein equations; preconditioned elliptic system; spherical harmonics
Categories
Ask authors/readers for more resources
We describe algorithms used in our construction of a fourth-order in time evolution for the full Einstein equations and assess the accuracy of some representative solutions. The scheme employs several novel geometric and numerical techniques, including a geometrically invariant coordinate gauge, which leads to a characteristic-transport formulation of the underlying hyperbolic system, combined with a method of lines evolution; convolution splines for radial interpolation, regridding, differentiation, and noise suppression; representations using spin-weighted spherical harmonics; and a spectral preconditioner for solving a class of first-order elliptic systems on S-2. Initial data for the evolution is unconstrained, subject only to a mild size condition. For sample initial data of intermediate strength (19% of the total mass in gravitational energy), the code is accurate to 1 part in 10(5), until null time z = 55m when the coordinate condition breaks down.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available