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 ...
28 Geometry of hypersurfaces The non vanishing of the Riemann tensor is reflected by the well-known property that the sum of angles of any triangle drawn at the surface of a sphere is larger than π (cf. Fig. 2.4). The unit vector n normal to Σ (and oriented towards the exterior of the sphere) has the following components with respect to the coordinates (X α ) = (x,y,z): n α x = x2 + y2 + z2 , y x 2 + y 2 + z 2 , It is then easy to compute ∇βn α = ∂n α /∂X β to get ∇βn α = (x 2 + y 2 + z 2 ) −3/2 ⎛ ⎝ z x2 + y2 + z2 . (2.54) y 2 + z 2 −xy −xz −xy x 2 + z 2 −yz −xz −yz x 2 + y 2 The natural basis associated with the coordinates (x i ) = (θ,ϕ) on Σ is ∂θ = (x 2 + y 2 ) −1/2 xz ∂x + yz ∂y − (x 2 + y 2 )∂z ⎞ ⎠. (2.55) (2.56) ∂ϕ = −y ∂x + x∂y. (2.57) The components of the extrinsic curvature tensor in this basis are obtained from Kij = K(∂i,∂j) = −∇βnα (∂i) α (∂j) β . We get Kij = Kθθ Kθϕ Kϕθ Kϕϕ = The trace of K with respect to γ is then −R 0 0 −Rsin 2 θ = − 1 R γij. (2.58) K = − 2 . (2.59) R With these examples, we have encountered hypersurfaces with intrinsic and extrinsic curvature both vanishing (the plane), the intrinsic curvature vanishing but not the extrinsic one (the cylinder), and with both curvatures non vanishing (the sphere). As we shall see in Sec. 2.5, the extrinsic curvature is not fully independent from the intrinsic one: they are related by the Gauss equation. 2.4 Spacelike hypersurface From now on, we focus on spacelike hypersurfaces, i.e. hypersurfaces Σ such that the induced metric γ is definite positive (Riemannian), or equivalently such that the unit normal vector n is timelike (cf. Secs. 2.3.1 and 2.3.2).
2.4.1 The orthogonal projector 2.4 Spacelike hypersurface 29 At each point p ∈ Σ, the space of all spacetime vectors can be orthogonally decomposed as Tp(M) = Tp(Σ) ⊕ Vect(n) , (2.60) where Vect(n) stands for the 1-dimensional subspace of Tp(M) generated by the vector n. Remark : The orthogonal decomposition (2.60) holds for spacelike and timelike hypersurfaces, but not for the null ones. Indeed for any normal n to a null hypersurface Σ, Vect(n) ⊂ Tp(Σ). The orthogonal projector onto Σ is the operator γ associated with the decomposition (2.60) according to γ : Tp(M) −→ Tp(Σ) v ↦−→ v + (n · v)n. In particular, as a direct consequence of n · n = −1, γ satisfies (2.61) γ(n) = 0. (2.62) Besides, it reduces to the identity operator for any vector tangent to Σ: ∀v ∈ Tp(Σ), γ(v) = v. (2.63) According to Eq. (2.61), the components of γ with respect to any basis (eα) of Tp(M) are γ α β = δα β + nα nβ . (2.64) We have noticed in Sec. 2.3.1 that the embedding Φ of Σ in M induces a mapping Tp(Σ) → Tp(M) (push-forward) and a mapping T ∗ p (M) → T ∗ p (Σ) (pull-back), but does not provide any mapping in the reverse ways, i.e. from Tp(M) to Tp(Σ) and from T ∗ p (Σ) to T ∗ p (M). The orthogonal projector naturally provides these reverse mappings: from its very definition, it is a mapping Tp(M) → Tp(Σ) and we can construct from it a mapping γ ∗ M : T ∗ p (Σ) → T ∗ p (M) by setting, for any linear form ω ∈ T ∗ p (Σ), γ ∗ Mω : Tp(M) −→ R v ↦−→ ω(γ(v)). (2.65) This clearly defines a linear form belonging to T ∗ p (M). Obviously, we can extend the operation γ ∗ M to any multilinear form A acting on Tp(Σ), by setting γ ∗ M A : Tp(M) n −→ R (v1,... ,vn) ↦−→ A(γ(v1),... ,γ(vn)) . (2.66) Let us apply this definition to the bilinear form on Σ constituted by the induced metric γ: γ ∗ M γ is then a bilinear form on M, which coincides with γ if its two arguments are vectors tangent to Σ and which gives zero if any of its argument is a vector orthogonal to Σ, i.e. parallel to n.
- 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 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 46 and 47: 46 Geometry of foliations Remark :
- 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
28 Geometry <strong>of</strong> hypersurfaces<br />
The non vanishing <strong>of</strong> the Riemann tensor is reflected by the well-known property that the<br />
sum <strong>of</strong> angles <strong>of</strong> any triangle drawn at the surface <strong>of</strong> a sphere is larger than π (cf. Fig. 2.4).<br />
The unit vector n normal to Σ (<strong>and</strong> oriented towards the exterior <strong>of</strong> the sphere) has the<br />
following components with respect to the coordinates (X α ) = (x,y,z):<br />
n α <br />
x<br />
= <br />
x2 + y2 + z2 ,<br />
y<br />
x 2 + y 2 + z 2 ,<br />
It is then easy to compute ∇βn α = ∂n α /∂X β to get<br />
∇βn α = (x 2 + y 2 + z 2 ) −3/2<br />
⎛<br />
⎝<br />
<br />
z<br />
<br />
x2 + y2 + z2 . (2.54)<br />
y 2 + z 2 −xy −xz<br />
−xy x 2 + z 2 −yz<br />
−xz −yz x 2 + y 2<br />
The natural basis associated with the coordinates (x i ) = (θ,ϕ) on Σ is<br />
∂θ = (x 2 + y 2 ) −1/2 xz ∂x + yz ∂y − (x 2 + y 2 <br />
)∂z<br />
⎞<br />
⎠. (2.55)<br />
(2.56)<br />
∂ϕ = −y ∂x + x∂y. (2.57)<br />
The components <strong>of</strong> the extrinsic curvature tensor in this basis are obtained from Kij =<br />
K(∂i,∂j) = −∇βnα (∂i) α (∂j) β . We get<br />
Kij =<br />
Kθθ Kθϕ<br />
Kϕθ Kϕϕ<br />
<br />
=<br />
The trace <strong>of</strong> K with respect to γ is then<br />
−R 0<br />
0 −Rsin 2 θ<br />
<br />
= − 1<br />
R γij. (2.58)<br />
K = − 2<br />
. (2.59)<br />
R<br />
With these examples, we have encountered hypersurfaces with intrinsic <strong>and</strong> extrinsic curvature<br />
both vanishing (the plane), the intrinsic curvature vanishing but not the extrinsic one (the<br />
cylinder), <strong>and</strong> with both curvatures non vanishing (the sphere). As we shall see in Sec. 2.5, the<br />
extrinsic curvature is not fully independent from the intrinsic one: they are related by the Gauss<br />
equation.<br />
2.4 Spacelike hypersurface<br />
From now on, we focus on spacelike hypersurfaces, i.e. hypersurfaces Σ such that the induced<br />
metric γ is definite positive (Riemannian), or equivalently such that the unit normal vector n<br />
is timelike (cf. Secs. 2.3.1 <strong>and</strong> 2.3.2).