OpenFOAM
v2312
The open source CFD toolbox
setMultiRegionDeltaT.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) 2011 OpenFOAM Foundation
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
setMultiRegionDeltaT
28
29
Description
30
Reset the timestep to maintain a constant maximum Courant numbers.
31
Reduction of time-step is immediate, but increase is damped to avoid
32
unstable oscillations.
33
34
\*---------------------------------------------------------------------------*/
35
36
if
(
adjustTimeStep
)
37
{
38
if
(
CoNum
== -GREAT)
39
{
40
CoNum
= SMALL;
41
}
42
43
if
(
DiNum
== -GREAT)
44
{
45
DiNum
= SMALL;
46
}
47
48
49
const
scalar TFactorFluid =
maxCo
/(
CoNum
+ SMALL);
50
const
scalar TFactorSolid =
maxDi
/(
DiNum
+ SMALL);
51
const
scalar TFactorFilm =
maxCo
/(
surfaceFilm
.CourantNumber() + SMALL);
52
53
const
scalar dt0 =
runTime
.deltaTValue();
54
55
runTime
.setDeltaT
56
(
57
min
58
(
59
dt0*
min
(
min
(TFactorFluid,
min
(TFactorFilm, TFactorSolid)), 1.2),
60
maxDeltaT
61
)
62
);
63
}
64
65
66
// ************************************************************************* //
DiNum
scalar DiNum
Definition:
solidRegionDiffusionNo.H:1
runTime
engineTime & runTime
Definition:
createEngineTime.H:13
maxDeltaT
maxDeltaT
Definition:
readEngineTimeControls.H:3
maxDi
scalar maxDi
Definition:
createRegionControls.H:23
Foam::min
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Definition:
hashSets.C:26
CoNum
scalar CoNum
Definition:
compressibleCourantNo.H:25
surfaceFilm
regionModels::surfaceFilmModel & surfaceFilm
Definition:
createFieldRefs.H:3
maxCo
scalar maxCo
Definition:
createTimeControls.H:29
adjustTimeStep
bool adjustTimeStep
Definition:
createTimeControls.H:26
applications
solvers
combustion
fireFoam
setMultiRegionDeltaT.H
Generated by
1.8.14