src/turbulenceModels/incompressible/LES/Smagorinsky2/Smagorinsky2.H

   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 Class
  25     Foam::incompressible::LESModels::Smagorinsky2
  26
  27 Description
  28     The Isochoric Smagorinsky Model for incompressible flows
  29
  30     Algebraic eddy viscosity SGS model founded on the assumption that
  31     local equilibrium prevails, hence
  32     @verbatim
  33         B = 2/3*k*I - 2*nuSgs*dev(D) - 2*cD2*delta*(D.dev(D));
  34         Beff = 2/3*k*I - 2*nuEff*dev(D) - 2*cD2*delta*(D.dev(D));
  35
  36     where
  37
  38         D = symm(grad(U));
  39         k = cI*delta^2*||D||^2
  40         nuSgs = ck*sqrt(k)*delta
  41         nuEff = nuSgs + nu
  42     @endverbatim
  43
  44 SourceFiles
  45     Smagorinsky2.C
  46
  47 \*---------------------------------------------------------------------------*/
  48
  49 #ifndef Smagorinsky2_H
  50 #define Smagorinsky2_H
  51
  52 #include "Smagorinsky.H"
  53
  54 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  55
  56 namespace Foam
  57 {
  58 namespace incompressible
  59 {
  60 namespace LESModels
  61 {
  62
  63 /*---------------------------------------------------------------------------*\
  64                            Class Smagorinsky2 Declaration
  65 \*---------------------------------------------------------------------------*/
  66
  67 class Smagorinsky2
  68 :
  69     public Smagorinsky
  70 {
  71     // Private data
  72
  73         dimensionedScalar cD2_;
  74
  75
  76     // Private Member Functions
  77
  78         // Disallow default bitwise copy construct and assignment
  79         Smagorinsky2(const Smagorinsky2&);
  80         Smagorinsky2& operator=(const Smagorinsky2&);
  81
  82
  83 public:
  84
  85     //- Runtime type information
  86     TypeName("Smagorinsky2");
  87
  88
  89     // Constructors
  90
  91         //- Construct from components
  92         Smagorinsky2
  93         (
  94             const volVectorField& U,
  95             const surfaceScalarField& phi,
  96             transportModel& transport
  97         );
  98
  99
 100     //- Destructor
 101     virtual ~Smagorinsky2()
 102     {}
 103
 104
 105     // Member Functions
 106
 107         //- Return B.
 108         virtual tmp<volSymmTensorField> B() const;
 109
 110         //- Returns div(B).
 111         // This is the additional term due to the filtering of the NSE.
 112         virtual tmp<fvVectorMatrix> divDevBeff(volVectorField& U) const;
 113
 114         //- Read LESProperties dictionary
 115         virtual bool read();
 116 };
 117
 118
 119 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 120
 121 } // End namespace LESModels
 122 } // End namespace incompressible
 123 } // End namespace Foam
 124
 125 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 126
 127 #endif
 128
 129 // ************************************************************************* //