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::MUSCLLimiter 00026 00027 Description 00028 Class with limiter function which returns the limiter for the 00029 van Leer's MUSCL differencing scheme. 00030 00031 Used in conjunction with the template class LimitedScheme. 00032 00033 SourceFiles 00034 MUSCL.C 00035 00036 \*---------------------------------------------------------------------------*/ 00037 00038 #ifndef MUSCL_H 00039 #define MUSCL_H 00040 00041 #include <OpenFOAM/vector.H> 00042 00043 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00044 00045 namespace Foam 00046 { 00047 00048 /*---------------------------------------------------------------------------*\ 00049 Class MUSCLLimiter Declaration 00050 \*---------------------------------------------------------------------------*/ 00051 00052 template<class LimiterFunc> 00053 class MUSCLLimiter 00054 : 00055 public LimiterFunc 00056 { 00057 00058 public: 00059 00060 MUSCLLimiter(Istream&) 00061 {} 00062 00063 scalar limiter 00064 ( 00065 const scalar cdWeight, 00066 const scalar faceFlux, 00067 const typename LimiterFunc::phiType& phiP, 00068 const typename LimiterFunc::phiType& phiN, 00069 const typename LimiterFunc::gradPhiType& gradcP, 00070 const typename LimiterFunc::gradPhiType& gradcN, 00071 const vector& d 00072 ) const 00073 { 00074 scalar r = LimiterFunc::r 00075 ( 00076 faceFlux, phiP, phiN, gradcP, gradcN, d 00077 ); 00078 00079 return max(min(min(2*r, 0.5*r + 0.5), 2), 0); 00080 } 00081 }; 00082 00083 00084 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00085 00086 } // End namespace Foam 00087 00088 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00089 00090 #endif 00091 00092 // ************************ vim: set sw=4 sts=4 et: ************************ //