OpenFOAM: API Guide
v2112
The open source CFD toolbox
sonicDyMFoam.C
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-2017 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
Application
27
sonicDyMFoam
28
29
Group
30
grpCompressibleSolvers grpMovingMeshSolvers
31
32
Description
33
Transient solver for trans-sonic/supersonic, turbulent flow of a
34
compressible gas, with optional mesh motion and mesh topology changes.
35
36
\*---------------------------------------------------------------------------*/
37
38
#include "
fvCFD.H
"
39
#include "
dynamicFvMesh.H
"
40
#include "
psiThermo.H
"
41
#include "
turbulentFluidThermoModel.H
"
42
#include "
pimpleControl.H
"
43
#include "
CorrectPhi.H
"
44
#include "
fvOptions.H
"
45
46
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
47
48
int
main(
int
argc,
char
*argv[])
49
{
50
argList::addNote
51
(
52
"Transient solver for trans-sonic/supersonic, turbulent flow"
53
" of a compressible gas.\n"
54
"With optional mesh motion and mesh topology changes."
55
);
56
57
#include "
postProcess.H
"
58
59
#include "
setRootCaseLists.H
"
60
#include "
createTime.H
"
61
#include "
createDynamicFvMesh.H
"
62
#include "createDyMControls.H"
63
#include "createFields.H"
64
#include "createFieldRefs.H"
65
#include "
createRhoUf.H
"
66
#include "compressibleCourantNo.H"
67
#include "setInitialDeltaT.H"
68
#include "initContinuityErrs.H"
69
70
turbulence
->validate();
71
72
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
73
74
Info
<<
"\nStarting time loop\n"
<<
endl
;
75
76
while
(
runTime
.run())
77
{
78
#include "
readDyMControls.H
"
79
80
{
81
// Store divrhoU from the previous mesh so that it can be mapped
82
// and used in correctPhi to ensure the corrected phi has the
83
// same divergence
84
volScalarField
divrhoU
85
(
86
"divrhoU"
,
87
fvc::div(fvc::absolute(
phi
,
rho
,
U
))
88
);
89
90
#include "compressibleCourantNo.H"
91
#include "setDeltaT.H"
92
93
++
runTime
;
94
95
Info
<<
"Time = "
<<
runTime
.timeName() <<
nl
<<
endl
;
96
97
// Store momentum to set rhoUf for introduced faces.
98
volVectorField
rhoU(
"rhoU"
,
rho
*
U
);
99
100
// Do any mesh changes
101
mesh
.update();
102
103
if
(
mesh
.changing())
104
{
105
MRF
.update();
106
107
if
(
correctPhi
)
108
{
109
// Calculate absolute flux from the mapped surface velocity
110
phi
=
mesh
.Sf() &
rhoUf
;
111
112
#include "correctPhi.H"
113
114
// Make the fluxes relative to the mesh-motion
115
fvc::makeRelative(
phi
,
rho
,
U
);
116
}
117
}
118
119
if
(
checkMeshCourantNo
)
120
{
121
#include "
meshCourantNo.H
"
122
}
123
}
124
125
#include "rhoEqn.H"
126
Info
<<
"rho min/max : "
<<
min
(
rho
).value() <<
" "
<<
max
(
rho
).value()
127
<<
endl
;
128
129
// --- Pressure-velocity PIMPLE corrector loop
130
while
(
pimple
.loop())
131
{
132
#include "UEqn.H"
133
#include "EEqn.H"
134
135
// --- Pressure corrector loop
136
while
(
pimple
.correct())
137
{
138
#include "pEqn.H"
139
}
140
141
if
(
pimple
.turbCorr())
142
{
143
turbulence
->correct();
144
}
145
}
146
147
runTime
.write();
148
149
runTime
.printExecutionTime(Info);
150
}
151
152
Info
<<
"End\n"
<<
endl
;
153
154
return
0;
155
}
156
157
158
// ************************************************************************* //
CorrectPhi.H
max
Y[inertIndex] max(0.0)
phi
surfaceScalarField & phi
Definition:
setRegionFluidFields.H:8
MRF
IOMRFZoneList & MRF
Definition:
setRegionFluidFields.H:22
pimple
pimpleControl & pimple
Definition:
setRegionFluidFields.H:56
U
U
Definition:
pEqn.H:72
createDynamicFvMesh.H
mesh
dynamicFvMesh & mesh
Definition:
createDynamicFvMesh.H:6
runTime
engineTime & runTime
Definition:
createEngineTime.H:13
rhoUf
autoPtr< surfaceVectorField > rhoUf
Definition:
createRhoUfIfPresent.H:33
createRhoUf.H
Creates and initialises the velocity velocity field rhoUf.
createTime.H
dynamicFvMesh.H
turbulence
compressible::turbulenceModel & turbulence
Definition:
setRegionFluidFields.H:30
fvCFD.H
fvOptions.H
meshCourantNo.H
Calculates and outputs the mean and maximum Courant Numbers.
Foam::volVectorField
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition:
volFieldsFwd.H:83
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition:
volFieldsFwd.H:82
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
Foam::min
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Definition:
hashSets.C:33
Foam::nl
constexpr char nl
The newline '\n' character (0x0a)
Definition:
Ostream.H:53
pimpleControl.H
postProcess.H
Execute application functionObjects to post-process existing results.
psiThermo.H
readDyMControls.H
checkMeshCourantNo
checkMeshCourantNo
Definition:
readDyMControls.H:9
correctPhi
correctPhi
Definition:
readDyMControls.H:3
rho
rho
Definition:
readInitialConditions.H:88
setRootCaseLists.H
turbulentFluidThermoModel.H
applications
solvers
compressible
sonicFoam
sonicDyMFoam
sonicDyMFoam.C
Generated by
1.9.5
OPENFOAM® is a registered
trademark
of OpenCFD Ltd.