src/turbulenceModels/compressible/LES/lowReOneEqEddy/lowReOneEqEddy.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 "lowReOneEqEddy.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33 namespace compressible
  34 {
  35 namespace LESModels
  36 {
  37
  38 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  39
  40 defineTypeNameAndDebug(lowReOneEqEddy, 0);
  41 addToRunTimeSelectionTable(LESModel, lowReOneEqEddy, dictionary);
  42
  43 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  44
  45 void lowReOneEqEddy::updateSubGridScaleFields()
  46 {
  47     // High Re eddy viscosity
  48     muSgs_ = ck_*rho()*sqrt(k_)*delta();
  49
  50     // low Re no corrected eddy viscosity
  51     muSgs_ -= (mu()/beta_)*(scalar(1) - exp(-beta_*muSgs_/mu()));
  52     muSgs_.correctBoundaryConditions();
  53
  54     alphaSgs_ = muSgs_/Prt_;
  55     alphaSgs_.correctBoundaryConditions();
  56 }
  57
  58
  59 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  60
  61 lowReOneEqEddy::lowReOneEqEddy
  62 (
  63     const volScalarField& rho,
  64     const volVectorField& U,
  65     const surfaceScalarField& phi,
  66     const basicThermo& thermoPhysicalModel,
  67     const word& turbulenceModelName,
  68     const word& modelName
  69 )
  70 :
  71     LESModel(modelName, rho, U, phi, thermoPhysicalModel, turbulenceModelName),
  72     GenEddyVisc(rho, U, phi, thermoPhysicalModel),
  73
  74     k_
  75     (
  76         IOobject
  77         (
  78             "k",
  79             runTime_.timeName(),
  80             mesh_,
  81             IOobject::MUST_READ,
  82             IOobject::AUTO_WRITE
  83         ),
  84         mesh_
  85     ),
  86
  87     ck_
  88     (
  89         dimensioned<scalar>::lookupOrAddToDict
  90         (
  91             "ck",
  92             coeffDict_,
  93             0.07
  94         )
  95     ),
  96     beta_
  97     (
  98         dimensioned<scalar>::lookupOrAddToDict
  99         (
 100             "beta",
 101             coeffDict_,
 102             0.01
 103         )
 104     )
 105 {
 106     updateSubGridScaleFields();
 107
 108     printCoeffs();
 109 }
 110
 111
 112 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 113
 114 void lowReOneEqEddy::correct(const tmp<volTensorField>& tgradU)
 115 {
 116     const volTensorField& gradU = tgradU();
 117
 118     GenEddyVisc::correct(gradU);
 119
 120     volScalarField divU(fvc::div(phi()/fvc::interpolate(rho())));
 121     volScalarField G(2*muSgs_*(gradU && dev(symm(gradU))));
 122
 123     tmp<fvScalarMatrix> kEqn
 124     (
 125         fvm::ddt(rho(), k_)
 126       + fvm::div(phi(), k_)
 127       - fvm::laplacian(DkEff(), k_)
 128      ==
 129         G
 130       - fvm::SuSp(2.0/3.0*rho()*divU, k_)
 131       - fvm::Sp(ce_*rho()*sqrt(k_)/delta(), k_)
 132     );
 133
 134     kEqn().relax();
 135     kEqn().solve();
 136
 137     bound(k_, kMin_);
 138
 139     updateSubGridScaleFields();
 140 }
 141
 142
 143 bool lowReOneEqEddy::read()
 144 {
 145     if (GenEddyVisc::read())
 146     {
 147         ck_.readIfPresent(coeffDict());
 148         beta_.readIfPresent(coeffDict());
 149
 150         return true;
 151     }
 152     else
 153     {
 154         return false;
 155     }
 156 }
 157
 158
 159 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 160
 161 } // End namespace LESModels
 162 } // End namespace compressible
 163 } // End namespace Foam
 164
 165 // ************************************************************************* //