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