00001 /*---------------------------------------------------------------------------*\ 00002 ========= | 00003 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox 00004 \\ / O peration | 00005 \\ / A nd | Copyright (C) 2008-2010 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 00013 the Free Software Foundation, either version 3 of the License, or 00014 (at your 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, see <http://www.gnu.org/licenses/>. 00023 00024 Class 00025 Foam::LESmodels::SpalartAllmarasIDDES 00026 00027 Description 00028 SpalartAllmarasIDDES LES turbulence model for incompressible flows 00029 00030 SourceFiles 00031 SpalartAllmarasIDDES.C 00032 00033 \*---------------------------------------------------------------------------*/ 00034 00035 #ifndef SpalartAllmarasIDDES_H 00036 #define SpalartAllmarasIDDES_H 00037 00038 #include <incompressibleLESModels/SpalartAllmaras.H> 00039 #include <incompressibleLESModels/IDDESDelta.H> 00040 00041 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00042 00043 namespace Foam 00044 { 00045 namespace incompressible 00046 { 00047 namespace LESModels 00048 { 00049 00050 /*---------------------------------------------------------------------------*\ 00051 Class SpalartAllmarasIDDES Declaration 00052 \*---------------------------------------------------------------------------*/ 00053 00054 class SpalartAllmarasIDDES 00055 : 00056 public SpalartAllmaras 00057 { 00058 // Private data 00059 00060 // Model constants 00061 00062 autoPtr<LESdelta> hmax_; 00063 autoPtr<LESdelta> IDDESDelta_; 00064 dimensionedScalar fwStar_; 00065 dimensionedScalar cl_; 00066 dimensionedScalar ct_; 00067 00068 00069 // Private member functions 00070 00071 tmp<volScalarField> alpha() const; 00072 tmp<volScalarField> ft(const volScalarField& S) const; 00073 tmp<volScalarField> fl(const volScalarField& S) const; 00074 00075 tmp<volScalarField> rd 00076 ( 00077 const volScalarField& visc, 00078 const volScalarField& S 00079 ) const; 00080 00081 //- Delay function 00082 tmp<volScalarField> fd(const volScalarField& S) const; 00083 00084 // Disallow default bitwise copy construct and assignment 00085 SpalartAllmarasIDDES(const SpalartAllmarasIDDES&); 00086 SpalartAllmarasIDDES& operator=(const SpalartAllmarasIDDES&); 00087 00088 00089 protected: 00090 00091 // Protected member functions 00092 00093 //- Length scale 00094 virtual tmp<volScalarField> dTilda(const volScalarField& S) const; 00095 00096 00097 public: 00098 00099 //- Runtime type information 00100 TypeName("SpalartAllmarasIDDES"); 00101 00102 00103 // Constructors 00104 00105 //- Construct from components 00106 SpalartAllmarasIDDES 00107 ( 00108 const volVectorField& U, 00109 const surfaceScalarField& phi, 00110 transportModel& transport 00111 ); 00112 00113 00114 //- Destructor 00115 virtual ~SpalartAllmarasIDDES() 00116 {} 00117 00118 00119 // Member Functions 00120 00121 //- Access function to filter width 00122 virtual const volScalarField& delta() const 00123 { 00124 return IDDESDelta_(); 00125 } 00126 00127 //- Read LESProperties dictionary 00128 virtual bool read(); 00129 }; 00130 00131 00132 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00133 00134 } // End namespace LESModels 00135 } // End namespace incompressible 00136 } // End namespace Foam 00137 00138 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00139 00140 #endif 00141 00142 // ************************ vim: set sw=4 sts=4 et: ************************ //