src/thermophysicalModels/specie/equationOfState/icoPolynomial/icoPolynomialI.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2008-2009 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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #include "icoPolynomial.H"
  28
  29 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  30
  31 template<int PolySize>
  32 inline Foam::icoPolynomial<PolySize>::icoPolynomial
  33 (
  34     const specie& sp,
  35     const Polynomial<PolySize>& rhoPoly
  36 )
  37 :
  38     specie(sp),
  39     rhoPolynomial_(rhoPoly)
  40 {}
  41
  42
  43 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  44
  45 template<int PolySize>
  46 inline Foam::icoPolynomial<PolySize>::icoPolynomial
  47 (
  48     const word& name,
  49     const icoPolynomial<PolySize>& ip
  50 )
  51 :
  52     specie(name, ip),
  53     rhoPolynomial_(ip.rhoPolynomial_)
  54 {}
  55
  56
  57 template<int PolySize>
  58 inline Foam::autoPtr<Foam::icoPolynomial<PolySize> >
  59 Foam::icoPolynomial<PolySize>::clone() const
  60 {
  61     return autoPtr<icoPolynomial<PolySize> >
  62     (
  63         new icoPolynomial<PolySize>(*this)
  64     );
  65 }
  66
  67
  68 template<int PolySize>
  69 inline Foam::autoPtr<Foam::icoPolynomial<PolySize> >
  70 Foam::icoPolynomial<PolySize>::New(Istream& is)
  71 {
  72     return autoPtr<icoPolynomial<PolySize> >(new icoPolynomial<PolySize>(is));
  73 }
  74
  75
  76 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  77
  78 template<int PolySize>
  79 inline Foam::scalar Foam::icoPolynomial<PolySize>::rho(scalar, scalar T) const
  80 {
  81     return rhoPolynomial_.evaluate(T);
  82 }
  83
  84
  85 template<int PolySize>
  86 inline Foam::scalar Foam::icoPolynomial<PolySize>::psi(scalar, scalar) const
  87 {
  88     return 0.0;
  89 }
  90
  91
  92 template<int PolySize>
  93 inline Foam::scalar Foam::icoPolynomial<PolySize>::Z(scalar, scalar) const
  94 {
  95     return 0.0;
  96 }
  97
  98
  99 // * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
 100
 101 template<int PolySize>
 102 inline void Foam::icoPolynomial<PolySize>::operator+=
 103 (
 104     const icoPolynomial<PolySize>& ip
 105 )
 106 {
 107     scalar molr1 = this->nMoles();
 108
 109     specie::operator+=(ip);
 110
 111     molr1 /= this->nMoles();
 112     scalar molr2 = ip.nMoles()/this->nMoles();
 113
 114     rhoPolynomial_ = molr1*rhoPolynomial_ + molr2*ip.rhoPolynomial_;
 115 }
 116
 117
 118 template<int PolySize>
 119 inline void Foam::icoPolynomial<PolySize>::operator-=
 120 (
 121     const icoPolynomial<PolySize>& ip
 122 )
 123 {
 124     scalar molr1 = this->nMoles();
 125     specie::operator-=(ip);
 126
 127     molr1 /= this->nMoles();
 128     scalar molr2 = ip.nMoles()/this->nMoles();
 129
 130     rhoPolynomial_ = molr1*rhoPolynomial_ - molr2*ip.rhoPolynomial_;
 131 }
 132
 133
 134 template<int PolySize>
 135 inline void Foam::icoPolynomial<PolySize>::operator*=(const scalar s)
 136 {
 137     specie::operator*=(s);
 138 }
 139
 140
 141 // * * * * * * * * * * * * * * * Friend Operators  * * * * * * * * * * * * * //
 142
 143 template<int PolySize>
 144 Foam::icoPolynomial<PolySize> Foam::operator+
 145 (
 146     const icoPolynomial<PolySize>& ip1,
 147     const icoPolynomial<PolySize>& ip2
 148 )
 149 {
 150     scalar mol1 = ip1.nMoles();
 151     scalar mol2 = ip2.nMoles();
 152     scalar nMoles = mol1 + mol2;
 153
 154     return icoPolynomial<PolySize>
 155     (
 156         static_cast<const specie&>(ip1)
 157       + static_cast<const specie&>(ip2),
 158         (mol1/nMoles)*ip1.rhoPolynomial_ + (mol2/nMoles)*ip2.rhoPolynomial_
 159     );
 160 }
 161
 162
 163 template<int PolySize>
 164 Foam::icoPolynomial<PolySize> Foam::operator-
 165 (
 166     const icoPolynomial<PolySize>& ip1,
 167     const icoPolynomial<PolySize>& ip2
 168 )
 169 {
 170     scalar mol1 = ip1.nMoles();
 171     scalar mol2 = ip2.nMoles();
 172     scalar nMoles = mol1 + mol2;
 173
 174     return icoPolynomial<PolySize>
 175     (
 176         static_cast<const specie&>(ip1)
 177       - static_cast<const specie&>(ip2),
 178         (mol1/nMoles)*ip1.rhoPolynomial_ - (mol2/nMoles)*ip2.rhoPolynomial_
 179     );
 180 }
 181
 182
 183 template<int PolySize>
 184 Foam::icoPolynomial<PolySize> Foam::operator*
 185 (
 186     const scalar s,
 187     const icoPolynomial<PolySize>& ip
 188 )
 189 {
 190     return icoPolynomial<PolySize>
 191     (
 192         s*static_cast<const specie&>(ip),
 193         ip.rhoPolynomial_
 194     );
 195 }
 196
 197
 198 template<int PolySize>
 199 Foam::icoPolynomial<PolySize> Foam::operator==
 200 (
 201     const icoPolynomial<PolySize>& ip1,
 202     const icoPolynomial<PolySize>& ip2
 203 )
 204 {
 205     return ip2 - ip1;
 206 }
 207
 208
 209 // ************************************************************************* //