src/finiteVolume/fields/fvPatchFields/basic/directionMixed/directionMixedFvPatchField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-2010 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
  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 "directionMixedFvPatchField.H"
  27 #include "symmTransformField.H"
  28
  29 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  30
  31 template<class Type>
  32 Foam::directionMixedFvPatchField<Type>::directionMixedFvPatchField
  33 (
  34     const fvPatch& p,
  35     const DimensionedField<Type, volMesh>& iF
  36 )
  37 :
  38     transformFvPatchField<Type>(p, iF),
  39     refValue_(p.size()),
  40     refGrad_(p.size()),
  41     valueFraction_(p.size())
  42 {}
  43
  44
  45 template<class Type>
  46 Foam::directionMixedFvPatchField<Type>::directionMixedFvPatchField
  47 (
  48     const directionMixedFvPatchField<Type>& ptf,
  49     const fvPatch& p,
  50     const DimensionedField<Type, volMesh>& iF,
  51     const fvPatchFieldMapper& mapper
  52 )
  53 :
  54     transformFvPatchField<Type>(ptf, p, iF, mapper),
  55     refValue_(ptf.refValue_, mapper),
  56     refGrad_(ptf.refGrad_, mapper),
  57     valueFraction_(ptf.valueFraction_, mapper)
  58 {}
  59
  60
  61 template<class Type>
  62 Foam::directionMixedFvPatchField<Type>::directionMixedFvPatchField
  63 (
  64     const fvPatch& p,
  65     const DimensionedField<Type, volMesh>& iF,
  66     const dictionary& dict
  67 )
  68 :
  69     transformFvPatchField<Type>(p, iF, dict),
  70     refValue_("refValue", dict, p.size()),
  71     refGrad_("refGradient", dict, p.size()),
  72     valueFraction_("valueFraction", dict, p.size())
  73 {
  74     evaluate();
  75 }
  76
  77
  78 template<class Type>
  79 Foam::directionMixedFvPatchField<Type>::directionMixedFvPatchField
  80 (
  81     const directionMixedFvPatchField<Type>& ptf,
  82     const DimensionedField<Type, volMesh>& iF
  83 )
  84 :
  85     transformFvPatchField<Type>(ptf, iF),
  86     refValue_(ptf.refValue_),
  87     refGrad_(ptf.refGrad_),
  88     valueFraction_(ptf.valueFraction_)
  89 {}
  90
  91
  92 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  93
  94 template<class Type>
  95 void Foam::directionMixedFvPatchField<Type>::autoMap
  96 (
  97     const fvPatchFieldMapper& m
  98 )
  99 {
 100     transformFvPatchField<Type>::autoMap(m);
 101     refValue_.autoMap(m);
 102     refGrad_.autoMap(m);
 103     valueFraction_.autoMap(m);
 104 }
 105
 106
 107 template<class Type>
 108 void Foam::directionMixedFvPatchField<Type>::rmap
 109 (
 110     const fvPatchField<Type>& ptf,
 111     const labelList& addr
 112 )
 113 {
 114     transformFvPatchField<Type>::rmap(ptf, addr);
 115
 116     const directionMixedFvPatchField<Type>& dmptf =
 117         refCast<const directionMixedFvPatchField<Type> >(ptf);
 118
 119     refValue_.rmap(dmptf.refValue_, addr);
 120     refGrad_.rmap(dmptf.refGrad_, addr);
 121     valueFraction_.rmap(dmptf.valueFraction_, addr);
 122 }
 123
 124
 125 template<class Type>
 126 Foam::tmp<Foam::Field<Type> >
 127 Foam::directionMixedFvPatchField<Type>::snGrad() const
 128 {
 129     const Field<Type> pif(this->patchInternalField());
 130
 131     tmp<Field<Type> > normalValue = transform(valueFraction_, refValue_);
 132
 133     tmp<Field<Type> > gradValue = pif + refGrad_/this->patch().deltaCoeffs();
 134
 135     tmp<Field<Type> > transformGradValue =
 136         transform(I - valueFraction_, gradValue);
 137
 138     return
 139         (normalValue + transformGradValue - pif)*
 140         this->patch().deltaCoeffs();
 141 }
 142
 143
 144 template<class Type>
 145 void Foam::directionMixedFvPatchField<Type>::evaluate(const Pstream::commsTypes)
 146 {
 147     if (!this->updated())
 148     {
 149         this->updateCoeffs();
 150     }
 151
 152     tmp<Field<Type> > normalValue = transform(valueFraction_, refValue_);
 153
 154     tmp<Field<Type> > gradValue =
 155         this->patchInternalField() + refGrad_/this->patch().deltaCoeffs();
 156
 157     tmp<Field<Type> > transformGradValue =
 158         transform(I - valueFraction_, gradValue);
 159
 160     Field<Type>::operator=(normalValue + transformGradValue);
 161
 162     transformFvPatchField<Type>::evaluate();
 163 }
 164
 165
 166 template<class Type>
 167 Foam::tmp<Foam::Field<Type> >
 168 Foam::directionMixedFvPatchField<Type>::snGradTransformDiag() const
 169 {
 170     vectorField diag(valueFraction_.size());
 171
 172     diag.replace
 173     (
 174         vector::X,
 175         sqrt(mag(valueFraction_.component(symmTensor::XX)))
 176     );
 177     diag.replace
 178     (
 179         vector::Y,
 180         sqrt(mag(valueFraction_.component(symmTensor::YY)))
 181     );
 182     diag.replace
 183     (
 184         vector::Z,
 185         sqrt(mag(valueFraction_.component(symmTensor::ZZ)))
 186     );
 187
 188     return transformFieldMask<Type>(pow<vector, pTraits<Type>::rank>(diag));
 189 }
 190
 191
 192 template<class Type>
 193 void Foam::directionMixedFvPatchField<Type>::write(Ostream& os) const
 194 {
 195     transformFvPatchField<Type>::write(os);
 196     refValue_.writeEntry("refValue", os);
 197     refGrad_.writeEntry("refGradient", os);
 198     valueFraction_.writeEntry("valueFraction", os);
 199     this->writeEntry("value", os);
 200 }
 201
 202
 203 // ************************************************************************* //