added cantera transport headers and removed the time loop

testApp/YEqn.H

@@ -1,12 +0,0 @@
-{
-    forAll(Y, specieI)
-    {
-        volScalarField& Yi = Y[specieI];
-
-        solve
-        (
-            fvm::ddt(rho, Yi) - chemistry.RR(specieI),
-            mesh.solver("Yi")
-        );
-    }
-}

testApp/chemFoam.C

@@ -31,6 +31,9 @@ Description
 
 \*---------------------------------------------------------------------------*/
 
+#include <cantera/transport.h>
+#include <cantera/IdealGasMix.h>
+
 #include "fvCFD.H"
 #include "psiReactionThermo.H"
 #include "psiChemistryModel.H"
@@ -62,27 +65,6 @@ int main(int argc, char *argv[])
 
     Info<< "\nStarting time loop\n" << endl;
 
-    while (runTime.run())
-    {
-        #include "readControls.H"
-
-        #include "setDeltaT.H"
-
-        runTime++;
-        Info<< "Time = " << runTime.timeName() << nl << endl;
-
-        #include "solveChemistry.H"
-        #include "YEqn.H"
-        #include "hEqn.H"
-        #include "pEqn.H"
-
-        #include "output.H"
-
-        Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
-            << "  ClockTime = " << runTime.elapsedClockTime() << " s"
-            << nl << endl;
-    }
-
     Info << "Number of steps = " << runTime.timeIndex() << endl;
     Info << "End" << nl << endl;
 

testApp/hEqn.H

@@ -1,14 +0,0 @@
-{
-    volScalarField& h = thermo.he();
-
-    if (constProp == "volume")
-    {
-        h[0] = u0 + p[0]/rho[0] + integratedHeat;
-    }
-    else
-    {
-        h[0] = h0 + integratedHeat;
-    }
-
-    thermo.correct();
-}

testApp/output.H

@@ -1,11 +0,0 @@
-    runTime.write();
-
-    Info<< "Sh = " << Sh
-        << ", T = " << thermo.T()[0]
-        << ", p = " << thermo.p()[0]
-        << ", " << Y[0].name() << " = " << Y[0][0]
-        << endl;
-
-    post<< runTime.value() << token::TAB << thermo.T()[0] << token::TAB
-        << thermo.p()[0] << endl;
-

testApp/pEqn.H

@@ -1,16 +0,0 @@
-{
-    rho = thermo.rho();
-    if (constProp == "volume")
-    {
-        scalar invW = 0.0;
-        forAll(Y, i)
-        {
-            invW += Y[i][0]/specieData[i].W();
-        }
-
-        Rspecific[0] = 1000.0*constant::physicoChemical::R.value()*invW;
-
-        p[0] = rho0*Rspecific[0]*thermo.T()[0];
-        rho[0] = rho0;
-    }
-}

testApp/readControls.H

@@ -1,3 +0,0 @@
-runTime.controlDict().lookup("adjustTimeStep") >> adjustTimeStep;
-
-maxDeltaT = readScalar(runTime.controlDict().lookup("maxDeltaT"));

testApp/setDeltaT.H

@@ -1,6 +0,0 @@
-if (adjustTimeStep)
-{
-    runTime.setDeltaT(min(dtChem, maxDeltaT));
-    Info<< "deltaT = " <<  runTime.deltaT().value() << endl;
-}
-

testApp/solveChemistry.H

@@ -1,3 +0,0 @@
-    dtChem = chemistry.solve(runTime.deltaT().value());
-    scalar Sh = chemistry.Sh()()[0]/rho[0];
-    integratedHeat += Sh*runTime.deltaT().value();