applications/solvers/compressible/rhoCentralFoam/BCs/T/smoluchowskiJumpTFvPatchScalarField.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 "smoluchowskiJumpTFvPatchScalarField.H"
  28 #include <OpenFOAM/addToRunTimeSelectionTable.H>
  29 #include <finiteVolume/fvPatchFieldMapper.H>
  30 #include <finiteVolume/volFields.H>
  31 #include <OpenFOAM/mathematicalConstants.H>
  32
  33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37
  38 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  39
  40 smoluchowskiJumpTFvPatchScalarField::smoluchowskiJumpTFvPatchScalarField
  41 (
  42     const fvPatch& p,
  43     const DimensionedField<scalar, volMesh>& iF
  44 )
  45 :
  46     mixedFvPatchScalarField(p, iF),
  47     accommodationCoeff_(1.0),
  48     Twall_(p.size(), 0.0),
  49     gamma_(1.4)
  50 {
  51     refValue() = 0.0;
  52     refGrad() = 0.0;
  53     valueFraction() = 0.0;
  54 }
  55
  56 smoluchowskiJumpTFvPatchScalarField::smoluchowskiJumpTFvPatchScalarField
  57 (
  58     const smoluchowskiJumpTFvPatchScalarField& ptf,
  59     const fvPatch& p,
  60     const DimensionedField<scalar, volMesh>& iF,
  61     const fvPatchFieldMapper& mapper
  62 )
  63 :
  64     mixedFvPatchScalarField(ptf, p, iF, mapper),
  65     accommodationCoeff_(ptf.accommodationCoeff_),
  66     Twall_(ptf.Twall_),
  67     gamma_(ptf.gamma_)
  68 {}
  69
  70
  71 smoluchowskiJumpTFvPatchScalarField::smoluchowskiJumpTFvPatchScalarField
  72 (
  73     const fvPatch& p,
  74     const DimensionedField<scalar, volMesh>& iF,
  75     const dictionary& dict
  76 )
  77 :
  78     mixedFvPatchScalarField(p, iF),
  79     accommodationCoeff_(readScalar(dict.lookup("accommodationCoeff"))),
  80     Twall_("Twall", dict, p.size()),
  81     gamma_(dict.lookupOrDefault<scalar>("gamma", 1.4))
  82 {
  83     if
  84     (
  85         mag(accommodationCoeff_) < SMALL
  86      || mag(accommodationCoeff_) > 2.0
  87     )
  88     {
  89         FatalIOErrorIn
  90         (
  91             "smoluchowskiJumpTFvPatchScalarField::"
  92             "smoluchowskiJumpTFvPatchScalarField"
  93             "("
  94             "    const fvPatch&,"
  95             "    const DimensionedField<scalar, volMesh>&,"
  96             "    const dictionary&"
  97             ")",
  98             dict
  99         )   << "unphysical accommodationCoeff specified"
 100             << "(0 < accommodationCoeff <= 1)" << endl
 101             << exit(FatalError);
 102     }
 103
 104     if (dict.found("value"))
 105     {
 106         fvPatchField<scalar>::operator=
 107         (
 108             scalarField("value", dict, p.size())
 109         );
 110     }
 111     else
 112     {
 113         fvPatchField<scalar>::operator=(patchInternalField());
 114     }
 115
 116     refValue() = *this;
 117     refGrad() = 0.0;
 118     valueFraction() = 0.0;
 119 }
 120
 121
 122 smoluchowskiJumpTFvPatchScalarField::smoluchowskiJumpTFvPatchScalarField
 123 (
 124     const smoluchowskiJumpTFvPatchScalarField& ptpsf,
 125     const DimensionedField<scalar, volMesh>& iF
 126 )
 127 :
 128     mixedFvPatchScalarField(ptpsf, iF),
 129     accommodationCoeff_(ptpsf.accommodationCoeff_),
 130     Twall_(ptpsf.Twall_),
 131     gamma_(ptpsf.gamma_)
 132 {}
 133
 134
 135 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 136
 137 // Map from self
 138 void smoluchowskiJumpTFvPatchScalarField::autoMap
 139 (
 140     const fvPatchFieldMapper& m
 141 )
 142 {
 143     mixedFvPatchScalarField::autoMap(m);
 144 }
 145
 146
 147 // Reverse-map the given fvPatchField onto this fvPatchField
 148 void smoluchowskiJumpTFvPatchScalarField::rmap
 149 (
 150     const fvPatchField<scalar>& ptf,
 151     const labelList& addr
 152 )
 153 {
 154     mixedFvPatchField<scalar>::rmap(ptf, addr);
 155 }
 156
 157
 158 // Update the coefficients associated with the patch field
 159 void smoluchowskiJumpTFvPatchScalarField::updateCoeffs()
 160 {
 161     if (updated())
 162     {
 163         return;
 164     }
 165
 166     const fvPatchScalarField& pmu =
 167         patch().lookupPatchField<volScalarField, scalar>("mu");
 168     const fvPatchScalarField& prho =
 169         patch().lookupPatchField<volScalarField, scalar>("rho");
 170     const fvPatchField<scalar>& ppsi =
 171         patch().lookupPatchField<volScalarField, scalar>("psi");
 172     const fvPatchVectorField& pU =
 173         patch().lookupPatchField<volVectorField, vector>("U");
 174
 175     // Prandtl number reading consistent with rhoCentralFoam
 176     const dictionary& thermophysicalProperties =
 177         db().lookupObject<IOdictionary>("thermophysicalProperties");
 178     dimensionedScalar Pr = dimensionedScalar("Pr", dimless, 1.0);
 179     if (thermophysicalProperties.found("Pr"))
 180     {
 181         Pr = thermophysicalProperties.lookup("Pr");
 182     }
 183
 184     Field<scalar> C2 = pmu/prho
 185         *sqrt(ppsi*mathematicalConstant::pi/2.0)
 186         *2.0*gamma_/Pr.value()/(gamma_ + 1.0)
 187         *(2.0 - accommodationCoeff_)/accommodationCoeff_;
 188
 189     Field<scalar> aCoeff = prho.snGrad() - prho/C2;
 190     Field<scalar> KEbyRho = 0.5*magSqr(pU);
 191
 192     valueFraction() = (1.0/(1.0 + patch().deltaCoeffs()*C2));
 193     refValue() = Twall_;
 194     refGrad() = 0.0;
 195
 196     mixedFvPatchScalarField::updateCoeffs();
 197 }
 198
 199
 200 // Write
 201 void smoluchowskiJumpTFvPatchScalarField::write(Ostream& os) const
 202 {
 203     fvPatchScalarField::write(os);
 204     os.writeKeyword("accommodationCoeff")
 205         << accommodationCoeff_ << token::END_STATEMENT << nl;
 206     Twall_.writeEntry("Twall", os);
 207     os.writeKeyword("gamma")
 208         << gamma_ << token::END_STATEMENT << nl;
 209     writeEntry("value", os);
 210 }
 211
 212
 213 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 214
 215 makePatchTypeField(fvPatchScalarField, smoluchowskiJumpTFvPatchScalarField);
 216
 217 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 218
 219 } // End namespace Foam
 220
 221 // ************************ vim: set sw=4 sts=4 et: ************************ //