PLUTO
|
Go to the source code of this file.
Macros | |
#define | C_IONS 3 /* in [1,5] */ |
#define | N_IONS 3 /* in [1,5] */ |
#define | O_IONS 3 /* in [1,5] */ |
#define | Ne_IONS 3 /* in [1,5] */ |
#define | S_IONS 3 /* in [1,5] */ |
#define | Fe_IONS 0 /* in [0,3] */ |
#define | C_EXPAND(a, b, c, d, e) ,a,b,c |
#define | N_EXPAND(a, b, c, d, e) ,a,b,c |
#define | O_EXPAND(a, b, c, d, e) ,a,b,c |
#define | Ne_EXPAND(a, b, c, d, e) ,a,b,c |
#define | S_EXPAND(a, b, c, d, e) ,a,b,c |
#define | Fe_EXPAND(a, b, c) |
#define | NIONS (3+C_IONS+N_IONS+O_IONS+Ne_IONS+S_IONS+Fe_IONS) |
Enumerations | |
enum | { X_HI = NFLX, X_HeI, Ne_EXPAND } |
Functions | |
double | GetMaxRate (double *, double *, double) |
double | H_MassFrac (void) |
double | CompEquil (double, double, double *) |
double | find_N_rho () |
void | Radiat (double *, double *) |
void | NormalizeIons (double *) |
anonymous enum |
Enumerator | |
---|---|
X_HI | |
X_HeI | |
Ne_EXPAND |
Definition at line 109 of file cooling.h.
double CompEquil | ( | double | N, |
double | T, | ||
double * | v | ||
) |
[in] | n | the particle number density (not needed, but kept for compatibility) |
[in] | T | the temperature (in K) for which equilibrium must be found. |
[in,out] | v | an array of primitive variables. On input, only density needs to be defined. On output, fractions will be updated to the equilibrium values. |
Compute the equilibrium ionization balance for (rho,T)
[in] | n | the particle number density (not needed, but kept for compatibility) |
[in] | T | the temperature (in K) for which equilibrium must be found. |
[in,out] | v | an array of primitive variables. On input, only density needs to be defined. On output, fractions will be updated to the equilibrium values. |
Compute the equilibrium ionization balance for (rho,T)
Definition at line 44 of file comp_equil.c.
double find_N_rho | ( | ) |
Definition at line 694 of file radiat.c.
Return an estimate of the maximum rate (dimension 1/time) in the chemical network. This will serve as a "stiffness" detector in the main ode integrator.
For integration to be carried explicitly all the time, return a small value (1.e-12).
double H_MassFrac | ( | void | ) |
Compute the mass fraction X of Hydrogen as function of the composition of the gas.
f_H A_H
X = -------------— f_k A_k
where
f_k : is the fractional abundance (by number), f_k = N_k/N_tot of atomic species (no matter ionization degree).
A_K : is the atomic weight
where N_{tot} is the total number density of particles
ARGUMENTS
none
Compute the mass fraction X of Hydrogen as function of the composition of the gas.
f_H A_H
X = -------------— f_k A_k
where
f_k : is the fractional abundance (by number), f_k = N_k/N_tot of atomic species (no matter ionization degree).
A_K : is the atomic weight
Note: In this module, f_H = 1.0
where N_{tot} is the total number density of particles
ARGUMENTS
none
Definition at line 721 of file radiat.c.
void NormalizeIons | ( | double * | ) |
void Radiat | ( | double * | v, |
double * | rhs | ||
) |
Cooling for optically thin plasma up to about 200,000 K Plasma composition: H, HeI-II, CI-V, NI-V, OI-V, NeI-V, SI-V Assumed abundances in elem_ab Uses S : Array = Variables vector x line points rhs : output for the system of ODE ibeg, iend : begin and end points of the current line
Cooling for neutral or singly ionized gas: good up to about 35,000 K in equilibrium or shocks in neutral gas up to about 80 km/s. Assumed abundances in ab Uses t : Kelvin dene : electron density cm*-3 fneut : hydrogen neutral fraction (adimensionale) ci,cr : H ionization and recombination rate coefficients
em(1) = TOTAL EMISSIVITY : (ergs cm**3 s**-1) em(2) = Ly alpha + two photon continuum: Aggarwal MNRAS 202, 10**4.3 K em(3) = H alpha: Aggarwal, Case B em(4) = He I 584 + two photon + 623 (all n=2 excitations): Berrington &Kingston,JPB 20 em(5) = C I 9850 + 9823: Mendoza, IAU 103, 5000 K em(6) = C II, 156 micron: Mendoza, 10,000 K em(7) = C II] 2325 A: Mendoza, 15,000 K em(8) = N I 5200 A: Mendoza, 7500 K em(9) = N II 6584 + 6548 A: Mendoza em(10) = O I 63 micron: Mendoza,2500 K em(11) = O I 6300 A + 6363 A: Mendoza, 7500 K em(12) = O II 3727: Mendoza em(13) = Mg II 2800: Mendoza em(14) = Si II 35 micron: Dufton&Kingston, MNRAS 248 em(15) = S II 6717+6727: Mendoza em(16) = Fe II 25 micron: Nussbaumer&Storey em(17) = Fe II 1.6 micron em(18) = thermal energy lost by ionization em(19) = thermal energy lost by recombination (2/3 kT per recombination. The ionization energy lost is not included here.
Provide r.h.s. for tabulated cooling.
Definition at line 94 of file radiat.c.