src/finiteVolume/fields/fvPatchFields/derived/timeVaryingUniformTotalPressure/timeVaryingUniformTotalPressureFvPatchScalarField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 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 "timeVaryingUniformTotalPressureFvPatchScalarField.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "fvPatchFieldMapper.H"
  30 #include "volFields.H"
  31 #include "surfaceFields.H"
  32
  33 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  34
  35 Foam::timeVaryingUniformTotalPressureFvPatchScalarField::
  36 timeVaryingUniformTotalPressureFvPatchScalarField
  37 (
  38     const fvPatch& p,
  39     const DimensionedField<scalar, volMesh>& iF
  40 )
  41 :
  42     fixedValueFvPatchScalarField(p, iF),
  43     UName_("undefined"),
  44     phiName_("undefined"),
  45     rhoName_("undefined"),
  46     psiName_("undefined"),
  47     gamma_(0.0),
  48     p0_(0.0),
  49     totalPressureTimeSeries_()
  50 {}
  51
  52
  53 Foam::timeVaryingUniformTotalPressureFvPatchScalarField::
  54 timeVaryingUniformTotalPressureFvPatchScalarField
  55 (
  56     const fvPatch& p,
  57     const DimensionedField<scalar, volMesh>& iF,
  58     const dictionary& dict
  59 )
  60 :
  61     fixedValueFvPatchScalarField(p, iF),
  62     UName_(dict.lookup("U")),
  63     phiName_(dict.lookup("phi")),
  64     rhoName_(dict.lookup("rho")),
  65     psiName_(dict.lookup("psi")),
  66     gamma_(readScalar(dict.lookup("gamma"))),
  67     p0_(readScalar(dict.lookup("p0"))),
  68     totalPressureTimeSeries_(dict)
  69 {
  70     if (dict.found("value"))
  71     {
  72         fvPatchField<scalar>::operator=
  73         (
  74             scalarField("value", dict, p.size())
  75         );
  76     }
  77     else
  78     {
  79         fvPatchField<scalar>::operator=(p0_);
  80     }
  81 }
  82
  83
  84 Foam::timeVaryingUniformTotalPressureFvPatchScalarField::
  85 timeVaryingUniformTotalPressureFvPatchScalarField
  86 (
  87     const timeVaryingUniformTotalPressureFvPatchScalarField& ptf,
  88     const fvPatch& p,
  89     const DimensionedField<scalar, volMesh>& iF,
  90     const fvPatchFieldMapper& mapper
  91 )
  92 :
  93     fixedValueFvPatchScalarField(ptf, p, iF, mapper),
  94     UName_(ptf.UName_),
  95     phiName_(ptf.phiName_),
  96     rhoName_(ptf.rhoName_),
  97     psiName_(ptf.psiName_),
  98     gamma_(ptf.gamma_),
  99     p0_(ptf.p0_),
 100     totalPressureTimeSeries_(ptf.totalPressureTimeSeries_)
 101 {}
 102
 103
 104 Foam::timeVaryingUniformTotalPressureFvPatchScalarField::
 105 timeVaryingUniformTotalPressureFvPatchScalarField
 106 (
 107     const timeVaryingUniformTotalPressureFvPatchScalarField& tppsf
 108 )
 109 :
 110     fixedValueFvPatchScalarField(tppsf),
 111     UName_(tppsf.UName_),
 112     phiName_(tppsf.phiName_),
 113     rhoName_(tppsf.rhoName_),
 114     psiName_(tppsf.psiName_),
 115     gamma_(tppsf.gamma_),
 116     p0_(tppsf.p0_),
 117     totalPressureTimeSeries_(tppsf.totalPressureTimeSeries_)
 118 {}
 119
 120
 121 Foam::timeVaryingUniformTotalPressureFvPatchScalarField::
 122 timeVaryingUniformTotalPressureFvPatchScalarField
 123 (
 124     const timeVaryingUniformTotalPressureFvPatchScalarField& tppsf,
 125     const DimensionedField<scalar, volMesh>& iF
 126 )
 127 :
 128     fixedValueFvPatchScalarField(tppsf, iF),
 129     UName_(tppsf.UName_),
 130     phiName_(tppsf.phiName_),
 131     rhoName_(tppsf.rhoName_),
 132     psiName_(tppsf.psiName_),
 133     gamma_(tppsf.gamma_),
 134     p0_(tppsf.p0_),
 135     totalPressureTimeSeries_(tppsf.totalPressureTimeSeries_)
 136 {}
 137
 138
 139 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 140
 141 void Foam::timeVaryingUniformTotalPressureFvPatchScalarField::updateCoeffs
 142 (
 143     const vectorField& Up
 144 )
 145 {
 146     if (updated())
 147     {
 148         return;
 149     }
 150
 151     p0_ = totalPressureTimeSeries_(this->db().time().timeOutputValue());
 152
 153     const fvsPatchField<scalar>& phip =
 154         patch().lookupPatchField<surfaceScalarField, scalar>(phiName_);
 155
 156     if (psiName_ == "none" && rhoName_ == "none")
 157     {
 158         operator==(p0_ - 0.5*(1.0 - pos(phip))*magSqr(Up));
 159     }
 160     else if (rhoName_ == "none")
 161     {
 162         const fvPatchField<scalar>& psip =
 163             patch().lookupPatchField<volScalarField, scalar>(psiName_);
 164
 165         if (gamma_ > 1.0)
 166         {
 167             scalar gM1ByG = (gamma_ - 1.0)/gamma_;
 168
 169             operator==
 170             (
 171                 p0_
 172                /pow
 173                 (
 174                     (1.0 + 0.5*psip*gM1ByG*(1.0 - pos(phip))*magSqr(Up)),
 175                     1.0/gM1ByG
 176                 )
 177             );
 178         }
 179         else
 180         {
 181             operator==(p0_/(1.0 + 0.5*psip*(1.0 - pos(phip))*magSqr(Up)));
 182         }
 183     }
 184     else if (psiName_ == "none")
 185     {
 186         const fvPatchField<scalar>& rho =
 187             patch().lookupPatchField<volScalarField, scalar>(rhoName_);
 188
 189         operator==(p0_ - 0.5*rho*(1.0 - pos(phip))*magSqr(Up));
 190     }
 191     else
 192     {
 193         FatalErrorIn
 194         (
 195             "timeVaryingUniformTotalPressureFvPatchScalarField::updateCoeffs()"
 196         )   << " rho or psi set inconsitently, rho = " << rhoName_
 197             << ", psi = " << psiName_ << '.' << nl
 198             << "    Set either rho or psi or neither depending on the "
 199                "definition of total pressure." << nl
 200             << "    Set the unused variables to 'none'."
 201             << "\n    on patch " << this->patch().name()
 202             << " of field " << this->dimensionedInternalField().name()
 203             << " in file " << this->dimensionedInternalField().objectPath()
 204             << exit(FatalError);
 205     }
 206
 207     fixedValueFvPatchScalarField::updateCoeffs();
 208 }
 209
 210
 211 void Foam::timeVaryingUniformTotalPressureFvPatchScalarField::updateCoeffs()
 212 {
 213     updateCoeffs(patch().lookupPatchField<volVectorField, vector>(UName_));
 214 }
 215
 216
 217 void Foam::timeVaryingUniformTotalPressureFvPatchScalarField::
 218 write(Ostream& os) const
 219 {
 220     fvPatchScalarField::write(os);
 221     os.writeKeyword("U") << UName_ << token::END_STATEMENT << nl;
 222     os.writeKeyword("phi") << phiName_ << token::END_STATEMENT << nl;
 223     os.writeKeyword("rho") << rhoName_ << token::END_STATEMENT << nl;
 224     os.writeKeyword("psi") << psiName_ << token::END_STATEMENT << nl;
 225     os.writeKeyword("gamma") << gamma_ << token::END_STATEMENT << nl;
 226     os.writeKeyword("p0") << p0_ << token::END_STATEMENT << endl;
 227     totalPressureTimeSeries_.write(os);
 228     writeEntry("value", os);
 229 }
 230
 231
 232 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 233
 234 namespace Foam
 235 {
 236     makePatchTypeField
 237     (
 238         fvPatchScalarField,
 239         timeVaryingUniformTotalPressureFvPatchScalarField
 240     );
 241 }
 242
 243 // ************************************************************************* //