src/thermophysicalModels/radiation/radiationModel/radiationModel/radiationModel.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 "radiationModel.H"
  28 #include "absorptionEmissionModel.H"
  29 #include "scatterModel.H"
  30 #include "fvm.H"
  31
  32 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  33
  34 namespace Foam
  35 {
  36     namespace radiation
  37     {
  38         defineTypeNameAndDebug(radiationModel, 0);
  39         defineRunTimeSelectionTable(radiationModel, dictionary);
  40     }
  41 }
  42
  43
  44 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  45
  46 Foam::radiation::radiationModel::radiationModel(const volScalarField& T)
  47 :
  48     IOdictionary
  49     (
  50         IOobject
  51         (
  52             "radiationProperties",
  53             T.time().constant(),
  54             T.mesh(),
  55             IOobject::MUST_READ,
  56             IOobject::NO_WRITE
  57         )
  58     ),
  59     mesh_(T.mesh()),
  60     time_(T.time()),
  61     T_(T),
  62     radiation_(false),
  63     coeffs_(dictionary::null),
  64     solverFreq_(0),
  65     absorptionEmission_(NULL),
  66     scatter_(NULL)
  67 {}
  68
  69
  70 Foam::radiation::radiationModel::radiationModel
  71 (
  72     const word& type,
  73     const volScalarField& T
  74 )
  75 :
  76     IOdictionary
  77     (
  78         IOobject
  79         (
  80             "radiationProperties",
  81             T.time().constant(),
  82             T.mesh(),
  83             IOobject::MUST_READ,
  84             IOobject::NO_WRITE
  85         )
  86     ),
  87     mesh_(T.mesh()),
  88     time_(T.time()),
  89     T_(T),
  90     radiation_(lookup("radiation")),
  91     coeffs_(subDict(type + "Coeffs")),
  92     solverFreq_(readLabel(lookup("solverFreq"))),
  93     absorptionEmission_(absorptionEmissionModel::New(*this, mesh_)),
  94     scatter_(scatterModel::New(*this, mesh_))
  95 {
  96     solverFreq_ = max(1, solverFreq_);
  97 }
  98
  99
 100 // * * * * * * * * * * * * * * * * Destructor    * * * * * * * * * * * * * * //
 101
 102 Foam::radiation::radiationModel::~radiationModel()
 103 {}
 104
 105
 106 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 107
 108 bool Foam::radiation::radiationModel::read()
 109 {
 110     if (regIOobject::read())
 111     {
 112         lookup("radiation") >> radiation_;
 113         coeffs_ = subDict(type() + "Coeffs");
 114
 115         return true;
 116     }
 117     else
 118     {
 119         return false;
 120     }
 121 }
 122
 123
 124 void Foam::radiation::radiationModel::correct()
 125 {
 126     if (!radiation_)
 127     {
 128         return;
 129     }
 130
 131     if (time_.timeIndex() % solverFreq_ == 0)
 132     {
 133         calculate();
 134     }
 135 }
 136
 137
 138 Foam::tmp<Foam::fvScalarMatrix> Foam::radiation::radiationModel::Sh
 139 (
 140     basicThermo& thermo
 141 ) const
 142 {
 143     volScalarField& h = thermo.h();
 144     const volScalarField cp = thermo.Cp();
 145     const volScalarField T3 = pow3(T_);
 146
 147     return
 148     (
 149         Ru()
 150       - fvm::Sp(4.0*Rp()*T3/cp, h)
 151       - Rp()*T3*(T_ - 4.0*h/cp)
 152     );
 153 }
 154
 155 Foam::tmp<Foam::fvScalarMatrix> Foam::radiation::radiationModel::Shs
 156 (
 157     basicThermo& thermo
 158 ) const
 159 {
 160     volScalarField& hs = thermo.hs();
 161     const volScalarField cp = thermo.Cp();
 162     const volScalarField T3 = pow3(T_);
 163
 164     return
 165     (
 166         Ru()
 167       - fvm::Sp(4.0*Rp()*T3/cp, hs)
 168       - Rp()*T3*(T_ - 4.0*hs/cp)
 169     );
 170 }
 171
 172 // ************************************************************************* //