alphaEqnsSubCycle.H
Go to the documentation of this file.
1 {
2  #include "alphaControls.H"
3 
5  phic = min(interface.cAlpha()*phic, max(phic));
6 
8 
9  if (nAlphaSubCycles > 1)
10  {
11  dimensionedScalar totalDeltaT = runTime.deltaT();
12  surfaceScalarField rhoPhiSum
13  (
14  IOobject
15  (
16  "rhoPhiSum",
17  runTime.timeName(),
18  mesh
19  ),
20  mesh,
21  dimensionedScalar("0", rhoPhi.dimensions(), 0)
22  );
23 
24  for
25  (
26  subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
27  !(++alphaSubCycle).end();
28  )
29  {
30  #include "alphaEqns.H"
31  rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
32  }
33 
34  rhoPhi = rhoPhiSum;
35  }
36  else
37  {
38  #include "alphaEqns.H"
39  }
40 }
label nAlphaSubCycles(readLabel(alphaControls.lookup("nAlphaSubCycles")))
volScalarField divU(fvc::div(fvc::absolute(phi, U)))
dimensioned< scalar > mag(const dimensioned< Type > &)
dynamicFvMesh & mesh
rhoPhi
Definition: rhoEqn.H:10
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
Definition: fvcMeshPhi.C:187
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
volScalarField & alpha1
Definition: createFields.H:15
surfaceScalarField & phi
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
interfaceProperties interface(alpha1, U, mixture())
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
U
Definition: pEqn.H:82