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

Compare commits

merge into: research:master
Branches Tags
research:master
research:mixture_cs_field_op
...
pull from: research:mixture_cs_field_op
Branches Tags
research:master
research:mixture_cs_field_op
Sign in to create a new pull request.

1 commit
master ... mixture_cs

Author SHA1 Message Date
ignis
d7cdb18820 substitute cell loop to field operations 2017-05-16 14:10:45 +09:00
2 changed files with 17 additions and 14 deletions
Show stats Expand all files Collapse all files

4
createFields.H
View file

@ -585,4 +585,6 @@ forAll (csSpecies, isp)
Info<< csSpeciesW << endl; Info<< csSpeciesW << endl;
scalarField sigma_m (Te); scalarField sigma_m (Te.internalField().size(), 0.0);
scalarField Xsum (Te.internalField().size(), 0.0);
scalarField meanEps (Te.internalField().size(), 0.0);

27
neEqn.H
View file

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