src/turbulenceModels/incompressible/LES/mixedSmagorinsky/mixedSmagorinsky.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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 "mixedSmagorinsky.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(mixedSmagorinsky, 0);
  41 addToRunTimeSelectionTable(LESModel, mixedSmagorinsky, dictionary);
  42
  43 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  44
  45 mixedSmagorinsky::mixedSmagorinsky
  46 (
  47     const volVectorField& U,
  48     const surfaceScalarField& phi,
  49     transportModel& transport,
  50     const word& turbulenceModelName,
  51     const word& modelName
  52 )
  53 :
  54     LESModel(modelName, U, phi, transport, turbulenceModelName),
  55     scaleSimilarity(U, phi, transport),
  56     Smagorinsky(U, phi, transport)
  57 {
  58     printCoeffs();
  59 }
  60
  61
  62 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  63
  64 tmp<volScalarField> mixedSmagorinsky::k() const
  65 {
  66     return
  67     (
  68         scaleSimilarity::k()
  69       + Smagorinsky::k()
  70     );
  71 }
  72
  73
  74 tmp<volScalarField> mixedSmagorinsky::epsilon() const
  75 {
  76     return
  77     (
  78         scaleSimilarity::epsilon()
  79       + Smagorinsky::epsilon()
  80     );
  81 }
  82
  83
  84 tmp<volSymmTensorField> mixedSmagorinsky::B() const
  85 {
  86     return
  87     (
  88         scaleSimilarity::B()
  89       + Smagorinsky::B()
  90     );
  91 }
  92
  93
  94 tmp<volSymmTensorField> mixedSmagorinsky::devBeff() const
  95 {
  96     return
  97     (
  98         scaleSimilarity::devBeff()
  99       + Smagorinsky::devBeff()
 100     );
 101 }
 102
 103
 104 tmp<fvVectorMatrix> mixedSmagorinsky::divDevBeff
 105 (
 106     volVectorField& U
 107 ) const
 108 {
 109     return
 110     (
 111         scaleSimilarity::divDevBeff(U)
 112       + Smagorinsky::divDevBeff(U)
 113     );
 114 }
 115
 116
 117 void mixedSmagorinsky::correct(const tmp<volTensorField>& gradU)
 118 {
 119     scaleSimilarity::correct(gradU);
 120     Smagorinsky::correct(gradU);
 121 }
 122
 123
 124 bool mixedSmagorinsky::read()
 125 {
 126     if (LESModel::read())
 127     {
 128         scaleSimilarity::read();
 129         Smagorinsky::read();
 130
 131         return true;
 132     }
 133     else
 134     {
 135         return false;
 136     }
 137 }
 138
 139
 140 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 141
 142 } // End namespace LESModels
 143 } // End namespace incompressible
 144 } // End namespace Foam
 145
 146 // ************************************************************************* //