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