github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-04-23 03:00:17 -04:00
Code Issues Packages Projects Releases Wiki Activity

If a non-liquid liquid solution is found, bump up the liquid density and try again

Browse Source 
Ok, this isn't perfect, but it seems to work pretty well.

Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian Bell
2014-09-15 16:18:40 +02:00
parent c0157c0814
commit 731bd0f345
1 changed files with 6 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
src/Backends/Helmholtz/HelmholtzEOSMixtureBackend.cpp
View File

@@ -1702,13 +1702,12 @@ long double HelmholtzEOSMixtureBackend::solver_rho_Tp(long double T, long double
if (!ValidNumber(rhomolar)){
throw ValueError();
}
// double rr = this->rhomolar();
// for (double rho = rr*0.1; rho < 1.1*rr; rho += 100){
// specify_phase(iphase_gas);
// this->update(DmolarT_INPUTS, rho, T);
// unspecify_phase();
// std::cout << format("%g %g\n", rho, this->p());
// }
if (phase == iphase_liquid && !is_pure_or_pseudopure && first_partial_deriv(iP, iDmolar, iT) < 0){
// Try again with a larger density in order to end up at the right solution
rhomolar = Newton(resid, rhomolar_guess*1.5, 1e-8, 100, errstring);
return rhomolar;
}
return rhomolar;
}
catch(std::exception &)