OpenFOAM logo
Open Source CFD Toolkit

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