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 \*---------------------------------------------------------------------------*/ 00025 00026 #include "bound.H" 00027 #include <finiteVolume/volFields.H> 00028 #include <finiteVolume/fvc.H> 00029 00030 // * * * * * * * * * * * * * * * Global Functions * * * * * * * * * * * * * // 00031 00032 void Foam::bound(volScalarField& vsf, const dimensionedScalar& vsf0) 00033 { 00034 scalar minVsf = min(vsf).value(); 00035 00036 if (minVsf < vsf0.value()) 00037 { 00038 Info<< "bounding " << vsf.name() 00039 << ", min: " << minVsf 00040 << " max: " << max(vsf).value() 00041 << " average: " << gAverage(vsf.internalField()) 00042 << endl; 00043 00044 vsf.internalField() = max 00045 ( 00046 max 00047 ( 00048 vsf.internalField(), 00049 fvc::average(max(vsf, vsf0))().internalField() 00050 *pos(-vsf.internalField()) 00051 ), 00052 vsf0.value() 00053 ); 00054 00055 vsf.boundaryField() = max(vsf.boundaryField(), vsf0.value()); 00056 } 00057 } 00058 00059 00060 // ************************ vim: set sw=4 sts=4 et: ************************ //