PLUTO
4.4-patch2
|
PVTE_LAW for a partially ionized gas. More...
#include "pluto.h"
Functions | |
static double | SahaXFrac (double T, double rho) |
double | InternalEnergyFunc (double *v, double T) |
void | GetMu (double T, double rho, double *mu) |
double | Gamma1 (double *v) |
Compute the internal energy for a partially ionized hydrogen gas.
Reference
double Gamma1 | ( | double * | v | ) |
Calculate the value of the first adiabatic index:
where p
and rho
are in c.g.s units. Note that if species are evolved explicitly (non-equilibrium chemistry), we set chi=1
.
The heat capacity at constant volume, cV
, is defined as the derivative of specific internal energy with respect to temperature:
and it is computed numerically using a centered derivative.
This function is needed (at present) only when computing the sound speed in the Riemann solver. Since this is only needed for an approximated value, 5/3 (upper bound) should be ok.
[in] | v | 1D array of primitive quantities |
void GetMu | ( | double | T, |
double | rho, | ||
double * | mu | ||
) |
Calculate the mean molecular weight for the case in which hydrogen fractions are estimated using Saha Equations.
[in] | T | Gas temperature in Kelvin. |
[in] | rho | Gas density (code units) |
[out] | mu | Mean molecular weight |
double InternalEnergyFunc | ( | double * | v, |
double | T | ||
) |
Compute the gas internal energy as a function of temperature and fractions (or density):
rhoe = rhoe(T,rho)
in LTE or CIE; rhoe = rhoe(T,X)
in non-equilibrium chemistry.[in] | v | 1D Array of primitive variables containing density and species. Other variables are ignored. |
[in] | T | Gas temperature |
rhoe
) in code units.
|
static |
Use Saha equation to compute the degree of ionization.