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 ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
Chapter 8<br />
The initial data problem<br />
Contents<br />
8.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125<br />
8.2 Conformal transverse-traceless method . . . . . . . . . . . . . . . . . 127<br />
8.3 Conformal thin s<strong>and</strong>wich method . . . . . . . . . . . . . . . . . . . . . 139<br />
8.4 Initial data for binary systems . . . . . . . . . . . . . . . . . . . . . . 145<br />
8.1 Introduction<br />
8.1.1 The initial data problem<br />
We have seen in Chap. 4 that thanks to the <strong>3+1</strong> decomposition, the resolution <strong>of</strong> Einstein<br />
equation amounts to solving a Cauchy problem, namely to evolve “forward in time” some initial<br />
data. However this is a Cauchy problem with constraints. This makes the set up <strong>of</strong> initial data<br />
a non trivial task, because these data must obey the constraints. Actually one may distinguish<br />
two problems:<br />
• The mathematical problem: given some hypersurface Σ0, find a Riemannian metric γ,<br />
a symmetric bilinear form K <strong>and</strong> some matter distribution (E,p) on Σ0 such that the<br />
Hamiltonian constraint (4.65) <strong>and</strong> the momentum constraint (4.66) are satisfied:<br />
R + K 2 − KijK ij = 16πE (8.1)<br />
DjK j<br />
i − DiK = 8πpi . (8.2)<br />
In addition, the matter distribution (E,p) may have some constraints from its own. We<br />
shall not discuss them here.<br />
• The astrophysical problem: make sure that the solution to the constraint equations has<br />
something to do with the physical system that one wish to study.