sigmaRadial.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 | Copyright (C) 2011-2015 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 
38 defineTypeNameAndDebug(sigmaRadial, 0);
39 
40 addToRunTimeSelectionTable(extrudeModel, sigmaRadial, dictionary);
41 
42 
43 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
44 
46 :
47  extrudeModel(typeName, dict),
48  RTbyg_(readScalar(coeffDict_.lookup("RTbyg"))),
49  pRef_(readScalar(coeffDict_.lookup("pRef"))),
50  pStrat_(readScalar(coeffDict_.lookup("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 // ************************************************************************* //
virtual ~sigmaRadial()
Destructor.
Definition: sigmaRadial.C:62
dictionary dict
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
dimensionedScalar log(const dimensionedScalar &ds)
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:137
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:253
sigmaRadial(const dictionary &dict)
Construct from dictionary.
Definition: sigmaRadial.C:45
defineTypeNameAndDebug(cyclicSector, 0)
Macros for easy insertion into run-time selection tables.
Top level extrusion model class.
Definition: extrudeModel.H:51
stressControl lookup("compactNormalStress") >> compactNormalStress
bool readScalar(const char *buf, doubleScalar &s)
Read whole of buf as a scalar. Return true if succesful.
Definition: doubleScalar.H:63
addToRunTimeSelectionTable(extrudeModel, cyclicSector, dictionary)
#define WarningInFunction
Report a warning using Foam::Warning.
scalar expansionRatio() const
Definition: extrudeModel.C:65
dimensioned< scalar > mag(const dimensioned< Type > &)
volScalarField & p
Namespace for OpenFOAM.