src/turbulenceModels/incompressible/LES/dynOneEqEddy/dynOneEqEddy.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::dynOneEqEddy
  27
  28 Description
  29     One Equation Eddy Viscosity Model for incompressible flows.
  30
  31     Eddy viscosity SGS model using a modeled balance equation to simulate
  32     the behaviour of k.
  33
  34     Thus
  35     @verbatim
  36         d/dt(k) + div(U*k) - div(nuSgs*grad(k))
  37         =
  38         -B*L - ce*k^3/2/delta
  39
  40     and
  41
  42         B = 2/3*k*I - 2*nuSgs*dev(D)
  43         Beff = 2/3*k*I - 2*nuEff*dev(D)
  44
  45     where
  46
  47         D = symm(grad(U));
  48         nuSgs = ck*sqrt(k)*delta
  49         nuEff = nuSgs + nu
  50     @endverbatim
  51
  52 SourceFiles
  53     dynOneEqEddy.C
  54
  55 \*---------------------------------------------------------------------------*/
  56
  57 #ifndef dynOneEqEddy_H
  58 #define dynOneEqEddy_H
  59
  60 #include "GenEddyVisc.H"
  61 #include "LESfilter.H"
  62
  63 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  64
  65 namespace Foam
  66 {
  67 namespace incompressible
  68 {
  69 namespace LESModels
  70 {
  71
  72 /*---------------------------------------------------------------------------*\
  73                            Class dynOneEqEddy Declaration
  74 \*---------------------------------------------------------------------------*/
  75
  76 class dynOneEqEddy
  77 :
  78     public GenEddyVisc
  79 {
  80     // Private data
  81
  82         volScalarField k_;
  83
  84         autoPtr<LESfilter> filterPtr_;
  85         LESfilter& filter_;
  86
  87
  88     // Private Member Functions
  89
  90         //- Update sub-grid scale fields
  91         void updateSubGridScaleFields(const volSymmTensorField& D);
  92
  93         //- Calculate ck, ce by filtering the velocity field U.
  94         dimensionedScalar ck(const volSymmTensorField& D) const;
  95         dimensionedScalar ce(const volSymmTensorField& D) const;
  96
  97         // Disallow default bitwise copy construct and assignment
  98         dynOneEqEddy(const dynOneEqEddy&);
  99         dynOneEqEddy& operator=(const dynOneEqEddy&);
 100
 101
 102 public:
 103
 104     //- Runtime type information
 105     TypeName("dynOneEqEddy");
 106
 107     // Constructors
 108
 109         //- Construct from components
 110         dynOneEqEddy
 111         (
 112             const volVectorField& U,
 113             const surfaceScalarField& phi,
 114             transportModel& transport
 115         );
 116
 117
 118     //- Destructor
 119     virtual ~dynOneEqEddy()
 120     {}
 121
 122
 123     // Member Functions
 124
 125         //- Return SGS kinetic energy
 126         virtual tmp<volScalarField> k() const
 127         {
 128             return k_;
 129         }
 130
 131         //- Return the effective diffusivity for k
 132         tmp<volScalarField> DkEff() const
 133         {
 134             return tmp<volScalarField>
 135             (
 136                 new volScalarField("DkEff", nuSgs_ + nu())
 137             );
 138         }
 139
 140         //- Correct Eddy-Viscosity and related properties
 141         virtual void correct(const tmp<volTensorField>& gradU);
 142
 143         //- Read LESProperties dictionary
 144         virtual bool read();
 145 };
 146
 147
 148 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 149
 150 } // End namespace LESModels
 151 } // End namespace incompressible
 152 } // End namespace Foam
 153
 154 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 155
 156 #endif
 157
 158 // ************************************************************************* //