All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
temperatureAndPressure.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 Global
25  temperatureAndPressure.H
26 
27 Description
28  Accumulates values for temperature and pressure measurement, and
29  calculates and outputs the average values at output times.
30  Requires temperatureAndPressureVariables.H to be declared before the
31  timeloop.
32 
33 \*---------------------------------------------------------------------------*/
34 
36 
38 
40 
42 
44 
46 
48 
50 
51 if (runTime.writeTime())
52 {
53  if (accumulatedNMols)
54  {
55  Info<< "calculating averages" << endl;
56 
58  (
60  *
61  (
63  -
65  )
66  );
67 
69  (
70  (
71  (accumulatedNMols/nAveragingSteps)
73  + accumulatedTotalrDotfSum/(6.0*nAveragingSteps)
74  )
75  /
77  );
78 
79  Info<< "----------------------------------------" << nl
80  << "Averaged properties" << nl
81  << "Average |velocity| = "
83  << "Average temperature = " << averageTemperature << nl
84  << "Average pressure = " << averagePressure << nl
85  << "----------------------------------------" << endl;
86  }
87  else
88  {
89  Info<< "Not averaging temperature and pressure: "
90  << "no molecules in system" << endl;
91  }
92 
94 
96 
98 
100 
101  accumulatedTotalPE = 0.0;
102 
104 
105  accumulatedNMols = 0;
106 
107  accumulatedDOFs = 0;
108 }
109 
110 
111 // ************************************************************************* //
scalar singleStepTotalrDotf
label singleStepNMols
scalar averagePressure
accumulatedTotalLinearKE
engineTime & runTime
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
accumulatedTotalPE
label k
Boltzmann constant.
scalar singleStepTotalMass
static const zero Zero
Definition: zero.H:97
accumulatedTotalLinearMomentum
dimensioned< scalar > magSqr(const dimensioned< Type > &)
static const char nl
Definition: Ostream.H:260
scalar averageTemperature
accumulatedTotalMass
accumulatedTotalAngularKE
label singleStepDOFs
messageStream Info
accumulatedTotalrDotfSum
dimensioned< scalar > mag(const dimensioned< Type > &)
scalar singleStepTotalLinearKE
vector singleStepTotalLinearMomentum(Zero)
scalar singleStepTotalPE
scalar singleStepTotalAngularKE