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 ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
52 <strong>3+1</strong> decomposition <strong>of</strong> Einstein equation<br />
Let us assume that the spacetime (M,g) is globally hyperbolic (cf. Sec. 3.2.1) <strong>and</strong> let be<br />
(Σt)t∈R by a foliation <strong>of</strong> M by a family <strong>of</strong> spacelike hypersurfaces. The foundation <strong>of</strong> the <strong>3+1</strong><br />
<strong>formalism</strong> amounts to projecting the Einstein equation (4.1) onto Σt <strong>and</strong> perpendicularly to Σt.<br />
To this purpose let us first consider the <strong>3+1</strong> decomposition <strong>of</strong> the stress-energy tensor.<br />
4.1.2 <strong>3+1</strong> decomposition <strong>of</strong> the stress-energy tensor<br />
From the very definition <strong>of</strong> a stress-energy tensor, the matter energy density as measured<br />
by the Eulerian observer introduced in Sec. 3.3.3 is<br />
E := T(n,n) . (4.3)<br />
This follows from the fact that the 4-velocity <strong>of</strong> the Eulerian observer in the unit normal vector<br />
n.<br />
Similarly, also from the very definition <strong>of</strong> a stress-energy tensor, the matter momentum<br />
density as measured by the Eulerian observer is the linear form<br />
i.e. the linear form defined by<br />
In components:<br />
p := −T(n,γ(.)) , (4.4)<br />
∀v ∈ Tp(M), 〈p,v〉 = −T(n,γ(v)). (4.5)<br />
pα = −Tµν n µ γ ν α. (4.6)<br />
Notice that, thanks to the projector γ, p is a linear form tangent to Σt.<br />
Remark : The momentum density p is <strong>of</strong>ten denoted j. Here we reserve the latter for electric<br />
current density.<br />
Finally, still from the very definition <strong>of</strong> a stress-energy tensor, the matter stress tensor<br />
as measured by the Eulerian observer is the bilinear form<br />
or, in components,<br />
S := γ ∗ T , (4.7)<br />
Sαβ = Tµνγ µ αγ ν β<br />
As for p, S is a tensor field tangent to Σt. Let us recall the physical interpretation <strong>of</strong> the stress<br />
tensor S: given two spacelike unit vectors e <strong>and</strong> e ′ (possibly equal) in the rest frame <strong>of</strong> the<br />
Eulerian observer (i.e. two unit vectors orthogonal to n), S(e,e ′ ) is the force in the direction e<br />
acting on the unit surface whose normal is e ′ . Let us denote by S the trace <strong>of</strong> S with respect<br />
to the metric γ (or equivalently with respect to the metric g):<br />
(4.8)<br />
S := γ ij Sij = g µν Sµν . (4.9)