OpenFOAM logo
Open Source CFD Toolkit

LaunderSharmaKE.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     LaunderSharmaKE
00027 
00028 Description
00029     Launder and Sharma low-Reynolds k-epsilon turbulence model for
00030     compressible and combusting flows.
00031 
00032 SourceFiles
00033     LaunderSharmaKE.C
00034     LaunderSharmaKECorrect.C
00035 
00036 \*---------------------------------------------------------------------------*/
00037 
00038 #ifndef compressibleLaunderSharmaKE_H
00039 #define compressibleLaunderSharmaKE_H
00040 
00041 #include "turbulenceModel.H"
00042 
00043 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00044 
00045 namespace Foam
00046 {
00047 namespace compressible
00048 {
00049 namespace turbulenceModels
00050 {
00051 
00052 /*---------------------------------------------------------------------------*\
00053                            Class LaunderSharmaKE Declaration
00054 \*---------------------------------------------------------------------------*/
00055 
00056 class LaunderSharmaKE
00057 :
00058     public turbulenceModel
00059 {
00060     // Private data
00061 
00062         dimensionedScalar Cmu;
00063         dimensionedScalar C1;
00064         dimensionedScalar C2;
00065         dimensionedScalar C3;
00066         dimensionedScalar alphak;
00067         dimensionedScalar alphaEps;
00068         dimensionedScalar alphah;
00069 
00070         volScalarField k_;
00071         volScalarField epsilon_;
00072         volScalarField mut_;
00073 
00074 
00075     // Private member functions
00076 
00077         tmp<volScalarField> fMu() const;
00078         tmp<volScalarField> f2() const;
00079 
00080 
00081 public:
00082 
00083     //- Runtime type information
00084     TypeName("LaunderSharmaKE");
00085 
00086     // Constructors
00087 
00088         //- from components
00089         LaunderSharmaKE
00090         (
00091             const volScalarField& rho,
00092             const volVectorField& U,
00093             const surfaceScalarField& phi,
00094             basicThermo& thermophysicalModel
00095         );
00096 
00097 
00098     // Destructor
00099 
00100         ~LaunderSharmaKE(){}
00101 
00102 
00103     // Member Functions
00104 
00105         tmp<volScalarField> mut() const
00106         {
00107             return mut_;
00108         }
00109 
00110         //- Return the effective diffusivity for k
00111         tmp<volScalarField> DkEff() const
00112         {
00113             return tmp<volScalarField> 
00114             (
00115                 new volScalarField("DkEff", alphak*mut_ + mu())
00116             );
00117         }
00118 
00119         //- Return the effective diffusivity for epsilon
00120         tmp<volScalarField> DepsilonEff() const
00121         {
00122             return tmp<volScalarField> 
00123             (
00124                 new volScalarField("DepsilonEff", alphaEps*mut_ + mu())
00125             );
00126         }
00127 
00128         //- Return the effective turbulent thermal diffusivity
00129         tmp<volScalarField> alphaEff() const
00130         {
00131             return tmp<volScalarField> 
00132             (
00133                 new volScalarField("alphaEff", alphah*mut_ + alpha())
00134             );
00135         }
00136 
00137         tmp<volScalarField> k() const
00138         {
00139             return k_;
00140         }
00141 
00142         tmp<volScalarField> epsilon() const
00143         {
00144             return epsilon_;
00145         }
00146 
00147         tmp<volTensorField> R() const;
00148 
00149         tmp<fvVectorMatrix> divRhoR(volVectorField& U) const;
00150 
00151         tmp<fvScalarMatrix> divRhoUh(volScalarField& h) const;
00152 
00153         void correct();
00154 
00155         //- Read turbulenceProperties dictionary
00156         bool read();
00157 };
00158 
00159 
00160 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00161 
00162 } // End namespace turbulenceModels
00163 } // End namespace compressible
00164 } // End namespace Foam
00165 
00166 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00167 
00168 #endif
00169 
00170 // ************************************************************************* //
For further information go to www.openfoam.org