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
.
.
.
This commit is contained in:
@ -973,8 +973,6 @@ scalar kinematicSingleLayer::CourantNumber() const
|
|||||||
const scalar CoNum =
|
const scalar CoNum =
|
||||||
0.5*gMax(sumPhi/regionMesh().V().field())*time_.deltaTValue();
|
0.5*gMax(sumPhi/regionMesh().V().field())*time_.deltaTValue();
|
||||||
|
|
||||||
Info<< "Film max Courant number: " << CoNum << endl;
|
|
||||||
|
|
||||||
return CoNum;
|
return CoNum;
|
||||||
}
|
}
|
||||||
|
|
||||||
@ -1018,6 +1016,7 @@ void kinematicSingleLayer::info()
|
|||||||
<< gSum((delta_*rho()*magSf())()) << nl
|
<< gSum((delta_*rho()*magSf())()) << nl
|
||||||
<< indent << "min/max(mag(U)) = " << gMin(mag(Uinternal)) << ", "
|
<< indent << "min/max(mag(U)) = " << gMin(mag(Uinternal)) << ", "
|
||||||
<< gMax(mag(Uinternal)) << nl
|
<< gMax(mag(Uinternal)) << nl
|
||||||
|
<< indent << "max Courant number = " << CourantNumber() << nl
|
||||||
<< indent << "min/max(delta) = " << gMin(deltaInternal) << ", "
|
<< indent << "min/max(delta) = " << gMin(deltaInternal) << ", "
|
||||||
<< gMax(deltaInternal) << nl
|
<< gMax(deltaInternal) << nl
|
||||||
<< indent << "coverage = "
|
<< indent << "coverage = "
|
||||||
|
|||||||
Reference in New Issue
Block a user