src/turbulenceModels/LES/incompressible/locDynOneEqEddy/locDynOneEqEddy.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2008 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 \*---------------------------------------------------------------------------*/
  26
  27 #include "locDynOneEqEddy.H"
  28 #include "addToRunTimeSelectionTable.H"
  29
  30 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  31
  32 namespace Foam
  33 {
  34 namespace incompressible
  35 {
  36 namespace LESModels
  37 {
  38
  39 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  40
  41 defineTypeNameAndDebug(locDynOneEqEddy, 0);
  42 addToRunTimeSelectionTable(LESModel, locDynOneEqEddy, dictionary);
  43
  44 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  45
  46 volScalarField locDynOneEqEddy::ck
  47 (
  48     const volSymmTensorField& D,
  49     const volScalarField& KK
  50 ) const
  51 {
  52     volSymmTensorField LL =
  53         simpleFilter_(dev(filter_(sqr(U())) - (sqr(filter_(U())))));
  54
  55     volSymmTensorField MM = simpleFilter_(-2.0*delta()*pow(KK, 0.5)*filter_(D));
  56
  57     volScalarField ck =
  58         simpleFilter_(0.5*(LL && MM))
  59        /(
  60             simpleFilter_(magSqr(MM))
  61           + dimensionedScalar("small", sqr(MM.dimensions()), VSMALL)
  62         );
  63
  64     return 0.5*(mag(ck) + ck);
  65 }
  66
  67
  68 volScalarField locDynOneEqEddy::ce
  69 (
  70     const volSymmTensorField& D,
  71     const volScalarField& KK
  72 ) const
  73 {
  74     volScalarField ce =
  75         simpleFilter_(nuEff()*(filter_(magSqr(D)) - magSqr(filter_(D))))
  76        /simpleFilter_(pow(KK, 1.5)/(2.0*delta()));
  77
  78     return 0.5*(mag(ce) + ce);
  79 }
  80
  81
  82 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  83
  84 locDynOneEqEddy::locDynOneEqEddy
  85 (
  86     const volVectorField& U,
  87     const surfaceScalarField& phi,
  88     transportModel& transport
  89 )
  90 :
  91     LESModel(typeName, U, phi, transport),
  92     GenEddyVisc(U, phi, transport),
  93
  94     k_
  95     (
  96         IOobject
  97         (
  98             "k",
  99             runTime_.timeName(),
 100             mesh_,
 101             IOobject::MUST_READ,
 102             IOobject::AUTO_WRITE
 103         ),
 104         mesh_
 105     ),
 106
 107     simpleFilter_(U.mesh()),
 108     filterPtr_(LESfilter::New(U.mesh(), coeffDict())),
 109     filter_(filterPtr_())
 110 {
 111     printCoeffs();
 112 }
 113
 114
 115 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 116
 117 locDynOneEqEddy::~locDynOneEqEddy()
 118 {}
 119
 120
 121 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 122
 123 void locDynOneEqEddy::correct(const tmp<volTensorField>& gradU)
 124 {
 125     LESModel::correct(gradU);
 126
 127     volSymmTensorField D = symm(gradU);
 128
 129     volScalarField KK = 0.5*(filter_(magSqr(U())) - magSqr(filter_(U())));
 130     KK.max(dimensionedScalar("small", KK.dimensions(), SMALL));
 131
 132     volScalarField P = 2.0*nuSgs_*magSqr(D);
 133
 134     solve
 135     (
 136        fvm::ddt(k_)
 137      + fvm::div(phi(), k_)
 138      - fvm::laplacian(DkEff(), k_)
 139     ==
 140        P
 141      - fvm::Sp(ce(D, KK)*sqrt(k_)/delta(), k_)
 142     );
 143
 144     bound(k_, k0());
 145
 146     nuSgs_ = ck(D, KK)*sqrt(k_)*delta();
 147     nuSgs_.correctBoundaryConditions();
 148 }
 149
 150
 151 bool locDynOneEqEddy::read()
 152 {
 153     if (GenEddyVisc::read())
 154     {
 155         filter_.read(coeffDict());
 156
 157         return true;
 158     }
 159     else
 160     {
 161         return false;
 162     }
 163 }
 164
 165
 166 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 167
 168 } // End namespace LESModels
 169 } // End namespace incompressible
 170 } // End namespace Foam
 171
 172 // ************************************************************************* //