src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions/nutWallFunctions/nutUTabulatedWallFunction/nutUTabulatedWallFunctionFvPatchScalarField.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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 "nutUTabulatedWallFunctionFvPatchScalarField.H"
  27 #include "RASModel.H"
  28 #include "fvPatchFieldMapper.H"
  29 #include "volFields.H"
  30 #include "addToRunTimeSelectionTable.H"
  31
  32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36 namespace incompressible
  37 {
  38 namespace RASModels
  39 {
  40
  41 // * * * * * * * * * * * * Protected Member Functions  * * * * * * * * * * * //
  42
  43 tmp<scalarField> nutUTabulatedWallFunctionFvPatchScalarField::calcNut() const
  44 {
  45     const label patchI = patch().index();
  46
  47     const RASModel& rasModel = db().lookupObject<RASModel>("RASProperties");
  48     const scalarField& y = rasModel.y()[patchI];
  49     const fvPatchVectorField& Uw = rasModel.U().boundaryField()[patchI];
  50     const scalarField magUp(mag(Uw.patchInternalField() - Uw));
  51     const scalarField magGradU(mag(Uw.snGrad()));
  52     const tmp<volScalarField> tnu = rasModel.nu();
  53     const volScalarField& nu = tnu();
  54     const scalarField& nuw = nu.boundaryField()[patchI];
  55
  56     return
  57         max
  58         (
  59             scalar(0),
  60             sqr(magUp/(calcUPlus(magUp*y/nuw) + ROOTVSMALL))
  61            /(magGradU + ROOTVSMALL)
  62           - nuw
  63         );
  64 }
  65
  66
  67 tmp<scalarField> nutUTabulatedWallFunctionFvPatchScalarField::calcUPlus
  68 (
  69     const scalarField& Rey
  70 ) const
  71 {
  72     tmp<scalarField> tuPlus(new scalarField(patch().size(), 0.0));
  73     scalarField& uPlus = tuPlus();
  74
  75     forAll(uPlus, faceI)
  76     {
  77         uPlus[faceI] = uPlusTable_.interpolateLog10(Rey[faceI]);
  78     }
  79
  80     return tuPlus;
  81 }
  82
  83
  84 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  85
  86 nutUTabulatedWallFunctionFvPatchScalarField::
  87 nutUTabulatedWallFunctionFvPatchScalarField
  88 (
  89     const fvPatch& p,
  90     const DimensionedField<scalar, volMesh>& iF
  91 )
  92 :
  93     nutkWallFunctionFvPatchScalarField(p, iF),
  94     uPlusTableName_("undefined-uPlusTableName"),
  95     uPlusTable_
  96     (
  97         IOobject
  98         (
  99             uPlusTableName_,
 100             patch().boundaryMesh().mesh().time().constant(),
 101             patch().boundaryMesh().mesh(),
 102             IOobject::NO_READ,
 103             IOobject::NO_WRITE,
 104             false
 105         ),
 106         false
 107     )
 108 {}
 109
 110
 111 nutUTabulatedWallFunctionFvPatchScalarField::
 112 nutUTabulatedWallFunctionFvPatchScalarField
 113 (
 114     const nutUTabulatedWallFunctionFvPatchScalarField& ptf,
 115     const fvPatch& p,
 116     const DimensionedField<scalar, volMesh>& iF,
 117     const fvPatchFieldMapper& mapper
 118 )
 119 :
 120     nutkWallFunctionFvPatchScalarField(ptf, p, iF, mapper),
 121     uPlusTableName_(ptf.uPlusTableName_),
 122     uPlusTable_(ptf.uPlusTable_)
 123 {}
 124
 125
 126 nutUTabulatedWallFunctionFvPatchScalarField::
 127 nutUTabulatedWallFunctionFvPatchScalarField
 128 (
 129     const fvPatch& p,
 130     const DimensionedField<scalar, volMesh>& iF,
 131     const dictionary& dict
 132 )
 133 :
 134     nutkWallFunctionFvPatchScalarField(p, iF, dict),
 135     uPlusTableName_(dict.lookup("uPlusTable")),
 136     uPlusTable_
 137     (
 138         IOobject
 139         (
 140             uPlusTableName_,
 141             patch().boundaryMesh().mesh().time().constant(),
 142             patch().boundaryMesh().mesh(),
 143             IOobject::MUST_READ_IF_MODIFIED,
 144             IOobject::NO_WRITE,
 145             false
 146         ),
 147         true
 148     )
 149 {}
 150
 151
 152 nutUTabulatedWallFunctionFvPatchScalarField::
 153 nutUTabulatedWallFunctionFvPatchScalarField
 154 (
 155     const nutUTabulatedWallFunctionFvPatchScalarField& wfpsf
 156 )
 157 :
 158     nutkWallFunctionFvPatchScalarField(wfpsf),
 159     uPlusTableName_(wfpsf.uPlusTableName_),
 160     uPlusTable_(wfpsf.uPlusTable_)
 161 {}
 162
 163
 164 nutUTabulatedWallFunctionFvPatchScalarField::
 165 nutUTabulatedWallFunctionFvPatchScalarField
 166 (
 167     const nutUTabulatedWallFunctionFvPatchScalarField& wfpsf,
 168     const DimensionedField<scalar, volMesh>& iF
 169 )
 170 :
 171     nutkWallFunctionFvPatchScalarField(wfpsf, iF),
 172     uPlusTableName_(wfpsf.uPlusTableName_),
 173     uPlusTable_(wfpsf.uPlusTable_)
 174 {}
 175
 176
 177 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 178
 179 tmp<scalarField> nutUTabulatedWallFunctionFvPatchScalarField::yPlus() const
 180 {
 181     const label patchI = patch().index();
 182
 183     const RASModel& rasModel = db().lookupObject<RASModel>("RASProperties");
 184     const scalarField& y = rasModel.y()[patchI];
 185     const fvPatchVectorField& Uw = rasModel.U().boundaryField()[patchI];
 186     const scalarField magUp(mag(Uw.patchInternalField() - Uw));
 187     const tmp<volScalarField> tnu = rasModel.nu();
 188     const volScalarField& nu = tnu();
 189     const scalarField& nuw = nu.boundaryField()[patchI];
 190     const scalarField Rey(magUp*y/nuw);
 191
 192     return Rey/(calcUPlus(Rey) + ROOTVSMALL);
 193 }
 194
 195
 196 void nutUTabulatedWallFunctionFvPatchScalarField::write(Ostream& os) const
 197 {
 198     fvPatchField<scalar>::write(os);
 199     os.writeKeyword("uPlusTable") << uPlusTableName_
 200         << token::END_STATEMENT << nl;
 201     writeEntry("value", os);
 202 }
 203
 204
 205 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 206
 207 makePatchTypeField
 208 (
 209     fvPatchScalarField,
 210     nutUTabulatedWallFunctionFvPatchScalarField
 211 );
 212
 213 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 214
 215 } // End namespace RASModels
 216 } // End namespace incompressible
 217 } // End namespace Foam
 218
 219 // ************************************************************************* //