src/turbulenceModels/incompressible/LES/scaleSimilarity/scaleSimilarity.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2010 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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your 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, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "scaleSimilarity.H"
  27
  28 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  29
  30 namespace Foam
  31 {
  32 namespace incompressible
  33 {
  34 namespace LESModels
  35 {
  36
  37 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  38
  39 defineTypeNameAndDebug(scaleSimilarity, 0);
  40
  41 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  42
  43 scaleSimilarity::scaleSimilarity
  44 (
  45     const volVectorField& U,
  46     const surfaceScalarField& phi,
  47     transportModel& transport
  48 )
  49 :
  50     LESModel(typeName, U, phi, transport),
  51     filterPtr_(LESfilter::New(U.mesh(), coeffDict())),
  52     filter_(filterPtr_())
  53 {
  54     printCoeffs();
  55 }
  56
  57
  58 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  59
  60 tmp<volScalarField> scaleSimilarity::k() const
  61 {
  62     return(0.5*(filter_(magSqr(U())) - magSqr(filter_(U()))));
  63 }
  64
  65
  66 tmp<volScalarField> scaleSimilarity::epsilon() const
  67 {
  68     volSymmTensorField D = symm(fvc::grad(U()));
  69
  70     return((filter_(sqr(U())) - sqr(filter_(U()))) && D);
  71 }
  72
  73
  74 tmp<volSymmTensorField> scaleSimilarity::B() const
  75 {
  76     return(filter_(sqr(U())) - sqr(filter_(U())));
  77 }
  78
  79
  80 tmp<volSymmTensorField> scaleSimilarity::devBeff() const
  81 {
  82     return dev(B());
  83 }
  84
  85
  86 tmp<fvVectorMatrix> scaleSimilarity::divDevBeff(volVectorField& U) const
  87 {
  88     return fvm::Su(fvc::div(devBeff()), U);
  89 }
  90
  91
  92 void scaleSimilarity::correct(const tmp<volTensorField>&)
  93 {}
  94
  95
  96 bool scaleSimilarity::read()
  97 {
  98     if (LESModel::read())
  99     {
 100         filter_.read(coeffDict());
 101
 102         return true;
 103     }
 104     else
 105     {
 106         return false;
 107     }
 108 }
 109
 110
 111 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 112
 113 } // End namespace LESModels
 114 } // End namespace incompressible
 115 } // End namespace Foam
 116
 117 // ************************************************************************* //