FreeFOAM The Cross-Platform CFD Toolkit
alphaEqnSubCycle.H
Go to the documentation of this file.
1 label nAlphaCorr
2 (
3  readLabel(piso.lookup("nAlphaCorr"))
4 );
5 
6 label nAlphaSubCycles
7 (
8  readLabel(piso.lookup("nAlphaSubCycles"))
9 );
10 
12 {
13  dimensionedScalar totalDeltaT = runTime.deltaT();
14  surfaceScalarField rhoPhiSum = 0.0*rhoPhi;
15 
16  for
17  (
18  subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
19  !(++alphaSubCycle).end();
20  )
21  {
22 # include "alphaEqn.H"
23  rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
24  }
25 
26  rhoPhi = rhoPhiSum;
27 }
28 else
29 {
30 # include "alphaEqn.H"
31 }
32 
33 interface.correct();
34 
35 rho == alpha1*rho1 + (scalar(1) - alpha1)*rho2;