FreeFOAM The Cross-Platform CFD Toolkit
spectEddyVisc.H
Go to the documentation of this file.
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::spectEddyVisc
26 
27 Description
28  The Isochoric spectral Eddy Viscosity Model for incompressible flows.
29 
30  Algebraic eddy viscosity SGS model founded on the assumption that
31  local equilibrium prevail and that viscous effects may be of importance.
32 
33  Thus,
34  @verbatim
35  B = 2/3*k*I - 2*nuSgs*dev(D)
36  Beff = 2/3*k*I - 2*nuEff*dev(D)
37 
38  where
39 
40  k = cK1*delta^(2/3)*eps^(2/3)*exp(-cK2*delta^(4/3)*nu*eps^(-1/3)) -
41  ck3*sqrt(eps*nu)*erfc(ck4*delta^(-2/3)*sqrt(nu)*eps^(-1/6)))
42 
43  nuEff = nu/(1 - exp(-cB*pow((nu/nuEff), 1.0/3.0)*pow(Re, -2.0/3.0)))
44  nuSgs = nuEff - nu
45 
46  Re = delta^2*mag(D)/nu
47  @endverbatim
48 
49 SourceFiles
50  spectEddyVisc.C
51 
52 \*---------------------------------------------------------------------------*/
53 
54 #ifndef spectEddyVisc_H
55 #define spectEddyVisc_H
56 
58 
59 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
60 
61 namespace Foam
62 {
63 namespace incompressible
64 {
65 namespace LESModels
66 {
67 
68 /*---------------------------------------------------------------------------*\
69  Class spectEddyVisc Declaration
70 \*---------------------------------------------------------------------------*/
71 
73 :
74  public GenEddyVisc
75 {
76  // Private data
77 
79  dimensionedScalar cK1_;
80  dimensionedScalar cK2_;
81  dimensionedScalar cK3_;
82  dimensionedScalar cK4_;
83 
84 
85  // Private Member Functions
86 
87  //- Update sub-grid scale fields
88  void updateSubGridScaleFields(const volTensorField& gradU);
89 
90  // Disallow default bitwise copy construct and assignment
92  spectEddyVisc& operator=(const spectEddyVisc&);
93 
94 
95 public:
96 
97  //- Runtime type information
98  TypeName("spectEddyVisc");
99 
100  // Constructors
101 
102  // Construct from components
104  (
105  const volVectorField& U,
106  const surfaceScalarField& phi,
108  );
109 
110 
111  //- Destructor
112  virtual ~spectEddyVisc()
113  {}
114 
115 
116  // Member Functions
117 
118  //- Return SGS kinetic energy
119  virtual tmp<volScalarField> k() const;
120 
121  //- Correct Eddy-Viscosity and related properties
122  virtual void correct(const tmp<volTensorField>&);
123 
124  //- Read LESProperties dictionary
125  virtual bool read();
126 };
127 
128 
129 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
130 
131 } // End namespace LESModels
132 } // End namespace incompressible
133 } // End namespace Foam
134 
135 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
136 
137 #endif
138 
139 // ************************ vim: set sw=4 sts=4 et: ************************ //