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 ...

luth.obspm.fr
from luth.obspm.fr More from this publisher
20.07.2013 Views

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.

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).

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!