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