src/thermophysicalModels/specie/reaction/reactionRate/LandauTellerReactionRate/LandauTellerReactionRateI.H

   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 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  28
  29 namespace Foam
  30 {
  31
  32 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  33
  34 // Construct from components
  35 inline LandauTellerReactionRate::LandauTellerReactionRate
  36 (
  37     const scalar A,
  38     const scalar beta,
  39     const scalar Ta,
  40     const scalar B,
  41     const scalar C
  42 )
  43 :
  44     A_(A),
  45     beta_(beta),
  46     Ta_(Ta),
  47     B_(B),
  48     C_(C)
  49 {}
  50
  51
  52 //- Construct from Istream
  53 inline LandauTellerReactionRate::LandauTellerReactionRate
  54 (
  55     const speciesTable&,
  56     Istream& is
  57 )
  58 :
  59     A_(readScalar(is.readBegin("LandauTellerReactionRate(Istream&)"))),
  60     beta_(readScalar(is)),
  61     Ta_(readScalar(is)),
  62     B_(readScalar(is)),
  63     C_(readScalar(is))
  64 {
  65     is.readEnd("LandauTellerReactionRate(Istream&)");
  66 }
  67
  68
  69 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  70
  71 inline scalar LandauTellerReactionRate::operator()
  72 (
  73     const scalar T,
  74     const scalar,
  75     const scalarField&
  76 ) const
  77 {
  78     scalar lta = A_;
  79
  80     if (mag(beta_) > VSMALL)
  81     {
  82         lta *= pow(T, beta_);
  83     }
  84
  85     scalar expArg = 0.0;
  86
  87     if (mag(Ta_) > VSMALL)
  88     {
  89         expArg -= Ta_/T;
  90     }
  91
  92     if (mag(B_) > VSMALL)
  93     {
  94         expArg += B_/cbrt(T);
  95     }
  96
  97     if (mag(C_) > VSMALL)
  98     {
  99         expArg += C_/pow(T, 2.0/3.0);
 100     }
 101
 102     if (mag(expArg) > VSMALL)
 103     {
 104         lta *= exp(expArg);
 105     }
 106
 107     return lta;
 108 }
 109
 110
 111 inline Ostream& operator<<(Ostream& os, const LandauTellerReactionRate& arr)
 112 {
 113     os  << token::BEGIN_LIST
 114         << arr.A_ << token::SPACE << arr.beta_ << token::SPACE << arr.Ta_
 115         << token::SPACE << arr.B_ << token::SPACE << arr.C_
 116         << token::END_LIST;
 117     return os;
 118 }
 119
 120
 121 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 122
 123 } // End namespace Foam
 124
 125 // ************************************************************************* //