zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

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

Browse Source 
Now consistent with the way the Courant number is calculated for other solvers
This commit is contained in:
Henry Weller
2015-06-17 20:20:49 +01:00
parent bc493180d9
commit 164448f8b1
1 changed files with 3 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
applications/solvers/compressible/rhoCentralFoam/centralCourantNo.H
View File

@ -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 = "