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

electron and ion drift relaxation

Browse source
This commit is contained in:
ignis 2017-04-04 10:14:18 +09:00
parent 63419630cc
commit 005aa102df
3 changed files with 14 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
applications/solvers/combustion/plasmaReactingFoam/YEqn.H
View file

@ -119,6 +119,11 @@ tmp<fv::convectionScheme<scalar> > mvConvection
{
phis[i] = phi;
phis[i] += fvc::interpolate((rho*Di/T*(eCharge*z/kB))*E) & mesh.Sf();
if (relaxDrift < 1.0 && relaxDrift > 0.0)
{
phis[i] *= relaxDrift;
}
}
if (ions.contains(Y[i].name()))

5
applications/solvers/combustion/plasmaReactingFoam/createFields.H
View file

@ -12,6 +12,11 @@ IOdictionary physicalProperties
)
);
scalar relaxDrift
(
physicalProperties.lookupOrDefault("relaxDrift", 1.0)
);
dimensionedScalar epsilon0
(
physicalProperties.lookup("epsilon0")

4
applications/solvers/combustion/plasmaReactingFoam/neEqn.H
View file

@ -33,6 +33,10 @@
Udrift = -(mue/ng)*E;
if (relaxDrift < 1.0 && relaxDrift > 0.0)
{
Udrift *= relaxDrift;
}
Uthermal = -((De/ng/Te)*fvc::grad(Te));
ve = (linearInterpolate(Udrift+Uthermal) & mesh.Sf()) + q;