Maximum temp for PY flash is now 1.5*Tmax

Closes https://github.com/CoolProp/CoolProp/issues/137

Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian Bell
2014-09-10 16:24:36 +02:00
parent a2a2c80323
commit 3844f281fd

View File

@@ -725,7 +725,7 @@ void FlashRoutines::HSU_P_flash(HelmholtzEOSMixtureBackend &HEOS, parameters oth
{
case iphase_gas:
{
Tmax = HEOS.Tmax()+1;
Tmax = 1.5*HEOS.Tmax();
if (saturation_called){ Tmin = HEOS.SatV->T();}else{Tmin = HEOS._TVanc.pt();}
break;
}
@@ -746,7 +746,7 @@ void FlashRoutines::HSU_P_flash(HelmholtzEOSMixtureBackend &HEOS, parameters oth
case iphase_supercritical_gas:
case iphase_supercritical:
{
Tmax = HEOS.Tmax()+1;
Tmax = 1.5*HEOS.Tmax();
// Sometimes the minimum pressure for the melting line is a bit above the triple point pressure
if (HEOS.has_melting_line() && HEOS._p > HEOS.calc_melting_line(iP_min, -1, -1)){
Tmin = HEOS.calc_melting_line(iT, iP, HEOS._p)-1e-3;