alphaEqnSubCycle.H
Go to the documentation of this file.
1 {
2  // Standard face-flux compression coefficient
3  surfaceScalarField phic(mixture.cAlpha()*mag(phi/mesh.magSf()));
4 
5  // Add the optional isotropic compression contribution
6  if (icAlpha > 0)
7  {
8  phic *= (1.0 - icAlpha);
9  phic += (mixture.cAlpha()*icAlpha)*fvc::interpolate(mag(U));
10  }
11 
12  if (nAlphaSubCycles > 1)
13  {
14  dimensionedScalar totalDeltaT = runTime.deltaT();
15  surfaceScalarField rhoPhiSum("rhoPhiSum", rhoPhi);
16 
17  for
18  (
19  subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
20  !(++alphaSubCycle).end();
21  )
22  {
23  #include "alphaEqn.H"
24  rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
25  }
26 
27  rhoPhi = rhoPhiSum;
28  }
29  else
30  {
31  #include "alphaEqn.H"
32  }
33 
34  rho == alpha1*rho1 + alpha2*rho2;
35 }
surfaceScalarField & phi
engineTime & runTime
scalar icAlpha(alphaControls.lookupOrDefault< scalar >("icAlpha", 0))
dynamicFvMesh & mesh
const volScalarField & alpha1
alpha2
Definition: alphaEqn.H:115
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
U
Definition: pEqn.H:72
surfaceScalarField phic(mixture.cAlpha() *mag(phi/mesh.magSf()))
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
rhoPhi
Definition: rhoEqn.H:10
dimensioned< scalar > mag(const dimensioned< Type > &)
const dimensionedScalar & rho2
Definition: createFields.H:44
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
phaseChangeTwoPhaseMixture & mixture
Definition: createFields.H:38
const dimensionedScalar & rho1
Definition: createFields.H:43
label nAlphaSubCycles(readLabel(alphaControls.lookup("nAlphaSubCycles")))