ENH: RASModel.C: clipping input to log
[OpenFOAM-1.7.x.git] / src / thermophysicalModels / thermophysicalFunctions / NSRDSfunctions / NSRDSfunc1 / NSRDSfunc1.H
blob291bdd5b0bc55fdd4cf354a070cb0f67dc3da6ea
1 /*---------------------------------------------------------------------------*\
2   =========                 |
3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
4    \\    /   O peration     |
5     \\  /    A nd           | Copyright (C) 1991-2010 OpenCFD Ltd.
6      \\/     M anipulation  |
7 -------------------------------------------------------------------------------
8 License
9     This file is part of OpenFOAM.
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.
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.
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/>.
24 Class
25     Foam::NSRDSfunc1
27 Description
28     NSRDS function number 101
30     Source:
31     @verbatim
32                       NSRDS - AICHE
33                  Data Compilation Tables
34                     of Properties of
35                      Pure Compounds
37         Design Institute for Physical Property Data
38           American Institute of Chemical Engineers
39                   345 East 47th Street
40                 New York, New York 10017
42          National Standard Reference Data System
43          American Institute of Chemical Engineers
45           T.E. Daubert       -       R.P. Danner
47             Department of Chemical Engineering
48             The Pennsylvania State University
49                 University Park, PA 16802
50     @endverbatim
52 \*---------------------------------------------------------------------------*/
54 #ifndef NSRDSfunc1_H
55 #define NSRDSfunc1_H
57 #include "thermophysicalFunction.H"
59 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
61 namespace Foam
64 /*---------------------------------------------------------------------------*\
65                            Class NSRDSfunc1 Declaration
66 \*---------------------------------------------------------------------------*/
68 class NSRDSfunc1
70     public thermophysicalFunction
72     // Private data
74         // NSRDS function 101 coefficients
75         scalar a_, b_, c_, d_, e_;
78 public:
80     //- Runtime type information
81     TypeName("NSRDSfunc1");
84     // Constructors
86         //- Construct from components
87         NSRDSfunc1(scalar a, scalar b, scalar c, scalar d, scalar e)
88         :
89             a_(a),
90             b_(b),
91             c_(c),
92             d_(d),
93             e_(e)
94         {}
96         //- Construct from Istream
97         NSRDSfunc1(Istream& is)
98         :
99             a_(readScalar(is)),
100             b_(readScalar(is)),
101             c_(readScalar(is)),
102             d_(readScalar(is)),
103             e_(readScalar(is))
104         {}
107     // Member Functions
109         //- Evaluate the function and return the result
110         scalar f(scalar, scalar T) const
111         {
112             return exp(a_ + b_/T + c_*log(T) + d_*pow(T, e_));
113         }
116         //- Write the function coefficients
117         void writeData(Ostream& os) const
118         {
119             os  << a_ << token::SPACE
120                 << b_ << token::SPACE
121                 << c_ << token::SPACE
122                 << d_ << token::SPACE
123                 << e_;
124         }
127     // Ostream Operator
129         friend Ostream& operator<<(Ostream& os, const NSRDSfunc1& f)
130         {
131             f.writeData(os);
132             return os;
133         }
137 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
139 } // End namespace Foam
141 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
143 #endif
145 // ************************************************************************* //