src/finiteVolume/fields/fvPatchFields/derived/buoyantPressure/buoyantPressureFvPatchScalarField.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 "buoyantPressureFvPatchScalarField.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "fvPatchFieldMapper.H"
  30 #include "volFields.H"
  31 #include "uniformDimensionedFields.H"
  32
  33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37
  38 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  39
  40 buoyantPressureFvPatchScalarField::
  41 buoyantPressureFvPatchScalarField
  42 (
  43     const fvPatch& p,
  44     const DimensionedField<scalar, volMesh>& iF
  45 )
  46 :
  47     fixedGradientFvPatchScalarField(p, iF),
  48     rhoName_("rho")
  49 {}
  50
  51
  52 buoyantPressureFvPatchScalarField::
  53 buoyantPressureFvPatchScalarField
  54 (
  55     const fvPatch& p,
  56     const DimensionedField<scalar, volMesh>& iF,
  57     const dictionary& dict
  58 )
  59 :
  60     fixedGradientFvPatchScalarField(p, iF),
  61     rhoName_(dict.lookupOrDefault<word>("rho", "rho"))
  62 {
  63     fvPatchField<scalar>::operator=(patchInternalField());
  64     gradient() = 0.0;
  65 }
  66
  67
  68 buoyantPressureFvPatchScalarField::
  69 buoyantPressureFvPatchScalarField
  70 (
  71     const buoyantPressureFvPatchScalarField& ptf,
  72     const fvPatch& p,
  73     const DimensionedField<scalar, volMesh>& iF,
  74     const fvPatchFieldMapper& mapper
  75 )
  76 :
  77     fixedGradientFvPatchScalarField(ptf, p, iF, mapper),
  78     rhoName_(ptf.rhoName_)
  79 {}
  80
  81
  82 buoyantPressureFvPatchScalarField::
  83 buoyantPressureFvPatchScalarField
  84 (
  85     const buoyantPressureFvPatchScalarField& ptf
  86 )
  87 :
  88     fixedGradientFvPatchScalarField(ptf),
  89     rhoName_(ptf.rhoName_)
  90 {}
  91
  92
  93 buoyantPressureFvPatchScalarField::
  94 buoyantPressureFvPatchScalarField
  95 (
  96     const buoyantPressureFvPatchScalarField& ptf,
  97     const DimensionedField<scalar, volMesh>& iF
  98 )
  99 :
 100     fixedGradientFvPatchScalarField(ptf, iF),
 101     rhoName_(ptf.rhoName_)
 102 {}
 103
 104
 105 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 106
 107 void buoyantPressureFvPatchScalarField::updateCoeffs()
 108 {
 109     if (updated())
 110     {
 111         return;
 112     }
 113
 114     const uniformDimensionedVectorField& g =
 115         db().lookupObject<uniformDimensionedVectorField>("g");
 116
 117     const fvPatchField<scalar>& rho =
 118         patch().lookupPatchField<volScalarField, scalar>(rhoName_);
 119
 120     // If the variable name is "p_rgh" or "pd" assume it is p - rho*g.h
 121     // and set the gradient appropriately.
 122     // Otherwise assume the variable is the static pressure.
 123     if
 124     (
 125         dimensionedInternalField().name() == "p_rgh"
 126      || dimensionedInternalField().name() == "pd"
 127     )
 128     {
 129         gradient() = -rho.snGrad()*(g.value() & patch().Cf());
 130     }
 131     else
 132     {
 133         gradient() = rho*(g.value() & patch().nf());
 134     }
 135
 136     fixedGradientFvPatchScalarField::updateCoeffs();
 137 }
 138
 139
 140 void buoyantPressureFvPatchScalarField::write(Ostream& os) const
 141 {
 142     fixedGradientFvPatchScalarField::write(os);
 143     os.writeKeyword("rho") << rhoName_ << token::END_STATEMENT << nl;
 144     writeEntry("value", os);
 145 }
 146
 147
 148 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 149
 150 makePatchTypeField
 151 (
 152     fvPatchScalarField,
 153     buoyantPressureFvPatchScalarField
 154 );
 155
 156 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 157
 158 } // End namespace Foam
 159
 160 // ************************************************************************* //