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