src/thermophysicalModels/specie/reaction/reactionRate/fallOffFunctions/TroeFallOffFunction/TroeFallOffFunctionI.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 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  28
  29 // Construct from components
  30 inline Foam::TroeFallOffFunction::TroeFallOffFunction
  31 (
  32     const scalar alpha,
  33     const scalar Tsss,
  34     const scalar Ts,
  35     const scalar Tss
  36 )
  37 :
  38     alpha_(alpha),
  39     Tsss_(Tsss),
  40     Ts_(Ts),
  41     Tss_(Tss)
  42 {}
  43
  44
  45 // Construct from Istream
  46 inline Foam::TroeFallOffFunction::TroeFallOffFunction(Istream& is)
  47 :
  48     alpha_(readScalar(is.readBegin("TroeFallOffFunction(Istream&)"))),
  49     Tsss_(readScalar(is)),
  50     Ts_(readScalar(is)),
  51     Tss_(readScalar(is))
  52 {
  53     is.readEnd("TroeFallOffFunction(Istream&)");
  54 }
  55
  56
  57 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  58
  59 inline Foam::scalar Foam::TroeFallOffFunction::operator()
  60 (
  61     const scalar T,
  62     const scalar Pr
  63 ) const
  64 {
  65     scalar logFcent = log10
  66     (
  67         max
  68         (
  69             (1 - alpha_)*exp(-T/Tsss_) + alpha_*exp(-T/Ts_) + exp(-Tss_/T),
  70             SMALL
  71         )
  72     );
  73
  74     scalar c = -0.4 - 0.67*logFcent;
  75     static const scalar d = 0.14;
  76     scalar n = 0.75 - 1.27*logFcent;
  77
  78     scalar logPr = log10(max(Pr, SMALL));
  79     return pow(10.0, logFcent/(1.0 + sqr((logPr + c)/(n - d*(logPr + c)))));
  80 }
  81
  82
  83 // * * * * * * * * * * * * * * * Ostream Operator  * * * * * * * * * * * * * //
  84
  85 inline Foam::Ostream& Foam::operator<<
  86 (
  87     Foam::Ostream& os,
  88     const Foam::TroeFallOffFunction& tfof
  89 )
  90 {
  91     os  << token::BEGIN_LIST
  92         << tfof.alpha_
  93         << token::SPACE << tfof.Tsss_
  94         << token::SPACE << tfof.Ts_
  95         << token::SPACE << tfof.Tss_
  96         << token::END_LIST;
  97
  98     return os;
  99 }
 100
 101
 102 // ************************************************************************* //