FreeFOAM The Cross-Platform CFD Toolkit
Hosted by SourceForge:
Get FreeFOAM at SourceForge.net.
            Fast, secure and Free Open Source software downloads

MUSCL.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-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: ************************ //
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines