This thorn is documented within two Jupyter notebooks, contained within the FishboneMoncriefID/doc/ directory:
Tutorial-FishboneMoncriefID.ipynb: Contains basic description, equations, and citations for this initial data type.
Tutorial-ETK_thorn-FishboneMoncriefID.ipynb: Constructs this Einstein Toolkit thorn.
random_max | Scope: private | REAL |
Description: Ceiling value of random perturbation to initial pressure, where perturbed pressure
= pressure*(1.0 + (random_min + (random_max-random_min)*RAND[0,1)))
| ||
| Range | Default: 0.02 | |
*:* | Any value
| |
random_min | Scope: private | REAL |
Description: Floor value of random perturbation to initial pressure, where perturbed pressure =
pressure*(1.0 + (random_min + (random_max-random_min)*RAND[0,1)))
| ||
| Range | Default: -0.02 | |
*:* | Any value
| |
a | Scope: restricted | REAL |
Description: The spin parameter of the black hole
| ||
| Range | Default: 0.9375 | |
0:1.0 | Positive values, up to 1. Negative disallowed, as certain roots are chosen in
the hydro fields setup. Check those before enabling negative spins!
| |
a_b | Scope: restricted | REAL |
Description: Scaling factor for the vector potential
| ||
| Range | Default: 1.0 | |
*:* | ||
gamma | Scope: restricted | REAL |
Description: Equation of state: P = kappa * rhoĝamma
| ||
| Range | Default: 1.3333333333333333333333333333 | |
0.0:* | Positive values
| |
kappa | Scope: restricted | REAL |
Description: Equation of state: P = kappa * rhoĝamma
| ||
| Range | Default: 1.0e-3 | |
0.0:* | Positive values
| |
m | Scope: restricted | REAL |
Description: Kerr-Schild BH mass. Probably should always set M=1.
| ||
| Range | Default: 1.0 | |
0.0:* | Must be positive
| |
r_at_max_density | Scope: restricted | REAL |
Description: Radius at maximum disk density. Needs to be > r_in
| ||
| Range | Default: 12.0 | |
0.0:* | Must be positive
| |
r_in | Scope: restricted | REAL |
Description: Fixes the inner edge of the disk
| ||
| Range | Default: 6.0 | |
0.0:* | Must be positive
| |
initial_hydro | Scope: shared from HYDROBASE | KEYWORD |
Extends ranges:
| ||
FishboneMoncriefID | Initial GRHD data from FishboneMoncriefID solution
| |
Implements:
fishbonemoncriefid
Inherits:
admbase
grid
hydrobase
This section lists all the variables which are assigned storage by thorn WVUThorns/FishboneMoncriefID. Storage can either last for the duration of the run (Always means that if this thorn is activated storage will be assigned, Conditional means that if this thorn is activated storage will be assigned for the duration of the run if some condition is met), or can be turned on for the duration of a schedule function.
| Always: | |
| ADMBase::metric[metric_timelevels] ADMBase::curv[metric_timelevels] ADMBase::lapse[lapse_timelevels] ADMBase::shift[shift_timelevels] | |
HydroBase_Initial
fishbonemoncrief_et_grhd_initial
set up general relativistic hydrodynamic (grhd) fields for fishbone-moncrief disk
| Language: | c | |
| Reads: | grid::x(everywhere) | |
| grid::y(everywhere) | ||
| grid::z(everywhere) | ||
| Type: | function | |
| Writes: | admbase::alp(everywhere) | |
| admbase::betax(everywhere) | ||
| admbase::betay(everywhere) | ||
| admbase::betaz(everywhere) | ||
| admbase::kxx(everywhere) | ||
| admbase::kxy(everywhere) | ||
| admbase::kxz(everywhere) | ||
| admbase::kyy(everywhere) | ||
| admbase::kyz(everywhere) | ||
| admbase::kzz(everywhere) | ||
| admbase::gxx(everywhere) | ||
| admbase::gxy(everywhere) | ||
| admbase::gxz(everywhere) | ||
| admbase::gyy(everywhere) | ||
| admbase::gyz(everywhere) | ||
| admbase::gzz(everywhere) | ||
| hydrobase::vel(everywhere) | ||
| hydrobase::rho(everywhere) | ||
| hydrobase::eps(everywhere) | ||
| hydrobase::press(everywhere) | ||
CCTK_INITIAL
fishbonemoncrief_et_grhd_initial__perturb_pressure
add random perturbation to initial pressure, after seed magnetic fields have been set up (in case we’d like the seed magnetic fields to depend on the pristine pressures)
| After: | seed_magnetic_fields | |
| Before: | illinoisgrmhd_id_converter | |
| Language: | c | |
| Type: | function | |