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

Fix Coverity CID 38566

Browse Source
This commit is contained in:
mikekaganski
2015-03-07 23:04:07 +10:00
parent 16b13f3f0f
commit c99805cc79
1 changed files with 8 additions and 11 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
src/HumidAirProp.cpp
View File

@@ -1013,19 +1013,16 @@ double DewpointTemperature(double T, double p, double psi_w)
class WetBulbSolver : public CoolProp::FuncWrapper1D
{
private:
double _T,_p,_W,LHS,RHS,v_bar_w,M_ha;
double _p,_W,LHS;
public:
WetBulbSolver(double T, double p, double psi_w){
_T = T;
_p = p;
_W = epsilon*psi_w/(1-psi_w);
WetBulbSolver(double T, double p, double psi_w)
: _p(p),_W(epsilon*psi_w/(1-psi_w))
{
//These things are all not a function of Twb
v_bar_w = MolarVolume(T,p,psi_w);
M_ha = MM_Water()*psi_w+(1-psi_w)*0.028966;
double v_bar_w = MolarVolume(T,p,psi_w),
M_ha = MM_Water()*psi_w+(1-psi_w)*0.028966;
LHS = MolarEnthalpy(T,p,psi_w,v_bar_w)*(1+_W)/M_ha;
};
~WetBulbSolver(){};
}
double call(double Twb)
{
double epsilon=0.621945;
@@ -1066,7 +1063,7 @@ public:
M_ha_wb = MM_Water()*psi_wb+(1-psi_wb)*0.028966;
v_bar_wb=MolarVolume(Twb,_p,psi_wb);
RHS = (MolarEnthalpy(Twb,_p,psi_wb,v_bar_wb)*(1+W_s_wb)/M_ha_wb+(_W-W_s_wb)*h_w);
double RHS = (MolarEnthalpy(Twb,_p,psi_wb,v_bar_wb)*(1+W_s_wb)/M_ha_wb+(_W-W_s_wb)*h_w);
if (!ValidNumber(LHS-RHS)){throw CoolProp::ValueError();}
return LHS - RHS;
}