github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-01-09 14:08:08 -05:00
Code Issues Packages Projects Releases Wiki Activity

Wetbulb solver T,B,P works better now; Closes #419 and #427

Browse Source 
Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian Bell
2015-02-09 22:03:35 -07:00
parent a02d973aad
commit f81ac8fec2
1 changed files with 3 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/HumidAirProp.cpp
View File

@@ -285,11 +285,13 @@ static double Secant_HAProps_W(const std::string &OutputName, const std::string
// Use a secant solve in order to yield a target output value for HAProps by altering humidity ratio
double x1=0,x2=0,x3=0,y1=0,y2=0,eps=1e-12,f=999,W=0.0001;
int iter=1;
if (OutputName == "B"){eps = 1e-7;}
while ((iter<=3 || std::abs(f)>eps) && iter<100)
{
if (iter == 1){x1 = W_guess; W = x1;}
if (iter == 2){x2 = W_guess*1.001; W = x2;}
if (iter == 2){x2 = W_guess*1.1; W = x2;}
if (iter > 2) {W = x2;}
f = HAPropsSI(OutputName,"W",W,Input1Name,Input1,Input2Name,Input2)-TargetVal;
if (iter == 1){y1 = f;}