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