research/eReactingFoam-4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

pure composition test for cantera mixture averaged diffusivity

Browse source
This commit is contained in:
Yeongdo Park 2018-11-02 09:48:49 -04:00
parent 0a89513335
commit fa0af403cf
1 changed files with 32 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

35
testApp/chemFoam.C
View file

@ -122,11 +122,40 @@ int main(int argc, char *argv[])
diff.correct();
// Cantera gas transport
gas_.setState_TPY(T[i], p0, XY.data());
gas_.setState_TPY(T[i], p[0], XY.data());
scalarField Dc (composition.species().size(), 0.0);
tr_->getMixDiffCoeffsMass(Dc.data());
scalarField Dij (composition.species().size()*composition.species().size(), 0.0);
const scalar relTol = 1.0e-15;
// check for singular case (Xi = 1)
label idxFracOne = -1;
scalar sumX = sum(XY);
forAll (XY, i)
{
if ((sumX-XY[i])/sumX < relTol)
{
idxFracOne = i;
break;
}
}
if (idxFracOne < 0)
{
tr_->getMixDiffCoeffsMass(Dc.data());
}
else
{
label nSp = nCanteraSp_;
tr_->getBinaryDiffCoeffs(nSp, Dij.data());
forAll(Dc, i)
{
Dc[i] = Dij[nSp*idxFracOne+i] / p[0];
}
}
forAll(thermo.composition().species(), k)
{