pttools.models.full
Full thermodynamics-based model
Classes
- class pttools.models.full.FullModel(thermo, V_s=0, V_b=0, T_crit_guess=None, allow_invalid=False, name=None, label_latex=None, label_unicode=None)
Bases:
ModelFull thermodynamics-based equation of state
Temperature limits should be set in the ThermoModel.
- Parameters:
thermo (ThermoModel) – model of the underlying thermodynamics. Some models don’t take this, but use their own approximations instead.
V_s (float) – the constant term in the expression of \(p\) in the symmetric phase
V_b (float) – the constant term in the expression of \(p\) in the broken phase
T_crit_guess (float)
allow_invalid (bool)
name (str)
label_latex (str)
label_unicode (str)
- DEFAULT_LABEL = 'Full model'
- DEFAULT_NAME: str = 'full'
- critical_temp_opt(temp)
Optimizer function for critical temperature
- Parameters:
temp (float)
- Return type:
float
- e_temp(temp, phase)
Energy density \(e(T,\phi)\), using Borsanyi et al., 2016, eq. S12 $$ e(T,\phi) = \frac{\pi^2}{30} g_e(T,\phi) T^4 $$ :param temp: temperature \(T\) :param phase: phase \(\phi\) :return: \(e(T,\phi)\)
- ge_temp(temp, phase)
- gen_cs2()
This function generates the Numba-jitted cs2 function to be used by the fluid integrator
- gp_temp(temp, phase)
- gs_temp(temp, phase)
- p_temp(temp, phase)
Pressure \(p(T,\phi)\) $$ p(T,\phi) = \frac{\pi^2}{90} g_p(T,\phi) T^4$$
- params_str()
Model parameters as a string
- Return type:
str
- s_temp(temp, phase)
Entropy density $s(T,phi), using Borsanyi et al., 2016, eq. S12\( \)$ s(T,\phi) = \frac{2\pi^2}{45} g_s(T) T^3$$ :param temp: temperature \(T\) :param phase: phase \(\phi\) :return: \(s(T,\phi)\)
- temp(w, phase)
Temperature \(T\)
- w(temp, phase)
Enthalpy density \(w\) $$ w = e + p = Ts = T \frac{dp}{dT} = \frac{2\pi^2}{45} g_s T^4 $$ For the steps please see Hindmarsh et al., 2021 page 23 and eq. 7.1. and :borsanyi_2016: eq. S12.