![]() |
|
00001 /*---------------------------------------------------------------------------*\ 00002 ========= | 00003 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox 00004 \\ / O peration | 00005 \\ / A nd | Copyright (C) 1991-2005 OpenCFD Ltd. 00006 \\/ M anipulation | 00007 ------------------------------------------------------------------------------- 00008 License 00009 This file is part of OpenFOAM. 00010 00011 OpenFOAM is free software; you can redistribute it and/or modify it 00012 under the terms of the GNU General Public License as published by the 00013 Free Software Foundation; either version 2 of the License, or (at your 00014 option) any later version. 00015 00016 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT 00017 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or 00018 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License 00019 for more details. 00020 00021 You should have received a copy of the GNU General Public License 00022 along with OpenFOAM; if not, write to the Free Software Foundation, 00023 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA 00024 00025 Description 00026 00027 \*---------------------------------------------------------------------------*/ 00028 00029 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00030 00031 namespace Foam 00032 { 00033 00034 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * // 00035 00036 // Construct from components 00037 inline ArrheniusReactionRate::ArrheniusReactionRate 00038 ( 00039 const scalar A, 00040 const scalar beta, 00041 const scalar Ta 00042 ) 00043 : 00044 A_(A), 00045 beta_(beta), 00046 Ta_(Ta) 00047 {} 00048 00049 00050 //- Construct from Istream 00051 inline ArrheniusReactionRate::ArrheniusReactionRate 00052 ( 00053 const speciesTable&, 00054 Istream& is 00055 ) 00056 : 00057 A_(readScalar(is.readBegin("ArrheniusReactionRate(Istream&)"))), 00058 beta_(readScalar(is)), 00059 Ta_(readScalar(is)) 00060 { 00061 is.readEnd("ArrheniusReactionRate(Istream&)"); 00062 } 00063 00064 00065 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // 00066 00067 inline scalar ArrheniusReactionRate::operator() 00068 ( 00069 const scalar T, 00070 const scalar, 00071 const scalarField& 00072 ) const 00073 { 00074 scalar ak = A_; 00075 00076 if (mag(beta_) > VSMALL) 00077 { 00078 ak *= pow(T, beta_); 00079 } 00080 00081 if (mag(Ta_) > VSMALL) 00082 { 00083 ak *= exp(-Ta_/T); 00084 } 00085 00086 return ak; 00087 } 00088 00089 00090 inline Ostream& operator<<(Ostream& os, const ArrheniusReactionRate& arr) 00091 { 00092 os << token::BEGIN_LIST 00093 << arr.A_ << token::SPACE << arr.beta_ << token::SPACE << arr.Ta_ 00094 << token::END_LIST; 00095 return os; 00096 } 00097 00098 00099 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 00100 00101 } // End namespace Foam 00102 00103 // ************************************************************************* //