alphaEqnSubCycle.H
Go to the documentation of this file.
1 {
2  // Standard face-flux compression coefficient
3  surfaceScalarField phic(interface.cAlpha()*mag(phi/mesh.magSf()));
4 
5  // Add the optional isotropic compression contribution
6  if (icAlpha > 0)
7  {
8  phic *= (1.0 - icAlpha);
9  phic += (interface.cAlpha()*icAlpha)*fvc::interpolate(mag(U));
10  }
11 
13 
14  if (nAlphaSubCycles > 1)
15  {
16  dimensionedScalar totalDeltaT = runTime.deltaT();
17  surfaceScalarField rhoPhiSum("rhoPhiSum", rhoPhi);
18 
19  for
20  (
21  subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
22  !(++alphaSubCycle).end();
23  )
24  {
25  #include "alphaEqn.H"
26  rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
27  }
28 
29  rhoPhi = rhoPhiSum;
30  }
31  else
32  {
33  #include "alphaEqn.H"
34  }
35 
36  rho == alpha1*rho1 + alpha2*rho2;
37 }
dimensioned< scalar > mag(const dimensioned< Type > &)
tmp< surfaceScalarField > interpolate(const RhoType &rho)
dynamicFvMesh & mesh
rhoPhi
Definition: rhoEqn.H:10
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
scalar icAlpha(alphaControls.lookupOrDefault< scalar >("icAlpha", 0))
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
rho1
Definition: pEqn.H:124
rho2
Definition: pEqn.H:125
volScalarField & alpha1
Definition: createFields.H:15
surfaceScalarField & phi
alpha2
Definition: alphaEqn.H:112
label nAlphaSubCycles(readLabel(alphaControls.lookup("nAlphaSubCycles")))
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
interfaceProperties interface(alpha1, U, mixture())
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
volScalarField divU(fvc::div(phi))
U
Definition: pEqn.H:82