next up previous contents index
Next: 7.1.13 Dust: dustscheme=dust_moment04_c2 Up: 7.1 Parameter file: rhd.par Previous: 7.1.11 Tensor-viscosity control   Contents   Index

Subsections


7.1.12 Dust/molecules/hydrogen ionization: General

CO5BOLD can handle a number of additional density arrays. They can be used to describe e.g. the mass density of dust distribution moments or number densities of molecules. These species are properly advected with the gas density. There is also already a small number of dust/molecule formation models available. These models have to be improved in the future and the influence on the radiation field (opacities, radiation pressure on dust) has to be taken into account.


7.1.12.1 character dustscheme

A scheme for dust or molecule formation and transport can be selected e.g. with
character dustscheme f=A80 b=80 n='Dust model' &
  c0='none (default), nosource, dust_simple_01, co_component01_01' &
  c1='dust_k3mon_01, dust_k3mon_02'
dust_k3mon_01
The following values are currently possible:


7.1.12.2 character dustreconstruction

While usually (i.e., when this parameter is not set) the reconstruction scheme for the additional density fields is the same as for the normal hydrodynamics quantities (see 7.1.8.4), it can make sense to choose a strictly monotonic scheme for these additional fields. This scheme can be more (e.g. PP) or less (e.g. PPmimo) aggressive. It can be chosen e.g. with
character dustreconstruction f=A80 b=80 n='Reconstruction method for qucs' &
  c0=Constant c1=Minmod/VanLeer/Superbee c2=PP/PPmimo &
  c3=FRmimo,FRmono
FRmono
Useful combinations are

7.1.12.3 real c_dust0X

There are ten parameters (real c_dust01 to c_dust10) to control each dust formation scheme in detail. A parameter can be given as in
real c_dust01 f=E15.8 b=4 n='Dust parameter 1'
0.0
The meaning (and unit) can vary from scheme to scheme. The default value is 0.0 in each case.
Important: The parameter real c_dust01 must be set to 1.0 in order to activate advection of particle densities for the CHEM (chemreacnet) and the HION (hiontd) module. The default value is 0.0, i.e. advection is switched off.

7.1.12.4 character chem_reacfile

The name of the input file containing the chemical reaction network.
character chem_reacfile f=A80 b=80 n='file name of reaction table'
chem.dat

7.1.12.5 character chem_reacpath

The path of the input file containing the chemical reaction network.
character chem_reacpath f=A80 b=80 n='path of reaction table'
/data/sven/cobold/dat/chem/

7.1.12.6 real chem_abumetal

Abundance of the representative metal ('M', if present) relative to hydrogen ( $\epsilon_\mathrm{M} = [M] / [H]$). A value of 1.0E-04 means there are $10^4$ hydrogen atoms for every metal atom. The metal number density for each grid cell is then derived via $n_\mathrm{H} = \epsilon_\mathrm{M}  \rho / m_\mathrm{H}$ (assuming a pure hydrogen gas). This parameter is ignored when a quc array for the metal is found in the input model.
real chem_abumetal f=E15.8 b=4 n='Chemical abundance of repres. metal'
1.0E-04

7.1.12.7 character hion_datapath

The path of all input files for HION.
character hion_datapath f=A80 b=80 n='HION data path'
/data/sven/cobold/dat/hion/

7.1.12.8 character hion_atomfile

The file name of the model atom for HION.
character hion_atomfile f=A80 b=80 n='HION atom file name'
H_6.atom

7.1.12.9 character hion_abufile

The name of the HION input file containing the chemical abundances.
character hion_abufile f=A80 b=80 n='HION abundance file name'
abundance.input

7.1.12.10 character hion_edensfile

The name of the HION input file containing the electron density table.
character hion_edensfile f=A80 b=80 n='HION electron density file name'
edens.dat

7.1.12.11 character hion_pffile

The name of the HION input file containing the partition functions.
character hion_pffile f=A80 b=80 n='HION partition function file name'
pf_kurucz.dat

7.1.12.12 real dtime_out_hion

Time increment for additional HION output. Positive values specify the time increment in seconds, negative values the increment in computational time steps. Setting the parameter to zero, suppresses the output.
real dtime_out_hion f=E15.8 b=4 n='Output file time step (HION)'
-10.0

7.1.12.13 integer hion_chunks

Number of chunks to use for HION in order to limit the required memory. Do not make it bigger than the number of points in the x2 direction.
integer hion_chunks f=I9 b=4 n='Number of HION chunks'
1


next up previous contents index
Next: 7.1.13 Dust: dustscheme=dust_moment04_c2 Up: 7.1 Parameter file: rhd.par Previous: 7.1.11 Tensor-viscosity control   Contents   Index