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

Cleaned up the implementation of the NASA9poly test

Browse source
This commit is contained in:
Ray Speth 2012-03-09 22:54:09 +00:00
parent 4534819d1d
commit 66ada1c423
1 changed files with 18 additions and 52 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

70
test_problems/NASA9poly_test/NASA9poly_test.cpp
View file

@ -2,32 +2,15 @@
* @file NASA9poly_test
* test problem for NASA 9 coefficient formulation
*/
#include <iostream>
#include <string>
#include <vector>
#include <string>
#include <iomanip>
using namespace std;
/*****************************************************************/
/*****************************************************************/
#include "cantera/transport.h"
#include "cantera/IdealGasMix.h"
#include "cantera/transport/TransportFactory.h"
using namespace Cantera;
#include <iostream>
#include <vector>
#include <memory>
void printDbl(double val)
{
if (fabs(val) < 5.0E-17) {
cout << " nil";
} else {
cout << val;
}
}
using namespace Cantera;
int main(int argc, char** argv)
{
@ -35,62 +18,45 @@ int main(int argc, char** argv)
_set_output_format(_TWO_DIGIT_EXPONENT);
#endif
try {
//throw ArraySizeError("MadeUpFunction", 12, 22);
IdealGasMix g("gasNASA9.xml", "gri30_mix");
int nsp = g.nSpecies();
size_t nsp = g.nSpecies();
double pres = 1.0E5;
vector_fp Xset(nsp, 0.0);
Xset[0] = 0.5 ;
Xset[1] = 0.5;
g.setState_TPX(1500.0, pres, DATA_PTR(Xset));
Xset[0] = 0.5 ;
Xset[1] = 0.5;
g.setState_TPX(1500.0, pres, &Xset[0]);
vector_fp cp_R(nsp, 0.0);
g.getCp_R(DATA_PTR(cp_R));
vector_fp H_RT(nsp, 0.0);
vector_fp S_R(nsp, 0.0);
vector_fp Gvalues(nsp, 0.0);
g.getEntropy_R(&S_R[0]);
g.getCp_R(&cp_R[0]);
g.getEnthalpy_RT(&H_RT[0]);
printf("Comparisons of H2 calculated via several equivalent classes:\n");
printf("1500 K and 1 atm:\n");
printf(" NasaThermo Nasa9 Nasa9_4reg \n");
printf(" cp/R: %11.6g %11.6g %11.6g\n", cp_R[0], cp_R[1], cp_R[2]);
vector_fp H_RT(nsp, 0.0);
g.getEnthalpy_RT(DATA_PTR(H_RT));
printf(" H/RT: %11.6g %11.6g %11.6g\n", H_RT[0], H_RT[1], H_RT[2]);
vector_fp S_R(nsp, 0.0);
g.getEntropy_R(DATA_PTR(S_R));
printf(" S/R: %11.6g %11.6g %11.6g\n", S_R[0], S_R[1], S_R[2]);
Transport* tran = newTransportMgr("Mix", &g);
std::auto_ptr<Transport> tran(newTransportMgr("Mix", &g));
vector_fp Gvalues(nsp, 0.0);
// MultiTransport * tranMix = dynamic_cast<MultiTransport *>(tran);
printf("Viscoscity and thermal Cond vs. T\n");
for (int k = 0; k < 40; k++) {
double T1 = 400. + 200. * k;
g.setState_TPX(T1, pres, DATA_PTR(Xset));
g.getPureGibbs(DATA_PTR(Gvalues));
//printf(" -- %13g %13.5g %13.5g %13.5g %13.5g \n",
// Gvalues[0], Gvalues[1], Gvalues[2],
// Gvalues[3], Gvalues[4]);
g.setState_TPX(T1, pres, &Xset[0]);
g.getPureGibbs(&Gvalues[0]);
double visc = tran->viscosity();
double cond = tran->thermalConductivity();
printf(" %13g %13.5g %13.5g\n", T1, visc, cond);
}
} catch (CanteraError& err) {
err.save();
// cout << err.what() << endl;
showErrors(cout);
std::cout << err.what() << std::endl;
return 1;
}
return 0;
}
/***********************************************************/