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

Added more cropping for better testing

Browse source
This commit is contained in:
Harry Moffat 2008-10-16 14:24:51 +00:00
parent 8b80bdff0e
commit f8461afd5f
2 changed files with 9 additions and 5 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
test_problems/surfSolverTest/results2_blessed.txt
View file

@ -32,13 +32,13 @@ Surface Phase: soot_interface (11)
Surface Temperature = 1.4e+03 Surface Temperature = 1.4e+03
Surface Pressure = 1.01e+05 Surface Pressure = 1.01e+05
Name Coverage SrcRate Name Coverage SrcRate
0 Csoot-* 0.0184677 -8.195044e-20 0 Csoot-* 0.0184677 0.000000e+00
1 Csoot-H 0.981532 -8.195044e-20 1 Csoot-H 0.981532 0.000000e+00
Sum of coverages = 1 Sum of coverages = 1
Surface Phase: soot_interface (11) Surface Phase: soot_interface (11)
Surface Temperature = 1.4e+03 Surface Temperature = 1.4e+03
Surface Pressure = 1.01e+05 Surface Pressure = 1.01e+05
Name Coverage SrcRate Name Coverage SrcRate
0 Csoot-* 0.0184677 -8.195044e-20 0 Csoot-* 0.0184677 0.000000e+00
1 Csoot-H 0.981532 -8.195044e-20 1 Csoot-H 0.981532 0.000000e+00
Sum of coverages = 1 Sum of coverages = 1

6
test_problems/surfSolverTest/surfaceSolver2.cpp
View file

@ -145,9 +145,13 @@ void printSurf(ostream &oooo, ThermoPhase *surfPhaseTP,
int nspSurf = surfPhaseTP->nSpecies(); int nspSurf = surfPhaseTP->nSpecies();
for (int k = 0; k < nspSurf; k++) { for (int k = 0; k < nspSurf; k++) {
kstart = iKin_ptr->kineticsSpeciesIndex(0, 2); kstart = iKin_ptr->kineticsSpeciesIndex(0, 2);
double srcK = src[kstart];
if (fabs(srcK) < 1.0E-8) {
srcK = 0.0;
}
sprintf(buf, "%4d %24s %14g %14e\n", sprintf(buf, "%4d %24s %14g %14e\n",
k, surfPhaseTP->speciesName(k).c_str(), k, surfPhaseTP->speciesName(k).c_str(),
x[k], src[kstart]); x[k], srcK);
oooo << buf; oooo << buf;
sum += x[k]; sum += x[k];
} }