61 lines
1.3 KiB
C++
61 lines
1.3 KiB
C++
/*
|
|
* $Author: hkmoffa $
|
|
* $Date: 2008/02/16 21:44:47 $
|
|
* $Revision: 1.5 $
|
|
*
|
|
*/
|
|
|
|
#include "cantera/Cantera.h"
|
|
#include "cantera/IdealGasMix.h"
|
|
#include "cantera/equilibrium.h"
|
|
|
|
using namespace std;
|
|
using namespace Cantera;
|
|
using namespace Cantera_CXX;
|
|
|
|
int main(int argc, char **argv) {
|
|
try {
|
|
IdealGasMix g("red1.xml", "gri30_mix");
|
|
|
|
#ifdef DEBUG_BASISOPTIMIZE
|
|
Cantera::BasisOptimize_print_lvl = 0;
|
|
#endif
|
|
#ifdef DEBUG_CHEMEQUIL
|
|
Cantera::ChemEquil_print_lvl = 0;
|
|
#endif
|
|
|
|
double pres = 1.0E5;
|
|
g.setState_TPX(2000.0, pres, "C2H2:0.9, CH:0.1");
|
|
|
|
MultiPhase mphase;
|
|
mphase.addPhase(&g, 10.0);
|
|
mphase.init();
|
|
int usedZeroedSpecies = 0;
|
|
vector_int orderVectorSpecies;
|
|
vector_int orderVectorElements;
|
|
|
|
bool doFormMatrix = true;
|
|
vector_fp formRxnMatrix;
|
|
|
|
int nc = BasisOptimize(&usedZeroedSpecies, doFormMatrix,
|
|
&mphase, orderVectorSpecies,
|
|
orderVectorElements,
|
|
formRxnMatrix);
|
|
|
|
cout << "number of components = " << nc << endl;
|
|
|
|
/*
|
|
* The ChemEquil solver throws an error for this case.
|
|
* The MultiPhaseEquil solver just gets the wrong result.
|
|
*/
|
|
equilibrate(g, "TP", -1);
|
|
cout << g;
|
|
|
|
return 0;
|
|
}
|
|
catch (CanteraError) {
|
|
showErrors(cerr);
|
|
cerr << "program terminating." << endl;
|
|
return -1;
|
|
}
|
|
}
|