src/thermophysicalModels/specie/reaction/reactionRate/FallOffReactionRate/FallOffReactionRateI.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-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 template<class ReactionRate, class FallOffFunction>
  29 inline Foam::FallOffReactionRate<ReactionRate, FallOffFunction>::
  30 FallOffReactionRate
  31 (
  32     const ReactionRate& k0,
  33     const ReactionRate& kInf,
  34     const FallOffFunction& F,
  35     const thirdBodyEfficiencies& tbes
  36 )
  37 :
  38     k0_(k0),
  39     kInf_(kInf),
  40     F_(F),
  41     thirdBodyEfficiencies_(tbes)
  42 {}
  43
  44
  45 template<class ReactionRate, class FallOffFunction>
  46 inline Foam::FallOffReactionRate<ReactionRate, FallOffFunction>::
  47 FallOffReactionRate
  48 (
  49     const speciesTable& species,
  50     Istream& is
  51 )
  52 :
  53     k0_(species, is.readBegin("FallOffReactionRate(Istream&)")),
  54     kInf_(species, is),
  55     F_(is),
  56     thirdBodyEfficiencies_(species, is)
  57 {
  58     is.readEnd("FallOffReactionRate(Istream&)");
  59 }
  60
  61
  62 template<class ReactionRate, class FallOffFunction>
  63 inline Foam::FallOffReactionRate<ReactionRate, FallOffFunction>::
  64 FallOffReactionRate
  65 (
  66     const speciesTable& species,
  67     const dictionary& dict
  68 )
  69 :
  70     k0_(species, dict),
  71     kInf_(species, dict),
  72     F_(dict),
  73     thirdBodyEfficiencies_(species, dict)
  74 {}
  75
  76
  77 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  78
  79 template<class ReactionRate, class FallOffFunction>
  80 inline Foam::scalar
  81 Foam::FallOffReactionRate<ReactionRate, FallOffFunction>::operator()
  82 (
  83     const scalar T,
  84     const scalar p,
  85     const scalarField& c
  86 ) const
  87 {
  88     scalar k0 = k0_(T, p, c);
  89     scalar kInf = kInf_(T, p, c);
  90
  91     scalar Pr = k0*thirdBodyEfficiencies_.M(c)/kInf;
  92
  93     return kInf*(Pr/(1 + Pr))*F_(T, Pr);
  94 }
  95
  96
  97 template<class ReactionRate, class FallOffFunction>
  98 inline void Foam::FallOffReactionRate<ReactionRate, FallOffFunction>::write
  99 (
 100     Ostream& os
 101 ) const
 102 {
 103     k0_.write(os);
 104     kInf_.write(os);
 105     F_.write(os);
 106     thirdBodyEfficiencies_.write(os);
 107 }
 108
 109
 110 template<class ReactionRate, class FallOffFunction>
 111 inline Foam::Ostream& Foam::operator<<
 112 (
 113     Ostream& os,
 114     const FallOffReactionRate<ReactionRate, FallOffFunction>& forr
 115 )
 116 {
 117     os  << token::BEGIN_LIST
 118         << forr.k0_ << token::SPACE
 119         << forr.kInf_ << token::SPACE
 120         << forr.F_ << token::SPACE
 121         << forr.thirdBodyEfficiencies_
 122         << token::END_LIST;
 123     return os;
 124 }
 125
 126
 127 // ************************************************************************* //