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

added constant mobility mode

Browse source
This commit is contained in:
ignis 2017-05-17 20:44:45 +09:00
parent bab311826e
commit de8600e4a7
2 changed files with 37 additions and 22 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
createFields.H
View file

@ -14,6 +14,11 @@ IOdictionary physicalProperties
)
);
Switch constETransport
(
physicalProperties.lookupOrDefault("constETransport", false)
);
scalar relaxDrift
(
physicalProperties.lookupOrDefault("relaxDrift", 1.0)
@ -586,3 +591,10 @@ Info<< csSpeciesW << endl;
scalarField sigma_m (Te);
if (constETransport)
{
mue.internalField() = mueNFac * mueN.value(0.0);
mue.correctBoundaryConditions();
}

3
neEqn.H
View file

@ -6,6 +6,8 @@
Te = T;
Te.correctBoundaryConditions();
if (!constETransport)
{
scalarField &W (csSpeciesW);
forAll (Te, cid)
@ -33,6 +35,7 @@
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);