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

removed mole fraction array

Browse source
This commit is contained in:
ignis 2017-05-16 14:23:06 +09:00
parent 6e8f90c595
commit bab311826e
1 changed files with 6 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

13
neEqn.H
View file

@ -6,22 +6,21 @@
Te = T;
Te.correctBoundaryConditions();
scalarField X (csSpecies.size(), 0.0);
scalarField &W (csSpeciesW);
forAll (Te, cid)
{
scalar eEps = (3./2.) * kB.value() * Te[cid]
/ eCharge.value();
scalar eEps = (3./2.) * kB.value() * Te[cid] / eCharge.value();
scalar sumX = 0.0;
sigma_m[cid] = 0;
forAll (csSpecies, isp)
{
X[isp] = Y[csSpeciesI[isp]][cid]/W[isp];
sigma_m[cid] +=
X[isp] * csList[isp].value(eEps);
scalar Xi = Y[csSpeciesI[isp]][cid]/W[isp];
sumX += Xi;
sigma_m[cid] += Xi * csList[isp].value(eEps);
}
sigma_m[cid] /= sum(X);
sigma_m[cid] /= sumX;
}
Info << min(sigma_m) << " / " << max(sigma_m) << endl;