next up previous contents index
Next: 7.1.4 Boundary conditions (general) Up: 7.1 Parameter file: rhd.par Previous: 7.1.2 Header   Contents   Index


7.1.3 Fundamental model parameters real teff

The effective temperature is one of the basic model parameters and is specified e.g. with
real teff  f=F13.3 b=4  n='Effective Temperature' u=K
(for a relatively cool star). Note that this parameter specifies the average output flux only for local models with closed radiative lower boundary and that in many cases the actual effective temperature can only be determined a posteriori. For other lower (or inner boundary conditions) the entropy of the instreaming material s_inflow (see below) is more important than teff itself. See Sect. for different heating modes controled with heat_mode. character grav_mode

Gravity is another characteristic of a stellar atmosphere. The type (or geometry) of the external gravity field has to be specified e.g. with
character grav_mode  f=A80 b=80  n='Type of gravity field' &
Three values are possible so far: real grav

In the case of a constant gravity the amount of the acceleration has to specified with
real grav  f=E15.8 b=4  n='Gravity' u=cm/s^2
Setting this value to zero switches off gravity. real mass_star

In the case of a central the mass (in cgs units) of the star has to be specified with
real mass_star f=E15.8 b=4 n='Stellar Mass'              u=g
 9.94500e+33 real r0_grav

To avoid the central singularity in a 1/r potential it is smoothed in the center to give a central potential of 1/r0_grav, specified with
real r0_grav   f=E15.8 b=4 n='Inner Smoothing Radius'    u=cm
This parameter should always be non-zero for a central potential of a global model. For local models, this parameter specifies the depth of the lower heating and damping region, see Sect. real r1_grav

The density in an atmosphere in hydrostatic equilibrium can decline to very low values. To artificially enlarge the pressure (and density) scale height in the outer layers of the star (the corners of the box) the gravity can be reduced by defining the potential at infinity to be 1/r1_grav, specified with
real r1_grav   f=E15.8 b=4 n='Outer Smoothing Radius'    u=cm &
  c0='0.0: Not used'
Setting this parameter to zero gives the usual $1/r$ behavior of the potential in the outer layers but also chooses another smoothing formula in the central part (where r0_grav is relevant). But a value somewhat larger than the remotest corner of the box effectively cancels this artificial smoothing in the outer layers. real r2_grav

Using r2_grav instead of r1_grav means that a potential function more appropriate for a polytropic stellar interior model is used. It can be set, e.g. with
real r2_grav   f=E15.8 b=4 n='Outer Smoothing Radius'    u=cm &
  c0='0.0: Not used'
Setting this parameter to zero means the potential parameterized with r1_grav is used. This parameter is similar to r1_grav, i.e., it enlarges the pressure scale height in the outer layers but does not change the smoothing formula in the center. This parameter is only effective if r1_grav is set to zero. real r0_core

To insert energy in a sphere different with a radius other than r0_grav (the default), the heating radius r0_core can be specified separately, e.g. with
real r0_core   f=E15.8 b=4 n='Core Radius'    u=cm
If this parameter is not specified or set to zero r0_grav is used as the radius of the core. character centrifugal_force

Usually, a centrifugal force is applied if nu_rotation$\ne$0. To switch it off even for a non-zero rotation rate, the parameter centrifugal_force can be used, e.g. by setting
character centrifugal_force  f=A80 b=80  n='Switch on/off the centrifugal force' &
  c0='on: default, off: no centrifugal force, even for non-zero nu_rotation'
Two alternative values are possible so far: By default, the centrifugal force is activated if the rotation rate nu_rotation is non-zero, but only for global models as the centrifugal force for local models would in effect only modify the value of the gravity constant. real nu_rotation

To transform onto a rotating coordinate-system a rotation rate can be specified with e.g.
real nu_rotation   f=E15.8 b=4 n='Rotation frequency'    u=1/s
The potential is modified by adding terms due to a centrifugal force unless centrifugal_force is set to off. In addition, a Coriolis force is applied during the hydrodynamics step. real ar_rotationaxis

The rotation axis can be specified with
real ar_rotationaxis f=E10.4 b=4 p=1 d=(1:3) n='Rotation axis' u=1
The default value is (0.0,0.0,1.0) and should be used for rotating global models. A different axis only makes sense for local models.

next up previous contents index
Next: 7.1.4 Boundary conditions (general) Up: 7.1 Parameter file: rhd.par Previous: 7.1.2 Header   Contents   Index