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

bound.C

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