1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 2004-2011 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
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
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/>.
28 NSRDS function number 114
33 Data Compilation Tables
37 Design Institute for Physical Property Data
38 American Institute of Chemical Engineers
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
52 \*---------------------------------------------------------------------------*/
57 #include "thermophysicalFunction.H"
59 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
64 /*---------------------------------------------------------------------------*\
65 Class NSRDSfunc14 Declaration
66 \*---------------------------------------------------------------------------*/
70 public thermophysicalFunction
74 // NSRDS function 114 coefficients
75 scalar Tc_, a_, b_, c_, d_;
80 //- Runtime type information
81 TypeName("NSRDSfunc14");
86 //- Construct from components
96 //- Construct from Istream
97 NSRDSfunc14(Istream& is);
99 //- Construct from dictionary
100 NSRDSfunc14(const dictionary& dict);
105 //- Evaluate the function and return the result
106 scalar f(scalar, scalar T) const
108 scalar Tdash = min(T, Tc_ - ROOTVSMALL);
110 scalar t = 1.0 - Tdash/Tc_;
112 a_*a_/(t + ROOTVSMALL) + b_ - t
115 + t*(a_*d_ + t*(c_*c_/3.0 + t*(0.5*c_*d_ + 0.2*d_*d_*t)))
119 //- Write the function coefficients
120 void writeData(Ostream& os) const
122 os << Tc_ << token::SPACE
123 << a_ << token::SPACE
124 << b_ << token::SPACE
125 << c_ << token::SPACE
132 friend Ostream& operator<<(Ostream& os, const NSRDSfunc14& f)
140 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
142 } // End namespace Foam
144 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
148 // ************************************************************************* //