src/finiteVolume/interpolation/surfaceInterpolation/schemes/midPoint/midPoint.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
   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::midPoint
  26
  27 Description
  28     Mid-point interpolation (weighting factors = 0.5) scheme class.
  29
  30 SourceFiles
  31     midPoint.C
  32
  33 \*---------------------------------------------------------------------------*/
  34
  35 #ifndef midPoint_H
  36 #define midPoint_H
  37
  38 #include "surfaceInterpolationScheme.H"
  39 #include "volFields.H"
  40
  41 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  42
  43 namespace Foam
  44 {
  45
  46 /*---------------------------------------------------------------------------*\
  47                            Class midPoint Declaration
  48 \*---------------------------------------------------------------------------*/
  49
  50 template<class Type>
  51 class midPoint
  52 :
  53     public surfaceInterpolationScheme<Type>
  54 {
  55     // Private Member Functions
  56
  57         //- Disallow default bitwise assignment
  58         void operator=(const midPoint&);
  59
  60
  61 public:
  62
  63     //- Runtime type information
  64     TypeName("midPoint");
  65
  66
  67     // Constructors
  68
  69         //- Construct from mesh
  70         midPoint(const fvMesh& mesh)
  71         :
  72             surfaceInterpolationScheme<Type>(mesh)
  73         {}
  74
  75         //- Construct from Istream
  76         midPoint(const fvMesh& mesh, Istream&)
  77         :
  78             surfaceInterpolationScheme<Type>(mesh)
  79         {}
  80
  81         //- Construct from faceFlux and Istream
  82         midPoint
  83         (
  84             const fvMesh& mesh,
  85             const surfaceScalarField&,
  86             Istream&
  87         )
  88         :
  89             surfaceInterpolationScheme<Type>(mesh)
  90         {}
  91
  92
  93     // Member Functions
  94
  95         //- Return the interpolation weighting factors
  96         tmp<surfaceScalarField> weights
  97         (
  98             const GeometricField<Type, fvPatchField, volMesh>&
  99         ) const
 100         {
 101             tmp<surfaceScalarField> taw
 102             (
 103                 new surfaceScalarField
 104                 (
 105                     IOobject
 106                     (
 107                         "midPointWeights",
 108                         this->mesh().time().timeName(),
 109                         this->mesh()
 110                     ),
 111                     this->mesh(),
 112                     dimensionedScalar("0.5", dimless, 0.5)
 113                 )
 114             );
 115
 116             surfaceScalarField::GeometricBoundaryField& awbf =
 117                 taw().boundaryField();
 118
 119             forAll(awbf, patchi)
 120             {
 121                 if (!awbf[patchi].coupled())
 122                 {
 123                     awbf[patchi] = 1.0;
 124                 }
 125             }
 126
 127             return taw;
 128         }
 129 };
 130
 131
 132 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 133
 134 } // End namespace Foam
 135
 136 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 137
 138 #endif
 139
 140 // ************************************************************************* //