Commit d28a72f1 authored by Henry Weller's avatar Henry Weller
Browse files

rhoCentralFoam: Updated Courant number calculation to be cell-based rather than face-based

Now consistent with the way the Courant number is calculated for other solvers
parent 0d3bfc0a
...@@ -31,16 +31,12 @@ Description ...@@ -31,16 +31,12 @@ Description
if (mesh.nInternalFaces()) if (mesh.nInternalFaces())
{ {
surfaceScalarField amaxSfbyDelta scalarField sumAmaxSf(fvc::surfaceSum(amaxSf)().internalField());
(
mesh.surfaceInterpolation::deltaCoeffs()*amaxSf
);
CoNum = max(amaxSfbyDelta/mesh.magSf()).value()*runTime.deltaTValue(); CoNum = 0.5*gMax(sumAmaxSf/mesh.V().field())*runTime.deltaTValue();
meanCoNum = meanCoNum =
(sum(amaxSfbyDelta)/sum(mesh.magSf())).value() 0.5*(gSum(sumAmaxSf)/gSum(mesh.V().field()))*runTime.deltaTValue();
*runTime.deltaTValue();
} }
Info<< "Mean and max Courant Numbers = " Info<< "Mean and max Courant Numbers = "
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment