src/turbulenceModels/LES/compressible/Smagorinsky/Smagorinsky.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 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::compressible::LESModels::Smagorinsky
  27
  28 Description
  29     The choric Smagorinsky Model for compressible flows.
  30
  31     Algebraic eddy viscosity SGS model founded on the assumption that
  32     local equilibrium prevails.
  33     Thus,
  34     @verbatim
  35
  36         B = 2/3*k*I - 2*nuSgs*dev(D)
  37
  38     where
  39
  40         D = symm(grad(U));
  41         k from rho*D:B + ce*rho*k^3/2/delta = 0
  42         muSgs = ck*rho*sqrt(k)*delta
  43     @endverbatim
  44
  45 SourceFiles
  46     Smagorinsky.C
  47
  48 \*---------------------------------------------------------------------------*/
  49
  50 #ifndef compressibleSmagorinsky_H
  51 #define compressibleSmagorinsky_H
  52
  53 #include "GenEddyVisc.H"
  54
  55 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  56
  57 namespace Foam
  58 {
  59 namespace compressible
  60 {
  61 namespace LESModels
  62 {
  63
  64 /*---------------------------------------------------------------------------*\
  65                            Class Smagorinsky Declaration
  66 \*---------------------------------------------------------------------------*/
  67
  68 class Smagorinsky
  69 :
  70     public GenEddyVisc
  71 {
  72     // Private data
  73
  74         dimensionedScalar ck_;
  75
  76
  77     // Private Member Functions
  78
  79         // Disallow default bitwise copy construct and assignment
  80         Smagorinsky(const Smagorinsky&);
  81         Smagorinsky& operator=(const Smagorinsky&);
  82
  83
  84 public:
  85
  86     //- Runtime type information
  87     TypeName("Smagorinsky");
  88
  89     // Constructors
  90
  91         //- Construct from components
  92         Smagorinsky
  93         (
  94             const volScalarField& rho,
  95             const volVectorField& U,
  96             const surfaceScalarField& phi,
  97             const basicThermo& thermoPhysicalModel
  98         );
  99
 100
 101     // Destructor
 102
 103         ~Smagorinsky()
 104         {}
 105
 106
 107     // Member Functions
 108
 109         //- Correct Eddy-Viscosity and related properties
 110         void correct(const tmp<volTensorField>& gradU);
 111
 112         //- Read turbulenceProperties dictionary
 113         bool read();
 114 };
 115
 116
 117 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 118
 119 } // End namespace LESModels
 120 } // End namespace compressible
 121 } // End namespace Foam
 122
 123 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 124
 125 #endif
 126
 127 // ************************************************************************* //