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::QUICKVLimiter 00026 00027 Description 00028 Class with limiter function which returns the limiter for the 00029 quadratic-upwind differencing scheme. 00030 00031 Note that the weighting factors are not bounded between upwind and 00032 central-differencing, some downwind contribution is possible although 00033 the interpolate is limited to be between the upwind and downwind cell 00034 values. 00035 00036 Used in conjunction with the template class LimitedScheme. 00037 00038 SourceFiles 00039 QUICKV.C 00040 00041 \*---------------------------------------------------------------------------*/ 00042 00043 #ifndef QUICKV_H 00044 #define QUICKV_H 00045 00046 #include <OpenFOAM/vector.H> 00047 00048 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00049 00050 namespace Foam 00051 { 00052 00053 /*---------------------------------------------------------------------------*\ 00054 Class QUICKVLimiter Declaration 00055 \*---------------------------------------------------------------------------*/ 00056 00057 template<class LimiterFunc> 00058 class QUICKVLimiter 00059 : 00060 public LimiterFunc 00061 { 00062 00063 public: 00064 00065 QUICKVLimiter(Istream&) 00066 {} 00067 00068 scalar limiter 00069 ( 00070 const scalar cdWeight, 00071 const scalar faceFlux, 00072 const typename LimiterFunc::phiType& phiP, 00073 const typename LimiterFunc::phiType& phiN, 00074 const typename LimiterFunc::gradPhiType& gradcP, 00075 const typename LimiterFunc::gradPhiType& gradcN, 00076 const vector& d 00077 ) const 00078 { 00079 vector gradfV = phiN - phiP; 00080 00081 scalar phiCD = gradfV & (cdWeight*phiP + (1 - cdWeight)*phiN); 00082 00083 scalar phiU, phif; 00084 00085 if (faceFlux > 0) 00086 { 00087 phiU = gradfV & phiP; 00088 phif = 0.5*(phiCD + phiU + (1 - cdWeight)*(gradfV & (d & gradcP))); 00089 } 00090 else 00091 { 00092 phiU = gradfV & phiN; 00093 phif = 0.5*(phiCD + phiU - cdWeight*(gradfV & (d & gradcN))); 00094 } 00095 00096 // Calculate the effective limiter for the QUICK interpolation 00097 scalar QLimiter = (phif - phiU)/stabilise(phiCD - phiU, SMALL); 00098 00099 // Limit the limiter between upwind and downwind 00100 return max(min(QLimiter, 2), 0); 00101 } 00102 }; 00103 00104 00105 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00106 00107 } // End namespace Foam 00108 00109 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00110 00111 #endif 00112 00113 // ************************ vim: set sw=4 sts=4 et: ************************ //
