00001 /*---------------------------------------------------------------------------*\ 00002 ========= | 00003 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox 00004 \\ / O peration | 00005 \\ / A nd | Copyright (C) 1991-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::OSPRELimiter 00026 00027 Description 00028 Class with limiter function which returns the limiter for the 00029 OSPRE differencing scheme based on r obtained from the LimiterFunc 00030 class. 00031 00032 Used in conjunction with the template class LimitedScheme. 00033 00034 SourceFiles 00035 OSPRE.C 00036 00037 \*---------------------------------------------------------------------------*/ 00038 00039 #ifndef OSPRE_H 00040 #define OSPRE_H 00041 00042 #include <OpenFOAM/vector.H> 00043 00044 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00045 00046 namespace Foam 00047 { 00048 00049 /*---------------------------------------------------------------------------*\ 00050 Class OSPRELimiter Declaration 00051 \*---------------------------------------------------------------------------*/ 00052 00053 template<class LimiterFunc> 00054 class OSPRELimiter 00055 : 00056 public LimiterFunc 00057 { 00058 00059 public: 00060 00061 OSPRELimiter(Istream&) 00062 {} 00063 00064 scalar limiter 00065 ( 00066 const scalar cdWeight, 00067 const scalar faceFlux, 00068 const typename LimiterFunc::phiType& phiP, 00069 const typename LimiterFunc::phiType& phiN, 00070 const typename LimiterFunc::gradPhiType& gradcP, 00071 const typename LimiterFunc::gradPhiType& gradcN, 00072 const vector& d 00073 ) const 00074 { 00075 scalar r = LimiterFunc::r 00076 ( 00077 faceFlux, phiP, phiN, gradcP, gradcN, d 00078 ); 00079 00080 scalar rrp1 = r*(r + 1); 00081 return 1.5*rrp1/(rrp1 + 1); 00082 } 00083 }; 00084 00085 00086 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00087 00088 } // End namespace Foam 00089 00090 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00091 00092 #endif 00093 00094 // ************************ vim: set sw=4 sts=4 et: ************************ //
