sigmaRadial.C
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-2019 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 \*---------------------------------------------------------------------------*/
25 
26 #include "sigmaRadial.H"
28 
29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
30 
31 namespace Foam
32 {
33 namespace extrudeModels
34 {
35 
36 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
37 
39 
41 
42 
43 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
44 
46 :
47  extrudeModel(typeName, dict),
48  RTbyg_(coeffDict_.lookup<scalar>("RTbyg")),
49  pRef_(coeffDict_.lookup<scalar>("pRef")),
50  pStrat_(coeffDict_.lookup<scalar>("pStrat"))
51 {
52  if (mag(expansionRatio() - 1.0) > small)
53  {
55  << "Ignoring expansionRatio setting." << endl;
56  }
57 }
58 
59 
60 // * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
61 
63 {}
64 
65 
66 // * * * * * * * * * * * * * * * * Operators * * * * * * * * * * * * * * * * //
67 
68 point sigmaRadial::operator()
69 (
70  const point& surfacePoint,
71  const vector& surfaceNormal,
72  const label layer
73 ) const
74 {
75  // radius of the surface
76  scalar rs = mag(surfacePoint);
77  vector rsHat = surfacePoint/rs;
78 
79  scalar p = pRef_ - layer*(pRef_ - pStrat_)/nLayers_;
80  scalar r = rs - RTbyg_*log(p/pRef_);
81 
82  return r*rsHat;
83 }
84 
85 
86 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
87 
88 } // End namespace extrudeModels
89 } // End namespace Foam
90 
91 // ************************************************************************* //
Macros for easy insertion into run-time selection tables.
A list of keyword definitions, which are a keyword followed by any number of values (e....
Definition: dictionary.H:162
Top level extrusion model class.
Definition: extrudeModel.H:52
scalar expansionRatio() const
Definition: extrudeModel.C:65
virtual ~sigmaRadial()
Destructor.
Definition: sigmaRadial.C:62
sigmaRadial(const dictionary &dict)
Construct from dictionary.
Definition: sigmaRadial.C:45
#define WarningInFunction
Report a warning using Foam::Warning.
defineTypeNameAndDebug(cyclicSector, 0)
addToRunTimeSelectionTable(extrudeModel, cyclicSector, dictionary)
Namespace for OpenFOAM.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:257
dimensionedScalar log(const dimensionedScalar &ds)
dimensioned< scalar > mag(const dimensioned< Type > &)
dictionary dict
volScalarField & p