src/dynamicMesh/meshMotion/fvMotionSolver/pointPatchFields/derived/angularOscillatingVelocity/angularOscillatingVelocityPointPatchVectorField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     4.0
   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 "angularOscillatingVelocityPointPatchVectorField.H"
  27 #include "pointPatchFields.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "foamTime.H"
  30 #include "polyMesh.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36
  37 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  38
  39 angularOscillatingVelocityPointPatchVectorField::
  40 angularOscillatingVelocityPointPatchVectorField
  41 (
  42     const pointPatch& p,
  43     const DimensionedField<vector, pointMesh>& iF
  44 )
  45 :
  46     fixedValuePointPatchVectorField(p, iF),
  47     axis_(vector::zero),
  48     origin_(vector::zero),
  49     angle0_(0.0),
  50     amplitude_(0.0),
  51     omega_(0.0),
  52     p0_(p.localPoints())
  53 {}
  54
  55
  56 angularOscillatingVelocityPointPatchVectorField::
  57 angularOscillatingVelocityPointPatchVectorField
  58 (
  59     const pointPatch& p,
  60     const DimensionedField<vector, pointMesh>& iF,
  61     const dictionary& dict
  62 )
  63 :
  64     fixedValuePointPatchVectorField(p, iF, dict),
  65     axis_(dict.lookup("axis")),
  66     origin_(dict.lookup("origin")),
  67     angle0_(readScalar(dict.lookup("angle0"))),
  68     amplitude_(readScalar(dict.lookup("amplitude"))),
  69     omega_(readScalar(dict.lookup("omega")))
  70 {
  71     if (!dict.found("value"))
  72     {
  73         updateCoeffs();
  74     }
  75
  76     if (dict.found("p0"))
  77     {
  78         p0_ = vectorField("p0", dict , p.size());
  79     }
  80     else
  81     {
  82         p0_ = p.localPoints();
  83     }
  84 }
  85
  86
  87 angularOscillatingVelocityPointPatchVectorField::
  88 angularOscillatingVelocityPointPatchVectorField
  89 (
  90     const angularOscillatingVelocityPointPatchVectorField& ptf,
  91     const pointPatch& p,
  92     const DimensionedField<vector, pointMesh>& iF,
  93     const PointPatchFieldMapper& mapper
  94 )
  95 :
  96     fixedValuePointPatchVectorField(ptf, p, iF, mapper),
  97     axis_(ptf.axis_),
  98     origin_(ptf.origin_),
  99     angle0_(ptf.angle0_),
 100     amplitude_(ptf.amplitude_),
 101     omega_(ptf.omega_),
 102     p0_(ptf.p0_, mapper)
 103 {}
 104
 105
 106 angularOscillatingVelocityPointPatchVectorField::
 107 angularOscillatingVelocityPointPatchVectorField
 108 (
 109     const angularOscillatingVelocityPointPatchVectorField& ptf,
 110     const DimensionedField<vector, pointMesh>& iF
 111 )
 112 :
 113     fixedValuePointPatchVectorField(ptf, iF),
 114     axis_(ptf.axis_),
 115     origin_(ptf.origin_),
 116     angle0_(ptf.angle0_),
 117     amplitude_(ptf.amplitude_),
 118     omega_(ptf.omega_),
 119     p0_(ptf.p0_)
 120 {}
 121
 122
 123 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 124
 125 void angularOscillatingVelocityPointPatchVectorField::autoMap
 126 (
 127     const PointPatchFieldMapper& m
 128 )
 129 {
 130     fixedValuePointPatchVectorField::autoMap(m);
 131     p0_.autoMap(m);
 132 }
 133
 134 void angularOscillatingVelocityPointPatchVectorField::rmap
 135 (
 136     const PointPatchField
 137         <pointPatchField, pointMesh, pointPatch, DummyMatrix, vector>& ptf,
 138     const labelList& addr
 139 )
 140 {
 141     const angularOscillatingVelocityPointPatchVectorField& aOVptf =
 142             refCast<const angularOscillatingVelocityPointPatchVectorField>(ptf);
 143
 144     fixedValuePointPatchVectorField::rmap(aOVptf, addr);
 145     p0_.rmap(aOVptf.p0_, addr);
 146 }
 147
 148 void angularOscillatingVelocityPointPatchVectorField::updateCoeffs()
 149 {
 150     if (this->updated())
 151     {
 152         return;
 153     }
 154
 155     const polyMesh& mesh = this->dimensionedInternalField().mesh()();
 156     const Time& t = mesh.time();
 157     const pointPatch& p = this->patch();
 158
 159     scalar angle = angle0_ + amplitude_*sin(omega_*t.value());
 160     vector axisHat = axis_/mag(axis_);
 161     vectorField p0Rel = p0_ - origin_;
 162
 163     vectorField::operator=
 164     (
 165         (
 166             p0_
 167           + p0Rel*(cos(angle) - 1)
 168           + (axisHat ^ p0Rel*sin(angle))
 169           + (axisHat & p0Rel)*(1 - cos(angle))*axisHat
 170           - p.localPoints()
 171         )/t.deltaT().value()
 172     );
 173
 174     fixedValuePointPatchVectorField::updateCoeffs();
 175 }
 176
 177
 178 void angularOscillatingVelocityPointPatchVectorField::write
 179 (
 180     Ostream& os
 181 ) const
 182 {
 183     pointPatchField<vector>::write(os);
 184     os.writeKeyword("axis")
 185         << axis_ << token::END_STATEMENT << nl;
 186     os.writeKeyword("origin")
 187         << origin_ << token::END_STATEMENT << nl;
 188     os.writeKeyword("angle0")
 189         << angle0_ << token::END_STATEMENT << nl;
 190     os.writeKeyword("amplitude")
 191         << amplitude_ << token::END_STATEMENT << nl;
 192     os.writeKeyword("omega")
 193         << omega_ << token::END_STATEMENT << nl;
 194     p0_.writeEntry("p0", os);
 195     writeEntry("value", os);
 196 }
 197
 198
 199 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 200
 201 makePointPatchTypeField
 202 (
 203     pointPatchVectorField,
 204     angularOscillatingVelocityPointPatchVectorField
 205 );
 206
 207 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 208
 209 } // End namespace Foam
 210
 211 // ************************************************************************* //