src/turbulenceModels/compressible/LES/derivedFvPatchFields/wallFunctions/muSgsWallFunctions/muSgsWallFunction/muSgsWallFunctionFvPatchScalarField.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 "muSgsWallFunctionFvPatchScalarField.H"
  28 #include "LESModel.H"
  29 #include "fvPatchFieldMapper.H"
  30 #include "volFields.H"
  31 #include "addToRunTimeSelectionTable.H"
  32
  33 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  34
  35 namespace Foam
  36 {
  37 namespace compressible
  38 {
  39 namespace LESModels
  40 {
  41
  42 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  43
  44 muSgsWallFunctionFvPatchScalarField::muSgsWallFunctionFvPatchScalarField
  45 (
  46     const fvPatch& p,
  47     const DimensionedField<scalar, volMesh>& iF
  48 )
  49 :
  50     fixedValueFvPatchScalarField(p, iF),
  51     UName_("U"),
  52     rhoName_("rho"),
  53     muName_("mu"),
  54     kappa_(0.41),
  55     E_(9.8)
  56 {}
  57
  58
  59 muSgsWallFunctionFvPatchScalarField::muSgsWallFunctionFvPatchScalarField
  60 (
  61     const muSgsWallFunctionFvPatchScalarField& ptf,
  62     const fvPatch& p,
  63     const DimensionedField<scalar, volMesh>& iF,
  64     const fvPatchFieldMapper& mapper
  65 )
  66 :
  67     fixedValueFvPatchScalarField(ptf, p, iF, mapper),
  68     UName_(ptf.UName_),
  69     rhoName_(ptf.rhoName_),
  70     muName_(ptf.muName_),
  71     kappa_(ptf.kappa_),
  72     E_(ptf.E_)
  73 {}
  74
  75
  76 muSgsWallFunctionFvPatchScalarField::muSgsWallFunctionFvPatchScalarField
  77 (
  78     const fvPatch& p,
  79     const DimensionedField<scalar, volMesh>& iF,
  80     const dictionary& dict
  81 )
  82 :
  83     fixedValueFvPatchScalarField(p, iF, dict),
  84     UName_(dict.lookupOrDefault<word>("U", "U")),
  85     rhoName_(dict.lookupOrDefault<word>("rho", "rho")),
  86     muName_(dict.lookupOrDefault<word>("mu", "mu")),
  87     kappa_(dict.lookupOrDefault<scalar>("kappa", 0.41)),
  88     E_(dict.lookupOrDefault<scalar>("E", 9.8))
  89 {}
  90
  91
  92 muSgsWallFunctionFvPatchScalarField::muSgsWallFunctionFvPatchScalarField
  93 (
  94     const muSgsWallFunctionFvPatchScalarField& mwfpsf
  95 )
  96 :
  97     fixedValueFvPatchScalarField(mwfpsf),
  98     UName_(mwfpsf.UName_),
  99     rhoName_(mwfpsf.rhoName_),
 100     muName_(mwfpsf.muName_),
 101     kappa_(mwfpsf.kappa_),
 102     E_(mwfpsf.E_)
 103 {}
 104
 105
 106 muSgsWallFunctionFvPatchScalarField::muSgsWallFunctionFvPatchScalarField
 107 (
 108     const muSgsWallFunctionFvPatchScalarField& mwfpsf,
 109     const DimensionedField<scalar, volMesh>& iF
 110 )
 111 :
 112     fixedValueFvPatchScalarField(mwfpsf, iF),
 113     UName_(mwfpsf.UName_),
 114     rhoName_(mwfpsf.rhoName_),
 115     muName_(mwfpsf.muName_),
 116     kappa_(mwfpsf.kappa_),
 117     E_(mwfpsf.E_)
 118 {}
 119
 120
 121 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 122
 123 void muSgsWallFunctionFvPatchScalarField::evaluate
 124 (
 125     const Pstream::commsTypes
 126 )
 127 {
 128     const scalarField& ry = patch().deltaCoeffs();
 129
 130     const fvPatchVectorField& U =
 131         patch().lookupPatchField<volVectorField, vector>(UName_);
 132
 133     scalarField magUp = mag(U.patchInternalField() - U);
 134
 135     const scalarField& muw =
 136         patch().lookupPatchField<volScalarField, scalar>(muName_);
 137
 138     const scalarField& rhow =
 139         patch().lookupPatchField<volScalarField, scalar>(rhoName_);
 140
 141     scalarField& muSgsw = *this;
 142
 143     scalarField magFaceGradU = mag(U.snGrad());
 144
 145     forAll(muSgsw, facei)
 146     {
 147         scalar magUpara = magUp[facei];
 148
 149         scalar utau =
 150             sqrt((muSgsw[facei] + muw[facei])*magFaceGradU[facei]/rhow[facei]);
 151
 152         if (utau > 0)
 153         {
 154             int iter = 0;
 155             scalar err = GREAT;
 156
 157             do
 158             {
 159                 scalar kUu = kappa_*magUpara/utau;
 160                 scalar fkUu = exp(kUu) - 1 - kUu*(1 + 0.5*kUu);
 161
 162                 scalar f =
 163                     - utau/(ry[facei]*muw[facei]/rhow[facei])
 164                     + magUpara/utau
 165                     + 1/E_*(fkUu - 1.0/6.0*kUu*sqr(kUu));
 166
 167                 scalar df =
 168                     - 1.0/(ry[facei]*muw[facei]/rhow[facei])
 169                     - magUpara/sqr(utau)
 170                     - 1/E_*kUu*fkUu/utau;
 171
 172                 scalar utauNew = utau - f/df;
 173                 err = mag((utau - utauNew)/utau);
 174                 utau = utauNew;
 175
 176             } while (utau > VSMALL && err > 0.01 && ++iter < 10);
 177
 178             muSgsw[facei] =
 179                 max
 180                 (
 181                     rhow[facei]*sqr(utau)/magFaceGradU[facei] - muw[facei],
 182                     0.0
 183                 );
 184         }
 185         else
 186         {
 187             muSgsw[facei] = 0;
 188         }
 189     }
 190 }
 191
 192
 193 void muSgsWallFunctionFvPatchScalarField::write(Ostream& os) const
 194 {
 195     fvPatchField<scalar>::write(os);
 196     writeEntryIfDifferent<word>(os, "U", "U", UName_);
 197     writeEntryIfDifferent<word>(os, "rho", "rho", rhoName_);
 198     writeEntryIfDifferent<word>(os, "mu", "mu", muName_);
 199     os.writeKeyword("kappa") << kappa_ << token::END_STATEMENT << nl;
 200     os.writeKeyword("E") << E_ << token::END_STATEMENT << nl;
 201     writeEntry("value", os);
 202 }
 203
 204
 205 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 206
 207 makePatchTypeField
 208 (
 209     fvPatchScalarField,
 210     muSgsWallFunctionFvPatchScalarField
 211 );
 212
 213 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 214
 215 } // End namespace LESModels
 216 } // End namespace compressible
 217 } // End namespace Foam
 218
 219 // ************************************************************************* //