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

Compare commits

base: research:master
Branches Tags
research:master
research:mixture_cs_field_op
..
compare: research:mixture_cs_field_op
Branches Tags
research:master
research:mixture_cs_field_op

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 27 additions and 38 deletions
Show stats Expand all files Collapse all files

16
createFields.H
View file

@ -14,11 +14,6 @@ IOdictionary physicalProperties
) )
); );
Switch constETransport
(
physicalProperties.lookupOrDefault("constETransport", false)
);
scalar relaxDrift scalar relaxDrift
( (
physicalProperties.lookupOrDefault("relaxDrift", 1.0) physicalProperties.lookupOrDefault("relaxDrift", 1.0)
@ -590,11 +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);
if (constETransport)
{
mue.internalField() = mueNFac * mueN.value(0.0);
mue.correctBoundaryConditions();
}

49
neEqn.H
View file

@ -6,37 +6,36 @@
Te = T; Te = T;
Te.correctBoundaryConditions(); Te.correctBoundaryConditions();
if (!constETransport) scalarField &W (csSpeciesW);
Xsum = 0.0;
sigma_m = 0.0;
meanEps = (3./2.) * kB.value() / eCharge.value() * Te.internalField();
forAll (csSpecies, isp)
{ {
scalarField &W (csSpeciesW); label k = csSpeciesI[isp];
forAll (Te, cid) scalarField X (Y[k].internalField() / W[isp]);
{
scalar eEps = (3./2.) * kB.value() * Te[cid] / eCharge.value();
scalar sumX = 0.0;
sigma_m[cid] = 0; Xsum += X;
forAll (csSpecies, isp)
{
scalar Xi = Y[csSpeciesI[isp]][cid]/W[isp];
sumX += Xi;
sigma_m[cid] += Xi * csList[isp].value(eEps);
}
sigma_m[cid] /= sumX;
}
Info << min(sigma_m) << " / " << max(sigma_m) << endl; sigma_m += X * csList[isp].value(meanEps);
/*
mue.internalField() = eCharge.value()
/ (sqrt(eMass.value() * kB.value() * Te) * sigma_m);
mue.correctBoundaryConditions();
*/
mue.internalField() = mueNFac * eCharge.value() / 3.0
/ (sqrt(pi * eMass.value() * kB.value() * Te / 8.0) * sigma_m);
mue.correctBoundaryConditions();
} }
sigma_m /= Xsum;
Info << min(sigma_m) << " / " << max(sigma_m) << endl;
/*
mue.internalField() = eCharge.value()
/ (sqrt(eMass.value() * kB.value() * Te) * sigma_m);
mue.correctBoundaryConditions();
*/
mue.internalField() = mueNFac * eCharge.value() / 3.0
/ (sqrt(pi * eMass.value() * kB.value() * Te / 8.0) * sigma_m);
mue.correctBoundaryConditions();
De = mue * Te * (kB / eCharge); De = mue * Te * (kB / eCharge);
De.correctBoundaryConditions(); De.correctBoundaryConditions();