src/thermophysicalModels/basic/derivedFvPatchFields/mixedEnthalpy/mixedEnthalpyFvPatchScalarField.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 "mixedEnthalpyFvPatchScalarField.H"
  27 #include "addToRunTimeSelectionTable.H"
  28 #include "fvPatchFieldMapper.H"
  29 #include "volFields.H"
  30 #include "basicThermo.H"
  31
  32 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  33
  34 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
  35 (
  36     const fvPatch& p,
  37     const DimensionedField<scalar, volMesh>& iF
  38 )
  39 :
  40     mixedFvPatchScalarField(p, iF)
  41 {
  42     valueFraction() = 0.0;
  43     refValue() = 0.0;
  44     refGrad() = 0.0;
  45 }
  46
  47
  48 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
  49 (
  50     const mixedEnthalpyFvPatchScalarField& ptf,
  51     const fvPatch& p,
  52     const DimensionedField<scalar, volMesh>& iF,
  53     const fvPatchFieldMapper& mapper
  54 )
  55 :
  56     mixedFvPatchScalarField(ptf, p, iF, mapper)
  57 {}
  58
  59
  60 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
  61 (
  62     const fvPatch& p,
  63     const DimensionedField<scalar, volMesh>& iF,
  64     const dictionary& dict
  65 )
  66 :
  67     mixedFvPatchScalarField(p, iF, dict)
  68 {}
  69
  70
  71 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
  72 (
  73     const mixedEnthalpyFvPatchScalarField& tppsf
  74 )
  75 :
  76     mixedFvPatchScalarField(tppsf)
  77 {}
  78
  79
  80 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
  81 (
  82     const mixedEnthalpyFvPatchScalarField& tppsf,
  83     const DimensionedField<scalar, volMesh>& iF
  84 )
  85 :
  86     mixedFvPatchScalarField(tppsf, iF)
  87 {}
  88
  89
  90 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  91
  92 void Foam::mixedEnthalpyFvPatchScalarField::updateCoeffs()
  93 {
  94     if (updated())
  95     {
  96         return;
  97     }
  98
  99     const basicThermo& thermo = db().lookupObject<basicThermo>
 100     (
 101         "thermophysicalProperties"
 102     );
 103
 104     const label patchi = patch().index();
 105
 106     mixedFvPatchScalarField& Tw = refCast<mixedFvPatchScalarField>
 107     (
 108         const_cast<fvPatchScalarField&>(thermo.T().boundaryField()[patchi])
 109     );
 110
 111     Tw.evaluate();
 112
 113     valueFraction() = Tw.valueFraction();
 114
 115     if
 116     (
 117         dimensionedInternalField().name() == db().mangleFileName("h")
 118     )
 119     {
 120         refValue() = thermo.h(Tw.refValue(), patchi);
 121         refGrad() = thermo.Cp(Tw, patchi)*Tw.refGrad()
 122           + patch().deltaCoeffs()*
 123             (
 124                 thermo.h(Tw, patchi)
 125               - thermo.h(Tw, patch().faceCells())
 126             );
 127     }
 128     else if
 129     (
 130         dimensionedInternalField().name() == db().mangleFileName("i")
 131     )
 132     {
 133     }
 134     else
 135     {
 136         refValue() = thermo.hs(Tw.refValue(), patchi);
 137         refGrad() = thermo.Cp(Tw, patchi)*Tw.refGrad()
 138           + patch().deltaCoeffs()*
 139             (
 140                 thermo.hs(Tw, patchi)
 141               - thermo.hs(Tw, patch().faceCells())
 142             );
 143     }
 144
 145     mixedFvPatchScalarField::updateCoeffs();
 146 }
 147
 148
 149 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 150
 151 namespace Foam
 152 {
 153     makePatchTypeField
 154     (
 155         fvPatchScalarField,
 156         mixedEnthalpyFvPatchScalarField
 157     );
 158 }
 159
 160
 161 // ************************************************************************* //