src/thermophysicalModels/solid/thermo/exponential/exponentialSolidThermo.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2010-2010 OpenCFD Ltd.
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  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.
  15
  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.
  20
  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/>.
  23
  24 Class
  25     Foam::exponentialSolidThermo
  26
  27 Description
  28     Exponential thermodynamics package for solids templated into the rhoType.
  29
  30 SourceFiles
  31     exponentialSolidThermoI.H
  32     exponentialSolidThermo.C
  33
  34 \*---------------------------------------------------------------------------*/
  35
  36 #ifndef exponentialSolidThermo_H
  37 #define exponentialSolidThermo_H
  38
  39 #include "scalar.H"
  40 #include "dictionary.H"
  41 #include "specie.H"
  42
  43 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  44
  45 namespace Foam
  46 {
  47
  48 template<class rhoType> class exponentialSolidThermo;
  49
  50 template<class rhoType>
  51 Ostream& operator<<
  52 (
  53     Ostream&,
  54     const exponentialSolidThermo<rhoType>&
  55 );
  56
  57
  58 /*---------------------------------------------------------------------------*\
  59                          Class exponentialSolidThermo Declaration
  60 \*---------------------------------------------------------------------------*/
  61
  62 template<class rhoType>
  63 class exponentialSolidThermo
  64 :
  65     public rhoType
  66 {
  67
  68     // Private data
  69
  70         scalar c0_;
  71         scalar n0_;
  72         scalar Tref_;
  73         scalar Hf_;
  74
  75
  76     // Private Member Functions
  77
  78         //- Check given temperature is within the range of the fitted coeffs
  79         inline void checkT(const scalar T) const;
  80
  81         //- Integrate Cp expression
  82         inline scalar integrateCp(const scalar T) const;
  83
  84
  85 public:
  86
  87     // Constructors
  88
  89         //- Construct from components
  90         inline exponentialSolidThermo
  91         (
  92             const rhoType& st,
  93             const scalar c0,
  94             const scalar n0,
  95             const scalar Tref,
  96             const scalar Hf
  97         );
  98
  99         //- Construct from dictionary
 100         exponentialSolidThermo(const dictionary&);
 101
 102         //- Construct as a named copy
 103         inline exponentialSolidThermo
 104         (
 105             const word&,
 106             const exponentialSolidThermo&
 107         );
 108
 109
 110     // Member Functions
 111
 112         //- Heat capacity at constant pressure [J/(kg K)]
 113         inline scalar Cp(const scalar T) const;
 114
 115         //- Enthalpy [J/kmol]
 116         inline scalar h(const scalar T) const;
 117
 118         //- Sensible enthalpy [J/kg]
 119         inline scalar hs(const scalar T) const;
 120
 121         //- Formation enthalpy [J/kg]
 122         inline scalar hf() const;
 123
 124
 125     // Ostream Operator
 126
 127         friend Ostream& operator<< <rhoType>
 128         (
 129             Ostream&,
 130             const exponentialSolidThermo&
 131         );
 132 };
 133
 134
 135 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 136
 137 } // End namespace Foam
 138
 139 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 140
 141
 142 #ifdef NoRepository
 143 #   include "exponentialSolidThermoI.H"
 144 #   include "exponentialSolidThermo.C"
 145 #endif
 146
 147 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 148
 149 #endif
 150
 151 // ************************************************************************* //