src/finiteVolume/fields/fvPatchFields/derived/fixedNormalSlip/fixedNormalSlipFvPatchField.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 "fixedNormalSlipFvPatchField.H"
  27 #include "symmTransformField.H"
  28
  29
  30 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  31
  32 template<class Type>
  33 Foam::fixedNormalSlipFvPatchField<Type>::fixedNormalSlipFvPatchField
  34 (
  35     const fvPatch& p,
  36     const DimensionedField<Type, volMesh>& iF
  37 )
  38 :
  39     transformFvPatchField<Type>(p, iF),
  40     fixedValue_(p.size(), pTraits<Type>::zero)
  41 {}
  42
  43
  44 template<class Type>
  45 Foam::fixedNormalSlipFvPatchField<Type>::fixedNormalSlipFvPatchField
  46 (
  47     const fixedNormalSlipFvPatchField<Type>& ptf,
  48     const fvPatch& p,
  49     const DimensionedField<Type, volMesh>& iF,
  50     const fvPatchFieldMapper& mapper
  51 )
  52 :
  53     transformFvPatchField<Type>(ptf, p, iF, mapper),
  54     fixedValue_(ptf.fixedValue_, mapper)
  55 {}
  56
  57
  58 template<class Type>
  59 Foam::fixedNormalSlipFvPatchField<Type>::fixedNormalSlipFvPatchField
  60 (
  61     const fvPatch& p,
  62     const DimensionedField<Type, volMesh>& iF,
  63     const dictionary& dict
  64 )
  65 :
  66     transformFvPatchField<Type>(p, iF),
  67     fixedValue_("fixedValue", dict, p.size())
  68 {
  69     evaluate();
  70 }
  71
  72
  73 template<class Type>
  74 Foam::fixedNormalSlipFvPatchField<Type>::fixedNormalSlipFvPatchField
  75 (
  76     const fixedNormalSlipFvPatchField<Type>& ptf
  77 )
  78 :
  79     transformFvPatchField<Type>(ptf),
  80     fixedValue_(ptf.fixedValue_)
  81 {}
  82
  83
  84 template<class Type>
  85 Foam::fixedNormalSlipFvPatchField<Type>::fixedNormalSlipFvPatchField
  86 (
  87     const fixedNormalSlipFvPatchField<Type>& ptf,
  88     const DimensionedField<Type, volMesh>& iF
  89 )
  90 :
  91     transformFvPatchField<Type>(ptf, iF),
  92     fixedValue_(ptf.fixedValue_)
  93 {}
  94
  95
  96 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  97
  98 template<class Type>
  99 void Foam::fixedNormalSlipFvPatchField<Type>::autoMap
 100 (
 101     const fvPatchFieldMapper& m
 102 )
 103 {
 104     transformFvPatchField<Type>::autoMap(m);
 105     fixedValue_.autoMap(m);
 106 }
 107
 108
 109 template<class Type>
 110 void Foam::fixedNormalSlipFvPatchField<Type>::rmap
 111 (
 112     const fvPatchField<Type>& ptf,
 113     const labelList& addr
 114 )
 115 {
 116     transformFvPatchField<Type>::rmap(ptf, addr);
 117
 118     const fixedNormalSlipFvPatchField<Type>& dmptf =
 119         refCast<const fixedNormalSlipFvPatchField<Type> >(ptf);
 120
 121     fixedValue_.rmap(dmptf.fixedValue_, addr);
 122 }
 123
 124
 125 template<class Type>
 126 Foam::tmp<Foam::Field<Type> >
 127 Foam::fixedNormalSlipFvPatchField<Type>::snGrad() const
 128 {
 129     const vectorField nHat(this->patch().nf());
 130     const Field<Type> pif(this->patchInternalField());
 131
 132     return
 133     (
 134         (nHat*(nHat & fixedValue_) + transform(I - sqr(nHat), pif)) - pif
 135     )*this->patch().deltaCoeffs();
 136 }
 137
 138
 139 template<class Type>
 140 void Foam::fixedNormalSlipFvPatchField<Type>::evaluate
 141 (
 142     const Pstream::commsTypes
 143 )
 144 {
 145     if (!this->updated())
 146     {
 147         this->updateCoeffs();
 148     }
 149
 150     const vectorField nHat(this->patch().nf());
 151
 152     Field<Type>::operator=
 153     (
 154         nHat*(nHat & fixedValue_)
 155       + transform(I - sqr(nHat), this->patchInternalField())
 156     );
 157
 158     transformFvPatchField<Type>::evaluate();
 159 }
 160
 161
 162 template<class Type>
 163 Foam::tmp<Foam::Field<Type> >
 164 Foam::fixedNormalSlipFvPatchField<Type>::snGradTransformDiag() const
 165 {
 166     const vectorField nHat(this->patch().nf());
 167     vectorField diag(nHat.size());
 168
 169     diag.replace(vector::X, mag(nHat.component(vector::X)));
 170     diag.replace(vector::Y, mag(nHat.component(vector::Y)));
 171     diag.replace(vector::Z, mag(nHat.component(vector::Z)));
 172
 173     return transformFieldMask<Type>(pow<vector, pTraits<Type>::rank>(diag));
 174 }
 175
 176
 177 template<class Type>
 178 void Foam::fixedNormalSlipFvPatchField<Type>::write(Ostream& os) const
 179 {
 180     transformFvPatchField<Type>::write(os);
 181     fixedValue_.writeEntry("fixedValue", os);
 182 }
 183
 184
 185 // ************************************************************************* //