src/turbulenceModels/incompressible/LES/dynMixedSmagorinsky/dynMixedSmagorinsky.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend 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 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #include "dynMixedSmagorinsky.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33 namespace incompressible
  34 {
  35 namespace LESModels
  36 {
  37
  38 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  39
  40 defineTypeNameAndDebug(dynMixedSmagorinsky, 0);
  41 addToRunTimeSelectionTable(LESModel, dynMixedSmagorinsky, dictionary);
  42
  43 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  44
  45 dynMixedSmagorinsky::dynMixedSmagorinsky
  46 (
  47     const volVectorField& U,
  48     const surfaceScalarField& phi,
  49     transportModel& transport
  50 )
  51 :
  52     LESModel(typeName, U, phi, transport),
  53     scaleSimilarity(U, phi, transport),
  54     dynSmagorinsky(U, phi, transport)
  55 {
  56     printCoeffs();
  57 }
  58
  59
  60 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  61
  62 tmp<volScalarField> dynMixedSmagorinsky::k() const
  63 {
  64     return
  65     (
  66         scaleSimilarity::k()
  67       + dynSmagorinsky::k()
  68     );
  69 }
  70
  71
  72 tmp<volScalarField> dynMixedSmagorinsky::epsilon() const
  73 {
  74     return
  75     (
  76         scaleSimilarity::epsilon()
  77       + dynSmagorinsky::epsilon()
  78     );
  79 }
  80
  81
  82 tmp<volSymmTensorField> dynMixedSmagorinsky::B() const
  83 {
  84     return
  85     (
  86         scaleSimilarity::B()
  87       + dynSmagorinsky::B()
  88     );
  89 }
  90
  91
  92 tmp<volSymmTensorField> dynMixedSmagorinsky::devBeff() const
  93 {
  94     return
  95     (
  96         scaleSimilarity::devBeff()
  97       + dynSmagorinsky::devBeff()
  98     );
  99 }
 100
 101
 102 tmp<fvVectorMatrix> dynMixedSmagorinsky::divDevBeff(volVectorField& U) const
 103 {
 104     return
 105     (
 106         scaleSimilarity::divDevBeff(U)
 107       + dynSmagorinsky::divDevBeff(U)
 108     );
 109 }
 110
 111
 112 void dynMixedSmagorinsky::correct(const tmp<volTensorField>& gradU)
 113 {
 114     scaleSimilarity::correct(gradU);
 115     dynSmagorinsky::correct(gradU());
 116 }
 117
 118
 119 bool dynMixedSmagorinsky::read()
 120 {
 121     if (LESModel::read())
 122     {
 123         scaleSimilarity::read();
 124         dynSmagorinsky::read();
 125
 126         return true;
 127     }
 128     else
 129     {
 130         return false;
 131     }
 132 }
 133
 134
 135 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 136
 137 } // End namespace LESModels
 138 } // namespace incompressible
 139 } // End namespace Foam
 140
 141 // ************************************************************************* //