OpenFOAM logo
Open Source CFD Toolkit

RNGkEpsilon.H

Go to the documentation of this file.
00001 /*---------------------------------------------------------------------------*\
00002   =========                 |
00003   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
00004    \\    /   O peration     |
00005     \\  /    A nd           | Copyright (C) 1991-2005 OpenCFD Ltd.
00006      \\/     M anipulation  |
00007 -------------------------------------------------------------------------------
00008 License
00009     This file is part of OpenFOAM.
00010 
00011     OpenFOAM is free software; you can redistribute it and/or modify it
00012     under the terms of the GNU General Public License as published by the
00013     Free Software Foundation; either version 2 of the License, or (at your
00014     option) any later version.
00015 
00016     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
00017     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
00018     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
00019     for more details.
00020 
00021     You should have received a copy of the GNU General Public License
00022     along with OpenFOAM; if not, write to the Free Software Foundation,
00023     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
00024 
00025 Class
00026     RNGkEpsilon
00027 
00028 Description
00029     Renormalisation group k-epsilon turbulence model.
00030 
00031 SourceFiles
00032     RNGkEpsilon.C
00033     RNGkEpsilonCorrect.C
00034 
00035 \*---------------------------------------------------------------------------*/
00036 
00037 #ifndef compressibleRNGkEpsilon_H
00038 #define compressibleRNGkEpsilon_H
00039 
00040 #include "turbulenceModel.H"
00041 
00042 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00043 
00044 namespace Foam
00045 {
00046 namespace compressible
00047 {
00048 namespace turbulenceModels
00049 {
00050 
00051 /*---------------------------------------------------------------------------*\
00052                            Class RNGkEpsilon Declaration
00053 \*---------------------------------------------------------------------------*/
00054 
00055 class RNGkEpsilon
00056 :
00057     public turbulenceModel
00058 {
00059     // Private data
00060 
00061         dimensionedScalar Cmu;
00062         dimensionedScalar C1;
00063         dimensionedScalar C2;
00064         dimensionedScalar C3;
00065         dimensionedScalar alphak;
00066         dimensionedScalar alphaEps;
00067         dimensionedScalar alphah;
00068         dimensionedScalar eta0;
00069         dimensionedScalar beta;
00070 
00071         volScalarField k_;
00072         volScalarField epsilon_;
00073         volScalarField mut_;
00074 
00075 
00076 public:
00077 
00078     //- Runtime type information
00079     TypeName("RNGkEpsilon");
00080 
00081     // Constructors
00082 
00083         //- from components
00084         RNGkEpsilon
00085         (
00086             const volScalarField& rho,
00087             const volVectorField& U,
00088             const surfaceScalarField& phi,
00089             basicThermo& thermophysicalModel
00090         );
00091 
00092 
00093     // Destructor
00094 
00095         ~RNGkEpsilon(){}
00096 
00097 
00098     // Member Functions
00099 
00100         tmp<volScalarField> mut() const
00101         {
00102             return mut_;
00103         }
00104 
00105         //- Return the effective diffusivity for k
00106         tmp<volScalarField> DkEff() const
00107         {
00108             return tmp<volScalarField> 
00109             (
00110                 new volScalarField("DkEff", alphak*mut_ + mu())
00111             );
00112         }
00113 
00114         //- Return the effective diffusivity for epsilon
00115         tmp<volScalarField> DepsilonEff() const
00116         {
00117             return tmp<volScalarField> 
00118             (
00119                 new volScalarField("DepsilonEff", alphaEps*mut_ + mu())
00120             );
00121         }
00122 
00123         //- Return the effective turbulent thermal diffusivity
00124         tmp<volScalarField> alphaEff() const
00125         {
00126             return tmp<volScalarField> 
00127             (
00128                 new volScalarField("alphaEff", alphah*mut_ + alpha())
00129             );
00130         }
00131 
00132         tmp<volScalarField> k() const
00133         {
00134             return k_;
00135         }
00136 
00137         tmp<volScalarField> epsilon() const
00138         {
00139             return epsilon_;
00140         }
00141 
00142         tmp<volTensorField> R() const;
00143 
00144         tmp<fvVectorMatrix> divRhoR(volVectorField& U) const;
00145 
00146         tmp<fvScalarMatrix> divRhoUh(volScalarField& h) const;
00147 
00148         void correct();
00149 
00150         //- Read turbulenceProperties dictionary
00151         bool read();
00152 };
00153 
00154 
00155 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00156 
00157 } // End namespace turbulenceModels
00158 } // End namespace compressible
00159 } // End namespace Foam
00160 
00161 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00162 
00163 #endif
00164 
00165 // ************************************************************************* //
For further information go to www.openfoam.org