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 ...
58 3+1 decomposition of Einstein equation Besides, since ∂i and ∂j are tangent to Σt, gij = g(∂i,∂j) = γ(∂i,∂j) = γij. (4.46) Collecting Eqs. (4.43), (4.45) and (4.46), we get the following expression of the metric components in terms of 3+1 quantities: gαβ = g00 g0j gi0 gij or, in terms of line elements [using Eq. (4.40)], −N2 + βkβ = k βj βi γij , (4.47) gµν dx µ dx ν = −N 2 dt 2 + γij(dx i + β i dt)(dx j + β j dt) . (4.48) The components of the inverse metric are given by the matrix inverse of (4.47): g αβ g00 g = 0j = g i0 g ij − 1 N2 βi N2 βj N2 γij − βiβj N2 . (4.49) Indeed, it is easily checked that the matrix product g αµ gµβ is equal to the identity matrix δ α β . Remark : Notice that gij = γij but that in general g ij = γ ij . One can deduce from the above formulæ a simple relation between the determinants of g and γ. Let us first define the latter ones by g := det(gαβ) , (4.50) γ := det(γij) . (4.51) Notice that g and γ depend upon the choice of the coordinates (x α ). They are not scalar quantities, but scalar densities. Using Cramer’s rule for expressing the inverse (g αβ ) of the matrix (gαβ), we have g 00 = C00 det(gαβ) C00 = , (4.52) g where C00 is the element (0,0) of the cofactor matrix associated with (gαβ). It is given by C00 = (−1) 0 M00 = M00, where M00 is the minor (0,0) of the matrix (gαβ), i.e. the determinant of the 3 × 3 matrix deduced from (gαβ) by suppressing the first line and the first column. From Eq. (4.47), we read M00 = det(γij) = γ. (4.53) Hence Eq. (4.52) becomes g 00 = γ . (4.54) g Expressing g 00 from Eq. (4.49) yields then g = −N 2 γ, or equivalently, √ −g = N √ γ . (4.55)
4.3 3+1 Einstein equation as a PDE system 59 4.2.4 Choice of coordinates via the lapse and the shift We have seen above that giving a coordinate system (x α ) on M such that the hypersurfaces x 0 = const. are spacelike determines uniquely a lapse function N and a shift vector β. The converse is true in the following sense: setting on some hypersurface Σ0 a scalar field N, a vector field β and a coordinate system (x i ) uniquely specifies a coordinate system (x α ) in some neighbourhood of Σ0, such that the hypersurface x 0 = 0 is Σ0. Indeed, the knowledge of the lapse function a each point of Σ0 determines a unique vector m = Nn and consequently the location of the “next” hypersurface Σδt by Lie transport along m (cf. Sec. 3.3.2). Graphically, we may also say that for each point of Σ0 the lapse function specifies how far is the point of Σδt located “above” it (“above” meaning perpendicularly to Σ0, cf. Fig. 3.2). Then the shift vector tells how to propagate the coordinate system (x i ) from Σ0 to Σδt (cf. Fig. 4.1). This way of choosing coordinates via the lapse function and the shift vector is one of the main topics in 3+1 numerical relativity and will be discussed in detail in Chap. 9. 4.3 3+1 Einstein equation as a PDE system 4.3.1 Lie derivatives along m as partial derivatives Let us consider the term LmK which occurs in the 3+1 Einstein equation (4.15). Thanks to Eq. (4.28), we can write LmK = L∂t K − Lβ K. (4.56) This implies that L∂t K is a tensor field tangent to Σt, since both LmK and Lβ K are tangent to Σt, the former by the property (3.32) and the latter because β and K are tangent to Σt. Moreover, if one uses tensor components with respect to a coordinate system (x α ) = (t,x i ) adapted to the foliation, the Lie derivative along ∂t reduces simply to the partial derivative with respect to t [cf. Eq. (A.3)]: L∂tKij = ∂Kij . (4.57) ∂t By means of formula (A.6), one can also express Lβ K in terms of partial derivatives: with k ∂Kij ∂β Lβ Kij = β + Kkj ∂xk k ∂β + Kik ∂xi k . (4.58) ∂xj Similarly, the relation (3.24) between Lmγ and K becomes L∂t γ − Lβ γ = −2NK, (4.59) L∂tγij = ∂γij . (4.60) ∂t and, evaluating the Lie derivative with the connection D instead of partial derivatives [cf. Eq. (A.8)]: i.e. Lβ γij = β k Dkγij +γkjDiβ k + γikDjβ k , (4.61) =0 Lβ γij = Diβj + Djβi. (4.62)
- 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 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 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
- Page 96 and 97: 96 Conformal decomposition hence Lm
- Page 98 and 99: 98 Conformal decomposition 6.5.2 Ha
- Page 100 and 101: 100 Conformal decomposition discuss
- Page 102 and 103: 102 Conformal decomposition Remark
- Page 104 and 105: 104 Asymptotic flatness and global
- Page 106 and 107: 106 Asymptotic flatness and global
58 <strong>3+1</strong> decomposition <strong>of</strong> Einstein equation<br />
Besides, since ∂i <strong>and</strong> ∂j are tangent to Σt,<br />
gij = g(∂i,∂j) = γ(∂i,∂j) = γij. (4.46)<br />
Collecting Eqs. (4.43), (4.45) <strong>and</strong> (4.46), we get the following expression <strong>of</strong> the metric components<br />
in terms <strong>of</strong> <strong>3+1</strong> quantities:<br />
gαβ =<br />
g00 g0j<br />
gi0 gij<br />
or, in terms <strong>of</strong> line elements [using Eq. (4.40)],<br />
<br />
−N2 + βkβ<br />
=<br />
k βj<br />
βi<br />
γij<br />
<br />
, (4.47)<br />
gµν dx µ dx ν = −N 2 dt 2 + γij(dx i + β i dt)(dx j + β j dt) . (4.48)<br />
The components <strong>of</strong> the inverse metric are given by the matrix inverse <strong>of</strong> (4.47):<br />
g αβ <br />
g00 g<br />
=<br />
0j<br />
<br />
=<br />
g i0 g ij<br />
− 1<br />
N2 βi N2 βj N2 γij − βiβj N2 <br />
. (4.49)<br />
Indeed, it is easily checked that the matrix product g αµ gµβ is equal to the identity matrix δ α β .<br />
Remark : Notice that gij = γij but that in general g ij = γ ij .<br />
One can deduce from the above formulæ a simple relation between the determinants <strong>of</strong> g<br />
<strong>and</strong> γ. Let us first define the latter ones by<br />
g := det(gαβ) , (4.50)<br />
γ := det(γij) . (4.51)<br />
Notice that g <strong>and</strong> γ depend upon the choice <strong>of</strong> the coordinates (x α ). They are not scalar<br />
quantities, but scalar densities. Using Cramer’s rule for expressing the inverse (g αβ ) <strong>of</strong> the<br />
matrix (gαβ), we have<br />
g 00 = C00<br />
det(gαβ)<br />
C00<br />
= , (4.52)<br />
g<br />
where C00 is the element (0,0) <strong>of</strong> the c<strong>of</strong>actor matrix associated with (gαβ). It is given by<br />
C00 = (−1) 0 M00 = M00, where M00 is the minor (0,0) <strong>of</strong> the matrix (gαβ), i.e. the determinant<br />
<strong>of</strong> the 3 × 3 matrix deduced from (gαβ) by suppressing the first line <strong>and</strong> the first column. From<br />
Eq. (4.47), we read<br />
M00 = det(γij) = γ. (4.53)<br />
Hence Eq. (4.52) becomes<br />
g 00 = γ<br />
. (4.54)<br />
g<br />
Expressing g 00 from Eq. (4.49) yields then g = −N 2 γ, or equivalently,<br />
√ −g = N √ γ . (4.55)