src/finiteVolume/fields/fvPatchFields/derived/inletOutletTotalTemperature/inletOutletTotalTemperatureFvPatchScalarField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-2011 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 "inletOutletTotalTemperatureFvPatchScalarField.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "fvPatchFieldMapper.H"
  29 #include "volFields.H"
  30 #include "surfaceFields.H"
  31
  32 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  33
  34 Foam::inletOutletTotalTemperatureFvPatchScalarField::
  35 inletOutletTotalTemperatureFvPatchScalarField
  36 (
  37     const fvPatch& p,
  38     const DimensionedField<scalar, volMesh>& iF
  39 )
  40 :
  41     inletOutletFvPatchScalarField(p, iF),
  42     UName_("U"),
  43     phiName_("phi"),
  44     psiName_("psi"),
  45     gamma_(0.0),
  46     T0_(p.size(), 0.0)
  47 {
  48     this->refValue() = pTraits<scalar>::zero;
  49     this->refGrad() = pTraits<scalar>::zero;
  50     this->valueFraction() = 0.0;
  51 }
  52
  53
  54 Foam::inletOutletTotalTemperatureFvPatchScalarField::
  55 inletOutletTotalTemperatureFvPatchScalarField
  56 (
  57     const inletOutletTotalTemperatureFvPatchScalarField& ptf,
  58     const fvPatch& p,
  59     const DimensionedField<scalar, volMesh>& iF,
  60     const fvPatchFieldMapper& mapper
  61 )
  62 :
  63     inletOutletFvPatchScalarField(ptf, p, iF, mapper),
  64     UName_(ptf.UName_),
  65     phiName_(ptf.phiName_),
  66     psiName_(ptf.psiName_),
  67     gamma_(ptf.gamma_),
  68     T0_(ptf.T0_, mapper)
  69 {}
  70
  71
  72 Foam::inletOutletTotalTemperatureFvPatchScalarField::
  73 inletOutletTotalTemperatureFvPatchScalarField
  74 (
  75     const fvPatch& p,
  76     const DimensionedField<scalar, volMesh>& iF,
  77     const dictionary& dict
  78 )
  79 :
  80     inletOutletFvPatchScalarField(p, iF),
  81     UName_(dict.lookupOrDefault<word>("U", "U")),
  82     phiName_(dict.lookupOrDefault<word>("phi", "phi")),
  83     psiName_(dict.lookupOrDefault<word>("psi", "psi")),
  84     gamma_(readScalar(dict.lookup("gamma"))),
  85     T0_("T0", dict, p.size())
  86 {
  87     this->refValue() = pTraits<scalar>::zero;
  88     if (dict.found("value"))
  89     {
  90         fvPatchField<scalar>::operator=
  91         (
  92             scalarField("value", dict, p.size())
  93         );
  94     }
  95     else
  96     {
  97         fvPatchField<scalar>::operator=(T0_);
  98     }
  99
 100     this->refGrad() = pTraits<scalar>::zero;
 101     this->valueFraction() = 0.0;
 102 }
 103
 104
 105 Foam::inletOutletTotalTemperatureFvPatchScalarField::
 106 inletOutletTotalTemperatureFvPatchScalarField
 107 (
 108     const inletOutletTotalTemperatureFvPatchScalarField& tppsf
 109 )
 110 :
 111     inletOutletFvPatchScalarField(tppsf),
 112     UName_(tppsf.UName_),
 113     phiName_(tppsf.phiName_),
 114     psiName_(tppsf.psiName_),
 115     gamma_(tppsf.gamma_),
 116     T0_(tppsf.T0_)
 117 {}
 118
 119
 120 Foam::inletOutletTotalTemperatureFvPatchScalarField::
 121 inletOutletTotalTemperatureFvPatchScalarField
 122 (
 123     const inletOutletTotalTemperatureFvPatchScalarField& tppsf,
 124     const DimensionedField<scalar, volMesh>& iF
 125 )
 126 :
 127     inletOutletFvPatchScalarField(tppsf, iF),
 128     UName_(tppsf.UName_),
 129     phiName_(tppsf.phiName_),
 130     psiName_(tppsf.psiName_),
 131     gamma_(tppsf.gamma_),
 132     T0_(tppsf.T0_)
 133 {}
 134
 135
 136 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 137
 138 void Foam::inletOutletTotalTemperatureFvPatchScalarField::autoMap
 139 (
 140     const fvPatchFieldMapper& m
 141 )
 142 {
 143     inletOutletFvPatchScalarField::autoMap(m);
 144     T0_.autoMap(m);
 145 }
 146
 147
 148 void Foam::inletOutletTotalTemperatureFvPatchScalarField::rmap
 149 (
 150     const fvPatchScalarField& ptf,
 151     const labelList& addr
 152 )
 153 {
 154     inletOutletFvPatchScalarField::rmap(ptf, addr);
 155
 156     const inletOutletTotalTemperatureFvPatchScalarField& tiptf =
 157         refCast<const inletOutletTotalTemperatureFvPatchScalarField>(ptf);
 158
 159     T0_.rmap(tiptf.T0_, addr);
 160 }
 161
 162
 163 void Foam::inletOutletTotalTemperatureFvPatchScalarField::updateCoeffs()
 164 {
 165     if (updated())
 166     {
 167         return;
 168     }
 169
 170     const fvPatchVectorField& Up =
 171         patch().lookupPatchField<volVectorField, vector>(UName_);
 172
 173     const fvsPatchField<scalar>& phip =
 174         patch().lookupPatchField<surfaceScalarField, scalar>(phiName_);
 175
 176     const fvPatchField<scalar>& psip =
 177         patch().lookupPatchField<volScalarField, scalar>(psiName_);
 178
 179     scalar gM1ByG = (gamma_ - 1.0)/gamma_;
 180
 181     this->refValue() =
 182         T0_/(1.0 + 0.5*psip*gM1ByG*(1.0 - pos(phip))*magSqr(Up));
 183     this->valueFraction() = 1.0 - pos(phip);
 184
 185     inletOutletFvPatchScalarField::updateCoeffs();
 186 }
 187
 188
 189 void Foam::inletOutletTotalTemperatureFvPatchScalarField::write(Ostream& os)
 190 const
 191 {
 192     fvPatchScalarField::write(os);
 193     writeEntryIfDifferent<word>(os, "U", "U", UName_);
 194     writeEntryIfDifferent<word>(os, "phi", "phi", phiName_);
 195     writeEntryIfDifferent<word>(os, "psi", "psi", psiName_);
 196     os.writeKeyword("gamma") << gamma_ << token::END_STATEMENT << nl;
 197     T0_.writeEntry("T0", os);
 198     writeEntry("value", os);
 199 }
 200
 201
 202 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 203
 204 namespace Foam
 205 {
 206     makePatchTypeField
 207     (
 208         fvPatchScalarField,
 209         inletOutletTotalTemperatureFvPatchScalarField
 210     );
 211 }
 212
 213 // ************************************************************************* //