OpenFOAM: API Guide
v2112
The open source CFD toolbox
surfaceCourantNo.H
Go to the documentation of this file.
1
/*---------------------------------------------------------------------------*\
2
========= |
3
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4
\\ / O peration |
5
\\ / A nd | www.openfoam.com
6
\\/ M anipulation |
7
-------------------------------------------------------------------------------
8
Copyright (C) 2016-2017 Wikki Ltd
9
-------------------------------------------------------------------------------
10
License
11
This file is part of OpenFOAM.
12
13
OpenFOAM is free software: you can redistribute it and/or modify it
14
under the terms of the GNU General Public License as published by
15
the Free Software Foundation, either version 3 of the License, or
16
(at your option) any later version.
17
18
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
19
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
20
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21
for more details.
22
23
You should have received a copy of the GNU General Public License
24
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
25
26
Global
27
surfaceCourantNo
28
29
Author
30
Hrvoje Jasak, Wikki Ltd.
31
32
Description
33
Calculates and outputs the mean and maximum Courant Numbers for the
34
Finite Area method.
35
36
\*---------------------------------------------------------------------------*/
37
38
scalar
CoNum
= 0.0;
39
scalar
meanCoNum
= 0.0;
40
scalar
velMag
= 0.0;
41
42
if
(
aMesh
.nInternalEdges())
43
{
44
edgeScalarField SfUfbyDelta
45
(
46
aMesh
.edgeInterpolation::deltaCoeffs()*mag(
phis
)
47
);
48
49
CoNum
=
max
(SfUfbyDelta/
aMesh
.magLe())
50
.value()*
runTime
.deltaT().value();
51
52
meanCoNum
= (sum(SfUfbyDelta)/sum(
aMesh
.magLe()))
53
.value()*
runTime
.deltaT().value();
54
55
velMag
=
max
(mag(
phis
)/
aMesh
.magLe()).value();
56
}
57
58
Info
<<
"Courant Number mean: "
<<
meanCoNum
59
<<
" max: "
<<
CoNum
60
<<
" velocity magnitude: "
<<
velMag
<<
endl
;
61
62
63
// ************************************************************************* //
max
Y[inertIndex] max(0.0)
runTime
engineTime & runTime
Definition:
createEngineTime.H:13
Foam::Info
messageStream Info
Information stream (stdout output on master, null elsewhere)
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition:
Ostream.H:372
phis
edgeScalarField phis(IOobject("phis", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), linearEdgeInterpolate(Us) &aMesh.Le())
aMesh
faMesh aMesh(mesh)
velMag
scalar velMag
Definition:
surfaceCourantNo.H:37
meanCoNum
scalar meanCoNum
Definition:
surfaceCourantNo.H:36
CoNum
scalar CoNum
Definition:
surfaceCourantNo.H:35
applications
solvers
finiteArea
liquidFilmFoam
surfaceCourantNo.H
Generated by
1.9.5
OPENFOAM® is a registered
trademark
of OpenCFD Ltd.