1volScalarField::Internal
Sp
10 dimensionedScalar(dgdt.dimensions(), Zero)
13volScalarField::Internal
Su
22 dimensionedScalar(dgdt.dimensions(), Zero)
27 if (dgdt[celli] > 0.0 &&
alpha1[celli] > 0.0)
29 Sp[celli] -= dgdt[celli]*
alpha1[celli];
30 Su[celli] += dgdt[celli]*
alpha1[celli];
32 else if (dgdt[celli] < 0.0 &&
alpha1[celli] < 1.0)
34 Sp[celli] += dgdt[celli]*(1.0 -
alpha1[celli]);
38volScalarField::Internal
divU
tmp< surfaceScalarField > phiCN(alphaPhic)
const volScalarField & alpha1
#define forAll(list, i)
Loop across all elements in list.