src/thermophysicalModels/solid/reaction/reactionRate/solidArrheniusReactionRate/solidArrheniusReactionRateI.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2010-2010 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 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  27
  28 inline Foam::solidArrheniusReactionRate::solidArrheniusReactionRate
  29 (
  30     const scalar A,
  31     const scalar Ta,
  32     const scalar Tcrit
  33 )
  34 :
  35     A_(A),
  36     Ta_(Ta),
  37     Tcrit_(Tcrit)
  38 {}
  39
  40
  41 inline Foam::solidArrheniusReactionRate::solidArrheniusReactionRate
  42 (
  43     const speciesTable&,
  44     Istream& is
  45 )
  46 :
  47     A_(readScalar(is.readBegin("solidArrheniusReaction(Istream&)"))),
  48     Ta_(readScalar(is)),
  49     Tcrit_(readScalar(is))
  50 {}
  51
  52
  53 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  54
  55 inline Foam::scalar Foam::solidArrheniusReactionRate::operator()
  56 (
  57     const scalar T,
  58     const scalar,
  59     const scalarField&
  60 ) const
  61 {
  62     scalar ak = A_;
  63
  64     if (T < Tcrit_)
  65     {
  66         ak *= 0.0;
  67     }
  68     else
  69     {
  70         ak *= exp(-Ta_/T);
  71     }
  72
  73     return ak;
  74 }
  75
  76
  77 inline Foam::Ostream& Foam::operator<<
  78 (
  79     Ostream& os,
  80     const solidArrheniusReactionRate& arr
  81 )
  82 {
  83     os  << token::BEGIN_LIST
  84         << arr.A_ << token::SPACE << arr.Ta_ << token::SPACE << arr.Tcrit_
  85         << token::END_LIST;
  86     return os;
  87 }
  88
  89
  90 // ************************************************************************* //