reverseLinear.H
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-2016 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 Class
25  Foam::reverseLinear
26 
27 Description
28  Inversed weight central-differencing interpolation scheme class.
29 
30  Useful for inverse weighted and harmonic interpolations.
31 
32 SourceFiles
33  reverseLinear.C
34 
35 \*---------------------------------------------------------------------------*/
36 
37 #ifndef reverseLinear_H
38 #define reverseLinear_H
39 
41 #include "volFields.H"
42 
43 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
44 
45 namespace Foam
46 {
47 
48 /*---------------------------------------------------------------------------*\
49  Class reverseLinear Declaration
50 \*---------------------------------------------------------------------------*/
51 
52 template<class Type>
53 class reverseLinear
54 :
55  public surfaceInterpolationScheme<Type>
56 {
57  // Private Member Functions
58 
59  //- Disallow default bitwise assignment
60  void operator=(const reverseLinear&);
61 
62 
63 public:
64 
65  //- Runtime type information
66  TypeName("reverseLinear");
67 
68 
69  // Constructors
70 
71  //- Construct from mesh
72  reverseLinear(const fvMesh& mesh)
73  :
74  surfaceInterpolationScheme<Type>(mesh)
75  {}
76 
77  //- Construct from Istream
79  :
80  surfaceInterpolationScheme<Type>(mesh)
81  {}
82 
83  //- Construct from faceFlux and Istream
85  (
86  const fvMesh& mesh,
87  const surfaceScalarField&,
88  Istream&
89  )
90  :
91  surfaceInterpolationScheme<Type>(mesh)
92  {}
93 
94 
95  // Member Functions
96 
97  //- Return the interpolation weighting factors
99  (
101  ) const
102  {
103  const fvMesh& mesh = this->mesh();
104 
105  tmp<surfaceScalarField> tcdWeights
106  (
107  mesh.surfaceInterpolation::weights()
108  );
109  const surfaceScalarField& cdWeights = tcdWeights();
110 
111  tmp<surfaceScalarField> treverseLinearWeights
112  (
114  (
115  IOobject
116  (
117  "reverseLinearWeights",
118  mesh.time().timeName(),
119  mesh
120  ),
121  mesh,
122  dimless
123  )
124  );
125  surfaceScalarField& reverseLinearWeights =
126  treverseLinearWeights.ref();
127 
128  reverseLinearWeights.primitiveFieldRef() =
129  1.0 - cdWeights.primitiveField();
130 
131  surfaceScalarField::Boundary& rlwbf =
132  reverseLinearWeights.boundaryFieldRef();
133 
134 
135  forAll(mesh.boundary(), patchi)
136  {
137  if (rlwbf[patchi].coupled())
138  {
139  rlwbf[patchi] = 1.0 - cdWeights.boundaryField()[patchi];
140  }
141  else
142  {
143  rlwbf[patchi] = cdWeights.boundaryField()[patchi];
144  }
145  }
146 
147  return treverseLinearWeights;
148  }
149 };
150 
151 
152 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
153 
154 } // End namespace Foam
155 
156 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
157 
158 #endif
159 
160 // ************************************************************************* //
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:428
TypeName("reverseLinear")
Runtime type information.
Inversed weight central-differencing interpolation scheme class.
Definition: reverseLinear.H:52
An Istream is an abstract base class for all input systems (streams, files, token lists etc)...
Definition: Istream.H:57
tmp< surfaceScalarField > weights(const GeometricField< Type, fvPatchField, volMesh > &) const
Return the interpolation weighting factors.
Definition: reverseLinear.H:98
static word timeName(const scalar, const int precision=precision_)
Return time name of given scalar time.
Definition: Time.C:644
const fvMesh & mesh() const
Return mesh reference.
const Time & time() const
Return the top-level database.
Definition: fvMesh.H:243
reverseLinear(const fvMesh &mesh)
Construct from mesh.
Definition: reverseLinear.H:71
Internal & ref()
Return a reference to the dimensioned internal field.
label patchi
const dimensionSet dimless(0, 0, 0, 0, 0, 0, 0)
Definition: dimensionSets.H:47
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:78
A class for managing temporary objects.
Definition: PtrList.H:53
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Definition: IOobject.H:92
Abstract base class for surface interpolation schemes.
Namespace for OpenFOAM.
const fvBoundaryMesh & boundary() const
Return reference to boundary mesh.
Definition: fvMesh.C:545