TEST_F(TransportFromScratch, thermalConductivityMulti)
{
Transport* trRef = newTransportMgr("Multi", ref.get());
MultiTransport trTest;
trTest.init(test.get());
for (int i = 0; i < 10; i++) {
double T = 300 + 111*i;
ref->setState_TPX(T, 5e5, "H2:0.5, O2:0.3, H2O:0.2");
test->setState_TPX(T, 5e5, "H2:0.5, O2:0.3, H2O:0.2");
EXPECT_DOUBLE_EQ(trRef->thermalConductivity(),
trTest.thermalConductivity()) << "T = " << T;
}
}
TEST_F(TransportFromScratch, multiDiffCoeffs)
{
Transport* trRef = newTransportMgr("Multi", ref.get());
MultiTransport trTest;
trTest.init(test.get());
size_t K = ref->nSpecies();
Array2D Dref(3,3);
Array2D Dtest(3,3);
ref->setState_TPX(400, 5e5, "H2:0.5, O2:0.3, H2O:0.2");
test->setState_TPX(400, 5e5, "H2:0.5, O2:0.3, H2O:0.2");
trRef->getMultiDiffCoeffs(K, &Dref(0,0));
trTest.getMultiDiffCoeffs(K, &Dtest(0,0));
for (size_t i=0; i < K; i++) {
for (size_t j=0; j < K; j++) {
EXPECT_DOUBLE_EQ(Dref(i,j), Dtest(i,j)) << "i = " << i << ", j = " << j;
}
}
}
