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

Fixed problems with HAPropsSI when wetbulb is an input by converting call to Brent. Slow but steady.

Browse Source 
Closes https://github.com/CoolProp/CoolProp/issues/311

Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian Bell
2014-12-09 00:08:27 -05:00
parent 1d1082d2be
commit e0e2591d73
1 changed files with 8 additions and 25 deletions
Download Patch File Download Diff File Expand all files Collapse all files

33
src/HumidAirProp.cpp
View File

@@ -1361,31 +1361,14 @@ double HAPropsSI(const std::string &OutputName, const std::string &Input1Name, d
double T_min = 210;
double T_max = 450;
T = -1;
// First try to use the secant solver to find T at a few different temperatures
for (T_guess = 210; T_guess < 450; T_guess += 60)
{
try{
T = Secant_HAProps_T(SecondaryInputName,(char *)"P",p,MainInputName,MainInputValue,SecondaryInputValue,T_guess);
double val = HAPropsSI(SecondaryInputName,(char *)"T",T,(char *)"P",p,MainInputName,MainInputValue);
if (!ValidNumber(T) || !ValidNumber(val) || !(T_min < T && T < T_max) || std::abs(val-SecondaryInputValue)>1e-3)
{
throw CoolProp::ValueError();
}
else
{
break;
}
}
catch (std::exception &){};
}
if (T < 0) // No solution found using secant
{
// Use the Brent's method solver to find T
T = Brent_HAProps_T(SecondaryInputName,(char *)"P",p,MainInputName,MainInputValue,SecondaryInputValue,T_min,T_max);
}
try{
// Use the Brent's method solver to find T. Slow but reliable
T = Brent_HAProps_T(SecondaryInputName, "P", p, MainInputName, MainInputValue, SecondaryInputValue, T_min,T_max);
}
catch(std::exception &e){
CoolProp::set_error_string(e.what());
return _HUGE;
}
// If you want the temperature, return it
if (Name2Type(OutputName)==GIVEN_T)