research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

[Test] Add an extra test of the equilibrium solvers

Browse source 
Testing the TP property pair helps identify whether issues arise from that
algorithm or the wrapper algorithms for other property pairs
This commit is contained in:
Ray Speth 2017-08-24 11:05:03 -04:00
parent 3820ca3147
commit b92bd24f92
1 changed files with 12 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
test/equil/equil_gas.cpp
View file

@ -224,6 +224,15 @@ TEST_F(GriMatrix, VcsNonideal_CH4_O2) { check_CH4_O2("vcs"); }
class PropertyPairs : public GriEquilibriumTest
{
public:
void check_TP(const std::string& solver) {
gas.setState_TPX(500, 1e5, "CH4:0.3, O2:0.3, N2:0.4");
save_elemental_mole_fractions();
gas.equilibrate("TP", solver);
EXPECT_NEAR(500, gas.temperature(), 1e-5);
EXPECT_NEAR(1e5, gas.pressure(), 1e-3);
check();
}
void check_HP(const std::string& solver) {
gas.setState_TPX(500, 1e5, "CH4:0.3, O2:0.3, N2:0.4");
double h0 = gas.enthalpy_mass();
@ -279,6 +288,9 @@ public:
}
};
TEST_F(PropertyPairs, ChemEquil_TP) { check_TP("element_potential"); }
TEST_F(PropertyPairs, MultiPhase_TP) { check_TP("gibbs"); }
TEST_F(PropertyPairs, VcsNonideal_TP) { check_TP("vcs"); }
TEST_F(PropertyPairs, ChemEquil_HP) { check_HP("element_potential"); }
TEST_F(PropertyPairs, MultiPhase_HP) { check_HP("gibbs"); }
TEST_F(PropertyPairs, VcsNonideal_HP) { check_HP("vcs"); }