alpha_g max(min(s2,1),0) "Bounded gas volume fraction"
alpha_l 1-alpha_g "Liquid volume fraction"
X_g alpha_g*rho_g/spf.rho "Gas mass fraction"
X_l alpha_l*rho_l/spf.rho "Liquid mass fraction"
phi_evap max(0,(alpha_geq-alpha_g)/t_evap) Evaporation
phi_cond min(0,(alpha_geq-alpha_g)/t_cond) Condensation
alpha_geq X_geq^(n_const)/(X_geq+(1-X_geq)*rho_g/rho_l) "Volume fraction of gas at equilibrium"
Delta_Tsub T_sat-T "Subcooled temperature difference"
Delta_Tsup T-T_sat "Superheated temperature difference"
t_evap K_evap*T_sat/max(Delta_Tsup,eps) "Time constant, evaporation"
t_cond K_cond*T_sat/max(Delta_Tsub,eps) "Time constant, condensation"
h_mix if(T>T_sat,h_l_sat+X_g*Delta_Hlg+ht.Cp*(Delta_Tsup),h_l) "Enthalpy of mixture"
h_l_sat cp_l*(T_sat-Tref) "Saturated liquid enthalpy"
X_geq max(h_mix-h_l_sat,0)/Delta_Hlg "Mass fraction at equilibrium"
h_l cp_l*(min(T,T_sat)-Tref) "Liquid enthalphy"
m_evap phi_evap*rho_g "Mass source term from evaporation"
m_cond phi_cond*rho_g "Mass source term from condensation"