nutkWallFunctionFvPatchScalarField.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-2020 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 
27 #include "momentumTransportModel.H"
28 #include "fvPatchFieldMapper.H"
29 #include "volFields.H"
30 #include "wallFvPatch.H"
32 
33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
34 
35 namespace Foam
36 {
37 
38 // * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
39 
41 {
42  const label patchi = patch().index();
43 
44  const momentumTransportModel& turbModel =
45  db().lookupObject<momentumTransportModel>
46  (
48  (
49  momentumTransportModel::typeName,
50  internalField().group()
51  )
52  );
53 
54  const scalarField& y = turbModel.y()[patchi];
55  const tmp<volScalarField> tk = turbModel.k();
56  const volScalarField& k = tk();
57  const tmp<scalarField> tnuw = turbModel.nu(patchi);
58  const scalarField& nuw = tnuw();
59 
60  const scalar Cmu25 = pow025(Cmu_);
61 
62  tmp<scalarField> tnutw(new scalarField(patch().size(), 0.0));
63  scalarField& nutw = tnutw.ref();
64 
65  forAll(nutw, facei)
66  {
67  label celli = patch().faceCells()[facei];
68 
69  scalar yPlus = Cmu25*y[facei]*sqrt(k[celli])/nuw[facei];
70 
71  if (yPlus > yPlusLam_)
72  {
73  nutw[facei] = nuw[facei]*(yPlus*kappa_/log(E_*yPlus) - 1.0);
74  }
75  }
76 
77  return tnutw;
78 }
79 
80 
81 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
82 
84 (
85  const fvPatch& p,
87 )
88 :
90 {}
91 
92 
94 (
96  const fvPatch& p,
98  const fvPatchFieldMapper& mapper
99 )
100 :
101  nutWallFunctionFvPatchScalarField(ptf, p, iF, mapper)
102 {}
103 
104 
106 (
107  const fvPatch& p,
109  const dictionary& dict
110 )
111 :
113 {}
114 
115 
117 (
119 )
120 :
122 {}
123 
124 
126 (
129 )
130 :
132 {}
133 
134 
135 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
136 
138 {
139  const label patchi = patch().index();
140 
141  const momentumTransportModel& turbModel =
142  db().lookupObject<momentumTransportModel>
143  (
145  (
146  momentumTransportModel::typeName,
147  internalField().group()
148  )
149  );
150 
151  const scalarField& y = turbModel.y()[patchi];
152 
153  const tmp<volScalarField> tk = turbModel.k();
154  const volScalarField& k = tk();
155  tmp<scalarField> kwc = k.boundaryField()[patchi].patchInternalField();
156  const tmp<scalarField> tnuw = turbModel.nu(patchi);
157  const scalarField& nuw = tnuw();
158 
159  return pow025(Cmu_)*y*sqrt(kwc)/nuw;
160 }
161 
162 
163 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
164 
166 (
169 );
170 
171 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
172 
173 } // End namespace Foam
174 
175 // ************************************************************************* //
virtual tmp< volScalarField > nu() const =0
Return the laminar viscosity.
const char *const group
Group name for atomic constants.
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
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:158
const Boundary & boundaryField() const
Return const-reference to the boundary field.
T & ref() const
Return non-const reference or generate a fatal error.
Definition: tmpI.H:174
dimensionedScalar sqrt(const dimensionedScalar &ds)
A finiteVolume patch using a polyPatch and a fvBoundaryMesh.
Definition: fvPatch.H:61
dimensionedScalar pow025(const dimensionedScalar &ds)
Abstract base class with a fat-interface to all derived classes covering all possible ways in which t...
Definition: fvPatchField.H:66
const nearWallDist & y() const
Return the near wall distances.
label k
Boltzmann constant.
Macros for easy insertion into run-time selection tables.
static const nutWallFunctionFvPatchScalarField & nutw(const momentumTransportModel &turbModel, const label patchi)
Return the nut patchField for the given wall patch.
scalar y
Turbulent viscosity wall-function boundary condition for high Reynolds number flows based on near-wal...
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
static word groupName(Name name, const word &group)
Foam::fvPatchFieldMapper.
virtual tmp< scalarField > yPlus() const
Calculate and return the yPlus at the boundary.
virtual tmp< volScalarField > k() const =0
Return the turbulence kinetic energy.
Abstract base class for turbulence models (RAS, LES and laminar).
label patchi
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
nutWallFunctionFvPatchScalarField(const fvPatch &, const DimensionedField< scalar, volMesh > &)
Construct from patch and internal field.
A class for managing temporary objects.
Definition: PtrList.H:53
scalar yPlusLam_
Y+ at the edge of the laminar sublayer.
makePatchTypeField(fvPatchScalarField, atmBoundaryLayerInletEpsilonFvPatchScalarField)
virtual tmp< scalarField > nut() const
Calculate the turbulence viscosity.
nutkWallFunctionFvPatchScalarField(const fvPatch &, const DimensionedField< scalar, volMesh > &)
Construct from patch and internal field.
Namespace for OpenFOAM.