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 ...
122 Asymptotic flatness and global quantities to the hypersurfaces Σt. Then n · φ = 0 and the integrand in the definition (7.100) is Accordingly Eq. (7.100) becomes ∇ µ φ ν dSµν = ∇µφν(s µ n ν − n µ s ν ) √ q d 2 y = 2∇µφν s µ n ν√ q d 2 y = −2s µ φν∇µn ν√ q d 2 y = 2Kijs i φ j√ q d 2 y. (7.103) JK = 1 Kijs 8π St i φ j√ q d 2 y . (7.104) Remark : Contrary to the 3+1 expression of the Komar mass which turned out to be very different from the expression of the ADM mass, the 3+1 expression of the Komar angular momentum as given by Eq. (7.104) is very similar to the expression of the angular momentum deduced from the Hamiltonian formalism, i.e. Eq. (7.63). The only differences are that it is no longer necessary to take the limit St → ∞ and that there is no trace term Kγijs i φ j in Eq. (7.104). Moreover, if one evaluates the Hamiltonian expression in the asymptotically maximal gauge (7.65) then K = O(r −3 ) and thanks to the asymptotic orthogonality of s and φ, γijs i φ j = O(1), so that Kγijs i φ j does not contribute to the integral and expressions (7.104) and (7.63) are then identical. Example : A trivial example is provided by Schwarzschild spacetime, which among other things is axisymmetric. For the 3+1 foliation associated with the Schwarzschild coordinates (7.16), the extrinsic curvature tensor K vanishes identically, so that Eq. (7.104) yields immediately JK = 0. For other foliations, like that associated with Eddington-Finkelstein coordinates, K is no longer zero but is such that Kijs i φ j = 0, yielding again JK = 0 (as it should be since the Komar angular momentum is independent of the foliation). Explicitely for Eddington-Finkelstein coordinates, Kijs i = − 2m r 2 1 + m r 1 + 2m r , 0, 0 , (7.105) (see e.g. Eq. (D.25) in Ref. [146]) and φ j = (0,0,1), so that obviously Kijs i φ j = 0. Example : The most natural non trivial example is certainly that of Kerr spacetime. Let us use the 3+1 foliation associated with the standard Boyer-Lindquist coordinates (t,r,θ,ϕ) and evaluate the integral (7.104) by choosing for St a sphere r = const. Then y a = (θ,ϕ). The Boyer-Lindquist components of φ are φ i = (0,0,1) and those of s are s i = (s r ,0,0) since γij is diagonal is these coordinates. The formula (7.104) then reduces to JK = 1 Krϕs 8π r=const r√ q dθ dϕ. (7.106)
7.6 Komar mass and angular momentum 123 The extrinsic curvature component Krϕ can be evaluated via formula (4.63), which reduces to 2NKij = Lβ γij since ∂γij/∂t = 0. From the Boyer-Lindquist line element (see e.g. Eq. (5.29) in Ref. [156]), we read the components of the shift: (β r ,β θ ,β ϕ 2amr ) = 0, 0, − (r2 + a2 )(r2 + a2 cos2 θ) + 2a2mr sin2 , (7.107) θ where m and a are the two parameters of the Kerr solution. Then, using Eq. (A.6), Krϕ = 1 2N Lβ γrϕ = 1 ϕ ∂γrϕ ∂β β +γϕϕ 2N ∂ϕ =0 ϕ ∂β + γrϕ ∂r ϕ = ∂ϕ =0 1 2N γϕϕ ∂βϕ . ∂r (7.108) Hence JK = 1 s 16π r=const r N γϕϕ ∂βϕ √ q dθ dϕ. (7.109) ∂r The values of s r , N, γϕϕ and √ q can all be read on the Boyer-Lindquist line element. However this is a bit tedious. To simplify things, let us evaluate JK only in the limit r → ∞. Then s r ∼ 1, N ∼ 1, γϕϕ ∼ r 2 sin 2 θ, √ q ∼ r 2 sin θ and, from Eq. (7.107), βϕ ∼ −2am/r3 , so that JK = 1 16π r 2 2 6am sin r4 r2 sin θ dθ dϕ = 3am π × 2π × sin 8π 0 3 θ dθ. (7.110) Hence, as expected, r=const JK = am. (7.111) Let us now find the 3+1 expression of the volume version (7.101) of the Komar angular momentum. We have n · φ = 0 and, from the 3+1 decomposition (4.10) of T: Hence formula (7.101) becomes with JK = T(n,φ) = −〈p,φ〉. (7.112) Σt 〈p,φ〉 √ γ d 3 x + J H K , (7.113) J H K = 1 Kijs 8π Ht i φ j√ q d 2 y . (7.114) Example : Let us consider a perfect fluid. Then p = (E + P)U [Eq. (5.61)], so that JK = (E + P)U · φ √ γ d 3 x + J H K . (7.115) Σt Taking φ = −y∂x +x∂y (symmetry axis = z-axis), the Newtonian limit of this expression is then JK = ρ(−yU x + xU y )dxdy dz, (7.116) Σt i.e. we recognize the standard expression for the angular momentum around the z-axis.
- 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
- Page 108 and 109: 108 Asymptotic flatness and global
- Page 110 and 111: 110 Asymptotic flatness and global
- Page 112 and 113: 112 Asymptotic flatness and global
- Page 114 and 115: 114 Asymptotic flatness and global
- Page 116 and 117: 116 Asymptotic flatness and global
- Page 118 and 119: 118 Asymptotic flatness and global
- Page 120 and 121: 120 Asymptotic flatness and global
- Page 124 and 125: 124 Asymptotic flatness and global
- Page 126 and 127: 126 The initial data problem Notice
- Page 128 and 129: 128 The initial data problem where
- Page 130 and 131: 130 The initial data problem 8.2.3
- Page 132 and 133: 132 The initial data problem where
- Page 134 and 135: 134 The initial data problem Figure
- Page 136 and 137: 136 The initial data problem Figure
- Page 138 and 139: 138 The initial data problem In par
- Page 140 and 141: 140 The initial data problem Accord
- Page 142 and 143: 142 The initial data problem Remark
- Page 144 and 145: 144 The initial data problem Remark
- 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
7.6 Komar mass <strong>and</strong> angular momentum 123<br />
The extrinsic curvature component Krϕ can be evaluated via formula (4.63), which reduces<br />
to 2NKij = Lβ γij since ∂γij/∂t = 0. From the Boyer-Lindquist line element (see e.g.<br />
Eq. (5.29) in Ref. [156]), we read the components <strong>of</strong> the shift:<br />
(β r ,β θ ,β ϕ <br />
2amr<br />
) = 0, 0, −<br />
(r2 + a2 )(r2 + a2 cos2 θ) + 2a2mr sin2 <br />
, (7.107)<br />
θ<br />
where m <strong>and</strong> a are the two parameters <strong>of</strong> the Kerr solution. Then, using Eq. (A.6),<br />
Krϕ = 1<br />
2N Lβ γrϕ = 1<br />
<br />
ϕ ∂γrϕ ∂β<br />
β +γϕϕ<br />
2N ∂ϕ<br />
<br />
=0<br />
ϕ ∂β<br />
+ γrϕ<br />
∂r ϕ <br />
=<br />
∂ϕ<br />
<br />
=0<br />
1<br />
2N γϕϕ<br />
∂βϕ .<br />
∂r<br />
(7.108)<br />
Hence<br />
JK = 1<br />
<br />
s<br />
16π r=const<br />
r<br />
N γϕϕ<br />
∂βϕ √<br />
q dθ dϕ. (7.109)<br />
∂r<br />
The values <strong>of</strong> s r , N, γϕϕ <strong>and</strong> √ q can all be read on the Boyer-Lindquist line element.<br />
However this is a bit tedious. To simplify things, let us evaluate JK only in the limit<br />
r → ∞. Then s r ∼ 1, N ∼ 1, γϕϕ ∼ r 2 sin 2 θ, √ q ∼ r 2 sin θ <strong>and</strong>, from Eq. (7.107),<br />
βϕ ∼ −2am/r3 , so that<br />
JK = 1<br />
<br />
16π<br />
r 2 2 6am<br />
sin<br />
r4 r2 sin θ dθ dϕ = 3am<br />
π<br />
× 2π × sin<br />
8π 0<br />
3 θ dθ. (7.110)<br />
Hence, as expected,<br />
r=const<br />
JK = am. (7.111)<br />
Let us now find the <strong>3+1</strong> expression <strong>of</strong> the volume version (7.101) <strong>of</strong> the Komar angular<br />
momentum. We have n · φ = 0 <strong>and</strong>, from the <strong>3+1</strong> decomposition (4.10) <strong>of</strong> T:<br />
Hence formula (7.101) becomes<br />
with<br />
<br />
JK =<br />
T(n,φ) = −〈p,φ〉. (7.112)<br />
Σt<br />
〈p,φ〉 √ γ d 3 x + J H K , (7.113)<br />
J H K = 1<br />
<br />
Kijs<br />
8π Ht<br />
i φ j√ q d 2 y . (7.114)<br />
Example : Let us consider a perfect fluid. Then p = (E + P)U [Eq. (5.61)], so that<br />
<br />
JK = (E + P)U · φ √ γ d 3 x + J H K . (7.115)<br />
Σt<br />
Taking φ = −y∂x +x∂y (symmetry axis = z-axis), the Newtonian limit <strong>of</strong> this expression<br />
is then<br />
<br />
JK = ρ(−yU x + xU y )dxdy dz, (7.116)<br />
Σt<br />
i.e. we recognize the st<strong>and</strong>ard expression for the angular momentum around the z-axis.