ArrheniusReactionRate.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2018 OpenFOAM Foundation
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 Class
25  Foam::ArrheniusReactionRate
26 
27 Description
28  Arrhenius reaction rate given by:
29 
30  k = A * T^beta * exp(-Ta/T)
31 
32 SourceFiles
33  ArrheniusReactionRateI.H
34 
35 \*---------------------------------------------------------------------------*/
36 
37 #ifndef ArrheniusReactionRate_H
38 #define ArrheniusReactionRate_H
39 
40 #include "scalarField.H"
41 #include "typeInfo.H"
42 
43 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
44 
45 namespace Foam
46 {
47 
48 // Forward declaration of friend functions and operators
49 
50 class ArrheniusReactionRate;
51 
52 Ostream& operator<<(Ostream&, const ArrheniusReactionRate&);
53 
54 
55 /*---------------------------------------------------------------------------*\
56  Class ArrheniusReactionRate Declaration
57 \*---------------------------------------------------------------------------*/
58 
60 {
61  // Private data
62 
63  scalar A_;
64  scalar beta_;
65  scalar Ta_;
66 
67 
68 public:
69 
70  // Constructors
71 
72  //- Construct from components
74  (
75  const scalar A,
76  const scalar beta,
77  const scalar Ta
78  );
79 
80  //- Construct from dictionary
82  (
83  const speciesTable& species,
84  const dictionary& dict
85  );
86 
87 
88  // Member Functions
89 
90  //- Return the type name
91  static word type()
92  {
93  return "Arrhenius";
94  }
95 
96  inline scalar operator()
97  (
98  const scalar p,
99  const scalar T,
100  const scalarField& c
101  ) const;
102 
103  inline scalar ddT
104  (
105  const scalar p,
106  const scalar T,
107  const scalarField& c
108  ) const;
109 
110  //- Third-body efficiencies (beta = 1-alpha)
111  // non-empty only for third-body reactions
112  // with enhanced molecularity (alpha != 1)
113  inline const List<Tuple2<label, scalar>>& beta() const;
114 
115  //- Species concentration derivative of the pressure dependent term
116  // By default this value is 1 as it multiplies the third-body term
117  inline void dcidc
118  (
119  const scalar p,
120  const scalar T,
121  const scalarField& c,
123  ) const;
124 
125  //- Temperature derivative of the pressure dependent term
126  // By default this value is 0 since ddT of molecularity is approx.0
127  inline scalar dcidT
128  (
129  const scalar p,
130  const scalar T,
131  const scalarField& c
132  ) const;
133 
134  //- Write to stream
135  inline void write(Ostream& os) const;
136 
137 
138  // Ostream Operator
139 
140  inline friend Ostream& operator<<
141  (
142  Ostream&,
143  const ArrheniusReactionRate&
144  );
145 };
146 
147 
148 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
149 
150 } // End namespace Foam
151 
152 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
153 
154 #include "ArrheniusReactionRateI.H"
155 
156 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
157 
158 #endif
159 
160 // ************************************************************************* //
dictionary dict
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:137
A 1D array of objects of type <T>, where the size of the vector is known and used for subscript bound...
Definition: HashTable.H:60
ArrheniusReactionRate(const scalar A, const scalar beta, const scalar Ta)
Construct from components.
A class for handling words, derived from string.
Definition: word.H:59
scalar ddT(const scalar p, const scalar T, const scalarField &c) const
An Ostream is an abstract base class for all output systems (streams, files, token lists...
Definition: Ostream.H:53
static word type()
Return the type name.
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
A wordList with hashed indices for faster lookup by name.
const dimensionedScalar c
Speed of light in a vacuum.
Ostream & operator<<(Ostream &, const ensightPart &)
Arrhenius reaction rate given by:
void dcidc(const scalar p, const scalar T, const scalarField &c, scalarField &dcidc) const
Species concentration derivative of the pressure dependent term.
volScalarField & p
void write(Ostream &os) const
Write to stream.
scalar dcidT(const scalar p, const scalar T, const scalarField &c) const
Temperature derivative of the pressure dependent term.
const List< Tuple2< label, scalar > > & beta() const
Third-body efficiencies (beta = 1-alpha)
Namespace for OpenFOAM.