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