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