00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037 #ifndef ReversibleReaction_H
00038 #define ReversibleReaction_H
00039
00040 #include "Reaction.H"
00041
00042
00043
00044 namespace Foam
00045 {
00046
00047
00048
00049
00050
00051 template<class ReactionThermo, class ReactionRate>
00052 class ReversibleReaction
00053 :
00054 public Reaction<ReactionThermo>
00055 {
00056
00057
00058 ReactionRate k_;
00059
00060
00061
00062
00063
00064 void operator=(const ReversibleReaction<ReactionThermo, ReactionRate>&);
00065
00066
00067 public:
00068
00069
00070 TypeName("reversible");
00071
00072
00073
00074
00075
00076 ReversibleReaction
00077 (
00078 const Reaction<ReactionThermo>& reaction,
00079 const ReactionRate& k
00080 );
00081
00082
00083 ReversibleReaction
00084 (
00085 const ReversibleReaction<ReactionThermo, ReactionRate>&,
00086 const speciesTable& species
00087 );
00088
00089
00090 ReversibleReaction
00091 (
00092 const speciesTable& species,
00093 const HashPtrTable<ReactionThermo>& thermoDatabase,
00094 Istream& is
00095 );
00096
00097
00098 virtual autoPtr<Reaction<ReactionThermo> > clone() const
00099 {
00100 return autoPtr<Reaction<ReactionThermo> >
00101 (
00102 new ReversibleReaction<ReactionThermo, ReactionRate>(*this)
00103 );
00104 }
00105
00106
00107 virtual autoPtr<Reaction<ReactionThermo> > clone
00108 (
00109 const speciesTable& species
00110 ) const
00111 {
00112 return autoPtr<Reaction<ReactionThermo> >
00113 (
00114 new ReversibleReaction<ReactionThermo, ReactionRate>
00115 (
00116 *this,
00117 species
00118 )
00119 );
00120 }
00121
00122
00123
00124
00125 virtual ~ReversibleReaction()
00126 {}
00127
00128
00129
00130
00131
00132
00133
00134 virtual scalar kf
00135 (
00136 const scalar T,
00137 const scalar p,
00138 const scalarField& c
00139 ) const;
00140
00141
00142 virtual scalar kr
00143 (
00144 const scalar kfwd,
00145 const scalar T,
00146 const scalar p,
00147 const scalarField& c
00148 ) const;
00149
00150
00151
00152
00153 virtual scalar kr
00154 (
00155 const scalar T,
00156 const scalar p,
00157 const scalarField& c
00158 ) const;
00159
00160
00161
00162 virtual void write(Ostream&) const;
00163 };
00164
00165
00166
00167
00168 }
00169
00170
00171
00172 #ifdef NoRepository
00173 # include "ReversibleReaction.C"
00174 #endif
00175
00176
00177
00178 #endif
00179
00180