src/thermophysicalModels/radiation/radiationModel/P1/P1.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 "P1.H"
  28 #include "addToRunTimeSelectionTable.H"
  29 #include "fvm.H"
  30
  31 #include "absorptionEmissionModel.H"
  32 #include "scatterModel.H"
  33 #include "mathematicalConstants.H"
  34 #include "radiationConstants.H"
  35
  36 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  37
  38 namespace Foam
  39 {
  40     namespace radiation
  41     {
  42         defineTypeNameAndDebug(P1, 0);
  43
  44         addToRunTimeSelectionTable
  45         (
  46             radiationModel,
  47             P1,
  48             dictionary
  49         );
  50     }
  51 }
  52
  53
  54 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  55
  56 Foam::radiation::P1::P1(const volScalarField& T)
  57 :
  58     radiationModel(typeName, T),
  59     G_
  60     (
  61         IOobject
  62         (
  63             "G",
  64             mesh_.time().timeName(),
  65             mesh_,
  66             IOobject::MUST_READ,
  67             IOobject::AUTO_WRITE
  68         ),
  69         mesh_
  70     ),
  71     a_
  72     (
  73         IOobject
  74         (
  75             "a",
  76             mesh_.time().timeName(),
  77             mesh_,
  78             IOobject::NO_READ,
  79             IOobject::NO_WRITE
  80         ),
  81         mesh_,
  82         dimensionedScalar("a", dimless/dimLength, 0.0)
  83     ),
  84     e_
  85     (
  86         IOobject
  87         (
  88             "e",
  89             mesh_.time().timeName(),
  90             mesh_,
  91             IOobject::NO_READ,
  92             IOobject::NO_WRITE
  93         ),
  94         mesh_,
  95         dimensionedScalar("a", dimless/dimLength, 0.0)
  96     ),
  97     E_
  98     (
  99         IOobject
 100         (
 101             "E",
 102             mesh_.time().timeName(),
 103             mesh_,
 104             IOobject::NO_READ,
 105             IOobject::NO_WRITE
 106         ),
 107         mesh_,
 108         dimensionedScalar("E", dimMass/dimLength/pow3(dimTime), 0.0)
 109     )
 110 {}
 111
 112
 113 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 114
 115 Foam::radiation::P1::~P1()
 116 {}
 117
 118
 119 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 120
 121 bool Foam::radiation::P1::read()
 122 {
 123     if (radiationModel::read())
 124     {
 125         // nothing to read
 126
 127         return true;
 128     }
 129     else
 130     {
 131         return false;
 132     }
 133 }
 134
 135
 136 void Foam::radiation::P1::calculate()
 137 {
 138     a_ = absorptionEmission_->a();
 139     e_ = absorptionEmission_->e();
 140     E_ = absorptionEmission_->E();
 141     const volScalarField sigmaEff = scatter_->sigmaEff();
 142
 143     // Construct diffusion
 144     const volScalarField gamma
 145     (
 146         IOobject
 147         (
 148             "gammaRad",
 149             G_.mesh().time().timeName(),
 150             G_.mesh(),
 151             IOobject::NO_READ,
 152             IOobject::NO_WRITE
 153         ),
 154         1.0/(3.0*a_ + sigmaEff)
 155     );
 156
 157     // Solve G transport equation
 158     solve
 159     (
 160         fvm::laplacian(gamma, G_)
 161       - fvm::Sp(a_, G_)
 162      ==
 163       - 4.0*(e_*radiation::sigmaSB*pow4(T_) + E_)
 164     );
 165 }
 166
 167
 168 Foam::tmp<Foam::volScalarField> Foam::radiation::P1::Rp() const
 169 {
 170     return tmp<volScalarField>
 171     (
 172         new volScalarField
 173         (
 174             IOobject
 175             (
 176                 "Rp",
 177                 mesh_.time().timeName(),
 178                 mesh_,
 179                 IOobject::NO_READ,
 180                 IOobject::NO_WRITE,
 181                 false
 182             ),
 183             4.0*absorptionEmission_->eCont()*radiation::sigmaSB
 184         )
 185     );
 186 }
 187
 188
 189 Foam::tmp<Foam::DimensionedField<Foam::scalar, Foam::volMesh> >
 190 Foam::radiation::P1::Ru() const
 191 {
 192     const DimensionedField<scalar, volMesh>& G =
 193         G_.dimensionedInternalField();
 194     const DimensionedField<scalar, volMesh> E =
 195         absorptionEmission_->ECont()().dimensionedInternalField();
 196     const DimensionedField<scalar, volMesh> a =
 197         absorptionEmission_->aCont()().dimensionedInternalField();
 198
 199     return  a*G - 4.0*E;
 200 }
 201
 202
 203 // ************************************************************************* //