src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.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 "Smagorinsky2.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33 namespace incompressible
  34 {
  35 namespace LESModels
  36 {
  37
  38 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  39
  40 defineTypeNameAndDebug(Smagorinsky2, 0);
  41 addToRunTimeSelectionTable(LESModel, Smagorinsky2, dictionary);
  42
  43 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  44
  45 Smagorinsky2::Smagorinsky2
  46 (
  47     const volVectorField& U,
  48     const surfaceScalarField& phi,
  49     transportModel& transport,
  50     const word& turbulenceModelName,
  51     const word& modelName
  52 )
  53 :
  54     LESModel(modelName, U, phi, transport, turbulenceModelName),
  55     Smagorinsky(U, phi, transport),
  56
  57     cD2_
  58     (
  59         dimensioned<scalar>::lookupOrAddToDict
  60         (
  61             "cD2",
  62             coeffDict_,
  63             0.02
  64         )
  65     )
  66 {
  67     printCoeffs();
  68 }
  69
  70
  71 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  72
  73 // Evaluate B (from the definition of an eddy viscosity model) and
  74 // return it.
  75
  76 tmp<volSymmTensorField> Smagorinsky2::B() const
  77 {
  78     volSymmTensorField D(dev(symm(fvc::grad(U()))));
  79
  80     return (((2.0/3.0)*I)*k() - 2.0*nuSgs_*D - (2.0*cD2_)*delta()*(D&D));
  81 }
  82
  83
  84 tmp<fvVectorMatrix> Smagorinsky2::divDevBeff
  85 (
  86     volVectorField& U
  87 ) const
  88 {
  89     volTensorField gradU(fvc::grad(U));
  90
  91     volSymmTensorField aniNuEff
  92     (
  93         "aniNuEff",
  94         I*nuEff() + cD2_*delta()*symm(gradU)
  95     );
  96
  97     return
  98     (
  99       - fvm::laplacian(aniNuEff, U) - fvc::div(nuEff()*dev(T(fvc::grad(U))))
 100     );
 101 }
 102
 103
 104 bool Smagorinsky2::read()
 105 {
 106     if (Smagorinsky::read())
 107     {
 108         cD2.readIfPresent(coeffDict());
 109
 110         return true;
 111     }
 112     else
 113     {
 114         return false;
 115     }
 116 }
 117
 118
 119 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 120
 121 } // End namespace LESModels
 122 } // End namespace incompressible
 123 } // End namespace Foam
 124
 125 // ************************************************************************* //