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

Make successive substitution a bit more conservative; closes #528

Browse Source
This commit is contained in:
Ian Bell
2015-03-05 23:08:14 -07:00
parent 01dfa4cdf0
commit 5832505ffc
1 changed files with 11 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
src/Backends/Helmholtz/VLERoutines.cpp
View File

@@ -1151,14 +1151,21 @@ void SaturationSolvers::successive_substitution(HelmholtzEOSMixtureBackend &HEOS
double dfdK = K[i]*z[i]/pow(1-beta+beta*K[i],(int)2);
df += dfdK*(deriv_liq-deriv_vap);
}
change = -f/df;
double omega = 1.0;
if (options.sstype == imposed_p){
T += change;
if (std::abs(change) > 0.05*T){
omega = 0.1;
}
T += omega*change;
}
else if (options.sstype == imposed_T){
p += change;
if (std::abs(change) > 0.05*p){
omega = 0.1;
}
p += omega*change;
}
x_and_y_from_K(beta, K, z, x, y);