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) 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 "icoPolynomial.H"
  27
  28 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  29
  30 template<int PolySize>
  31 inline Foam::icoPolynomial<PolySize>::icoPolynomial
  32 (
  33     const specie& sp,
  34     const Polynomial<PolySize>& rhoCoeffs
  35 )
  36 :
  37     specie(sp),
  38     rhoCoeffs_(rhoCoeffs)
  39 {}
  40
  41
  42 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  43
  44 template<int PolySize>
  45 inline Foam::icoPolynomial<PolySize>::icoPolynomial
  46 (
  47     const icoPolynomial<PolySize>& ip
  48 )
  49 :
  50     specie(ip),
  51     rhoCoeffs_(ip.rhoCoeffs_)
  52 {}
  53
  54
  55 template<int PolySize>
  56 inline Foam::icoPolynomial<PolySize>::icoPolynomial
  57 (
  58     const word& name,
  59     const icoPolynomial<PolySize>& ip
  60 )
  61 :
  62     specie(name, ip),
  63     rhoCoeffs_(ip.rhoCoeffs_)
  64 {}
  65
  66
  67 template<int PolySize>
  68 inline Foam::autoPtr<Foam::icoPolynomial<PolySize> >
  69 Foam::icoPolynomial<PolySize>::clone() const
  70 {
  71     return autoPtr<icoPolynomial<PolySize> >
  72     (
  73         new icoPolynomial<PolySize>(*this)
  74     );
  75 }
  76
  77
  78 template<int PolySize>
  79 inline Foam::autoPtr<Foam::icoPolynomial<PolySize> >
  80 Foam::icoPolynomial<PolySize>::New(Istream& is)
  81 {
  82     return autoPtr<icoPolynomial<PolySize> >(new icoPolynomial<PolySize>(is));
  83 }
  84
  85
  86 template<int PolySize>
  87 inline Foam::autoPtr<Foam::icoPolynomial<PolySize> >
  88 Foam::icoPolynomial<PolySize>::New(const dictionary& dict)
  89 {
  90     return autoPtr<icoPolynomial<PolySize> >
  91     (
  92         new icoPolynomial<PolySize>(dict)
  93     );
  94 }
  95
  96
  97 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  98
  99 template<int PolySize>
 100 inline Foam::scalar Foam::icoPolynomial<PolySize>::rho(scalar, scalar T) const
 101 {
 102     return rhoCoeffs_.value(T)/this->W();
 103 }
 104
 105
 106 template<int PolySize>
 107 inline Foam::scalar Foam::icoPolynomial<PolySize>::psi(scalar, scalar) const
 108 {
 109     return 0.0;
 110 }
 111
 112
 113 template<int PolySize>
 114 inline Foam::scalar Foam::icoPolynomial<PolySize>::Z(scalar, scalar) const
 115 {
 116     return 0.0;
 117 }
 118
 119
 120 // * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
 121
 122 template<int PolySize>
 123 inline Foam::icoPolynomial<PolySize>& Foam::icoPolynomial<PolySize>::operator=
 124 (
 125     const icoPolynomial<PolySize>& ip
 126 )
 127 {
 128     specie::operator=(ip);
 129
 130     rhoCoeffs_ = ip.rhoCoeffs_;
 131
 132     return *this;
 133 }
 134
 135
 136 template<int PolySize>
 137 inline void Foam::icoPolynomial<PolySize>::operator+=
 138 (
 139     const icoPolynomial<PolySize>& ip
 140 )
 141 {
 142     scalar molr1 = this->nMoles();
 143
 144     specie::operator+=(ip);
 145
 146     molr1 /= this->nMoles();
 147     scalar molr2 = ip.nMoles()/this->nMoles();
 148
 149     rhoCoeffs_ = molr1*rhoCoeffs_ + molr2*ip.rhoCoeffs_;
 150 }
 151
 152
 153 template<int PolySize>
 154 inline void Foam::icoPolynomial<PolySize>::operator-=
 155 (
 156     const icoPolynomial<PolySize>& ip
 157 )
 158 {
 159     scalar molr1 = this->nMoles();
 160
 161     specie::operator-=(ip);
 162
 163     molr1 /= this->nMoles();
 164     scalar molr2 = ip.nMoles()/this->nMoles();
 165
 166     rhoCoeffs_ = molr1*rhoCoeffs_ - molr2*ip.rhoCoeffs_;
 167 }
 168
 169
 170 template<int PolySize>
 171 inline void Foam::icoPolynomial<PolySize>::operator*=(const scalar s)
 172 {
 173     specie::operator*=(s);
 174 }
 175
 176
 177 // * * * * * * * * * * * * * * * Friend Operators  * * * * * * * * * * * * * //
 178
 179 template<int PolySize>
 180 Foam::icoPolynomial<PolySize> Foam::operator+
 181 (
 182     const icoPolynomial<PolySize>& ip1,
 183     const icoPolynomial<PolySize>& ip2
 184 )
 185 {
 186     scalar nMoles = ip1.nMoles() + ip2.nMoles();
 187     scalar molr1 = ip1.nMoles()/nMoles;
 188     scalar molr2 = ip2.nMoles()/nMoles;
 189
 190     return icoPolynomial<PolySize>
 191     (
 192         static_cast<const specie&>(ip1)
 193       + static_cast<const specie&>(ip2),
 194         molr1*ip1.rhoCoeffs_ + molr2*ip2.rhoCoeffs_
 195     );
 196 }
 197
 198
 199 template<int PolySize>
 200 Foam::icoPolynomial<PolySize> Foam::operator-
 201 (
 202     const icoPolynomial<PolySize>& ip1,
 203     const icoPolynomial<PolySize>& ip2
 204 )
 205 {
 206     scalar nMoles = ip1.nMoles() + ip2.nMoles();
 207     scalar molr1 = ip1.nMoles()/nMoles;
 208     scalar molr2 = ip2.nMoles()/nMoles;
 209
 210     return icoPolynomial<PolySize>
 211     (
 212         static_cast<const specie&>(ip1)
 213       - static_cast<const specie&>(ip2),
 214         molr1*ip1.rhoCoeffs_ - molr2*ip2.rhoCoeffs_
 215     );
 216 }
 217
 218
 219 template<int PolySize>
 220 Foam::icoPolynomial<PolySize> Foam::operator*
 221 (
 222     const scalar s,
 223     const icoPolynomial<PolySize>& ip
 224 )
 225 {
 226     return icoPolynomial<PolySize>
 227     (
 228         s*static_cast<const specie&>(ip),
 229         ip.rhoCoeffs_
 230     );
 231 }
 232
 233
 234 template<int PolySize>
 235 Foam::icoPolynomial<PolySize> Foam::operator==
 236 (
 237     const icoPolynomial<PolySize>& ip1,
 238     const icoPolynomial<PolySize>& ip2
 239 )
 240 {
 241     return ip2 - ip1;
 242 }
 243
 244
 245 // ************************************************************************* //