diff --git a/src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp b/src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp index 837354c7..c3b90c89 100644 --- a/src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp +++ b/src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp @@ -1368,14 +1368,14 @@ void HelmholtzEOSMixtureBackend::T_phase_determination_pure_or_pseudopure(int ot if (other == iP) { - if (value > HEOS.SatL->p()*(1e-10 + 1)){ + if (value > HEOS.SatL->p()*(1e-8 + 1)){ this->_phase = iphase_liquid; _Q = -1000; return; } - else if (value < HEOS.SatV->p()*(1 - 1e-10)){ + else if (value < HEOS.SatV->p()*(1 - 1e-8)){ this->_phase = iphase_gas; _Q = 1000; return; } else{ - throw ValueError(format("Saturation pressure [%g Pa] corresponding to T [%g K] is within 1e-8 %% of given p [%Lg Pa]", HEOS.SatL->p(), _T, value)); + throw ValueError(format("Saturation pressure [%g Pa] corresponding to T [%g K] is within 1e-6 %% of given p [%Lg Pa]", HEOS.SatL->p(), _T, value)); } } diff --git a/src/Tests/CoolProp-Tests.cpp b/src/Tests/CoolProp-Tests.cpp index 6accd0cc..1c36dfe0 100644 --- a/src/Tests/CoolProp-Tests.cpp +++ b/src/Tests/CoolProp-Tests.cpp @@ -704,7 +704,7 @@ TEST_CASE("Tests for solvers in P,T flash using Water", "[flash],[PT]") } SECTION("Subcritical p slightly subcooled should be ok") { - double Ts, ps, rho, dT = 1e-8; + double Ts, ps, rho, dT = 1e-4; CHECK_NOTHROW(Ts = PropsSI("T","P",101325,"Q",0,"Water")); CAPTURE(Ts); CHECK(ValidNumber(Ts)); @@ -715,7 +715,7 @@ TEST_CASE("Tests for solvers in P,T flash using Water", "[flash],[PT]") } SECTION("Subcritical p slightly superheated should be ok") { - double Ts, ps, rho, dT = 1e-8; + double Ts, ps, rho, dT = 1e-4; CHECK_NOTHROW(Ts = PropsSI("T","P",101325,"Q",0,"Water")); CAPTURE(Ts); CHECK(ValidNumber(Ts)); @@ -755,7 +755,7 @@ TEST_CASE("Tests for solvers in P,Y flash using Water", "[flash],[PH],[PS],[PU]" ss2 << "Subcritical barely superheated P," << ykey; SECTION(ss2.str(), "") { - double dT = 1e-6; + double dT = 1e-3; CHECK_NOTHROW(Ts=PropsSI("T","P",101325,"Q",0,"Water")); CHECK(ValidNumber(Ts)); CAPTURE(Ts); @@ -789,7 +789,7 @@ TEST_CASE("Tests for solvers in P,Y flash using Water", "[flash],[PH],[PS],[PU]" ss4 << "Subcritical barely subcooled P," << ykey; SECTION(ss4.str(), "") { - double dT = -1e-6; + double dT = -1e-3; CHECK_NOTHROW(Ts=PropsSI("T","P",101325,"Q",0,"Water")); CHECK(ValidNumber(Ts)); CAPTURE(Ts);