alphaEqnSubCycle.H
Go to the documentation of this file.
2 {
3  dimensionedScalar totalDeltaT = runTime.deltaT();
4  surfaceScalarField rhoPhiSum
5  (
6  IOobject
7  (
8  "rhoPhiSum",
9  runTime.timeName(),
10  mesh
11  ),
12  mesh,
13  dimensionedScalar(rhoPhi.dimensions(), 0)
14  );
15 
16  tmp<volScalarField> trSubDeltaT;
17 
18  if (LTS)
19  {
20  trSubDeltaT =
21  fv::localEulerDdt::localRSubDeltaT(mesh, nAlphaSubCycles);
22  }
23 
24  for
25  (
26  subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
27  !(++alphaSubCycle).end();
28  )
29  {
30  #include "alphaEqn.H"
31  rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
32  }
33 
34  rhoPhi = rhoPhiSum;
35 }
36 else
37 {
38  #include "alphaEqn.H"
39 }
40 
41 rho == mixture.rho();
volScalarField & alpha1(mixture.alpha1())
fvMesh & mesh
bool LTS
Definition: createRDeltaT.H:1
const label nAlphaSubCycles(alphaControls.lookup< label >("nAlphaSubCycles"))
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
rhoPhi
Definition: rhoEqn.H:10
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField