Abstract
This thorn creates time symmetric 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.
The purpose of this thorn is to create (time symmetric) initial data for a Brill wave spacetime. It does so by starting from a three–metric of the form originally considered by Brill
(1) |
where is a free function subject to certain regularity and fall-off conditions, and is a conformal factor to be solved for.
Thorn IDBrillData provides three choices for the function: an exponential form, (IDBrillData::q_function = "exp")
(2) |
a generalized form of the function first written down by Eppley (IDBrillData::q_function = "eppley")
(3) |
and the (default) Gundlach function which includes the Holz form (IDBrillData::q_function = "gundlach")
(4) |
Substituting the metric into the Hamiltonian constraint gives an elliptic equation for the conformal factor which is then numerically solved for a given function :
(5) |
where the conformal Ricci scalar is found to be
Assuming the initial data to be time symmetric means that the momentum constraints are trivially satisfied.
In the case of axisymmetry (that is in the above expressions for ), the Hamiltonian constraint can be written as an elliptic equation for with just the flat space Laplacian,
(7) |
If the initial data is chosen to be ADMBase::initial_data = "brilldata2D" then this elliptic equation is solved rather than the equation above.
Brill initial data is activated by choosing the CactusEinstein/ADMBase parameter initial_data to be brilldata, or for the case of axisymmetry brilldata2D can also be used.
The parameter IDBrillData::q_function chooses the form of the function to be used, defaulting to the Gundlach expression.
Additional IDBrillData parameters for each form of fix the remaining freedom:
(exp_a, exp_b, exp_sigmaz,exp_rho0)
(eppley_a, eppley_b, eppley_sigmarho, eppley_r0, eppley_sigmar, eppley_c)
(gundlach_a, gundlach_b, gundlach_sigmarho, gundlach_r0, gundlach_sigmar, gundlach_c)
(brill3d_d, brill3d_m, brill3d_e, brill3d_n, brill3d_phi0)
Note that the default expression is
IDBrillData can use the elliptic solvers (type LinMetric) provided by CactusEinstein/EllSOR,
AEIThorns/BAM_Elliptic, or CactusElliptic/EllPETSc to solve the equation resulting from the Hamiltonian
constraint. In all cases the parameter thresh sets the threshold for the elliptic solve. The choice of elliptic solver
is made through the parameter brill_solver: