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-2024 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 "speciesTable.H"
41 #include "scalarField.H"
42 #include "typeInfo.H"
43 
44 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
45 
46 namespace Foam
47 {
48 
49 // Forward declaration of friend functions and operators
50 
51 class ArrheniusReactionRate;
52 
53 Ostream& operator<<(Ostream&, const ArrheniusReactionRate&);
54 
55 
56 /*---------------------------------------------------------------------------*\
57  Class ArrheniusReactionRate Declaration
58 \*---------------------------------------------------------------------------*/
59 
61 {
62  // Private Data
63 
64  scalar beta_;
65  scalar A_;
66  scalar Ta_;
67 
68 
69 public:
70 
71  // Constructors
72 
73  //- Construct from components
75  (
76  const scalar A,
77  const scalar beta,
78  const scalar Ta
79  );
80 
81  //- Construct from dictionary
83  (
84  const speciesTable& species,
85  const dimensionSet& dims,
86  const dictionary& dict
87  );
88 
89 
90  // Member Functions
91 
92  //- Return the type name
93  static word type()
94  {
95  return "Arrhenius";
96  }
97 
98  //- Pre-evaluation hook
99  inline void preEvaluate() const;
100 
101  //- Post-evaluation hook
102  inline void postEvaluate() const;
103 
104  //- Return the rate
105  inline scalar operator()
106  (
107  const scalar p,
108  const scalar T,
109  const scalarField& c,
110  const label li
111  ) const;
112 
113  //- The derivative of the rate w.r.t. temperature
114  inline scalar ddT
115  (
116  const scalar p,
117  const scalar T,
118  const scalarField& c,
119  const label li
120  ) const;
121 
122  //- Is the rate a function of concentration?
123  inline bool hasDdc() const;
124 
125  //- The derivative of the rate w.r.t. concentration
126  inline void ddc
127  (
128  const scalar p,
129  const scalar T,
130  const scalarField& c,
131  const label li,
133  ) const;
134 
135  //- Write to stream
136  inline void write(Ostream& os) const;
137 
138 
139  // Ostream Operator
140 
141  inline friend Ostream& operator<<
142  (
143  Ostream&,
144  const ArrheniusReactionRate&
145  );
146 };
147 
148 
149 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
150 
151 } // End namespace Foam
152 
153 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
154 
155 #include "ArrheniusReactionRateI.H"
156 
157 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
158 
159 #endif
160 
161 // ************************************************************************* //
static const Foam::dimensionedScalar A("A", Foam::dimPressure, 611.21)
Arrhenius reaction rate given by:
void preEvaluate() const
Pre-evaluation hook.
ArrheniusReactionRate(const scalar A, const scalar beta, const scalar Ta)
Construct from components.
void postEvaluate() const
Post-evaluation hook.
static word type()
Return the type name.
void write(Ostream &os) const
Write to stream.
void ddc(const scalar p, const scalar T, const scalarField &c, const label li, scalarField &ddc) const
The derivative of the rate w.r.t. concentration.
scalar ddT(const scalar p, const scalar T, const scalarField &c, const label li) const
The derivative of the rate w.r.t. temperature.
bool hasDdc() const
Is the rate a function of concentration?
An Ostream is an abstract base class for all output systems (streams, files, token lists,...
Definition: Ostream.H:57
A list of keyword definitions, which are a keyword followed by any number of values (e....
Definition: dictionary.H:162
Dimension set for the base types.
Definition: dimensionSet.H:125
A wordList with hashed indices for faster lookup by name.
A class for handling words, derived from string.
Definition: word.H:62
const dimensionedScalar c
Speed of light in a vacuum.
Namespace for OpenFOAM.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
Ostream & operator<<(Ostream &os, const fvConstraints &constraints)
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
dictionary dict
volScalarField & p
Basic run-time type information using word as the type's name. Used to enhance the standard RTTI to c...