Bypassing stability conditions and curing logarithmic singularities
By Jörg Frauendiener and Jörg Hennig
Assume you want to model a general relativistic spacetime. Due to the annoying limitations of conventional computers, like finite memory and processing speed, it is tempting to focus on a finite portion of the spacetime. Then, without waiting endlessly, one can obtain an approximate description of this portion. One just has to choose a suitable numerical method and solve the field equations for the metric at some set of grid-points. While this approach is standard, it introduces unpleasant problems. Firstly, the set of equations needs to be complemented with boundary conditions at the outer edges of this finite portion, in order to obtain a complete mathematical problem. This, however, is quite unphysical as usually no information about the actual behaviour at such an artificial boundary is available. Consequently, spurious gravitational radiation enters the numerical domain. Secondly, if one is interested in accurately describing gravitational waves, one should recall that these are only well-defined at infinity. Hence it is desirable to extend the simulation up to infinity.