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
12 volScalarField divU
13 (
14 mesh.moving()
15 ? fvc::div(phi + mesh.phi())
16 : fvc::div(phi)
17 );
18
20 {
21 dimensionedScalar totalDeltaT = runTime.deltaT();
22 surfaceScalarField rhoPhiSum
23 (
24 IOobject
25 (
26 "rhoPhiSum",
27 runTime.timeName(),
28 mesh
29 ),
30 mesh,
31 dimensionedScalar(rhoPhi.dimensions(), Zero)
32 );
33
34 for
35 (
36 subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
37 !(++alphaSubCycle).end();
38 )
39 {
40 #include "alphaEqn.H"
41 rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
42 }
43
44 rhoPhi = rhoPhiSum;
45 }
46 else
47 {
48 #include "alphaEqn.H"
49 }
50
52}
surfaceScalarField phic(mixture.cAlpha() *mag(alphaPhic/mesh.magSf()))
rhoPhi
Definition: rhoEqn.H:10
surfaceScalarField & phi
const volScalarField & alpha1
volScalarField & rho2
const volScalarField & alpha2
volScalarField & rho1
U
Definition: pEqn.H:72
dynamicFvMesh & mesh
engineTime & runTime
zeroField divU
Definition: alphaSuSp.H:3
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
interfaceProperties interface(alpha1, U, thermo->transportPropertiesDict())
label nAlphaSubCycles(alphaControls.get< label >("nAlphaSubCycles"))
scalar icAlpha(alphaControls.getOrDefault< scalar >("icAlpha", 0))