src/finiteVolume/fields/fvPatchFields/derived/movingWallVelocity/movingWallVelocityFvPatchVectorField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "movingWallVelocityFvPatchVectorField.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "volFields.H"
  29 #include "surfaceFields.H"
  30 #include "fvcMeshPhi.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36
  37 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  38
  39 movingWallVelocityFvPatchVectorField::movingWallVelocityFvPatchVectorField
  40 (
  41     const fvPatch& p,
  42     const DimensionedField<vector, volMesh>& iF
  43 )
  44 :
  45     fixedValueFvPatchVectorField(p, iF)
  46 {}
  47
  48
  49 movingWallVelocityFvPatchVectorField::movingWallVelocityFvPatchVectorField
  50 (
  51     const movingWallVelocityFvPatchVectorField& ptf,
  52     const fvPatch& p,
  53     const DimensionedField<vector, volMesh>& iF,
  54     const fvPatchFieldMapper& mapper
  55 )
  56 :
  57     fixedValueFvPatchVectorField(ptf, p, iF, mapper)
  58 {}
  59
  60
  61 movingWallVelocityFvPatchVectorField::movingWallVelocityFvPatchVectorField
  62 (
  63     const fvPatch& p,
  64     const DimensionedField<vector, volMesh>& iF,
  65     const dictionary& dict
  66 )
  67 :
  68     fixedValueFvPatchVectorField(p, iF)
  69 {
  70     fvPatchVectorField::operator=(vectorField("value", dict, p.size()));
  71 }
  72
  73
  74 movingWallVelocityFvPatchVectorField::movingWallVelocityFvPatchVectorField
  75 (
  76     const movingWallVelocityFvPatchVectorField& pivpvf
  77 )
  78 :
  79     fixedValueFvPatchVectorField(pivpvf)
  80 {}
  81
  82
  83 movingWallVelocityFvPatchVectorField::movingWallVelocityFvPatchVectorField
  84 (
  85     const movingWallVelocityFvPatchVectorField& pivpvf,
  86     const DimensionedField<vector, volMesh>& iF
  87 )
  88 :
  89     fixedValueFvPatchVectorField(pivpvf, iF)
  90 {}
  91
  92
  93 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  94
  95 void movingWallVelocityFvPatchVectorField::updateCoeffs()
  96 {
  97     if (updated())
  98     {
  99         return;
 100     }
 101
 102     const fvMesh& mesh = dimensionedInternalField().mesh();
 103
 104     if (mesh.changing())
 105     {
 106         const fvPatch& p = patch();
 107         const polyPatch& pp = p.patch();
 108         const pointField& oldPoints = mesh.oldPoints();
 109
 110         vectorField oldFc(pp.size());
 111
 112         forAll(oldFc, i)
 113         {
 114             oldFc[i] = pp[i].centre(oldPoints);
 115         }
 116
 117         // Get wall-parallel mesh motion velocity from geometry
 118         vectorField Up =
 119             (pp.faceCentres() - oldFc)/mesh.time().deltaT().value();
 120
 121         const volVectorField& U =
 122             mesh.lookupObject<volVectorField>
 123             (
 124                 dimensionedInternalField().name()
 125             );
 126
 127         scalarField phip =
 128             p.patchField<surfaceScalarField, scalar>(fvc::meshPhi(U));
 129
 130         vectorField n = p.nf();
 131         const scalarField& magSf = p.magSf();
 132         scalarField Un = phip/(magSf + VSMALL);
 133
 134         // Adjust for surface-normal mesh motion flux
 135         vectorField::operator=(Up + n*(Un - (n & Up)));
 136     }
 137
 138     fixedValueFvPatchVectorField::updateCoeffs();
 139 }
 140
 141
 142 void movingWallVelocityFvPatchVectorField::write(Ostream& os) const
 143 {
 144     fvPatchVectorField::write(os);
 145     writeEntry("value", os);
 146 }
 147
 148
 149 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 150
 151 makePatchTypeField
 152 (
 153     fvPatchVectorField,
 154     movingWallVelocityFvPatchVectorField
 155 );
 156
 157 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 158
 159 } // End namespace Foam
 160
 161 // ************************************************************************* //