applications/solvers/combustion/PDRFoam/PDRModels/XiEqModels/basicXiSubXiEq/basicXiSubXiEq.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-2011 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
  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 "basicXiSubXiEq.H"
  27 #include "addToRunTimeSelectionTable.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 namespace Foam
  32 {
  33 namespace XiEqModels
  34 {
  35     defineTypeNameAndDebug(basicSubGrid, 0);
  36     addToRunTimeSelectionTable(XiEqModel, basicSubGrid, dictionary);
  37 }
  38 }
  39
  40
  41 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  42
  43 Foam::XiEqModels::basicSubGrid::basicSubGrid
  44 (
  45     const dictionary& XiEqProperties,
  46     const hhuCombustionThermo& thermo,
  47     const compressible::RASModel& turbulence,
  48     const volScalarField& Su
  49 )
  50 :
  51     XiEqModel(XiEqProperties, thermo, turbulence, Su),
  52
  53     B_
  54     (
  55         IOobject
  56         (
  57             "B",
  58             Su.mesh().facesInstance(),
  59             Su.mesh(),
  60             IOobject::MUST_READ,
  61             IOobject::NO_WRITE
  62         ),
  63         Su.mesh()
  64     ),
  65
  66     XiEqModel_(XiEqModel::New(XiEqModelCoeffs_, thermo, turbulence, Su))
  67 {}
  68
  69
  70 // * * * * * * * * * * * * * * * * Destructors * * * * * * * * * * * * * * * //
  71
  72 Foam::XiEqModels::basicSubGrid::~basicSubGrid()
  73 {}
  74
  75
  76 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
  77
  78 Foam::tmp<Foam::volScalarField> Foam::XiEqModels::basicSubGrid::XiEq() const
  79 {
  80     const fvMesh& mesh = Su_.mesh();
  81     const volVectorField& U = mesh.lookupObject<volVectorField>("U");
  82
  83     const volScalarField& Nv = mesh.lookupObject<volScalarField>("Nv");
  84     const volSymmTensorField& nsv =
  85         mesh.lookupObject<volSymmTensorField>("nsv");
  86
  87     volScalarField magU(mag(U));
  88     volVectorField Uhat
  89     (
  90         U/(mag(U) + dimensionedScalar("Usmall", U.dimensions(), 1e-4))
  91     );
  92
  93     const scalarField Cw = pow(mesh.V(), 2.0/3.0);
  94
  95     tmp<volScalarField> tN
  96     (
  97         new volScalarField
  98         (
  99             IOobject
 100             (
 101                 "tN",
 102                 mesh.time().constant(),
 103                 mesh,
 104                 IOobject::NO_READ,
 105                 IOobject::NO_WRITE
 106             ),
 107             mesh,
 108             dimensionedScalar("zero", Nv.dimensions(), 0.0),
 109             zeroGradientFvPatchVectorField::typeName
 110         )
 111     );
 112
 113     volScalarField& N = tN();
 114
 115     N.internalField() = Nv.internalField()*Cw;
 116
 117     tmp<volSymmTensorField> tns
 118     (
 119         new volSymmTensorField
 120         (
 121             IOobject
 122             (
 123                 "tns",
 124                 U.mesh().time().timeName(),
 125                 U.mesh(),
 126                 IOobject::NO_READ,
 127                 IOobject::NO_WRITE
 128             ),
 129             U.mesh(),
 130             dimensionedSymmTensor
 131             (
 132                 "zero",
 133                 nsv.dimensions(),
 134                 pTraits<symmTensor>::zero
 135             ),
 136              zeroGradientFvPatchSymmTensorField::typeName
 137         )
 138     );
 139
 140     volSymmTensorField& ns = tns();
 141
 142     ns.internalField() = nsv.internalField()*Cw;
 143
 144     volScalarField n(max(N - (Uhat & ns & Uhat), scalar(1e-4)));
 145
 146     volScalarField b((Uhat & B_ & Uhat)/sqrt(n));
 147
 148     volScalarField up(sqrt((2.0/3.0)*turbulence_.k()));
 149
 150     volScalarField XiSubEq
 151     (
 152         scalar(1)
 153       + max(2.2*sqrt(b), min(0.34*magU/up*sqrt(b), scalar(1.6)))
 154       * min(n, scalar(1))
 155     );
 156
 157     return (XiSubEq*XiEqModel_->XiEq());
 158 }
 159
 160
 161 bool Foam::XiEqModels::basicSubGrid::read(const dictionary& XiEqProperties)
 162 {
 163     XiEqModel::read(XiEqProperties);
 164
 165     return XiEqModel_->read(XiEqModelCoeffs_);
 166 }
 167
 168
 169 // ************************************************************************* //