src/turbulenceModels/incompressible/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
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36 namespace incompressible
  37 {
  38
  39 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  40
  41 turbulentHeatFluxTemperatureFvPatchScalarField::
  42 turbulentHeatFluxTemperatureFvPatchScalarField
  43 (
  44     const fvPatch& p,
  45     const DimensionedField<scalar, volMesh>& iF
  46 )
  47 :
  48     fixedGradientFvPatchScalarField(p, iF),
  49     q_(p.size(), 0.0),
  50     alphaEffName_("undefinedAlphaEff"),
  51     CpName_("undefinedCp")
  52 {}
  53
  54
  55 turbulentHeatFluxTemperatureFvPatchScalarField::
  56 turbulentHeatFluxTemperatureFvPatchScalarField
  57 (
  58     const turbulentHeatFluxTemperatureFvPatchScalarField& ptf,
  59     const fvPatch& p,
  60     const DimensionedField<scalar, volMesh>& iF,
  61     const fvPatchFieldMapper& mapper
  62 )
  63 :
  64     fixedGradientFvPatchScalarField(ptf, p, iF, mapper),
  65     q_(ptf.q_, mapper),
  66     alphaEffName_(ptf.alphaEffName_),
  67     CpName_(ptf.CpName_)
  68 {}
  69
  70
  71 turbulentHeatFluxTemperatureFvPatchScalarField::
  72 turbulentHeatFluxTemperatureFvPatchScalarField
  73 (
  74     const fvPatch& p,
  75     const DimensionedField<scalar, volMesh>& iF,
  76     const dictionary& dict
  77 )
  78 :
  79     fixedGradientFvPatchScalarField(p, iF),
  80     q_("q", dict, p.size()),
  81     alphaEffName_(dict.lookup("alphaEff")),
  82     CpName_(dict.lookup("Cp"))
  83 {
  84     fvPatchField<scalar>::operator=(patchInternalField());
  85     gradient() = 0.0;
  86 }
  87
  88
  89 turbulentHeatFluxTemperatureFvPatchScalarField::
  90 turbulentHeatFluxTemperatureFvPatchScalarField
  91 (
  92     const turbulentHeatFluxTemperatureFvPatchScalarField& thftpsf
  93 )
  94 :
  95     fixedGradientFvPatchScalarField(thftpsf),
  96     q_(thftpsf.q_),
  97     alphaEffName_(thftpsf.alphaEffName_),
  98     CpName_(thftpsf.CpName_)
  99 {}
 100
 101
 102 turbulentHeatFluxTemperatureFvPatchScalarField::
 103 turbulentHeatFluxTemperatureFvPatchScalarField
 104 (
 105     const turbulentHeatFluxTemperatureFvPatchScalarField& thftpsf,
 106     const DimensionedField<scalar, volMesh>& iF
 107 )
 108 :
 109     fixedGradientFvPatchScalarField(thftpsf, iF),
 110     q_(thftpsf.q_),
 111     alphaEffName_(thftpsf.alphaEffName_),
 112     CpName_(thftpsf.CpName_)
 113 {}
 114
 115
 116 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 117
 118 void turbulentHeatFluxTemperatureFvPatchScalarField::autoMap
 119 (
 120     const fvPatchFieldMapper& m
 121 )
 122 {
 123     scalarField::autoMap(m);
 124     q_.autoMap(m);
 125 }
 126
 127
 128 void turbulentHeatFluxTemperatureFvPatchScalarField::rmap
 129 (
 130     const fvPatchScalarField& ptf,
 131     const labelList& addr
 132 )
 133 {
 134     fixedGradientFvPatchScalarField::rmap(ptf, addr);
 135
 136     const turbulentHeatFluxTemperatureFvPatchScalarField& thftptf =
 137         refCast<const turbulentHeatFluxTemperatureFvPatchScalarField>
 138         (
 139             ptf
 140         );
 141
 142     q_.rmap(thftptf.q_, addr);
 143 }
 144
 145
 146 void turbulentHeatFluxTemperatureFvPatchScalarField::updateCoeffs()
 147 {
 148     if (updated())
 149     {
 150         return;
 151     }
 152
 153     const scalarField& alphaEffp =
 154         patch().lookupPatchField<volScalarField, scalar>(alphaEffName_);
 155
 156     const scalarField& Cpp =
 157         patch().lookupPatchField<volScalarField, scalar>(CpName_);
 158
 159     gradient() = q_/(Cpp*alphaEffp);
 160
 161     fixedGradientFvPatchScalarField::updateCoeffs();
 162 }
 163
 164
 165 void turbulentHeatFluxTemperatureFvPatchScalarField::write(Ostream& os) const
 166 {
 167     fixedGradientFvPatchScalarField::write(os);
 168     q_.writeEntry("q", os);
 169     os.writeKeyword("alphaEff") << alphaEffName_ << token::END_STATEMENT << nl;
 170     os.writeKeyword("Cp") << CpName_ << token::END_STATEMENT << nl;
 171     writeEntry("value", os);
 172 }
 173
 174
 175 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 176
 177 makePatchTypeField
 178 (
 179     fvPatchScalarField,
 180     turbulentHeatFluxTemperatureFvPatchScalarField
 181 );
 182
 183
 184 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 185
 186 } // End namespace incompressible
 187 } // End namespace Foam
 188
 189
 190 // ************************************************************************* //
 191