EOS_Polytrope

Ian Hawke

22/4/2002

Abstract

EOS_Polytrope

1 The equations

This equation provides a polytropic equation of state to thorns using the CactusEOS interface found in EOS_Base. As such it’s a fake, as EOS_Base assumes that, e.g., the pressure is a function of both density and specific internal energy. Here the pressure is just a function of the density, and is set appropriately (the specific internal energy is always ignored).

The two fluid constants are K (eos_k) and Γ (eos_gamma), which default to 100 and 2 respectively. The formulas that are applied under the appropriate EOS_Base function calls are

P = KρΓ (1) 𝜖 = KρΓ1 Γ 1 (2) ρ = P (Γ 1)𝜖 (3) P ρ = KΓρΓ1 (4) P 𝜖 = 0. (5)

To calculate the units of the Cactus quantities and back, remember that G = c = M = 1 in Cactus.
Here is one example how to calculate densities:

ρCactus = G3M 2 c6 ρ 1.6167 1021m3 kg ρ = 1.6167 1018cm3 g ρ (6)

and one example for calculating K (for Γ = 2):

KCactus = c4 G3M2 K 6.8824 1011 m5 kg s2 K = 6.8824 104 cm5 g s2 K (7)

2 Parameters




eos_gamma
Scope: restricted  REAL



Description: Adiabatic Index for Ideal Fluid



Range   Default: 2.0
:






eos_k
Scope: restricted  REAL



Description: Polytropic constant



Range   Default: 80.0
:






gamma_ini
Scope: restricted  REAL



Description: Polytropic Gamma used for the initial model (e.g., by RNSID)



Range   Default: 2.0
:






use_cgs
Scope: restricted  BOOLEAN



Description: Use the CGS units



  Default: no



3 Interfaces

General

Implements:

eos_2d_polytrope

Inherits:

eos_base

Uses header:

EOS_Base.h

EOS_Base.inc

4 Schedule

This section lists all the variables which are assigned storage by thorn EinsteinEOS/EOS_Polytrope. 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

NONE

Scheduled Functions

CCTK_STARTUP

  eos_polytrope_startup

  setup the polytropic eos

 

 Language:fortran
  Type: function