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 ...
46 Geometry of foliations Remark : In many numerical relativity articles, Eq. (3.27) is used to define the extrinsic curvature tensor of the hypersurface Σt. It is worth to keep in mind that this equation has a meaning only because Σt is member of a foliation. Indeed the right-hand side is the derivative of the induced metric in a direction which is not parallel to the hypersurface and therefore this quantity could not be defined for a single hypersurface, as considered in Chap. 2. 3.3.6 Evolution of the orthogonal projector Let us now evaluate the Lie derivative of the orthogonal projector onto Σt along the normal evolution vector. Using Eqs. (A.8) and (3.22), we have i.e. Lmγ α β = mµ ∇µγ α β − γµ β∇µm α + γ α µ = Nn µ ∇µ(n α nβ) + γ µ β −γ α µ = N( n µ ∇µn α =N −1 D α N µ ∇βm α NK µ + D α N nµ − n α ∇µN NK µ β + Dµ N nβ − n µ ∇βN nβ + n α n µ ∇µnβ =N −1 ) + NK DβN α β − nαDβN − NK α β − DαN nβ = 0, (3.30) Lmγ = 0 . (3.31) An important consequence of this is that the Lie derivative along m of any tensor field T tangent to Σt is a tensor field tangent to Σt: T tangent to Σt =⇒ LmT tangent to Σt . (3.32) Indeed a distinctive feature of a tensor field tangent to Σt is γ ∗ T = T. (3.33) Assume for instance that T is a tensor field of type 1 1 . Then the above equation writes [cf. Eq. (2.71)] γ α µ γνβT µ ν = T α β . (3.34) Taking the Lie derivative along m of this relation, employing the Leibniz rule and making use of Eq. (3.31), leads to α Lm γ µγ ν βT µ ν = LmT α β Lmγ α µ γ =0 ν βT µ ν + γ α µ Lmγ ν β T =0 µ ν + γ α µγ ν β LmT µ ν = LmT α β γ ∗ LmT = LmT. (3.35) This shows that LmT is tangent to Σt. The proof is readily extended to any type of tensor field tangent to Σt. Notice that the property (3.32) generalizes that obtained for vectors in Sec. 3.3.2 [cf. Eq. (3.13)].
3.4 Last part of the 3+1 decomposition of the Riemann tensor 47 Remark : An illustration of property (3.32) is provided by Eq. (3.24), which says that Lmγ is −2NK: K being tangent to Σt, we have immediately that Lmγ is tangent to Σt. Remark : Contrary to Ln γ and Lmγ, which are related by Eq. (3.26), Ln γ and Lmγ are not proportional. Indeed a calculation similar to that which lead to Eq. (3.26) gives Therefore the property Lmγ = 0 implies Ln γ = 1 N Lmγ + n ⊗ D ln N. (3.36) Ln γ = n ⊗ D ln N = 0. (3.37) Hence the privileged role played by m regarding the evolution of the hypersurfaces Σt is not shared by n; this merely reflects that the hypersurfaces are Lie dragged by m, not by n. 3.4 Last part of the 3+1 decomposition of the Riemann tensor 3.4.1 Last non trivial projection of the spacetime Riemann tensor In Chap. 2, we have formed the fully projected part of the spacetime Riemann tensor, i.e. γ ∗ 4 Riem, yielding the Gauss equation [Eq. (2.92)], as well as the part projected three times onto Σt and once along the normal n, yielding the Codazzi equation [Eq. (2.101)]. These two decompositions involve only fields tangents to Σt and their derivatives in directions parallel to Σt, namely γ, K, Riem and DK. This is why they could be defined for a single hypersurface. In the present section, we form the projection of the spacetime Riemann tensor twice onto Σt and twice along n. As we shall see, this involves a derivative in the direction normal to the hypersurface. As for the Codazzi equation, the starting point of the calculation is the Ricci identity applied to the vector n, i.e. Eq. (2.98). But instead of projecting it totally onto Σt, let us project it only twice onto Σt and once along n: γαµn σ γ ν β (∇ν∇σn µ − ∇σ∇νn µ ) = γαµn σ γ ν β 4 R µ ρνσn ρ . (3.38) By means of Eq. (3.20), we get successively γαµ n ρ γ ν β nσ 4 R µ ρνσ = γαµn σ γ ν β [−∇ν(K µ σ + D µ ln N nσ) + ∇σ(K µ ν + D µ ln N nν)] = γαµn σ γ ν β [ − ∇νK µ σ − ∇νnσ D µ ln N − nσ∇νD µ ln N +∇σK µ ν + ∇σnν D µ ln N + nν∇σD µ ln N ] = γαµγ ν β [Kµ σ ∇νn σ + ∇νD µ ln N + n σ ∇σK µ ν + Dν ln N D µ ln N] = −KασK σ β + DβDα ln N + γ µ αγ ν β nσ ∇σKµν + Dα ln NDβ ln N = −KασK σ β + 1 N DβDαN + γ µ α γν β nσ ∇σKµν. (3.39)
- Page 1: arXiv:gr-qc/0703035v1 6 Mar 2007 3+
- Page 4 and 5: 4 CONTENTS 3.3.6 Evolution of the o
- Page 6 and 7: 6 CONTENTS 8 The initial data probl
- Page 8 and 9: 8 CONTENTS
- Page 10 and 11: 10 CONTENTS
- Page 12 and 13: 12 Introduction 3D (no symmetry at
- Page 14 and 15: 14 Introduction
- Page 16 and 17: 16 Geometry of hypersurfaces in two
- Page 18 and 19: 18 Geometry of hypersurfaces 2.2.3
- Page 20 and 21: 20 Geometry of hypersurfaces Figure
- Page 22 and 23: 22 Geometry of hypersurfaces •
- Page 24 and 25: 24 Geometry of hypersurfaces The ei
- Page 26 and 27: 26 Geometry of hypersurfaces Figure
- Page 28 and 29: 28 Geometry of hypersurfaces The no
- Page 30 and 31: 30 Geometry of hypersurfaces Since
- Page 32 and 33: 32 Geometry of hypersurfaces 2.4.3
- Page 34 and 35: 34 Geometry of hypersurfaces 2.5 Ga
- Page 36 and 37: 36 Geometry of hypersurfaces Exampl
- Page 38 and 39: 38 Geometry of hypersurfaces
- Page 40 and 41: 40 Geometry of foliations Figure 3.
- Page 42 and 43: 42 Geometry of foliations 3.3.2 Nor
- Page 44 and 45: 44 Geometry of foliations means Eq.
- Page 48 and 49: 48 Geometry of foliations Note that
- Page 50 and 51: 50 Geometry of foliations
- Page 52 and 53: 52 3+1 decomposition of Einstein eq
- Page 54 and 55: 54 3+1 decomposition of Einstein eq
- Page 56 and 57: 56 3+1 decomposition of Einstein eq
- Page 58 and 59: 58 3+1 decomposition of Einstein eq
- Page 60 and 61: 60 3+1 decomposition of Einstein eq
- Page 62 and 63: 62 3+1 decomposition of Einstein eq
- Page 64 and 65: 64 3+1 decomposition of Einstein eq
- Page 66 and 67: 66 3+1 decomposition of Einstein eq
- Page 68 and 69: 68 3+1 decomposition of Einstein eq
- Page 70 and 71: 70 3+1 decomposition of Einstein eq
- Page 72 and 73: 72 3+1 equations for matter and ele
- Page 74 and 75: 74 3+1 equations for matter and ele
- Page 76 and 77: 76 3+1 equations for matter and ele
- Page 78 and 79: 78 3+1 equations for matter and ele
- Page 80 and 81: 80 3+1 equations for matter and ele
- Page 82 and 83: 82 3+1 equations for matter and ele
- Page 84 and 85: 84 Conformal decomposition equivale
- Page 86 and 87: 86 Conformal decomposition As an ex
- Page 88 and 89: 88 Conformal decomposition 6.2.4 Co
- Page 90 and 91: 90 Conformal decomposition 6.3.1 Ge
- Page 92 and 93: 92 Conformal decomposition where K
- Page 94 and 95: 94 Conformal decomposition to write
46 Geometry <strong>of</strong> foliations<br />
Remark : In many <strong>numerical</strong> <strong>relativity</strong> articles, Eq. (3.27) is used to define the extrinsic curvature<br />
tensor <strong>of</strong> the hypersurface Σt. It is worth to keep in mind that this equation has<br />
a meaning only because Σt is member <strong>of</strong> a foliation. Indeed the right-h<strong>and</strong> side is the<br />
derivative <strong>of</strong> the induced metric in a direction which is not parallel to the hypersurface<br />
<strong>and</strong> therefore this quantity could not be defined for a single hypersurface, as considered in<br />
Chap. 2.<br />
3.3.6 Evolution <strong>of</strong> the orthogonal projector<br />
Let us now evaluate the Lie derivative <strong>of</strong> the orthogonal projector onto Σt along the normal<br />
evolution vector. Using Eqs. (A.8) <strong>and</strong> (3.22), we have<br />
i.e.<br />
Lmγ α β = mµ ∇µγ α β − γµ<br />
β∇µm α + γ α µ<br />
= Nn µ ∇µ(n α nβ) + γ µ<br />
β<br />
−γ α µ<br />
= N( n µ ∇µn α<br />
<br />
=N −1 D α N<br />
µ<br />
∇βm<br />
<br />
α<br />
NK µ + D α N nµ − n α ∇µN <br />
<br />
NK µ<br />
β + Dµ N nβ − n µ <br />
∇βN<br />
nβ + n α n µ ∇µnβ<br />
<br />
=N −1 ) + NK<br />
DβN<br />
α β − nαDβN − NK α β − DαN nβ<br />
= 0, (3.30)<br />
Lmγ = 0 . (3.31)<br />
An important consequence <strong>of</strong> this is that the Lie derivative along m <strong>of</strong> any tensor field T tangent<br />
to Σt is a tensor field tangent to Σt:<br />
T tangent to Σt =⇒ LmT tangent to Σt . (3.32)<br />
Indeed a distinctive feature <strong>of</strong> a tensor field tangent to Σt is<br />
γ ∗ T = T. (3.33)<br />
Assume for instance that T is a tensor field <strong>of</strong> type 1 1 . Then the above equation writes [cf.<br />
Eq. (2.71)]<br />
γ α µ γνβT µ ν = T α β . (3.34)<br />
Taking the Lie derivative along m <strong>of</strong> this relation, employing the Leibniz rule <strong>and</strong> making use<br />
<strong>of</strong> Eq. (3.31), leads to<br />
α<br />
Lm γ µγ ν βT µ <br />
ν = LmT α β<br />
Lmγ α µ γ<br />
<br />
=0<br />
ν βT µ ν + γ α µ Lmγ ν β T<br />
<br />
=0<br />
µ ν + γ α µγ ν β LmT µ ν = LmT α β<br />
γ ∗ LmT = LmT. (3.35)<br />
This shows that LmT is tangent to Σt. The pro<strong>of</strong> is readily extended to any type <strong>of</strong> tensor field<br />
tangent to Σt. Notice that the property (3.32) generalizes that obtained for vectors in Sec. 3.3.2<br />
[cf. Eq. (3.13)].