Go to the documentation of this file.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 inline Foam::SRIFallOffFunction::SRIFallOffFunction
00030 (
00031 const scalar a,
00032 const scalar b,
00033 const scalar c,
00034 const scalar d,
00035 const scalar e
00036 )
00037 :
00038 a_(a),
00039 b_(b),
00040 c_(c),
00041 d_(d),
00042 e_(e)
00043 {}
00044
00045
00046
00047 inline Foam::SRIFallOffFunction::SRIFallOffFunction(Istream& is)
00048 :
00049 a_(readScalar(is.readBegin("SRIFallOffFunction(Istream&)"))),
00050 b_(readScalar(is)),
00051 c_(readScalar(is)),
00052 d_(readScalar(is)),
00053 e_(readScalar(is))
00054 {
00055 is.readEnd("SRIFallOffFunction(Istream&)");
00056 }
00057
00058
00059
00060
00061 inline Foam::scalar Foam::SRIFallOffFunction::operator()
00062 (
00063 const scalar T,
00064 const scalar Pr
00065 ) const
00066 {
00067 scalar X = 1.0/(1.0 + sqr(log10(max(Pr, SMALL))));
00068 return d_*pow(a_*exp(-b_/T) + exp(-T/c_), X)*pow(T, e_);
00069 }
00070
00071
00072
00073
00074 inline Foam::Ostream& Foam::operator<<
00075 (
00076 Foam::Ostream& os,
00077 const Foam::SRIFallOffFunction& srifof
00078 )
00079 {
00080 os << token::BEGIN_LIST
00081 << srifof.a_
00082 << token::SPACE << srifof.b_
00083 << token::SPACE << srifof.c_
00084 << token::SPACE << srifof.d_
00085 << token::SPACE << srifof.e_
00086 << token::END_LIST;
00087
00088 return os;
00089 }
00090
00091
00092