src/finiteVolume/finiteVolume/snGradSchemes/limitedSnGrad/limitedSnGrad.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  11     OpenFOAM is free software; you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     option) any later version.
  15
  16     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  17     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  18     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM; if not, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 Description
  26     snGrad scheme with limited non-orthogonal correction.
  27
  28     The limiter is controlled by a coefficient with a value between 0 and 1
  29     which when 0 switches the correction off and the scheme behaves as
  30     uncorrectedSnGrad, when set to 1 the full correction is applied and the
  31     scheme behaves as correctedSnGrad and when set to 0.5 the limiter is
  32     calculated such that the non-orthogonal contribution does not exceed the
  33     orthogonal part.
  34
  35 \*---------------------------------------------------------------------------*/
  36
  37 #include "fv.H"
  38 #include "limitedSnGrad.H"
  39 #include "volFields.H"
  40 #include "surfaceFields.H"
  41 #include "correctedSnGrad.H"
  42 #include "localMax.H"
  43
  44 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  45
  46 namespace Foam
  47 {
  48
  49 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  50
  51 namespace fv
  52 {
  53
  54 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
  55
  56 template<class Type>
  57 limitedSnGrad<Type>::~limitedSnGrad()
  58 {}
  59
  60
  61 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  62
  63 template<class Type>
  64 tmp<GeometricField<Type, fvsPatchField, surfaceMesh> >
  65 limitedSnGrad<Type>::correction
  66 (
  67     const GeometricField<Type, fvPatchField, volMesh>& vf
  68 ) const
  69 {
  70     GeometricField<Type, fvsPatchField, surfaceMesh> corr =
  71         correctedSnGrad<Type>(this->mesh()).correction(vf);
  72
  73     surfaceScalarField limiter
  74     (
  75         min
  76         (
  77             limitCoeff_
  78            *mag(snGradScheme<Type>::snGrad(vf, deltaCoeffs(vf), "orthSnGrad"))
  79            /(
  80                 (1 - limitCoeff_)*mag(corr)
  81               + dimensionedScalar("small", corr.dimensions(), SMALL)
  82             ),
  83             dimensionedScalar("one", dimless, 1.0)
  84         )
  85     );
  86
  87     if (fv::debug)
  88     {
  89         Info<< "limitedSnGrad :: limiter min: " << min(limiter.internalField())
  90             << " max: "<< max(limiter.internalField())
  91             << " avg: " << average(limiter.internalField()) << endl;
  92     }
  93
  94     return limiter*corr;
  95 }
  96
  97
  98 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  99
 100 } // End namespace fv
 101
 102 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 103
 104 } // End namespace Foam
 105
 106 // ************************************************************************* //