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