src/finiteVolume/fvMesh/fvPatches/fvPatch/fvPatch.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "fvPatch.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "fvBoundaryMesh.H"
  30 #include "fvMesh.H"
  31 #include "primitiveMesh.H"
  32 #include "volFields.H"
  33 #include "surfaceFields.H"
  34
  35 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  36
  37 namespace Foam
  38 {
  39
  40 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  41
  42 defineTypeNameAndDebug(fvPatch, 0);
  43 defineRunTimeSelectionTable(fvPatch, polyPatch);
  44 addToRunTimeSelectionTable(fvPatch, fvPatch, polyPatch);
  45
  46
  47 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  48
  49 fvPatch::fvPatch(const polyPatch& p, const fvBoundaryMesh& bm)
  50 :
  51     polyPatch_(p),
  52     boundaryMesh_(bm)
  53 {}
  54
  55
  56 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  57
  58 fvPatch::~fvPatch()
  59 {}
  60
  61
  62 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  63
  64 bool fvPatch::constraintType(const word& pt)
  65 {
  66     return fvPatchField<scalar>::patchConstructorTablePtr_->found(pt);
  67 }
  68
  69
  70 wordList fvPatch::constraintTypes()
  71 {
  72     wordList cTypes(polyPatchConstructorTablePtr_->size());
  73
  74     label i = 0;
  75
  76     for
  77     (
  78         polyPatchConstructorTable::iterator cstrIter =
  79             polyPatchConstructorTablePtr_->begin();
  80         cstrIter != polyPatchConstructorTablePtr_->end();
  81         ++cstrIter
  82     )
  83     {
  84         if (constraintType(cstrIter.key()))
  85         {
  86             cTypes[i++] = cstrIter.key();
  87         }
  88     }
  89
  90     cTypes.setSize(i);
  91
  92     return cTypes;
  93 }
  94
  95
  96 const unallocLabelList& fvPatch::faceCells() const
  97 {
  98     return polyPatch_.faceCells();
  99 }
 100
 101
 102 const vectorField& fvPatch::Cf() const
 103 {
 104     return boundaryMesh().mesh().Cf().boundaryField()[index()];
 105 }
 106
 107
 108 tmp<vectorField> fvPatch::Cn() const
 109 {
 110     tmp<vectorField> tcc(new vectorField(size()));
 111     vectorField& cc = tcc();
 112
 113     const unallocLabelList& faceCells = this->faceCells();
 114
 115     // get reference to global cell centres
 116     const vectorField& gcc = boundaryMesh().mesh().cellCentres();
 117
 118     forAll (faceCells, faceI)
 119     {
 120         cc[faceI] = gcc[faceCells[faceI]];
 121     }
 122
 123     return tcc;
 124 }
 125
 126
 127 tmp<vectorField> fvPatch::nf() const
 128 {
 129     return Sf()/magSf();
 130 }
 131
 132
 133 const vectorField& fvPatch::Sf() const
 134 {
 135     return boundaryMesh().mesh().Sf().boundaryField()[index()];
 136 }
 137
 138
 139 const scalarField& fvPatch::magSf() const
 140 {
 141     return boundaryMesh().mesh().magSf().boundaryField()[index()];
 142 }
 143
 144
 145 tmp<vectorField> fvPatch::delta() const
 146 {
 147     return Cf() - Cn();
 148 }
 149
 150
 151 void fvPatch::makeWeights(scalarField& w) const
 152 {
 153     w = 1.0;
 154 }
 155
 156
 157 void fvPatch::makeDeltaCoeffs(scalarField& dc) const
 158 {
 159     dc = 1.0/(nf() & delta());
 160 }
 161
 162
 163 void fvPatch::initMovePoints()
 164 {}
 165
 166
 167 void fvPatch::movePoints()
 168 {}
 169
 170
 171 const scalarField& fvPatch::deltaCoeffs() const
 172 {
 173     return boundaryMesh().mesh().deltaCoeffs().boundaryField()[index()];
 174 }
 175
 176
 177 const scalarField& fvPatch::weights() const
 178 {
 179     return boundaryMesh().mesh().weights().boundaryField()[index()];
 180 }
 181
 182
 183 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 184
 185 } // End namespace Foam
 186
 187 // ************************************************************************* //