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

196 Conformal Killing operator and conformal vector Laplacian Replace each occurrence of ∂/∂x j by the component ξj of the linear form ξ, thereby obtaining a mapping which is no longer differential, i.e. that involves only values of the fields at the point p; this is the principal symbol of ˜∆L at p with respect to ξ: P (p,ξ) : Tp(Σ) −→ Tp(Σ) v = (v i ) ↦−→ P (p,ξ)(v) = ˜γ jk (p)ξjξk v i + 1 3 ˜γik (p)ξkξjv j , (B.13) The differential operator ˜∆L is said to be elliptic on Σ iff the principal symbol P (p,ξ) is an isomorphism for every p ∈ Σ and every non-vanishing linear form ξ ∈ T ∗ p (Σ). It is said to be strongly elliptic if all the eigenvalues of P (p,ξ) are non-vanishing and have the same sign. To check whether it is the case, let us consider the bilinear form ˜P (p,ξ) associated to the endomorphism P (p,ξ) by the conformal metric: ∀(v,w) ∈ Tp(Σ) 2 , ˜P (p,ξ)(v,w) = ˜γ v, P (p,ξ)(w) . (B.14) Its matrix ˜ Pij is deduced from the matrix P i j of P (p,ξ) by lowering the index i with ˜γ(p). We get ˜Pij = ˜γ kl (p)ξkξl ˜γij(p) + 1 3 ξiξj. (B.15) Hence ˜P (p,ξ) is clearly a symmetric bilinear form. Moreover it is positive definite for ξ = 0: for any vector v ∈ Tp(Σ) such that v = 0, we have ˜P (p,ξ)(v,v) = ˜γ kl (p)ξkξl ˜γij(p)v i v j + 1 3 (ξiv i ) 2 > 0, (B.16) where the > 0 follows from the positive definite character of ˜γ. ˜P (p,ξ) being positive definite symmetric bilinear form, we conclude that P (p,ξ) is an isomorphism and that all its eigenvalues are real and strictly positive. Therefore ˜∆L is a strongly elliptic operator. B.2.3 Kernel Let us now determine the kernel of ˜∆L. Clearly this kernel contains the kernel of the conformal Killing operator ˜L. Actually it is not larger than that kernel: ker ˜∆L = ker ˜L . (B.17) Let us establish this property. For any vector field v ∈ T (Σ), we have ˜γijv Σ i ∆L ˜ v j ˜γ d 3 x = = = ˜γijv Σ i Dl( ˜ ˜ Lv) jl ˜γ d 3 x Σ ∂Σ ˜Dl ˜γijv i ( ˜ Lv) jl ˜γijv i ( ˜ Lv) jl 2 ˜sl ˜q d y − − ˜γij ˜ Dlv i ( ˜ Lv) jl ˜γ d 3 x ˜γij Σ ˜ Dlv i ( ˜ Lv) jl ˜γ d 3 x, (B.18) where the Gauss-Ostrogradsky theorem has been used to get the last line. We shall consider two situations for (Σ,γ):

B.2 Conformal vector Laplacian 197 • Σ is a closed manifold, i.e. is compact without boundary; • (Σ, ˜γ) is an asymptotically flat manifold, in the sense made precise in Sec. 7.2. In the former case the lack of boundary of Σ implies that the first integral in the right-hand side of Eq. (B.18) is zero. In the latter case, we will restrict our attention to vectors v which decay at spatial infinity according to (cf. Sec. 7.2) v i = O(r −1 ) (B.19) ∂v i ∂x j = O(r−2 ), (B.20) where the components are to be taken with respect to the asymptotically Cartesian coordinate system (x i ) introduced in Sec. 7.2. The behavior (B.19)-(B.20) implies v i ( ˜ Lv) jl = O(r −3 ), (B.21) so that the surface integral in Eq. (B.18) vanishes. So for both cases of Σ closed or asymptotically flat, Eq. (B.18) reduces to ˜γijv Σ i ∆L ˜ v j ˜γ d 3 x = − In view of the right-hand side integrand, let us evaluate ˜γij Σ ˜ Dlv i ( ˜ Lv) jl ˜γ d 3 x. (B.22) ˜γij˜γkl( ˜ Lv) ik ( ˜ Lv) jl = ˜γij˜γkl( ˜ D i v k + ˜ D k v i )( ˜ Lv) jl − 2 3 ˜ Dmv m ˜γ ik ˜γij = =δ k j ˜γkl ˜ Djv k + ˜γij ˜ Dlv i ( ˜ Lv) jl − 2 3 ˜ Dmv m ˜γjl( ˜ Lv) jl =0 ˜γkl( ˜ Lv) jl = 2˜γij ˜ Dlv i ( ˜ Lv) jl , (B.23) where we have used the symmetry and the traceless property of ( ˜ Lv) jl to get the last line. Hence Eq. (B.22) becomes Σ ˜γijv i ∆L ˜ v j ˜γ d 3 x = − 1 ˜γij˜γkl( 2 Σ ˜ Lv) ik ( ˜ Lv) jl ˜γ d 3 x. (B.24) Let us assume now that v ∈ ker ˜∆L: ˜ ∆L v j = 0. Then the left-hand side of the above equation vanishes, leaving ˜γij˜γkl( Σ ˜ Lv) ik ( ˜ Lv) jl ˜γ d 3 x = 0. (B.25) Since ˜γ is a positive definite metric, we conclude that ( ˜ Lv) ij = 0, i.e. that v ∈ ker ˜L. This demonstrates property (B.17). Hence the “harmonic functions” of the conformal vector Laplacian ˜∆L are nothing but the conformal Killing vectors (one should add “which vanish at spatial infinity as (B.19)-(B.20)” in the case of an asymptotically flat space).

196 Conformal Killing operator <strong>and</strong> conformal vector Laplacian<br />

Replace each occurrence <strong>of</strong> ∂/∂x j by the component ξj <strong>of</strong> the linear form ξ, thereby obtaining<br />

a mapping which is no longer differential, i.e. that involves only values <strong>of</strong> the fields at the point<br />

p; this is the principal symbol <strong>of</strong> ˜∆L at p with respect to ξ:<br />

P (p,ξ) : Tp(Σ) −→ Tp(Σ)<br />

v = (v i ) ↦−→ P (p,ξ)(v) =<br />

<br />

˜γ jk (p)ξjξk v i + 1<br />

3 ˜γik (p)ξkξjv j<br />

<br />

,<br />

(B.13)<br />

The differential operator ˜∆L is said to be elliptic on Σ iff the principal symbol P (p,ξ) is an<br />

isomorphism for every p ∈ Σ <strong>and</strong> every non-vanishing linear form ξ ∈ T ∗<br />

p (Σ). It is said to be<br />

strongly elliptic if all the eigenvalues <strong>of</strong> P (p,ξ) are non-vanishing <strong>and</strong> have the same sign. To<br />

check whether it is the case, let us consider the bilinear form ˜P (p,ξ) associated to the endomorphism<br />

P (p,ξ) by the conformal metric:<br />

∀(v,w) ∈ Tp(Σ) 2 , ˜P (p,ξ)(v,w) = ˜γ v, P (p,ξ)(w) . (B.14)<br />

Its matrix ˜ Pij is deduced from the matrix P i j <strong>of</strong> P (p,ξ) by lowering the index i with ˜γ(p). We<br />

get<br />

˜Pij = ˜γ kl (p)ξkξl ˜γij(p) + 1<br />

3 ξiξj. (B.15)<br />

Hence ˜P (p,ξ) is clearly a symmetric bilinear form. Moreover it is positive definite for ξ = 0: for<br />

any vector v ∈ Tp(Σ) such that v = 0, we have<br />

˜P (p,ξ)(v,v) = ˜γ kl (p)ξkξl ˜γij(p)v i v j + 1<br />

3 (ξiv i ) 2 > 0, (B.16)<br />

where the > 0 follows from the positive definite character <strong>of</strong> ˜γ. ˜P (p,ξ) being positive definite<br />

symmetric bilinear form, we conclude that P (p,ξ) is an isomorphism <strong>and</strong> that all its eigenvalues<br />

are real <strong>and</strong> strictly positive. Therefore ˜∆L is a strongly elliptic operator.<br />

B.2.3 Kernel<br />

Let us now determine the kernel <strong>of</strong> ˜∆L. Clearly this kernel contains the kernel <strong>of</strong> the conformal<br />

Killing operator ˜L. Actually it is not larger than that kernel:<br />

ker ˜∆L = ker ˜L . (B.17)<br />

Let us establish this property. For any vector field v ∈ T (Σ), we have<br />

<br />

<br />

˜γijv<br />

Σ<br />

i ∆L<br />

˜ v j ˜γ d 3 x =<br />

=<br />

=<br />

<br />

<br />

˜γijv<br />

Σ<br />

i Dl( ˜ ˜ Lv) jl ˜γ d 3 x<br />

Σ<br />

∂Σ<br />

<br />

˜Dl ˜γijv i ( ˜ Lv) jl<br />

˜γijv i ( ˜ Lv) jl <br />

2<br />

˜sl ˜q d y −<br />

− ˜γij ˜ Dlv i ( ˜ Lv) jl ˜γ d 3 x<br />

<br />

˜γij<br />

Σ<br />

˜ Dlv i ( ˜ Lv) jl ˜γ d 3 x,<br />

(B.18)<br />

where the Gauss-Ostrogradsky theorem has been used to get the last line. We shall consider<br />

two situations for (Σ,γ):

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

Saved successfully!

Ooh no, something went wrong!