src/thermophysicalModels/specie/transport/const/constTransportI.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2004-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 \*---------------------------------------------------------------------------*/
  25
  26 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  27
  28 template<class Thermo>
  29 inline Foam::constTransport<Thermo>::constTransport
  30 (
  31     const Thermo& t,
  32     const scalar mu,
  33     const scalar Pr
  34 )
  35 :
  36     Thermo(t),
  37     mu_(mu),
  38     rPr_(1.0/Pr)
  39 {}
  40
  41
  42 template<class Thermo>
  43 inline Foam::constTransport<Thermo>::constTransport
  44 (
  45     const word& name,
  46     const constTransport& ct
  47 )
  48 :
  49     Thermo(name, ct),
  50     mu_(ct.mu_),
  51     rPr_(ct.rPr_)
  52 {}
  53
  54
  55 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  56
  57 template<class Thermo>
  58 inline Foam::scalar Foam::constTransport<Thermo>::mu(const scalar) const
  59 {
  60     return mu_;
  61 }
  62
  63
  64 template<class Thermo>
  65 inline Foam::scalar Foam::constTransport<Thermo>::kappa(const scalar T) const
  66 {
  67     return this->Cp(T)*mu(T)*rPr_;
  68 }
  69
  70
  71 template<class Thermo>
  72 inline Foam::scalar Foam::constTransport<Thermo>::alpha(const scalar T) const
  73 {
  74     scalar Cp_ = this->Cp(T);
  75
  76     scalar deltaT = T - specie::Tstd;
  77     scalar CpBar =
  78         (deltaT*(this->H(T) - this->H(specie::Tstd)) + Cp_)/(sqr(deltaT) + 1);
  79
  80     return Cp_*mu(T)*rPr_/CpBar;
  81 }
  82
  83
  84 // * * * * * * * * * * * * * * * Member Operators  * * * * * * * * * * * * * //
  85
  86 template<class Thermo>
  87 inline Foam::constTransport<Thermo>& Foam::constTransport<Thermo>::operator=
  88 (
  89     const constTransport<Thermo>& ct
  90 )
  91 {
  92     Thermo::operator=(ct);
  93
  94     mu_ = ct.mu_;
  95     rPr_ = ct.rPr_;
  96
  97     return *this;
  98 }
  99
 100
 101 // * * * * * * * * * * * * * * * Friend Operators  * * * * * * * * * * * * * //
 102
 103 template<class Thermo>
 104 inline Foam::constTransport<Thermo> Foam::operator+
 105 (
 106     const constTransport<Thermo>& ct1,
 107     const constTransport<Thermo>& ct2
 108 )
 109 {
 110     Thermo t
 111     (
 112         static_cast<const Thermo&>(ct1) + static_cast<const Thermo&>(ct2)
 113     );
 114
 115     scalar molr1 = ct1.nMoles()/t.nMoles();
 116     scalar molr2 = ct2.nMoles()/t.nMoles();
 117
 118     return constTransport<Thermo>
 119     (
 120         t,
 121         molr1*ct1.mu_ + molr2*ct2.mu_,
 122         molr1*ct1.rPr_ + molr2*ct2.rPr_
 123     );
 124 }
 125
 126
 127 template<class Thermo>
 128 inline Foam::constTransport<Thermo> Foam::operator-
 129 (
 130     const constTransport<Thermo>& ct1,
 131     const constTransport<Thermo>& ct2
 132 )
 133 {
 134     Thermo t
 135     (
 136         static_cast<const Thermo&>(ct1) - static_cast<const Thermo&>(ct2)
 137     );
 138
 139     scalar molr1 = ct1.nMoles()/t.nMoles();
 140     scalar molr2 = ct2.nMoles()/t.nMoles();
 141
 142     return constTransport<Thermo>
 143     (
 144         t,
 145         molr1*ct1.mu_ - molr2*ct2.mu_,
 146         molr1*ct1.rPr_ - molr2*ct2.rPr_
 147     );
 148 }
 149
 150
 151 template<class Thermo>
 152 inline Foam::constTransport<Thermo> Foam::operator*
 153 (
 154     const scalar s,
 155     const constTransport<Thermo>& ct
 156 )
 157 {
 158     return constTransport<Thermo>
 159     (
 160         s*static_cast<const Thermo&>(ct),
 161         ct.mu_,
 162         ct.rPr_
 163     );
 164 }
 165
 166
 167 template<class Thermo>
 168 inline Foam::constTransport<Thermo> Foam::operator==
 169 (
 170     const constTransport<Thermo>& ct1,
 171     const constTransport<Thermo>& ct2
 172 )
 173 {
 174     return ct2 - ct1;
 175 }
 176
 177
 178 // ************************************************************************* //