Relativistic slim disk model by Alexander Sadowski.
This model gives the structure of the polytropic relativistic slim disk model as described in Sadowski et al. (2011).
The model interpolates through a set of tabulated numerical solutions. These solution are computed for a reference mass M=10 Msun, and although the model can rescale the disk structure to any mass, one should not go very much off the stellar mass range. Let's say masses 5-20 M_sun should give a reasonably good results.
The model reads and interpolates tables with numerical slim disk solutions. Those tables are included in data subfolder in a compressed archive and have to be untared before using the model
cd models/slimdisk/data
tar -xvf slimdisk-polytropic.tar.xzThe model is initialized using diskmodel_init() function that takes three arguments: black-hole mass, black-hole spin and options. The options parameter specifies all other settings, where the most imporatnt is the path to the folder with the slim disk solutions. An example code to intialize the model is
float bh_mass = 10.0;
float bh_spin = 0.9;
float mdot = 0.1;
char options[256];
sprintf(options, "path=./data,mdot=%f", mdot);
int res = diskmodel_init(bh_mass, bh_spin, options);
if (res != 0) return res;
float F = diskmodel_flux(10.0);
diskmodel_done();See the main() function in disk-sd.c.
The initialization function diskmodel_init() takes the following model-specific options through the params parameter:
| Option | Description |
|---|---|
| path | disk path to the directory with slim disk solutions (default current directory) |
| alpha | value of alpha-viscosity paramater (range 0.005-0.1; default 0.1) |
| mdot | mass accretion rate in units of Mdot_Edd (range 0.01-500; default 0.1) |
| lumi | total disk luminosity in units of L_Edd; takes precedence over mdot option |
$ makeAfter compilation, two files are produced:
- a shared object library
disk-sd.so - a test executable
disk-sd-test
The shared object library can be used in custom programs. The test executable demonstrates a simple dump of quantities for a given choice of spin and mass accretion rate.