src/turbulenceModels/LES/incompressible/dynOneEqEddy/dynOneEqEddy.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 "dynOneEqEddy.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(dynOneEqEddy, 0);
  42 addToRunTimeSelectionTable(LESModel, dynOneEqEddy, dictionary);
  43
  44 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  45
  46 dimensionedScalar dynOneEqEddy::ck(const volSymmTensorField& D) const
  47 {
  48     volScalarField KK = 0.5*(filter_(magSqr(U())) - magSqr(filter_(U())));
  49
  50     volSymmTensorField LL = dev(filter_(sqr(U())) - sqr(filter_(U())));
  51
  52     volSymmTensorField MM =
  53         delta()*(filter_(sqrt(k_)*D) - 2*sqrt(KK + filter_(k_))*filter_(D));
  54
  55     dimensionedScalar MMMM = average(magSqr(MM));
  56
  57     if (MMMM.value() > VSMALL)
  58     {
  59         return average(LL && MM)/MMMM;
  60     }
  61     else
  62     {
  63         return 0.0;
  64     }
  65 }
  66
  67
  68 dimensionedScalar dynOneEqEddy::ce(const volSymmTensorField& D) const
  69 {
  70     volScalarField KK = 0.5*(filter_(magSqr(U())) - magSqr(filter_(U())));
  71
  72     volScalarField mm =
  73         pow(KK + filter_(k_), 1.5)/(2*delta()) - filter_(pow(k_, 1.5))/delta();
  74
  75     volScalarField ee =
  76         2*delta()*ck(D)
  77        *(
  78            filter_(sqrt(k_)*magSqr(D))
  79          - 2*sqrt(KK + filter_(k_))*magSqr(filter_(D))
  80         );
  81
  82     dimensionedScalar mmmm = average(magSqr(mm));
  83
  84     if (mmmm.value() > VSMALL)
  85     {
  86         return average(ee*mm)/mmmm;
  87     }
  88     else
  89     {
  90         return 0.0;
  91     }
  92 }
  93
  94
  95 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  96
  97 dynOneEqEddy::dynOneEqEddy
  98 (
  99     const volVectorField& U,
 100     const surfaceScalarField& phi,
 101     transportModel& transport
 102 )
 103 :
 104     LESModel(typeName, U, phi, transport),
 105     GenEddyVisc(U, phi, transport),
 106
 107     k_
 108     (
 109         IOobject
 110         (
 111             "k",
 112             runTime_.timeName(),
 113             mesh_,
 114             IOobject::MUST_READ,
 115             IOobject::AUTO_WRITE
 116         ),
 117         mesh_
 118     ),
 119
 120     filterPtr_(LESfilter::New(U.mesh(), coeffDict())),
 121     filter_(filterPtr_())
 122 {
 123     printCoeffs();
 124 }
 125
 126
 127 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 128
 129 dynOneEqEddy::~dynOneEqEddy()
 130 {}
 131
 132
 133 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 134
 135 void dynOneEqEddy::correct(const tmp<volTensorField>& gradU)
 136 {
 137     GenEddyVisc::correct(gradU);
 138
 139     volSymmTensorField D = symm(gradU);
 140
 141     volScalarField P = 2.0*nuSgs_*magSqr(D);
 142
 143     solve
 144     (
 145        fvm::ddt(k_)
 146      + fvm::div(phi(), k_)
 147      - fvm::laplacian(DkEff(), k_)
 148     ==
 149        P
 150      - fvm::Sp(ce(D)*sqrt(k_)/delta(), k_)
 151     );
 152
 153     bound(k_, k0());
 154
 155     nuSgs_ = ck(D)*sqrt(k_)*delta();
 156     nuSgs_.correctBoundaryConditions();
 157 }
 158
 159
 160 bool dynOneEqEddy::read()
 161 {
 162     if (GenEddyVisc::read())
 163     {
 164         filter_.read(coeffDict());
 165
 166         return true;
 167     }
 168     else
 169     {
 170         return false;
 171     }
 172 }
 173
 174
 175 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 176
 177 } // End namespace LESModels
 178 } // End namespace incompressible
 179 } // End namespace Foam
 180
 181 // ************************************************************************* //