alphaEqnSubCycle.H
Go to the documentation of this file.
1 #include "alphaControls.H"
2 
4 {
5  dimensionedScalar totalDeltaT = runTime.deltaT();
6  surfaceScalarField rhoPhiSum
7  (
8  IOobject
9  (
10  "rhoPhiSum",
11  runTime.timeName(),
12  mesh
13  ),
14  mesh,
15  dimensionedScalar("0", rhoPhi.dimensions(), 0)
16  );
17 
18  for
19  (
20  subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
21  !(++alphaSubCycle).end();
22  )
23  {
24  #include "alphaEqn.H"
25  rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
26  }
27 
28  rhoPhi = rhoPhiSum;
29 }
30 else
31 {
32  #include "alphaEqn.H"
33 }
34 
engineTime & runTime
dynamicFvMesh & mesh
const volScalarField & alpha1
alpha2
Definition: alphaEqn.H:115
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
rhoPhi
Definition: rhoEqn.H:10
const dimensionedScalar & rho2
Definition: createFields.H:40
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
const dimensionedScalar & rho1
Definition: createFields.H:39
label nAlphaSubCycles(readLabel(alphaControls.lookup("nAlphaSubCycles")))