src/finiteVolume/fields/fvPatchFields/derived/pressureInletOutletVelocity/pressureInletOutletVelocityFvPatchVectorField.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 "pressureInletOutletVelocityFvPatchVectorField.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "volFields.H"
  30 #include "surfaceFields.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36
  37 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  38
  39 pressureInletOutletVelocityFvPatchVectorField::
  40 pressureInletOutletVelocityFvPatchVectorField
  41 (
  42     const fvPatch& p,
  43     const DimensionedField<vector, volMesh>& iF
  44 )
  45 :
  46     directionMixedFvPatchVectorField(p, iF),
  47     phiName_("phi")
  48 {
  49     refValue() = vector::zero;
  50     refGrad() = vector::zero;
  51     valueFraction() = symmTensor::zero;
  52 }
  53
  54
  55 pressureInletOutletVelocityFvPatchVectorField::
  56 pressureInletOutletVelocityFvPatchVectorField
  57 (
  58     const pressureInletOutletVelocityFvPatchVectorField& ptf,
  59     const fvPatch& p,
  60     const DimensionedField<vector, volMesh>& iF,
  61     const fvPatchFieldMapper& mapper
  62 )
  63 :
  64     directionMixedFvPatchVectorField(ptf, p, iF, mapper),
  65     phiName_(ptf.phiName_)
  66 {
  67     if (ptf.tangentialVelocity_.size())
  68     {
  69         tangentialVelocity_ = mapper(ptf.tangentialVelocity_);
  70     }
  71 }
  72
  73
  74 pressureInletOutletVelocityFvPatchVectorField::
  75 pressureInletOutletVelocityFvPatchVectorField
  76 (
  77     const fvPatch& p,
  78     const DimensionedField<vector, volMesh>& iF,
  79     const dictionary& dict
  80 )
  81 :
  82     directionMixedFvPatchVectorField(p, iF),
  83     phiName_(dict.lookupOrDefault<word>("phi", "phi"))
  84 {
  85     fvPatchVectorField::operator=(vectorField("value", dict, p.size()));
  86
  87     if (dict.found("tangentialVelocity"))
  88     {
  89         setTangentialVelocity
  90         (
  91             vectorField("tangentialVelocity", dict, p.size())
  92         );
  93     }
  94     else
  95     {
  96         refValue() = vector::zero;
  97     }
  98
  99     refGrad() = vector::zero;
 100     valueFraction() = symmTensor::zero;
 101 }
 102
 103
 104 pressureInletOutletVelocityFvPatchVectorField::
 105 pressureInletOutletVelocityFvPatchVectorField
 106 (
 107     const pressureInletOutletVelocityFvPatchVectorField& pivpvf
 108 )
 109 :
 110     directionMixedFvPatchVectorField(pivpvf),
 111     phiName_(pivpvf.phiName_),
 112     tangentialVelocity_(pivpvf.tangentialVelocity_)
 113 {}
 114
 115
 116 pressureInletOutletVelocityFvPatchVectorField::
 117 pressureInletOutletVelocityFvPatchVectorField
 118 (
 119     const pressureInletOutletVelocityFvPatchVectorField& pivpvf,
 120     const DimensionedField<vector, volMesh>& iF
 121 )
 122 :
 123     directionMixedFvPatchVectorField(pivpvf, iF),
 124     phiName_(pivpvf.phiName_),
 125     tangentialVelocity_(pivpvf.tangentialVelocity_)
 126 {}
 127
 128
 129 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 130
 131 void pressureInletOutletVelocityFvPatchVectorField::
 132 setTangentialVelocity(const vectorField& tangentialVelocity)
 133 {
 134     tangentialVelocity_ = tangentialVelocity;
 135     vectorField n = patch().nf();
 136     refValue() = tangentialVelocity_ - n*(n & tangentialVelocity_);
 137 }
 138
 139
 140 void pressureInletOutletVelocityFvPatchVectorField::autoMap
 141 (
 142     const fvPatchFieldMapper& m
 143 )
 144 {
 145     directionMixedFvPatchVectorField::autoMap(m);
 146     if (tangentialVelocity_.size())
 147     {
 148         tangentialVelocity_.autoMap(m);
 149     }
 150 }
 151
 152
 153 void pressureInletOutletVelocityFvPatchVectorField::rmap
 154 (
 155     const fvPatchVectorField& ptf,
 156     const labelList& addr
 157 )
 158 {
 159     directionMixedFvPatchVectorField::rmap(ptf, addr);
 160
 161     if (tangentialVelocity_.size())
 162     {
 163         const pressureInletOutletVelocityFvPatchVectorField& tiptf =
 164             refCast<const pressureInletOutletVelocityFvPatchVectorField>(ptf);
 165
 166         tangentialVelocity_.rmap(tiptf.tangentialVelocity_, addr);
 167     }
 168 }
 169
 170
 171 void pressureInletOutletVelocityFvPatchVectorField::updateCoeffs()
 172 {
 173     if (updated())
 174     {
 175         return;
 176     }
 177
 178     const fvsPatchField<scalar>& phip =
 179         patch().lookupPatchField<surfaceScalarField, scalar>(phiName_);
 180
 181     valueFraction() = neg(phip)*(I - sqr(patch().nf()));
 182
 183     directionMixedFvPatchVectorField::updateCoeffs();
 184     directionMixedFvPatchVectorField::evaluate();
 185 }
 186
 187
 188 void pressureInletOutletVelocityFvPatchVectorField::write(Ostream& os) const
 189 {
 190     fvPatchVectorField::write(os);
 191     if (phiName_ != "phi")
 192     {
 193         os.writeKeyword("phi") << phiName_ << token::END_STATEMENT << nl;
 194     }
 195     if (tangentialVelocity_.size())
 196     {
 197         tangentialVelocity_.writeEntry("tangentialVelocity", os);
 198     }
 199     writeEntry("value", os);
 200 }
 201
 202
 203 // * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
 204
 205 void pressureInletOutletVelocityFvPatchVectorField::operator=
 206 (
 207     const fvPatchField<vector>& pvf
 208 )
 209 {
 210     vectorField normalValue = transform(valueFraction(), refValue());
 211     vectorField transformGradValue = transform(I - valueFraction(), pvf);
 212     fvPatchField<vector>::operator=(normalValue + transformGradValue);
 213 }
 214
 215
 216 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 217
 218 makePatchTypeField
 219 (
 220     fvPatchVectorField,
 221     pressureInletOutletVelocityFvPatchVectorField
 222 );
 223
 224 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 225
 226 } // End namespace Foam
 227
 228 // ************************************************************************* //