FreeFOAM The Cross-Platform CFD Toolkit
alphaEqnSubCycle.H
Go to the documentation of this file.
2 (
3  IOobject
4  (
5  "rhoPhi",
6  runTime.timeName(),
7  mesh
8  ),
9  mesh,
10  dimensionedScalar("0", dimensionSet(1, 0, -1, 0, 0), 0)
11 );
12 
13 {
14  label nAlphaCorr
15  (
16  readLabel(piso.lookup("nAlphaCorr"))
17  );
18 
19  label nAlphaSubCycles
20  (
21  readLabel(piso.lookup("nAlphaSubCycles"))
22  );
23 
24  surfaceScalarField phic = mag(phi/mesh.magSf());
25  phic = min(interface.cAlpha()*phic, max(phic));
26 
27  volScalarField divU = fvc::div(phi);
28 
29  dimensionedScalar totalDeltaT = runTime.deltaT();
30 
31  if (nAlphaSubCycles > 1)
32  {
33  for
34  (
35  subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
36  !(++alphaSubCycle).end();
37  )
38  {
39  #include "alphaEqn.H"
40  }
41  }
42  else
43  {
44  #include "alphaEqn.H"
45  }
46 
47  if (nOuterCorr == 1)
48  {
49  interface.correct();
50  }
51 
52  rho == alpha1*rho1 + (scalar(1) - alpha1)*rho2;
53 }
54 
55 // ************************ vim: set sw=4 sts=4 et: ************************ //