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

A bit more documentation for the DQ solver; see #822

Browse Source
This commit is contained in:
Ian Bell
2015-10-23 13:39:29 +02:00
parent 37f75ddc23
commit 66228a4f2f
1 changed files with 3 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/Backends/Helmholtz/FlashRoutines.cpp
View File

@@ -264,6 +264,7 @@ public:
HEOS.update(QT_INPUTS, 0, T); // Doesn't matter whether liquid or vapor, we are just doing a full VLE call for given T
double rhoL = HEOS.saturated_liquid_keyed_output(iDmolar);
double rhoV = HEOS.saturated_vapor_keyed_output(iDmolar);
/// Error between calculated and target vapor quality based on densities
return (1/rhomolar - 1/rhoL)/(1/rhoV - 1/rhoL) - Q_target;
}
double deriv(double T){
@@ -280,8 +281,9 @@ void FlashRoutines::DQ_flash(HelmholtzEOSMixtureBackend &HEOS)
options.use_logdelta = false;
HEOS.specify_phase(iphase_twophase);
if (HEOS.is_pure_or_pseudopure){
// Bump the temperatures to hopefully yield more reliable results
double Tmax = HEOS.T_critical() - 0.1;
double Tmin = HEOS.Tmin();
double Tmin = HEOS.Tmin() + 0.1;
std::string errstr;
DQ_flash_residual resid(HEOS, HEOS._rhomolar, HEOS._Q);
Brent(resid, Tmin, Tmax, DBL_EPSILON, 1e-10, 100, errstr);