src/finiteVolume/fields/fvPatchFields/derived/cylindricalInletVelocity/cylindricalInletVelocityFvPatchVectorField.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 the
  13     Free Software Foundation; either version 3 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 "cylindricalInletVelocityFvPatchVectorField.H"
  28 #include "volFields.H"
  29 #include "addToRunTimeSelectionTable.H"
  30 #include "fvPatchFieldMapper.H"
  31 #include "surfaceFields.H"
  32 #include "mathematicalConstants.H"
  33
  34 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  35
  36 Foam::
  37 cylindricalInletVelocityFvPatchVectorField::
  38 cylindricalInletVelocityFvPatchVectorField
  39 (
  40     const fvPatch& p,
  41     const DimensionedField<vector, volMesh>& iF
  42 )
  43 :
  44     fixedValueFvPatchField<vector>(p, iF),
  45     axialVelocity_(0),
  46     centre_(pTraits<vector>::zero),
  47     axis_(pTraits<vector>::zero),
  48     rpm_(0),
  49     radialVelocity_(0)
  50 {}
  51
  52
  53 Foam::
  54 cylindricalInletVelocityFvPatchVectorField::
  55 cylindricalInletVelocityFvPatchVectorField
  56 (
  57     const cylindricalInletVelocityFvPatchVectorField& ptf,
  58     const fvPatch& p,
  59     const DimensionedField<vector, volMesh>& iF,
  60     const fvPatchFieldMapper& mapper
  61 )
  62 :
  63     fixedValueFvPatchField<vector>(ptf, p, iF, mapper),
  64     axialVelocity_(ptf.axialVelocity_),
  65     centre_(ptf.centre_),
  66     axis_(ptf.axis_),
  67     rpm_(ptf.rpm_),
  68     radialVelocity_(ptf.radialVelocity_)
  69 {}
  70
  71
  72 Foam::
  73 cylindricalInletVelocityFvPatchVectorField::
  74 cylindricalInletVelocityFvPatchVectorField
  75 (
  76     const fvPatch& p,
  77     const DimensionedField<vector, volMesh>& iF,
  78     const dictionary& dict
  79 )
  80 :
  81     fixedValueFvPatchField<vector>(p, iF, dict),
  82     axialVelocity_(readScalar(dict.lookup("axialVelocity"))),
  83     centre_(dict.lookup("centre")),
  84     axis_(dict.lookup("axis")),
  85     rpm_(readScalar(dict.lookup("rpm"))),
  86     radialVelocity_(readScalar(dict.lookup("radialVelocity")))
  87 {}
  88
  89
  90 Foam::
  91 cylindricalInletVelocityFvPatchVectorField::
  92 cylindricalInletVelocityFvPatchVectorField
  93 (
  94     const cylindricalInletVelocityFvPatchVectorField& ptf
  95 )
  96 :
  97     fixedValueFvPatchField<vector>(ptf),
  98     axialVelocity_(ptf.axialVelocity_),
  99     centre_(ptf.centre_),
 100     axis_(ptf.axis_),
 101     rpm_(ptf.rpm_),
 102     radialVelocity_(ptf.radialVelocity_)
 103 {}
 104
 105
 106 Foam::
 107 cylindricalInletVelocityFvPatchVectorField::
 108 cylindricalInletVelocityFvPatchVectorField
 109 (
 110     const cylindricalInletVelocityFvPatchVectorField& ptf,
 111     const DimensionedField<vector, volMesh>& iF
 112 )
 113 :
 114     fixedValueFvPatchField<vector>(ptf, iF),
 115     axialVelocity_(ptf.axialVelocity_),
 116     centre_(ptf.centre_),
 117     axis_(ptf.axis_),
 118     rpm_(ptf.rpm_),
 119     radialVelocity_(ptf.radialVelocity_)
 120 {}
 121
 122
 123 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 124
 125 void Foam::cylindricalInletVelocityFvPatchVectorField::updateCoeffs()
 126 {
 127     if (updated())
 128     {
 129         return;
 130     }
 131
 132     vector hatAxis = axis_/mag(axis_);
 133
 134     const vectorField r(patch().Cf() - centre_);
 135     tmp<vectorField> d =  r - (hatAxis & r)*hatAxis;
 136
 137     tmp<vectorField> tangVel =
 138         (rpm_*constant::mathematical::pi/30.0)*(hatAxis) ^ d;
 139
 140     operator==(tangVel + axis_*axialVelocity_ + radialVelocity_*d);
 141
 142     fixedValueFvPatchField<vector>::updateCoeffs();
 143 }
 144
 145
 146 void Foam::cylindricalInletVelocityFvPatchVectorField::write(Ostream& os) const
 147 {
 148     fvPatchField<vector>::write(os);
 149     os.writeKeyword("axialVelocity") << axialVelocity_ <<
 150         token::END_STATEMENT << nl;
 151     os.writeKeyword("centre") << centre_ << token::END_STATEMENT << nl;
 152     os.writeKeyword("axis") << axis_ << token::END_STATEMENT << nl;
 153     os.writeKeyword("rpm") << rpm_ << token::END_STATEMENT << nl;
 154     os.writeKeyword("radialVelocity") << radialVelocity_ <<
 155         token::END_STATEMENT << nl;
 156     writeEntry("value", os);
 157 }
 158
 159
 160 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 161
 162 namespace Foam
 163 {
 164    makePatchTypeField
 165    (
 166        fvPatchVectorField,
 167        cylindricalInletVelocityFvPatchVectorField
 168    );
 169 }
 170
 171
 172 // ************************************************************************* //