34 readScalar(runTime.controlDict().lookup(
"maxAlphaCo"))
40 if (
mesh.nInternalFaces())
45 pos(alphaf - 0.01)*
pos(0.99 - alphaf)
46 *
mesh.surfaceInterpolation::deltaCoeffs()*
mag(
phi);
48 alphaCoNum =
max(SfUfbyDelta/
mesh.magSf())
49 .value()*runTime.deltaT().value();
51 meanAlphaCoNum = (
sum(SfUfbyDelta)/
sum(
mesh.magSf()))
52 .value()*runTime.deltaT().value();
55 Info<<
"Interface Courant Number mean: " << meanAlphaCoNum
56 <<
" max: " << alphaCoNum <<
endl;