src/thermophysicalModels/thermophysicalFunctions/NSRDSfunctions/NSRDSfunc6/NSRDSfunc6.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation, either version 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Class
  25     Foam::NSRDSfunc6
  26
  27 Description
  28     NSRDS function number 106
  29
  30     Source:
  31     @verbatim
  32                       NSRDS - AICHE
  33                  Data Compilation Tables
  34                     of Properties of
  35                      Pure Compounds
  36
  37         Design Institute for Physical Property Data
  38           American Institute of Chemical Engineers
  39                   345 East 47th Street
  40                 New York, New York 10017
  41
  42          National Standard Reference Data System
  43          American Institute of Chemical Engineers
  44
  45           T.E. Daubert       -       R.P. Danner
  46
  47             Department of Chemical Engineering
  48             The Pennsylvania State University
  49                 University Park, PA 16802
  50     @endverbatim
  51
  52 \*---------------------------------------------------------------------------*/
  53
  54 #ifndef NSRDSfunc6_H
  55 #define NSRDSfunc6_H
  56
  57 #include "thermophysicalFunction.H"
  58
  59 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  60
  61 namespace Foam
  62 {
  63
  64 /*---------------------------------------------------------------------------*\
  65                            Class NSRDSfunc6 Declaration
  66 \*---------------------------------------------------------------------------*/
  67
  68 class NSRDSfunc6
  69 :
  70     public thermophysicalFunction
  71 {
  72     // Private data
  73
  74         // NSRDS function 106 coefficients
  75         scalar Tc_, a_, b_, c_, d_, e_;
  76
  77
  78 public:
  79
  80     //- Runtime type information
  81     TypeName("NSRDSfunc6");
  82
  83
  84     // Constructors
  85
  86         //- Construct from components
  87         NSRDSfunc6(scalar Tc, scalar a, scalar b, scalar c, scalar d, scalar e)
  88         :
  89             Tc_(Tc),
  90             a_(a),
  91             b_(b),
  92             c_(c),
  93             d_(d),
  94             e_(e)
  95         {}
  96
  97         //- Construct from Istream
  98         NSRDSfunc6(Istream& is)
  99         :
 100             Tc_(readScalar(is)),
 101             a_(readScalar(is)),
 102             b_(readScalar(is)),
 103             c_(readScalar(is)),
 104             d_(readScalar(is)),
 105             e_(readScalar(is))
 106         {}
 107
 108
 109     // Member Functions
 110
 111         //- Evaluate the function and return the result
 112         scalar f(scalar, scalar T) const
 113         {
 114             scalar Tr = T/Tc_;
 115             return a_*pow(1 - Tr, ((e_*Tr + d_)*Tr + c_)*Tr + b_);
 116         }
 117
 118
 119         //- Write the function coefficients
 120         void writeData(Ostream& os) const
 121         {
 122             os  << Tc_ << token::SPACE
 123                 << a_ << token::SPACE
 124                 << b_ << token::SPACE
 125                 << c_ << token::SPACE
 126                 << d_ << token::SPACE
 127                 << e_;
 128         }
 129
 130
 131     // Ostream Operator
 132
 133         friend Ostream& operator<<(Ostream& os, const NSRDSfunc6& f)
 134         {
 135             f.writeData(os);
 136             return os;
 137         }
 138 };
 139
 140
 141 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 142
 143 } // End namespace Foam
 144
 145 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 146
 147 #endif
 148
 149 // ************************************************************************* //