3+1 formalism and bases of numerical relativity - LUTh ...
3+1 formalism and bases of numerical relativity - LUTh ...
3+1 formalism and bases of numerical relativity - LUTh ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
156 Choice <strong>of</strong> foliation <strong>and</strong> spatial coordinates<br />
Figure 9.3: Kruskal-Szekeres diagram showing the maximal slicing <strong>of</strong> Schwarzschild spacetime defined by the<br />
st<strong>and</strong>ard Schwarzschild time coordinate t. As for Fig. 9.1, R st<strong>and</strong>s for Schwarzschild radial coordinate (areal<br />
radius), so that R = 0 is the singularity <strong>and</strong> R = 2m is the event horizon, whereas r st<strong>and</strong>s for the isotropic<br />
radial coordinate [cf. Eq. (8.118)].<br />
Besides its nice geometrical definition, an interesting property <strong>of</strong> maximal slicing is the<br />
singularity avoidance. This is related to the fact that the set <strong>of</strong> the Eulerian observers <strong>of</strong><br />
a maximal foliation define an incompressible flow: indeed, thanks to Eq. (2.77), the condition<br />
K = 0 is equivalent to the incompressibility condition<br />
∇ · n = 0 (9.17)<br />
for the 4-velocity field n <strong>of</strong> the Eulerian observers. If we compare with the Eulerian observers<br />
<strong>of</strong> geodesic slicings (Sec. 9.2.1), who have the tendency to squeeze, we may say that maximalslicing<br />
Eulerian observers do not converge because they are accelerating (DN = 0) in order<br />
to balance the focusing effect <strong>of</strong> gravity. Loosely speaking, the incompressibility prevents the<br />
Eulerian observers from converging towards the central singularity if the latter forms during the<br />
time evolution. This is illustrated by the following example in Schwarzschild spacetime.<br />
Example : Let us consider the time development <strong>of</strong> the initial data constructed in Secs. 8.2.5<br />
<strong>and</strong> 8.3.3, namely the momentarily static slice tS = 0 <strong>of</strong> Schwarzschild spacetime (with the<br />
Einstein-Rosen bridge). A first maximal slicing development <strong>of</strong> these initial data is that<br />
based on Schwarzschild time coordinate tS <strong>and</strong> discussed above (Fig. 9.3). The corresponding<br />
lapse function is given by Eq. (9.16) <strong>and</strong> is antisymmetric about the minimal surface<br />
r = m/2 (throat). There exists a second maximal-slicing development <strong>of</strong> the same initial<br />
data but with a lapse which is symmetric about the throat. It has been found in 1973 by<br />
Estabrook, Wahlquist, Christensen, DeWitt, Smarr <strong>and</strong> Tsiang [118], as well as Reinhart