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

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   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 \*---------------------------------------------------------------------------*/
  25
  26 #include "exponentialSolidThermo.H"
  27 #include "specie.H"
  28
  29 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  30
  31
  32 template<class rhoType>
  33 inline void Foam::exponentialSolidThermo<rhoType>::checkT(const scalar T) const
  34 {
  35     if (T < 0.0)
  36     {
  37         FatalErrorIn
  38         (
  39             "exponentialSolidThermo<rhoType>::checkT(const scalar T) const"
  40         )   << "attempt to use exponentialSolidThermo<rhoType>"
  41                " out of temperature range "
  42             << T
  43             << abort(FatalError);
  44     }
  45 }
  46
  47
  48 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  49
  50 template<class rhoType>
  51 inline Foam::exponentialSolidThermo<rhoType>::exponentialSolidThermo
  52 (
  53     const word& name,
  54     const exponentialSolidThermo& jt
  55 )
  56 :
  57     rhoType(name, jt),
  58     c0_(jt.c0_),
  59     n0_(jt.n0_),
  60     Tref_(jt.Tref_),
  61     Hf_(jt.Hf_)
  62 {}
  63
  64
  65 template<class rhoType>
  66 inline Foam::exponentialSolidThermo<rhoType>::exponentialSolidThermo
  67 (
  68     const rhoType& st,
  69     const scalar c0,
  70     const scalar n0,
  71     const scalar Tref,
  72     const scalar Hf
  73 )
  74 :
  75     rhoType(st),
  76     c0_(c0),
  77     n0_(n0),
  78     Tref_(Tref),
  79     Hf_(Hf)
  80 {}
  81
  82
  83 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  84
  85 template<class rhoType>
  86 inline Foam::scalar Foam::exponentialSolidThermo<rhoType>::integrateCp
  87 (
  88     const scalar T
  89 ) const
  90 {
  91     return
  92     (
  93         c0_*pow(T, n0_ + 1.0)
  94        /(pow(Tref_, n0_)*(n0_ + 1.0))
  95     );
  96 }
  97
  98
  99 template<class rhoType>
 100 inline Foam::scalar Foam::exponentialSolidThermo<rhoType>::Cp
 101 (
 102     const scalar T
 103 ) const
 104 {
 105     return c0_*pow(T/Tref_, n0_);
 106 }
 107
 108
 109 template<class rhoType>
 110 inline Foam::scalar Foam::exponentialSolidThermo<rhoType>::h
 111 (
 112     const scalar T
 113 ) const
 114 {
 115     scalar hOffset = integrateCp(specie::Tstd);
 116
 117     return
 118     (
 119         integrateCp(T) + Hf_ - hOffset
 120     );
 121 }
 122
 123
 124 template<class rhoType>
 125 inline Foam::scalar Foam::exponentialSolidThermo<rhoType>::hs
 126 (
 127     const scalar T
 128 ) const
 129 {
 130     return h(T) - Hf_;
 131 }
 132
 133
 134 template<class rhoType>
 135 inline Foam::scalar Foam::exponentialSolidThermo<rhoType>::hf() const
 136 {
 137     return Hf_;
 138 }
 139
 140
 141 // ************************************************************************* //