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 ...
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).
- Page 146 and 147: 146 The initial data problem 8.4.1
- Page 148 and 149: 148 The initial data problem Since
- Page 150 and 151: 150 The initial data problem • fo
- Page 152 and 153: 152 Choice of foliation and spatial
- Page 154 and 155: 154 Choice of foliation and spatial
- Page 156 and 157: 156 Choice of foliation and spatial
- Page 158 and 159: 158 Choice of foliation and spatial
- Page 160 and 161: 160 Choice of foliation and spatial
- Page 162 and 163: 162 Choice of foliation and spatial
- Page 164 and 165: 164 Choice of foliation and spatial
- Page 166 and 167: 166 Choice of foliation and spatial
- Page 168 and 169: 168 Choice of foliation and spatial
- Page 170 and 171: 170 Choice of foliation and spatial
- Page 172 and 173: 172 Choice of foliation and spatial
- Page 174 and 175: 174 Choice of foliation and spatial
- Page 176 and 177: 176 Evolution schemes been used by
- Page 178 and 179: 178 Evolution schemes Comparing wit
- Page 180 and 181: 180 Evolution schemes Now the ∇-d
- Page 182 and 183: 182 Evolution schemes coordinates (
- Page 184 and 185: 184 Evolution schemes = 1 2 −∆
- Page 186 and 187: 186 Evolution schemes corresponds t
- Page 188 and 189: 188 Evolution schemes
- Page 190 and 191: 190 Lie derivative Figure A.1: Geom
- Page 192 and 193: 192 Lie derivative
- Page 194 and 195: 194 Conformal Killing operator and
- Page 198 and 199: 198 Conformal Killing operator and
- Page 200 and 201: 200 BIBLIOGRAPHY [13] A. Anderson a
- Page 202 and 203: 202 BIBLIOGRAPHY [43] T.W. Baumgart
- Page 204 and 205: 204 BIBLIOGRAPHY [75] M. Campanelli
- Page 206 and 207: 206 BIBLIOGRAPHY [105] G. Darmois :
- Page 208 and 209: 208 BIBLIOGRAPHY [135] H. Friedrich
- Page 210 and 211: 210 BIBLIOGRAPHY [163] J. Isenberg
- Page 212 and 213: 212 BIBLIOGRAPHY [195] S. Nissanke
- Page 214 and 215: 214 BIBLIOGRAPHY [226] M. Shibata :
- Page 216 and 217: 216 BIBLIOGRAPHY [255] K. Taniguchi
- Page 218 and 219: Index 1+log slicing, 161 3+1 formal
- Page 220: 220 INDEX Ricci identity, 18 Ricci
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 (Σ,γ):