1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2009 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
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.
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
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
25 \*---------------------------------------------------------------------------*/
27 #include "mixedEnthalpyFvPatchScalarField.H"
28 #include "addToRunTimeSelectionTable.H"
29 #include "fvPatchFieldMapper.H"
30 #include "volFields.H"
31 #include "basicThermo.H"
33 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
35 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
38 const DimensionedField<scalar, volMesh>& iF
41 mixedFvPatchScalarField(p, iF)
43 valueFraction() = 0.0;
49 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
51 const mixedEnthalpyFvPatchScalarField& ptf,
53 const DimensionedField<scalar, volMesh>& iF,
54 const fvPatchFieldMapper& mapper
57 mixedFvPatchScalarField(ptf, p, iF, mapper)
61 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
64 const DimensionedField<scalar, volMesh>& iF,
65 const dictionary& dict
68 mixedFvPatchScalarField(p, iF, dict)
72 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
74 const mixedEnthalpyFvPatchScalarField& tppsf
77 mixedFvPatchScalarField(tppsf)
81 Foam::mixedEnthalpyFvPatchScalarField::mixedEnthalpyFvPatchScalarField
83 const mixedEnthalpyFvPatchScalarField& tppsf,
84 const DimensionedField<scalar, volMesh>& iF
87 mixedFvPatchScalarField(tppsf, iF)
91 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
93 void Foam::mixedEnthalpyFvPatchScalarField::updateCoeffs()
100 const basicThermo& thermo = db().lookupObject<basicThermo>
102 "thermophysicalProperties"
105 const label patchi = patch().index();
107 mixedFvPatchScalarField& Tw = refCast<mixedFvPatchScalarField>
109 const_cast<fvPatchScalarField&>(thermo.T().boundaryField()[patchi])
114 valueFraction() = Tw.valueFraction();
116 if (dimensionedInternalField().name() == "h")
118 refValue() = thermo.h(Tw.refValue(), patchi);
119 refGrad() = thermo.Cp(Tw, patchi)*Tw.refGrad()
120 + patch().deltaCoeffs()*
123 - thermo.h(Tw, patch().faceCells())
128 refValue() = thermo.hs(Tw.refValue(), patchi);
129 refGrad() = thermo.Cp(Tw, patchi)*Tw.refGrad()
130 + patch().deltaCoeffs()*
132 thermo.hs(Tw, patchi)
133 - thermo.hs(Tw, patch().faceCells())
137 mixedFvPatchScalarField::updateCoeffs();
141 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
148 mixedEnthalpyFvPatchScalarField
153 // ************************************************************************* //