applications/solvers/compressible/rhopSonicFoam/BCs/rhoE/mixedRhoEFvPatchScalarField.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 "mixedRhoEFvPatchScalarField.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "fvPatchFieldMapper.H"
  30 #include "volFields.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36
  37 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  38
  39 mixedRhoEFvPatchScalarField::mixedRhoEFvPatchScalarField
  40 (
  41     const fvPatch& p,
  42     const DimensionedField<scalar, volMesh>& iF
  43 )
  44 :
  45     mixedFvPatchScalarField(p, iF)
  46 {
  47     refValue() = 0.0;
  48     refGrad() = 0.0;
  49     valueFraction() = 0.0;
  50 }
  51
  52
  53 mixedRhoEFvPatchScalarField::mixedRhoEFvPatchScalarField
  54 (
  55     const mixedRhoEFvPatchScalarField& 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 mixedRhoEFvPatchScalarField::mixedRhoEFvPatchScalarField
  66 (
  67     const fvPatch& p,
  68     const DimensionedField<scalar, volMesh>& iF,
  69     const dictionary& dict
  70 )
  71 :
  72     mixedFvPatchScalarField(p, iF)
  73 {
  74     if (dict.found("value"))
  75     {
  76         fvPatchField<scalar>::operator=
  77         (
  78             scalarField("value", dict, p.size())
  79         );
  80     }
  81     else
  82     {
  83         fvPatchField<scalar>::operator=(patchInternalField());
  84     }
  85
  86     refValue() = *this;
  87     refGrad() = 0.0;
  88     valueFraction() = 0.0;
  89 }
  90
  91
  92 mixedRhoEFvPatchScalarField::mixedRhoEFvPatchScalarField
  93 (
  94     const mixedRhoEFvPatchScalarField& ptpsf
  95 )
  96 :
  97     mixedFvPatchScalarField(ptpsf)
  98 {}
  99
 100
 101 mixedRhoEFvPatchScalarField::mixedRhoEFvPatchScalarField
 102 (
 103     const mixedRhoEFvPatchScalarField& ptpsf,
 104     const DimensionedField<scalar, volMesh>& iF
 105 )
 106 :
 107     mixedFvPatchScalarField(ptpsf, iF)
 108 {}
 109
 110
 111 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 112
 113 void mixedRhoEFvPatchScalarField::autoMap
 114 (
 115     const fvPatchFieldMapper& m
 116 )
 117 {
 118     mixedFvPatchScalarField::autoMap(m);
 119 }
 120
 121
 122 void mixedRhoEFvPatchScalarField::rmap
 123 (
 124     const fvPatchField<scalar>& ptf,
 125     const labelList& addr
 126 )
 127 {
 128     mixedFvPatchField<scalar>::rmap(ptf, addr);
 129 }
 130
 131
 132 void mixedRhoEFvPatchScalarField::updateCoeffs()
 133 {
 134     if (updated())
 135     {
 136         return;
 137     }
 138
 139     const fvPatchField<scalar>& rhop =
 140         patch().lookupPatchField<volScalarField, scalar>("rho");
 141     const fvPatchField<vector>& rhoUp =
 142         patch().lookupPatchField<volVectorField, vector>("rhoU");
 143 //    fvPatchField<scalar>& Tp =
 144 //        patch().lookupPatchField<volScalarField, scalar>("T");
 145
 146     const volScalarField& T = db().lookupObject<volScalarField>("T");
 147     const label patchi = patch().index();
 148     fvPatchScalarField& Tp =
 149         const_cast<fvPatchScalarField&>(T.boundaryField()[patchi]);
 150
 151     Tp.evaluate();
 152
 153     const dictionary& thermodynamicProperties = db().lookupObject<IOdictionary>
 154     (
 155         "thermodynamicProperties"
 156     );
 157
 158     dimensionedScalar Cv(thermodynamicProperties.lookup("Cv"));
 159
 160     valueFraction() = rhop.snGrad()/
 161         (rhop.snGrad() - rhop*this->patch().deltaCoeffs());
 162
 163     refValue() = 0.5*rhop*magSqr(rhoUp/rhop);
 164     refGrad() =
 165         rhop*Cv.value()*Tp.snGrad()
 166       + (
 167             refValue()
 168           - (0.5*rhop.patchInternalField()*
 169                 magSqr(rhoUp.patchInternalField()/rhop.patchInternalField()))
 170         )*patch().deltaCoeffs();
 171
 172     mixedFvPatchScalarField::updateCoeffs();
 173 }
 174
 175
 176 void mixedRhoEFvPatchScalarField::write(Ostream& os) const
 177 {
 178     fvPatchScalarField::write(os);
 179     os.writeKeyword("valueFraction") << valueFraction()
 180         << token::END_STATEMENT << endl;
 181     os.writeKeyword("refValue") << refValue() << token::END_STATEMENT << endl;
 182     os.writeKeyword("refGrad") << refGrad() << token::END_STATEMENT << endl;
 183     writeEntry("value", os);
 184 }
 185
 186
 187 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 188
 189 makePatchTypeField(fvPatchScalarField, mixedRhoEFvPatchScalarField);
 190
 191 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 192
 193 } // End namespace Foam
 194
 195 // ************************************************************************* //