github/CoolProp
1
1
Fork 0
mirror of https://github.com/CoolProp/CoolProp.git synced 2026-02-09 13:25:12 -05:00
Code Issues Packages Projects Releases Wiki Activity

Added failing tests for same reference state between REFPROP and Coolprop for enthalpy and entropy

Browse Source 
Signed-off-by: Ian Bell <ian.h.bell@gmail.com>
This commit is contained in:
Ian Bell
2014-06-03 12:54:57 +02:00
parent b78439b72c
commit 17bffeb91f
1 changed files with 43 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

43
src/Backends/REFPROP/REFPROPMixtureBackend.cpp
View File

@@ -1206,3 +1206,46 @@ void REFPROPMixtureBackend::update(long input_pair, double value1, double value2
}
} /* namespace CoolProp */
#ifdef ENABLE_CATCH
#include "CoolProp.h"
#include "catch.hpp"
TEST_CASE("Check REFPROP H,S reference states equal to CoolProp","[REFPROP]")
{
std::vector<std::string> ss = strsplit(CoolProp::get_global_param_string("FluidsList"),',');
for (std::vector<std::string>::iterator it = ss.begin(); it != ss.end(); ++it)
{
std::string Name = (*it);
std::string RPName = CoolProp::get_fluid_param_string((*it),"REFPROP_name");
// Skip fluids not in REFPROP
if (RPName.find("N/A") == 0){continue;}
std::tr1::shared_ptr<CoolProp::AbstractState> S1(CoolProp::AbstractState::factory("HEOS", (*it)));
double Tr = S1->T_critical();
S1->update(CoolProp::QT_INPUTS, 0, Tr*0.9);
double h_CP = S1->hmass();
double s_CP = S1->smass();
std::tr1::shared_ptr<CoolProp::AbstractState> S2(CoolProp::AbstractState::factory("REFPROP", RPName));
S2->update(CoolProp::QT_INPUTS, 0, Tr*0.9);
double h_RP = S2->hmass();
double s_RP = S2->smass();
CAPTURE(Name);
CAPTURE(RPName);
CAPTURE(h_CP);
CAPTURE(h_RP);
CAPTURE(s_CP);
CAPTURE(s_CP);
double DH = (S1->hmass()-S2->hmass())/S1->cpmass();
double DS = (S1->smass()-S2->smass())/S1->cpmass();
CHECK(fabs(DH) < 1e-3);
CHECK(fabs(DS) < 1e-3);
}
}
#endif