src/turbulenceModels/compressible/turbulenceModel/derivedFvPatchFields/turbulentHeatFluxTemperature/turbulentHeatFluxTemperatureFvPatchScalarField.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
  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 "turbulentHeatFluxTemperatureFvPatchScalarField.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "fvPatchFieldMapper.H"
  29 #include "volFields.H"
  30
  31
  32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36     // declare specialization within 'Foam' namespace
  37     template<>
  38     const char* NamedEnum
  39     <
  40         Foam::compressible::
  41         turbulentHeatFluxTemperatureFvPatchScalarField::heatSourceType,
  42         2
  43     >::names[] =
  44     {
  45         "power",
  46         "flux"
  47     };
  48
  49 }
  50
  51
  52 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  53
  54
  55 namespace Foam
  56 {
  57
  58 namespace compressible
  59 {
  60
  61 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  62
  63 const NamedEnum
  64 <
  65     turbulentHeatFluxTemperatureFvPatchScalarField::heatSourceType,
  66     2
  67 > turbulentHeatFluxTemperatureFvPatchScalarField::heatSourceTypeNames_;
  68
  69
  70 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  71
  72 turbulentHeatFluxTemperatureFvPatchScalarField::
  73 turbulentHeatFluxTemperatureFvPatchScalarField
  74 (
  75     const fvPatch& p,
  76     const DimensionedField<scalar, volMesh>& iF
  77 )
  78 :
  79     fixedGradientFvPatchScalarField(p, iF),
  80     temperatureCoupledBase(patch(), "undefined", "undefined-K"),
  81     heatSource_(hsPower),
  82     q_(p.size(), 0.0)
  83 {}
  84
  85
  86 turbulentHeatFluxTemperatureFvPatchScalarField::
  87 turbulentHeatFluxTemperatureFvPatchScalarField
  88 (
  89     const turbulentHeatFluxTemperatureFvPatchScalarField& ptf,
  90     const fvPatch& p,
  91     const DimensionedField<scalar, volMesh>& iF,
  92     const fvPatchFieldMapper& mapper
  93 )
  94 :
  95     fixedGradientFvPatchScalarField(ptf, p, iF, mapper),
  96     temperatureCoupledBase(patch(), ptf.KMethod(), ptf.KName()),
  97     heatSource_(ptf.heatSource_),
  98     q_(ptf.q_, mapper)
  99 {}
 100
 101
 102 turbulentHeatFluxTemperatureFvPatchScalarField::
 103 turbulentHeatFluxTemperatureFvPatchScalarField
 104 (
 105     const fvPatch& p,
 106     const DimensionedField<scalar, volMesh>& iF,
 107     const dictionary& dict
 108 )
 109 :
 110     fixedGradientFvPatchScalarField(p, iF),
 111     temperatureCoupledBase(patch(), dict),
 112     heatSource_(heatSourceTypeNames_.read(dict.lookup("heatSource"))),
 113     q_("q", dict, p.size())
 114 {
 115     fvPatchField<scalar>::operator=(patchInternalField());
 116     gradient() = 0.0;
 117 }
 118
 119
 120 turbulentHeatFluxTemperatureFvPatchScalarField::
 121 turbulentHeatFluxTemperatureFvPatchScalarField
 122 (
 123     const turbulentHeatFluxTemperatureFvPatchScalarField& thftpsf
 124 )
 125 :
 126     fixedGradientFvPatchScalarField(thftpsf),
 127     temperatureCoupledBase(patch(), thftpsf.KMethod(), thftpsf.KName()),
 128     heatSource_(thftpsf.heatSource_),
 129     q_(thftpsf.q_)
 130 {}
 131
 132
 133 turbulentHeatFluxTemperatureFvPatchScalarField::
 134 turbulentHeatFluxTemperatureFvPatchScalarField
 135 (
 136     const turbulentHeatFluxTemperatureFvPatchScalarField& thftpsf,
 137     const DimensionedField<scalar, volMesh>& iF
 138 )
 139 :
 140     fixedGradientFvPatchScalarField(thftpsf, iF),
 141     temperatureCoupledBase(patch(), thftpsf.KMethod(), thftpsf.KName()),
 142     heatSource_(thftpsf.heatSource_),
 143     q_(thftpsf.q_)
 144 {}
 145
 146
 147 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 148
 149 void turbulentHeatFluxTemperatureFvPatchScalarField::autoMap
 150 (
 151     const fvPatchFieldMapper& m
 152 )
 153 {
 154     fixedGradientFvPatchScalarField::autoMap(m);
 155     q_.autoMap(m);
 156 }
 157
 158
 159 void turbulentHeatFluxTemperatureFvPatchScalarField::rmap
 160 (
 161     const fvPatchScalarField& ptf,
 162     const labelList& addr
 163 )
 164 {
 165     fixedGradientFvPatchScalarField::rmap(ptf, addr);
 166
 167     const turbulentHeatFluxTemperatureFvPatchScalarField& thftptf =
 168         refCast<const turbulentHeatFluxTemperatureFvPatchScalarField>
 169         (
 170             ptf
 171         );
 172
 173     q_.rmap(thftptf.q_, addr);
 174 }
 175
 176
 177 void turbulentHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
 178 {
 179     if (updated())
 180     {
 181         return;
 182     }
 183
 184     const scalarField& Tp = *this;
 185
 186     switch (heatSource_)
 187     {
 188         case hsPower:
 189         {
 190             const scalar Ap = gSum(patch().magSf());
 191             gradient() = q_/(Ap*K(Tp));
 192             break;
 193         }
 194         case hsFlux:
 195         {
 196             gradient() = q_/K(Tp);
 197             break;
 198         }
 199         default:
 200         {
 201             FatalErrorIn
 202             (
 203                 "turbulentHeatFluxTemperatureFvPatchScalarField::updateCoeffs()"
 204             )   << "Unknown heat source type. Valid types are: "
 205                 << heatSourceTypeNames_ << nl << exit(FatalError);
 206         }
 207     }
 208
 209     fixedGradientFvPatchScalarField::updateCoeffs();
 210 }
 211
 212
 213 void turbulentHeatFluxTemperatureFvPatchScalarField::write
 214 (
 215     Ostream& os
 216 ) const
 217 {
 218     fixedGradientFvPatchScalarField::write(os);
 219     os.writeKeyword("heatSource") << heatSourceTypeNames_[heatSource_]
 220         << token::END_STATEMENT << nl;
 221     temperatureCoupledBase::write(os);
 222     q_.writeEntry("q", os);
 223     writeEntry("value", os);
 224 }
 225
 226
 227 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 228
 229 makePatchTypeField
 230 (
 231     fvPatchScalarField,
 232     turbulentHeatFluxTemperatureFvPatchScalarField
 233 );
 234
 235
 236 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 237
 238 } // End namespace compressible
 239 } // End namespace Foam
 240
 241
 242 // ************************************************************************* //
 243