OpenFOAM logo
Open Source CFD Toolkit

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-2005 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 the
00013     Free Software Foundation; either version 2 of the License, or (at your
00014     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, write to the Free Software Foundation,
00023     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
00024 
00025 Class
00026     MUSCLLimiter
00027 
00028 Description
00029     Class with limiter function which returns the limiter for the
00030     van Leer's MUSCL differencing scheme.
00031 
00032     Used in conjunction with the template class LimitedScheme.
00033 
00034 SourceFiles
00035     MUSCL.C
00036 
00037 \*---------------------------------------------------------------------------*/
00038 
00039 #ifndef MUSCL_H
00040 #define MUSCL_H
00041 
00042 #include "vector.H"
00043 
00044 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00045 
00046 namespace Foam
00047 {
00048 
00049 /*---------------------------------------------------------------------------*\
00050                            Class MUSCLWeight Declaration
00051 \*---------------------------------------------------------------------------*/
00052 
00053 template<class LimiterFunc>
00054 class MUSCLLimiter
00055 :
00056     public LimiterFunc
00057 {
00058 
00059 public:
00060 
00061     MUSCLLimiter(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             cdWeight, faceFlux, phiP, phiN, gradcP, gradcN, d
00078         );
00079 
00080         return max(min(min(2*r, 0.5*r + 0.5), 2), 0);
00081     }
00082 };
00083 
00084 
00085 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00086 
00087 } // End namespace Foam
00088 
00089 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
00090 
00091 #endif
00092 
00093 // ************************************************************************* //
For further information go to www.openfoam.org