mirror of
https://develop.openfoam.com/Development/openfoam.git
synced 2025-12-28 03:37:59 +00:00
dbe5e08213
- added pdFinal option
- labelled muf - grad(muf) used, so can now manipulate in fvSolution etc
- added absolute/relative around the calculation of divU
46 lines
942 B
C
46 lines
942 B
C
{
|
|
label nAlphaCorr
|
|
(
|
|
readLabel(piso.lookup("nAlphaCorr"))
|
|
);
|
|
|
|
label nAlphaSubCycles
|
|
(
|
|
readLabel(piso.lookup("nAlphaSubCycles"))
|
|
);
|
|
|
|
surfaceScalarField phic = mag(phi/mesh.magSf());
|
|
phic = min(interface.cAlpha()*phic, max(phic));
|
|
|
|
fvc::makeAbsolute(phi, U);
|
|
volScalarField divU = fvc::div(phi);
|
|
fvc::makeRelative(phi, U);
|
|
|
|
if (nAlphaSubCycles > 1)
|
|
{
|
|
dimensionedScalar totalDeltaT = runTime.deltaT();
|
|
surfaceScalarField rhoPhiSum = 0.0*rhoPhi;
|
|
|
|
for
|
|
(
|
|
subCycle<volScalarField> alphaSubCycle(alpha1, nAlphaSubCycles);
|
|
!(++alphaSubCycle).end();
|
|
)
|
|
{
|
|
# include "alphaEqns.H"
|
|
rhoPhiSum += (runTime.deltaT()/totalDeltaT)*rhoPhi;
|
|
}
|
|
|
|
rhoPhi = rhoPhiSum;
|
|
}
|
|
else
|
|
{
|
|
# include "alphaEqns.H"
|
|
}
|
|
|
|
if (oCorr == 0)
|
|
{
|
|
interface.correct();
|
|
}
|
|
}
|