src/turbulenceModels/incompressible/LES/oneEqEddy/oneEqEddy.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 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 "oneEqEddy.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(oneEqEddy, 0);
  42 addToRunTimeSelectionTable(LESModel, oneEqEddy, dictionary);
  43
  44
  45 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  46
  47 void oneEqEddy::updateSubGridScaleFields()
  48 {
  49     nuSgs_ = ck_*sqrt(k_)*delta();
  50     nuSgs_.correctBoundaryConditions();
  51 }
  52
  53
  54 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  55
  56 oneEqEddy::oneEqEddy
  57 (
  58     const volVectorField& U,
  59     const surfaceScalarField& phi,
  60     transportModel& transport
  61 )
  62 :
  63     LESModel(typeName, U, phi, transport),
  64     GenEddyVisc(U, phi, transport),
  65
  66     k_
  67     (
  68         IOobject
  69         (
  70             "k",
  71             runTime_.timeName(),
  72             mesh_,
  73             IOobject::MUST_READ,
  74             IOobject::AUTO_WRITE
  75         ),
  76         mesh_
  77     ),
  78
  79     ck_
  80     (
  81         dimensioned<scalar>::lookupOrAddToDict
  82         (
  83             "ck",
  84             coeffDict_,
  85             0.094
  86         )
  87     )
  88 {
  89     updateSubGridScaleFields();
  90
  91     printCoeffs();
  92 }
  93
  94
  95 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  96
  97 void oneEqEddy::correct(const tmp<volTensorField>& gradU)
  98 {
  99     GenEddyVisc::correct(gradU);
 100
 101     volScalarField G = 2.0*nuSgs_*magSqr(symm(gradU));
 102
 103     fvScalarMatrix kEqn
 104     (
 105        fvm::ddt(k_)
 106      + fvm::div(phi(), k_)
 107      - fvm::laplacian(DkEff(), k_)
 108     ==
 109        G
 110      - fvm::Sp(ce_*sqrt(k_)/delta(), k_)
 111     );
 112
 113     kEqn.relax();
 114     kEqn.solve();
 115
 116     bound(k_, k0());
 117
 118     updateSubGridScaleFields();
 119 }
 120
 121
 122 bool oneEqEddy::read()
 123 {
 124     if (GenEddyVisc::read())
 125     {
 126         ck_.readIfPresent(coeffDict());
 127
 128         return true;
 129     }
 130     else
 131     {
 132         return false;
 133     }
 134 }
 135
 136
 137 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 138
 139 } // End namespace LESModels
 140 } // End namespace incompressible
 141 } // End namespace Foam
 142
 143 // ************************************************************************* //