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) 1991-2010 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
  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 )
  51 :
  52     LESModel(typeName, U, phi, transport),
  53     Smagorinsky(U, phi, transport),
  54
  55     cD2_
  56     (
  57         dimensioned<scalar>::lookupOrAddToDict
  58         (
  59             "cD2",
  60             coeffDict_,
  61             0.02
  62         )
  63     )
  64 {
  65     printCoeffs();
  66 }
  67
  68
  69 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  70
  71 // Evaluate B (from the definition of an eddy viscosity model) and
  72 // return it.
  73
  74 tmp<volSymmTensorField> Smagorinsky2::B() const
  75 {
  76     volSymmTensorField D = dev(symm(fvc::grad(U())));
  77
  78     return (((2.0/3.0)*I)*k() - 2.0*nuSgs_*D - (2.0*cD2_)*delta()*(D&D));
  79 }
  80
  81
  82 tmp<fvVectorMatrix> Smagorinsky2::divDevBeff
  83 (
  84     volVectorField& U
  85 ) const
  86 {
  87     volTensorField gradU = fvc::grad(U);
  88
  89     volSymmTensorField aniNuEff
  90     (
  91         "aniNuEff",
  92         I*nuEff() + cD2_*delta()*symm(gradU)
  93     );
  94
  95     return
  96     (
  97       - fvm::laplacian(aniNuEff, U) - fvc::div(nuEff()*dev(fvc::grad(U)().T()))
  98     );
  99 }
 100
 101
 102 bool Smagorinsky2::read()
 103 {
 104     if (Smagorinsky::read())
 105     {
 106         cD2.readIfPresent(coeffDict());
 107
 108         return true;
 109     }
 110     else
 111     {
 112         return false;
 113     }
 114 }
 115
 116
 117 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 118
 119 } // End namespace LESModels
 120 } // End namespace incompressible
 121 } // End namespace Foam
 122
 123 // ************************************************************************* //