From ebe0b6f0cb77be4e1caa859f7c88c23dff592186 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Wed, 23 Jun 2021 10:14:32 +0100 Subject: [PATCH] surfaceFilmModels::kinematicSingleLayer: Added maximum Courant number to info e.g. in the tutorials/combustion/buoyantReactingFoam/Lagrangian/cylinder log: . . . Surface film: thermoSingleLayer added mass = 0.043 current mass = 0.04300000007 min/max(mag(U)) = 4.103129368e-18, 0.3516130906 max Courant number = 0.1622700892 min/max(delta) = 0, 0.0009353648867 coverage = 0.01952279011 ejected mass = 0 - patch: region0_to_wallFilmRegion_wallFilmFaces: 0 transferred mass = 0 - patch: region0_to_wallFilmRegion_wallFilmFaces: 0 min/mean/max(T) = 299.9999757, 300, 300.0000115 mass phase change = 0 vapourisation rate = 0 . . . --- .../kinematicSingleLayer/kinematicSingleLayer.C | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C b/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C index 3d540f4274..15f55c0b7d 100644 --- a/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C +++ b/src/regionModels/surfaceFilmModels/kinematicSingleLayer/kinematicSingleLayer.C @@ -973,8 +973,6 @@ scalar kinematicSingleLayer::CourantNumber() const const scalar CoNum = 0.5*gMax(sumPhi/regionMesh().V().field())*time_.deltaTValue(); - Info<< "Film max Courant number: " << CoNum << endl; - return CoNum; } @@ -1018,6 +1016,7 @@ void kinematicSingleLayer::info() << gSum((delta_*rho()*magSf())()) << nl << indent << "min/max(mag(U)) = " << gMin(mag(Uinternal)) << ", " << gMax(mag(Uinternal)) << nl + << indent << "max Courant number = " << CourantNumber() << nl << indent << "min/max(delta) = " << gMin(deltaInternal) << ", " << gMax(deltaInternal) << nl << indent << "coverage = "