src/turbulenceModels/compressible/LES/homogeneousDynOneEqEddy/homogeneousDynOneEqEddy.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 "homogeneousDynOneEqEddy.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33 namespace compressible
  34 {
  35 namespace LESModels
  36 {
  37
  38 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  39
  40 defineTypeNameAndDebug(homogeneousDynOneEqEddy, 0);
  41 addToRunTimeSelectionTable(LESModel, homogeneousDynOneEqEddy, dictionary);
  42
  43 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  44
  45 void homogeneousDynOneEqEddy::updateSubGridScaleFields(const volSymmTensorField& D)
  46 {
  47     muSgs_ = ck_(D)*rho()*sqrt(k_)*delta();
  48     muSgs_.correctBoundaryConditions();
  49
  50     alphaSgs_ = muSgs_/Prt_;
  51     alphaSgs_.correctBoundaryConditions();
  52 }
  53
  54
  55 dimensionedScalar homogeneousDynOneEqEddy::ck_(const volSymmTensorField& D) const
  56 {
  57     volScalarField KK(0.5*(filter_(magSqr(U())) - magSqr(filter_(U()))));
  58
  59     volSymmTensorField LL(dev(filter_(sqr(U())) - (sqr(filter_(U())))));
  60
  61     volSymmTensorField MM
  62     (
  63         delta()*(filter_(sqrt(k_)*D) - 2*sqrt(KK + filter_(k_))*filter_(D))
  64     );
  65
  66     return average(LL && MM)/average(magSqr(MM));
  67 }
  68
  69
  70 dimensionedScalar homogeneousDynOneEqEddy::ce_(const volSymmTensorField& D) const
  71 {
  72     volScalarField KK(0.5*(filter_(magSqr(U())) - magSqr(filter_(U()))));
  73
  74     volScalarField mm
  75     (
  76         pow(KK + filter_(k_), 1.5)/(2*delta()) - filter_(pow(k_, 1.5))/delta()
  77     );
  78
  79     volScalarField ee
  80     (
  81         2*delta()*ck_(D)
  82        *(
  83             filter_(sqrt(k_)*magSqr(D))
  84           - 2*sqrt(KK + filter_(k_))*magSqr(filter_(D))
  85         )
  86     );
  87
  88     return average(ee*mm)/average(mm*mm);
  89 }
  90
  91
  92 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  93
  94 homogeneousDynOneEqEddy::homogeneousDynOneEqEddy
  95 (
  96     const volScalarField& rho,
  97     const volVectorField& U,
  98     const surfaceScalarField& phi,
  99     const basicThermo& thermoPhysicalModel,
 100     const word& turbulenceModelName,
 101     const word& modelName
 102 )
 103 :
 104     LESModel(modelName, rho, U, phi, thermoPhysicalModel, turbulenceModelName),
 105     GenEddyVisc(rho, U, phi, thermoPhysicalModel),
 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     updateSubGridScaleFields(dev(symm(fvc::grad(U))));
 124
 125     printCoeffs();
 126 }
 127
 128
 129 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 130
 131 void homogeneousDynOneEqEddy::correct(const tmp<volTensorField>& tgradU)
 132 {
 133     const volTensorField& gradU = tgradU();
 134
 135     GenEddyVisc::correct(gradU);
 136
 137     volSymmTensorField D(dev(symm(gradU)));
 138     volScalarField divU(fvc::div(phi()/fvc::interpolate(rho())));
 139     volScalarField G(2*muSgs_*(gradU && D));
 140
 141     tmp<fvScalarMatrix> kEqn
 142     (
 143         fvm::ddt(rho(), k_)
 144       + fvm::div(phi(), k_)
 145       - fvm::laplacian(DkEff(), k_)
 146      ==
 147         G
 148       - fvm::SuSp(2.0/3.0*rho()*divU, k_)
 149       - fvm::Sp(ce_(D)*rho()*sqrt(k_)/delta(), k_)
 150     );
 151
 152     kEqn().relax();
 153     kEqn().solve();
 154
 155     bound(k_, kMin_);
 156
 157     updateSubGridScaleFields(D);
 158 }
 159
 160
 161 bool homogeneousDynOneEqEddy::read()
 162 {
 163     if (GenEddyVisc::read())
 164     {
 165         filter_.read(coeffDict());
 166
 167         return true;
 168     }
 169     else
 170     {
 171         return false;
 172     }
 173 }
 174
 175
 176 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 177
 178 } // End namespace LESModels
 179 } // End namespace compressible
 180 } // End namespace Foam
 181
 182 // ************************************************************************* //