src/thermophysicalModels/reactionThermo/mixtures/homogeneousMixture/homogeneousMixture.C

   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 "homogeneousMixture.H"
  27 #include "fvMesh.H"
  28
  29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
  30
  31 template<class ThermoType>
  32 const char* Foam::homogeneousMixture<ThermoType>::specieNames_[1] = {"b"};
  33
  34
  35 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  36
  37 template<class ThermoType>
  38 Foam::homogeneousMixture<ThermoType>::homogeneousMixture
  39 (
  40     const dictionary& thermoDict,
  41     const fvMesh& mesh
  42 )
  43 :
  44     basicMultiComponentMixture
  45     (
  46         thermoDict,
  47         speciesTable(nSpecies_, specieNames_),
  48         mesh
  49     ),
  50
  51     reactants_(thermoDict.subDict("reactants")),
  52     products_(thermoDict.subDict("products")),
  53     mixture_("mixture", reactants_),
  54     b_(Y("b"))
  55 {}
  56
  57
  58 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
  59
  60 template<class ThermoType>
  61 const ThermoType& Foam::homogeneousMixture<ThermoType>::mixture
  62 (
  63     const scalar b
  64 ) const
  65 {
  66     if (b > 0.999)
  67     {
  68         return reactants_;
  69     }
  70     else if (b < 0.001)
  71     {
  72         return products_;
  73     }
  74     else
  75     {
  76         mixture_ = b/reactants_.W()*reactants_;
  77         mixture_ += (1 - b)/products_.W()*products_;
  78
  79         return mixture_;
  80     }
  81 }
  82
  83
  84 template<class ThermoType>
  85 void Foam::homogeneousMixture<ThermoType>::read(const dictionary& thermoDict)
  86 {
  87     reactants_ = ThermoType(thermoDict.subDict("reactants"));
  88     products_ = ThermoType(thermoDict.subDict("products"));
  89 }
  90
  91
  92 template<class ThermoType>
  93 const ThermoType& Foam::homogeneousMixture<ThermoType>::getLocalThermo
  94 (
  95     const label specieI
  96 ) const
  97 {
  98     if (specieI == 0)
  99     {
 100         return reactants_;
 101     }
 102     else if (specieI == 1)
 103     {
 104         return products_;
 105     }
 106     else
 107     {
 108         FatalErrorIn
 109         (
 110             "const ThermoType& Foam::homogeneousMixture<ThermoType>::"
 111             "getLocalThermo"
 112             "("
 113                 "const label "
 114             ") const"
 115         )   << "Unknown specie index " << specieI << ". Valid indices are 0..1"
 116             << abort(FatalError);
 117
 118         return reactants_;
 119     }
 120 }
 121
 122
 123 template<class ThermoType>
 124 Foam::scalar Foam::homogeneousMixture<ThermoType>::nMoles
 125 (
 126     const label specieI
 127 ) const
 128 {
 129     return getLocalThermo(specieI).nMoles();
 130 }
 131
 132
 133 template<class ThermoType>
 134 Foam::scalar Foam::homogeneousMixture<ThermoType>::W
 135 (
 136     const label specieI
 137 ) const
 138 {
 139     return getLocalThermo(specieI).W();
 140 }
 141
 142
 143 template<class ThermoType>
 144 Foam::scalar Foam::homogeneousMixture<ThermoType>::Cp
 145 (
 146     const label specieI,
 147     const scalar T
 148 ) const
 149 {
 150     return getLocalThermo(specieI).Cp(T);
 151 }
 152
 153
 154 template<class ThermoType>
 155 Foam::scalar Foam::homogeneousMixture<ThermoType>::Cv
 156 (
 157     const label specieI,
 158     const scalar T
 159 ) const
 160 {
 161     return getLocalThermo(specieI).Cv(T);
 162 }
 163
 164
 165 template<class ThermoType>
 166 Foam::scalar Foam::homogeneousMixture<ThermoType>::H
 167 (
 168     const label specieI,
 169     const scalar T
 170 ) const
 171 {
 172     return getLocalThermo(specieI).H(T);
 173 }
 174
 175
 176 template<class ThermoType>
 177 Foam::scalar Foam::homogeneousMixture<ThermoType>::Hs
 178 (
 179     const label specieI,
 180     const scalar T
 181 ) const
 182 {
 183     return getLocalThermo(specieI).Hs(T);
 184 }
 185
 186
 187 template<class ThermoType>
 188 Foam::scalar Foam::homogeneousMixture<ThermoType>::Hc
 189 (
 190     const label specieI
 191 ) const
 192 {
 193     return getLocalThermo(specieI).Hc();
 194 }
 195
 196
 197 template<class ThermoType>
 198 Foam::scalar Foam::homogeneousMixture<ThermoType>::S
 199 (
 200     const label specieI,
 201     const scalar T
 202 ) const
 203 {
 204     return getLocalThermo(specieI).S(T);
 205 }
 206
 207
 208 template<class ThermoType>
 209 Foam::scalar Foam::homogeneousMixture<ThermoType>::E
 210 (
 211     const label specieI,
 212     const scalar T
 213 ) const
 214 {
 215     return getLocalThermo(specieI).E(T);
 216 }
 217
 218
 219 template<class ThermoType>
 220 Foam::scalar Foam::homogeneousMixture<ThermoType>::G
 221 (
 222     const label specieI,
 223     const scalar T
 224 ) const
 225 {
 226     return getLocalThermo(specieI).G(T);
 227 }
 228
 229
 230 template<class ThermoType>
 231 Foam::scalar Foam::homogeneousMixture<ThermoType>::A
 232 (
 233     const label specieI,
 234     const scalar T
 235 ) const
 236 {
 237     return getLocalThermo(specieI).A(T);
 238 }
 239
 240
 241 template<class ThermoType>
 242 Foam::scalar Foam::homogeneousMixture<ThermoType>::mu
 243 (
 244     const label specieI,
 245     const scalar T
 246 ) const
 247 {
 248     return getLocalThermo(specieI).mu(T);
 249 }
 250
 251
 252 template<class ThermoType>
 253 Foam::scalar Foam::homogeneousMixture<ThermoType>::kappa
 254 (
 255     const label specieI,
 256     const scalar T
 257 ) const
 258 {
 259     return getLocalThermo(specieI).kappa(T);
 260 }
 261
 262
 263 template<class ThermoType>
 264 Foam::scalar Foam::homogeneousMixture<ThermoType>::alpha
 265 (
 266     const label specieI,
 267     const scalar T
 268 ) const
 269 {
 270     return getLocalThermo(specieI).alpha(T);
 271 }
 272
 273
 274 // ************************************************************************* //