blob: f8203823b00c40525f1f5659b2b307d842bae092 (
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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
|
% Thorn documentation template
\documentclass{article}
\begin{document}
\title{IDBrillData}
\author{Carsten Gundlach}
\date{6 September 1999}
\maketitle
\abstract{This thorn creates initial data for Brill wave spacetimes.
It can create both axisymmetric data (in a 3D cartesian grid), as
well as data with an angular dependency.}
\section{Purpose}
The purpose of this thorn is to create initial data for a Brill wave
spacetime. It does so by starting from a three--metric of the form
originally considered by Brill
\begin{equation}
ds^2 = \Psi^4 \left[ e^{2q} \left( d\rho^2 + dz^2 \right)
+ \rho^2 d\phi^2 \right] =\Psi^4 \hat{ds}^{2},
\label{eqn:brillmetric}
\end{equation}
where $q$ is a free function subject to certain regularity and
fall-off conditions and $\Psi$ is a conformal factor to be solved for.
The thorn considers several different forms of the function $q$
depending on certain parameters that will be described below.
Substituting the metric above into the Hamiltonian constraint results
in an elliptic equation for the conformal factor $\Psi$ that can be
solved numerically once the function $q$ has been specified. The
initial data is also assumed to be time-symmetric, so the momentum
constraints are trivially satisfied.
The thorn is activated by choosing the CactusEinstein/ADMBase parameter
``initial\_data'' in one of the following two ways:
\begin{itemize}
\item initial\_data = ``brilldata'': Axisymmetric Brill wave initial data
(but calculated in a cartesian grid!).
\item initial\_data = ``brilldata3D'': Brill wave initial data with an
angular dependency.
\end{itemize}
\section{Parameters for the thorn}
The thorn is controlled by the following parameters:
\begin{itemize}
\item brill\_q (INT): Form of the function $q$ [0,1,2] (default 2):
\begin{itemize}
\item brill\_q = 0:
\[
q = a \; \frac{\rho^{2+b}}{r^2} \left( \frac{z}{\sigma_z} \right)^2
e^{-(\rho - \rho_0^2)}
\]
\item brill\_q = 1:
\[
q = a \left( \frac{\rho}{\sigma_\rho} \right)^b \frac{1}{1 + \left[
\left( r^2 - r_0^2 \right) / \sigma_r^2 \right]^{c/2}}
\]
\item brill\_q = 2:
\[
q = a \left( \frac{\rho}{\sigma_\rho} \right)^b e^{-\left[
\left( r^2 - r_0^2 \right) / \sigma_r^2 \right]^{c/2}}
\]
\item If one specifies 3D data (see above), the function $q$ is multiplied
by an additional factor with an angular dependency:
\[
q \rightarrow q \left[ 1 + d \frac{\rho^m}{1 + e \rho^m}
\cos^2 \left( n \phi + \phi_0 \right) \right]
\]
\end{itemize}
\item brill\_a (REAL): Amplitude (default 0.0).
\item brill\_b (REAL): $b$ in above expressions (default 2.0).
\item brill\_c (REAL): $c$ in above expressions (default 2.0).
\item brill\_d (REAL): $d$ in above expressions (default 0.0).
\item brill\_e (REAL): $e$ in above expressions (default 1.0).
\item brill\_m (REAL): $m$ in above expressions (default 2.0).
\item brill\_n (REAL): $n$ in above expressions (default 2.0).
\item brill\_r0 (REAL): $r_0$ in above expressions (default 0.0).
\item brill\_rho0 (REAL): $\rho_0$ in above expressions
(default 0.0).
\item brill\_phi0 (REAL): $\phi_0$ in above expressions
(default 0.0).
\item brill\_sr (REAL): $\sigma_r$ in above expressions
(default 1.0).
\item brill\_srho (REAL): $\sigma_\rho$ in above
expressions (default 1.0).
\end{itemize}
The elliptic solver is controlled by the additional parameters:
\begin{itemize}
\item solver (KEYWORD): Elliptic solver used to solve the
hamiltonian constraint [sor/petsc/bam] (default ``sor'').
\item thresh (REAL): Threshold for elliptic solver (default
0.00001).
\end{itemize}
% Automatically created from the ccl files
% Do not worry for now.
\include{interface}
\include{param}
\include{schedule}
\end{document}
|
