src/finiteVolume/fields/fvPatchFields/derived/fixedFluxPressure/fixedFluxPressureFvPatchScalarField.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 "fixedFluxPressureFvPatchScalarField.H"
  28 #include "fvPatchFieldMapper.H"
  29 #include "volFields.H"
  30 #include "surfaceFields.H"
  31 #include "addToRunTimeSelectionTable.H"
  32
  33 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  34
  35 Foam::fixedFluxPressureFvPatchScalarField::fixedFluxPressureFvPatchScalarField
  36 (
  37     const fvPatch& p,
  38     const DimensionedField<scalar, volMesh>& iF
  39 )
  40 :
  41     fixedGradientFvPatchScalarField(p, iF),
  42     UName_("U"),
  43     phiName_("phi"),
  44     rhoName_("rho"),
  45     adjoint_(false)
  46 {}
  47
  48
  49 Foam::fixedFluxPressureFvPatchScalarField::fixedFluxPressureFvPatchScalarField
  50 (
  51     const fixedFluxPressureFvPatchScalarField& ptf,
  52     const fvPatch& p,
  53     const DimensionedField<scalar, volMesh>& iF,
  54     const fvPatchFieldMapper& mapper
  55 )
  56 :
  57     fixedGradientFvPatchScalarField(ptf, p, iF, mapper),
  58     UName_(ptf.UName_),
  59     phiName_(ptf.phiName_),
  60     rhoName_(ptf.rhoName_),
  61     adjoint_(ptf.adjoint_)
  62 {}
  63
  64
  65 Foam::fixedFluxPressureFvPatchScalarField::fixedFluxPressureFvPatchScalarField
  66 (
  67     const fvPatch& p,
  68     const DimensionedField<scalar, volMesh>& iF,
  69     const dictionary& dict
  70 )
  71 :
  72     fixedGradientFvPatchScalarField(p, iF),
  73     UName_(dict.lookupOrDefault<word>("U", "U")),
  74     phiName_(dict.lookupOrDefault<word>("phi", "phi")),
  75     rhoName_(dict.lookupOrDefault<word>("rho", "rho")),
  76     adjoint_(dict.lookup("adjoint"))
  77 {
  78     if (dict.found("gradient"))
  79     {
  80         gradient() = scalarField("gradient", dict, p.size());
  81         fixedGradientFvPatchScalarField::updateCoeffs();
  82         fixedGradientFvPatchScalarField::evaluate();
  83     }
  84     else
  85     {
  86         fvPatchField<scalar>::operator=(patchInternalField());
  87         gradient() = 0.0;
  88     }
  89 }
  90
  91
  92 Foam::fixedFluxPressureFvPatchScalarField::fixedFluxPressureFvPatchScalarField
  93 (
  94     const fixedFluxPressureFvPatchScalarField& wbppsf
  95 )
  96 :
  97     fixedGradientFvPatchScalarField(wbppsf),
  98     UName_(wbppsf.UName_),
  99     phiName_(wbppsf.phiName_),
 100     rhoName_(wbppsf.rhoName_),
 101     adjoint_(wbppsf.adjoint_)
 102 {}
 103
 104
 105 Foam::fixedFluxPressureFvPatchScalarField::fixedFluxPressureFvPatchScalarField
 106 (
 107     const fixedFluxPressureFvPatchScalarField& wbppsf,
 108     const DimensionedField<scalar, volMesh>& iF
 109 )
 110 :
 111     fixedGradientFvPatchScalarField(wbppsf, iF),
 112     UName_(wbppsf.UName_),
 113     phiName_(wbppsf.phiName_),
 114     rhoName_(wbppsf.rhoName_),
 115     adjoint_(wbppsf.adjoint_)
 116 {}
 117
 118
 119 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 120
 121 void Foam::fixedFluxPressureFvPatchScalarField::updateCoeffs()
 122 {
 123     if (updated())
 124     {
 125         return;
 126     }
 127
 128     const fvPatchField<vector>& Up =
 129         patch().lookupPatchField<volVectorField, vector>(UName_);
 130
 131     const surfaceScalarField& phi =
 132         db().lookupObject<surfaceScalarField>(phiName_);
 133
 134     fvsPatchField<scalar> phip =
 135         patch().patchField<surfaceScalarField, scalar>(phi);
 136
 137     if (phi.dimensions() == dimDensity*dimVelocity*dimArea)
 138     {
 139         const fvPatchField<scalar>& rhop =
 140             patch().lookupPatchField<volScalarField, scalar>(rhoName_);
 141
 142         phip /= rhop;
 143     }
 144
 145     const fvPatchField<scalar>& rAp =
 146         patch().lookupPatchField<volScalarField, scalar>("(1|A("+UName_+"))");
 147
 148     if (adjoint_)
 149     {
 150         gradient() = ((patch().Sf() & Up) - phip)/patch().magSf()/rAp;
 151     }
 152     else
 153     {
 154         gradient() = (phip - (patch().Sf() & Up))/patch().magSf()/rAp;
 155     }
 156
 157     fixedGradientFvPatchScalarField::updateCoeffs();
 158 }
 159
 160
 161 void Foam::fixedFluxPressureFvPatchScalarField::write(Ostream& os) const
 162 {
 163     fvPatchScalarField::write(os);
 164     if (UName_ != "U")
 165     {
 166         os.writeKeyword("U") << UName_ << token::END_STATEMENT << nl;
 167     }
 168     if (phiName_ != "phi")
 169     {
 170         os.writeKeyword("phi") << phiName_ << token::END_STATEMENT << nl;
 171     }
 172     if (rhoName_ != "rho")
 173     {
 174         os.writeKeyword("rho") << rhoName_ << token::END_STATEMENT << nl;
 175     }
 176     os.writeKeyword("adjoint") << adjoint_ << token::END_STATEMENT << nl;
 177     gradient().writeEntry("gradient", os);
 178 }
 179
 180
 181 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 182
 183 namespace Foam
 184 {
 185     makePatchTypeField
 186     (
 187         fvPatchScalarField,
 188         fixedFluxPressureFvPatchScalarField
 189     );
 190 }
 191
 192 // ************************************************************************* //