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