alphaSuSp.H
Go to the documentation of this file.
1volScalarField::Internal Sp
2(
3 IOobject
4 (
5 "Sp",
6 runTime.timeName(),
7 mesh
8 ),
9 mesh,
10 dimensionedScalar(dgdt.dimensions(), Zero)
11);
12
13volScalarField::Internal Su
14(
15 IOobject
16 (
17 "Su",
18 runTime.timeName(),
19 mesh
20 ),
21 mesh,
22 dimensionedScalar(dgdt.dimensions(), Zero)
23);
24
25forAll(dgdt, celli)
26{
27 if (dgdt[celli] > 0.0)
28 {
29 Sp[celli] -= dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
30 Su[celli] += dgdt[celli]/max(1.0 - alpha1[celli], 1e-4);
31 }
32 else if (dgdt[celli] < 0.0)
33 {
34 Sp[celli] += dgdt[celli]/max(alpha1[celli], 1e-4);
35 }
36}
37
38volScalarField::Internal divU
39(
40 mesh.moving()
41 ? fvc::div(phiCN() + mesh.phi())
42 : fvc::div(phiCN())
43);
tmp< surfaceScalarField > phiCN(alphaPhic)
Y[inertIndex] max(0.0)
const volScalarField & alpha1
dynamicFvMesh & mesh
engineTime & runTime
zeroField divU
Definition: alphaSuSp.H:3
zeroField Su
Definition: alphaSuSp.H:1
zeroField Sp
Definition: alphaSuSp.H:2
volScalarField & e
Definition: createFields.H:11
#define forAll(list, i)
Loop across all elements in list.
Definition: stdFoam.H:333