FishboneMoncriefID: Black Hole Accretion Disk Initial Data

Author: Zachariah B. Etienne

\( \)Date: 2019-09-28 (Sat, 28 Sep 2019) \( \)

This thorn is documented within two Jupyter notebooks, contained within the FishboneMoncriefID/doc/ directory:

  1. Tutorial-FishboneMoncriefID.ipynb: Contains basic description, equations, and citations for this initial data type.

  2. Tutorial-ETK_thorn-FishboneMoncriefID.ipynb: Constructs this Einstein Toolkit thorn.

1 Parameters

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

2 Interfaces

General

Implements:

fishbonemoncriefid

Inherits:

admbase

grid

hydrobase

3 Schedule

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.

Storage

 

Always:  
ADMBase::metric[metric_timelevels] ADMBase::curv[metric_timelevels] ADMBase::lapse[lapse_timelevels] ADMBase::shift[shift_timelevels]  
   

Scheduled Functions

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