src/thermophysicalModels/basic/derivedFvPatchFields/mixedInternalEnergy/mixedInternalEnergyFvPatchScalarField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "mixedInternalEnergyFvPatchScalarField.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "fvPatchFieldMapper.H"
  29 #include "volFields.H"
  30 #include "basicThermo.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36
  37 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  38
  39 mixedInternalEnergyFvPatchScalarField::mixedInternalEnergyFvPatchScalarField
  40 (
  41     const fvPatch& p,
  42     const DimensionedField<scalar, volMesh>& iF
  43 )
  44 :
  45     mixedFvPatchScalarField(p, iF)
  46 {
  47     valueFraction() = 0.0;
  48     refValue() = 0.0;
  49     refGrad() = 0.0;
  50 }
  51
  52
  53 mixedInternalEnergyFvPatchScalarField::mixedInternalEnergyFvPatchScalarField
  54 (
  55     const mixedInternalEnergyFvPatchScalarField& ptf,
  56     const fvPatch& p,
  57     const DimensionedField<scalar, volMesh>& iF,
  58     const fvPatchFieldMapper& mapper
  59 )
  60 :
  61     mixedFvPatchScalarField(ptf, p, iF, mapper)
  62 {}
  63
  64
  65 mixedInternalEnergyFvPatchScalarField::mixedInternalEnergyFvPatchScalarField
  66 (
  67     const fvPatch& p,
  68     const DimensionedField<scalar, volMesh>& iF,
  69     const dictionary& dict
  70 )
  71 :
  72     mixedFvPatchScalarField(p, iF, dict)
  73 {}
  74
  75
  76 mixedInternalEnergyFvPatchScalarField::mixedInternalEnergyFvPatchScalarField
  77 (
  78     const mixedInternalEnergyFvPatchScalarField& tppsf
  79 )
  80 :
  81     mixedFvPatchScalarField(tppsf)
  82 {}
  83
  84
  85 mixedInternalEnergyFvPatchScalarField::mixedInternalEnergyFvPatchScalarField
  86 (
  87     const mixedInternalEnergyFvPatchScalarField& tppsf,
  88     const DimensionedField<scalar, volMesh>& iF
  89 )
  90 :
  91     mixedFvPatchScalarField(tppsf, iF)
  92 {}
  93
  94
  95 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  96
  97 void mixedInternalEnergyFvPatchScalarField::updateCoeffs()
  98 {
  99     if (updated())
 100     {
 101         return;
 102     }
 103
 104     const basicThermo& thermo = db().lookupObject<basicThermo>
 105     (
 106         "thermophysicalProperties"
 107     );
 108
 109     const label patchi = patch().index();
 110
 111     mixedFvPatchScalarField& Tw = refCast<mixedFvPatchScalarField>
 112     (
 113         const_cast<fvPatchScalarField&>(thermo.T().boundaryField()[patchi])
 114     );
 115
 116     Tw.evaluate();
 117
 118     valueFraction() = Tw.valueFraction();
 119     refValue() = thermo.e(Tw.refValue(), patchi);
 120     refGrad() = thermo.Cv(Tw, patchi)*Tw.refGrad()
 121       + patch().deltaCoeffs()*
 122         (
 123             thermo.e(Tw, patchi)
 124           - thermo.e(Tw, patch().faceCells())
 125         );
 126
 127     mixedFvPatchScalarField::updateCoeffs();
 128 }
 129
 130
 131 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 132
 133 makePatchTypeField(fvPatchScalarField, mixedInternalEnergyFvPatchScalarField);
 134
 135 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 136
 137 } // End namespace Foam
 138
 139 // ************************************************************************* //