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

plasmaReactingFoam write electron velocities and fvm::div instead of mvConvection->fvmDiv

Browse source
This commit is contained in:
ignis 2016-11-14 03:25:07 +09:00
parent 5e2f512f84
commit 3d72ced0d4
3 changed files with 33 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

2
applications/solvers/combustion/plasmaReactingFoam/EEqn.H
View file

@ -3,7 +3,7 @@
fvScalarMatrix EEqn
(
fvm::ddt(rho, he) + mvConvection->fvmDiv(phi, he)
fvm::ddt(rho, he) + fvm::div(phi, he, "div(phi,Yi_h)")
+ fvc::ddt(rho, K) + fvc::div(phi, K)
+ (
he.name() == "e"

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

@ -132,9 +132,9 @@ tmp<fv::convectionScheme<scalar> > mvConvection
if (Y[i].name() == electronSpecie)
{
Udrift = - linearInterpolate
((mue/ng)*E + ((De/ng/Te)*fvc::grad(Te)));
ve = (Udrift & mesh.Sf()) + q;
Udrift = -(mue/ng)*E;
Uthermal = -((De/ng/Te)*fvc::grad(Te));
ve = (linearInterpolate(Udrift+Uthermal) & mesh.Sf()) + q;
// Wall electron flux correction
forAll (wallPatcheIDs, pidx)
@ -216,10 +216,10 @@ tmp<fv::convectionScheme<scalar> > mvConvection
fvm::ddt(rho, Yi)
+
( nCharge != 0
? mvConvection->fvmDiv(phi_drift, Yi)
? fvm::div(phi_drift, Yi, "div(phi,Yi_h)")
: ( neutrals.contains(Y[i].name())
? mvConvection->fvmDiv(phi_neutral, Yi)
: mvConvection->fvmDiv(phi, Yi)
? fvm::div(phi_neutral, Yi, "div(phi,Yi_h)")
: fvm::div(phi, Yi, "div(phi,Yi_h)")
)
)
// - fvm::laplacian(turbulence->muEff(), Yi)

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

@ -330,7 +330,31 @@ bolos.presolve();
*/
Info<< "Calculating face flux field ve\n" << endl;
surfaceVectorField Udrift ("Udrift", - linearInterpolate(mue*E/ng));
volVectorField Udrift
(
IOobject
(
"Udrift",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
- (mue*E/ng)
);
volVectorField Uthermal
(
IOobject
(
"Uthermal",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
- ((De/ng/Te)*fvc::grad(Te))
);
surfaceScalarField ve
(
@ -342,7 +366,7 @@ surfaceScalarField ve
IOobject::NO_READ,
IOobject::NO_WRITE
),
Udrift & mesh.Sf()
linearInterpolate(Udrift+Uthermal) & mesh.Sf()
);