src/thermophysicalModels/radiationModels/radiationModel/fvDOM/absorptionCoeffs/absorptionCoeffs.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2008-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 "absorptionCoeffs.H"
  27 #include "IOstreams.H"
  28
  29 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  30
  31 Foam::radiation::absorptionCoeffs::absorptionCoeffs(Istream& is)
  32 :
  33    Tcommon_(readScalar(is)),
  34    Tlow_(readScalar(is)),
  35    Thigh_(readScalar(is)),
  36    invTemp_(readBool(is))
  37 {
  38     for (label coefLabel=0; absorptionCoeffs::nCoeffs_; coefLabel++)
  39     {
  40         is >> highACoeffs_[coefLabel];
  41     }
  42
  43     for (label coefLabel=0; absorptionCoeffs::nCoeffs_; coefLabel++)
  44     {
  45         is >> lowACoeffs_[coefLabel];
  46     }
  47 }
  48
  49
  50 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * //
  51
  52 Foam::radiation::absorptionCoeffs::~absorptionCoeffs()
  53 {}
  54
  55
  56 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  57
  58 void Foam::radiation::absorptionCoeffs::checkT(const scalar T) const
  59 {
  60     if (T < Tlow_ || T > Thigh_)
  61     {
  62         WarningIn
  63         (
  64             "absorptionCoeffs::checkT(const scalar T) const"
  65         )   << "usinf absCoeff out of temperature range:" << nl
  66             << "    " << Tlow_ << " -> " << Thigh_ << ";  T = " << T
  67             << nl << endl;
  68     }
  69 }
  70
  71
  72 const Foam::radiation::absorptionCoeffs::coeffArray&
  73 Foam::radiation::absorptionCoeffs::coeffs
  74 (
  75     const scalar T
  76 ) const
  77 {
  78     checkT(T);
  79
  80     if (T < Tcommon_)
  81     {
  82         return lowACoeffs_;
  83     }
  84     else
  85     {
  86         return highACoeffs_;
  87     }
  88 }
  89
  90
  91 void Foam::radiation::absorptionCoeffs::initialise(Istream&)
  92 {
  93     absorptionCoeffs(Istream);
  94 }
  95
  96
  97 void Foam::radiation::absorptionCoeffs::initialise(const dictionary& dict)
  98 {
  99     dict.lookup("Tcommon") >> Tcommon_;
 100     dict.lookup("Tlow") >> Tlow_;
 101     dict.lookup("Thigh") >> Thigh_;
 102     dict.lookup("invTemp") >> invTemp_;
 103
 104     dict.lookup("loTcoeffs") >> lowACoeffs_;
 105     dict.lookup("hiTcoeffs") >> highACoeffs_;
 106 }
 107
 108
 109 // ************************************************************************* //