src/finiteVolume/fields/fvPatchFields/constraint/wedge/wedgeFvPatchField.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 "wedgeFvPatch.H"
  28 #include "wedgeFvPatchField.H"
  29 #include "transformField.H"
  30 #include "symmTransform.H"
  31 #include "diagTensor.H"
  32
  33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37
  38 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  39
  40 template<class Type>
  41 wedgeFvPatchField<Type>::wedgeFvPatchField
  42 (
  43     const fvPatch& p,
  44     const DimensionedField<Type, volMesh>& iF
  45 )
  46 :
  47     transformFvPatchField<Type>(p, iF)
  48 {}
  49
  50
  51 template<class Type>
  52 wedgeFvPatchField<Type>::wedgeFvPatchField
  53 (
  54     const wedgeFvPatchField<Type>& ptf,
  55     const fvPatch& p,
  56     const DimensionedField<Type, volMesh>& iF,
  57     const fvPatchFieldMapper& mapper
  58 )
  59 :
  60     transformFvPatchField<Type>(ptf, p, iF, mapper)
  61 {
  62     if (!isType<wedgeFvPatch>(this->patch()))
  63     {
  64         FatalErrorIn
  65         (
  66             "wedgeFvPatchField<Type>::wedgeFvPatchField\n"
  67             "(\n"
  68             "    const wedgeFvPatchField<Type>& ptf,\n"
  69             "    const fvPatch& p,\n"
  70             "    const DimensionedField<Type, volMesh>& iF,\n"
  71             "    const fvPatchFieldMapper& mapper\n"
  72             ")\n"
  73         )   << "\n    patch type '" << p.type()
  74             << "' not constraint type '" << typeName << "'"
  75             << "\n    for patch " << p.name()
  76             << " of field " << this->dimensionedInternalField().name()
  77             << " in file " << this->dimensionedInternalField().objectPath()
  78             << exit(FatalIOError);
  79     }
  80 }
  81
  82
  83 template<class Type>
  84 wedgeFvPatchField<Type>::wedgeFvPatchField
  85 (
  86     const fvPatch& p,
  87     const DimensionedField<Type, volMesh>& iF,
  88     const dictionary& dict
  89 )
  90 :
  91     transformFvPatchField<Type>(p, iF, dict)
  92 {
  93     if (!isType<wedgeFvPatch>(p))
  94     {
  95         FatalIOErrorIn
  96         (
  97             "wedgeFvPatchField<Type>::wedgeFvPatchField\n"
  98             "(\n"
  99             "    const fvPatch& p,\n"
 100             "    const Field<Type>& field,\n"
 101             "    dictionary& dict\n"
 102             ")\n",
 103             dict
 104         )   << "\n    patch type '" << p.type()
 105             << "' not constraint type '" << typeName << "'"
 106             << "\n    for patch " << p.name()
 107             << " of field " << this->dimensionedInternalField().name()
 108             << " in file " << this->dimensionedInternalField().objectPath()
 109             << exit(FatalIOError);
 110     }
 111
 112     evaluate();
 113 }
 114
 115
 116 template<class Type>
 117 wedgeFvPatchField<Type>::wedgeFvPatchField
 118 (
 119     const wedgeFvPatchField<Type>& ptf
 120 )
 121 :
 122     transformFvPatchField<Type>(ptf)
 123 {}
 124
 125
 126 template<class Type>
 127 wedgeFvPatchField<Type>::wedgeFvPatchField
 128 (
 129     const wedgeFvPatchField<Type>& ptf,
 130     const DimensionedField<Type, volMesh>& iF
 131 )
 132 :
 133     transformFvPatchField<Type>(ptf, iF)
 134 {}
 135
 136
 137 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 138
 139 template<class Type>
 140 tmp<Field<Type> > wedgeFvPatchField<Type>::snGrad() const
 141 {
 142     Field<Type> pif = this->patchInternalField();
 143     return
 144     (
 145         transform(refCast<const wedgeFvPatch>(this->patch()).cellT(), pif) - pif
 146     )*(0.5*this->patch().deltaCoeffs());
 147 }
 148
 149
 150 template<class Type>
 151 void wedgeFvPatchField<Type>::evaluate(const Pstream::commsTypes)
 152 {
 153     if (!this->updated())
 154     {
 155         this->updateCoeffs();
 156     }
 157
 158     fvPatchField<Type>::operator==
 159     (
 160         transform
 161         (
 162             refCast<const wedgeFvPatch>(this->patch()).faceT(),
 163             this->patchInternalField()
 164         )
 165     );
 166 }
 167
 168
 169 template<class Type>
 170 tmp<Field<Type> > wedgeFvPatchField<Type>::snGradTransformDiag() const
 171 {
 172     diagTensor diagT =
 173         0.5*diag(I - refCast<const wedgeFvPatch>(this->patch()).cellT());
 174
 175     return tmp<Field<Type> >
 176     (
 177         new Field<Type>
 178         (
 179             this->size(),
 180             transformMask<Type>
 181             (
 182                 pow
 183                 (
 184                     reinterpret_cast<const vector&>(diagT),
 185                     pTraits<typename powProduct<vector, pTraits<Type>::rank>
 186                     ::type>::zero
 187                 )
 188             )
 189         )
 190     );
 191 }
 192
 193
 194 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 195
 196 } // End namespace Foam
 197
 198 // ************************************************************************* //