src/turbulenceModels/compressible/RAS/laminar/laminar.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 "laminar.H"
  28 #include "addToRunTimeSelectionTable.H"
  29
  30 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  31
  32 namespace Foam
  33 {
  34 namespace compressible
  35 {
  36 namespace RASModels
  37 {
  38
  39 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  40
  41 defineTypeNameAndDebug(laminar, 0);
  42 addToRunTimeSelectionTable(RASModel, laminar, dictionary);
  43
  44 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  45
  46 laminar::laminar
  47 (
  48     const volScalarField& rho,
  49     const volVectorField& U,
  50     const surfaceScalarField& phi,
  51     const basicThermo& thermophysicalModel
  52 )
  53 :
  54     RASModel(typeName, rho, U, phi, thermophysicalModel)
  55 {}
  56
  57
  58 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  59
  60 tmp<volScalarField> laminar::mut() const
  61 {
  62     return tmp<volScalarField>
  63     (
  64         new volScalarField
  65         (
  66             IOobject
  67             (
  68                 "mut",
  69                 runTime_.timeName(),
  70                 mesh_,
  71                 IOobject::NO_READ,
  72                 IOobject::NO_WRITE
  73             ),
  74             mesh_,
  75             dimensionedScalar("mut", mu().dimensions(), 0.0)
  76         )
  77     );
  78 }
  79
  80
  81 tmp<volScalarField> laminar::k() const
  82 {
  83     return tmp<volScalarField>
  84     (
  85         new volScalarField
  86         (
  87             IOobject
  88             (
  89                 "k",
  90                 runTime_.timeName(),
  91                 mesh_,
  92                 IOobject::NO_READ,
  93                 IOobject::NO_WRITE
  94             ),
  95             mesh_,
  96             dimensionedScalar("k", sqr(U_.dimensions()), 0.0)
  97         )
  98     );
  99 }
 100
 101
 102 tmp<volScalarField> laminar::epsilon() const
 103 {
 104     return tmp<volScalarField>
 105     (
 106         new volScalarField
 107         (
 108             IOobject
 109             (
 110                 "epsilon",
 111                 runTime_.timeName(),
 112                 mesh_,
 113                 IOobject::NO_READ,
 114                 IOobject::NO_WRITE
 115             ),
 116             mesh_,
 117             dimensionedScalar
 118             (
 119                 "epsilon", sqr(U_.dimensions())/dimTime, 0.0
 120             )
 121         )
 122     );
 123 }
 124
 125
 126 tmp<volSymmTensorField> laminar::R() const
 127 {
 128     return tmp<volSymmTensorField>
 129     (
 130         new volSymmTensorField
 131         (
 132             IOobject
 133             (
 134                 "R",
 135                 runTime_.timeName(),
 136                 mesh_,
 137                 IOobject::NO_READ,
 138                 IOobject::NO_WRITE
 139             ),
 140             mesh_,
 141             dimensionedSymmTensor
 142             (
 143                 "R", sqr(U_.dimensions()), symmTensor::zero
 144             )
 145         )
 146     );
 147 }
 148
 149
 150 tmp<volSymmTensorField> laminar::devRhoReff() const
 151 {
 152     return tmp<volSymmTensorField>
 153     (
 154         new volSymmTensorField
 155         (
 156             IOobject
 157             (
 158                 "devRhoReff",
 159                 runTime_.timeName(),
 160                 mesh_,
 161                 IOobject::NO_READ,
 162                 IOobject::NO_WRITE
 163             ),
 164            -mu()*dev(twoSymm(fvc::grad(U_)))
 165         )
 166     );
 167 }
 168
 169
 170 tmp<fvVectorMatrix> laminar::divDevRhoReff(volVectorField& U) const
 171 {
 172     return
 173     (
 174       - fvm::laplacian(muEff(), U)
 175       - fvc::div(muEff()*dev2(fvc::grad(U)().T()))
 176     );
 177 }
 178
 179
 180 bool laminar::read()
 181 {
 182     return RASModel::read();
 183 }
 184
 185
 186 void laminar::correct()
 187 {}
 188
 189
 190 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 191
 192 } // End namespace RASModels
 193 } // End namespace compressible
 194 } // End namespace Foam
 195
 196 // ************************************************************************* //