blob: 04d7f8f2164dcb9c10fdf3030e0dee26a9baa362 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
|
\documentclass{article}
\begin{document}
\title{ADMConstraints}
\author{Gabrielle Allen}
\date{April 2002}
\maketitle
\abstract{Calculate the values of the Hamiltonian and Momentum Constraints on the computational grid for output.}
\section{Physics}
This analysis thorn calculates the values of the Hamiltonian and
momentum constraint equations on the computational grid during a spacetime
evolution. The Hamiltonian Constraint (grid function: ham) is given by
\begin{equation}
H = R - K^{i}_{\;\;j} K^{j}_{\;\;i} + trK^2 - 16 \pi \rho
\end{equation}
The Momentum Constraints (grid functions momx, momy, momz) are given by
\begin{equation}
M_i = \nabla_j K_i^{\;\;j} - \nabla_i trK - 8 \pi S_i
\end{equation}
where
\begin{equation}
S_i = -g_{ia} n_b T^{ab} = - g_i^a n^b T_{ab} = -\frac{1}{\alpha} (T_{i0} - \beta^j T_{ij})
\end{equation}
\section{Comments}
\begin{itemize}
\item Symmetry boundary conditions ({\tt CactusBase/CartGrid3D}) are implemented.
\item Default behaviour is to apply so called flat boundary conditions to the calculated values of the constraint equations.
\item Excision ({\tt Excision/LegoExcision}) is included.
\item The 3-metric is only included through the Macros in {\tt CactusEinstein/Einstein} and as such the physical metric is used.
\item Matter is included using the {\tt CalcTmuNu} mechanism.
\item The constraint calculation is default performed at the ANALYSIS time bin, if the constraint values are needed at every iteration, the parameter {\tt constraints\_persist} should be used to schedule the calculation instead at POSTSTEP.
\end{itemize}
% Automatically created from the ccl files by using gmake thorndoc
\include{interface}
\include{param}
\include{schedule}
\end{document}
|
