src/thermophysicalModels/reactionThermo/mixtures/multiComponentMixture/multiComponentMixture.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 "multiComponentMixture.H"
  27
  28 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  29
  30 template<class ThermoType>
  31 const ThermoType& Foam::multiComponentMixture<ThermoType>::constructSpeciesData
  32 (
  33     const dictionary& thermoDict
  34 )
  35 {
  36     forAll(species_, i)
  37     {
  38         speciesData_.set
  39         (
  40             i,
  41             new ThermoType(thermoDict.subDict(species_[i]))
  42         );
  43     }
  44
  45     return speciesData_[0];
  46 }
  47
  48
  49 template<class ThermoType>
  50 void Foam::multiComponentMixture<ThermoType>::correctMassFractions()
  51 {
  52     // It changes Yt patches to "calculated"
  53     volScalarField Yt("Yt", 1.0*Y_[0]);
  54
  55     for (label n=1; n<Y_.size(); n++)
  56     {
  57         Yt += Y_[n];
  58     }
  59
  60     forAll(Y_, n)
  61     {
  62         Y_[n] /= Yt;
  63     }
  64 }
  65
  66
  67 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  68
  69 template<class ThermoType>
  70 Foam::multiComponentMixture<ThermoType>::multiComponentMixture
  71 (
  72     const dictionary& thermoDict,
  73     const wordList& specieNames,
  74     const HashPtrTable<ThermoType>& specieThermoData,
  75     const fvMesh& mesh
  76 )
  77 :
  78     basicMultiComponentMixture(thermoDict, specieNames, mesh),
  79     speciesData_(species_.size()),
  80     mixture_("mixture", *specieThermoData[specieNames[0]])
  81 {
  82     forAll(species_, i)
  83     {
  84         speciesData_.set
  85         (
  86             i,
  87             new ThermoType(*specieThermoData[species_[i]])
  88         );
  89     }
  90
  91     correctMassFractions();
  92 }
  93
  94
  95 template<class ThermoType>
  96 Foam::multiComponentMixture<ThermoType>::multiComponentMixture
  97 (
  98     const dictionary& thermoDict,
  99     const fvMesh& mesh
 100 )
 101 :
 102     basicMultiComponentMixture(thermoDict, thermoDict.lookup("species"), mesh),
 103     speciesData_(species_.size()),
 104     mixture_("mixture", constructSpeciesData(thermoDict))
 105 {
 106     correctMassFractions();
 107 }
 108
 109
 110 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 111
 112 template<class ThermoType>
 113 const ThermoType& Foam::multiComponentMixture<ThermoType>::cellMixture
 114 (
 115     const label celli
 116 ) const
 117 {
 118     mixture_ = Y_[0][celli]/speciesData_[0].W()*speciesData_[0];
 119
 120     for (label n=1; n<Y_.size(); n++)
 121     {
 122         mixture_ += Y_[n][celli]/speciesData_[n].W()*speciesData_[n];
 123     }
 124
 125     return mixture_;
 126 }
 127
 128
 129 template<class ThermoType>
 130 const ThermoType& Foam::multiComponentMixture<ThermoType>::patchFaceMixture
 131 (
 132     const label patchi,
 133     const label facei
 134 ) const
 135 {
 136     mixture_ =
 137         Y_[0].boundaryField()[patchi][facei]
 138        /speciesData_[0].W()*speciesData_[0];
 139
 140     for (label n=1; n<Y_.size(); n++)
 141     {
 142         mixture_ +=
 143             Y_[n].boundaryField()[patchi][facei]
 144            /speciesData_[n].W()*speciesData_[n];
 145     }
 146
 147     return mixture_;
 148 }
 149
 150
 151 template<class ThermoType>
 152 void Foam::multiComponentMixture<ThermoType>::read
 153 (
 154     const dictionary& thermoDict
 155 )
 156 {
 157     forAll(species_, i)
 158     {
 159         speciesData_[i] = ThermoType(thermoDict.subDict(species_[i]));
 160     }
 161 }
 162
 163
 164 template<class ThermoType>
 165 Foam::scalar Foam::multiComponentMixture<ThermoType>::nMoles
 166 (
 167     const label specieI
 168 ) const
 169 {
 170     return speciesData_[specieI].nMoles();
 171 }
 172
 173
 174 template<class ThermoType>
 175 Foam::scalar Foam::multiComponentMixture<ThermoType>::W
 176 (
 177     const label specieI
 178 ) const
 179 {
 180     return speciesData_[specieI].W();
 181 }
 182
 183
 184 template<class ThermoType>
 185 Foam::scalar Foam::multiComponentMixture<ThermoType>::Cp
 186 (
 187     const label specieI,
 188     const scalar T
 189 ) const
 190 {
 191     return speciesData_[specieI].Cp(T);
 192 }
 193
 194
 195 template<class ThermoType>
 196 Foam::scalar Foam::multiComponentMixture<ThermoType>::Cv
 197 (
 198     const label specieI,
 199     const scalar T
 200 ) const
 201 {
 202     return speciesData_[specieI].Cv(T);
 203 }
 204
 205
 206 template<class ThermoType>
 207 Foam::scalar Foam::multiComponentMixture<ThermoType>::H
 208 (
 209     const label specieI,
 210     const scalar T
 211 ) const
 212 {
 213     return speciesData_[specieI].H(T);
 214 }
 215
 216
 217 template<class ThermoType>
 218 Foam::scalar Foam::multiComponentMixture<ThermoType>::Hs
 219 (
 220     const label specieI,
 221     const scalar T
 222 ) const
 223 {
 224     return speciesData_[specieI].Hs(T);
 225 }
 226
 227
 228 template<class ThermoType>
 229 Foam::scalar Foam::multiComponentMixture<ThermoType>::Hc
 230 (
 231     const label specieI
 232 ) const
 233 {
 234     return speciesData_[specieI].Hc();
 235 }
 236
 237
 238 template<class ThermoType>
 239 Foam::scalar Foam::multiComponentMixture<ThermoType>::S
 240 (
 241     const label specieI,
 242     const scalar T
 243 ) const
 244 {
 245     return speciesData_[specieI].S(T);
 246 }
 247
 248
 249 template<class ThermoType>
 250 Foam::scalar Foam::multiComponentMixture<ThermoType>::E
 251 (
 252     const label specieI,
 253     const scalar T
 254 ) const
 255 {
 256     return speciesData_[specieI].E(T);
 257 }
 258
 259
 260 template<class ThermoType>
 261 Foam::scalar Foam::multiComponentMixture<ThermoType>::G
 262 (
 263     const label specieI,
 264     const scalar T
 265 ) const
 266 {
 267     return speciesData_[specieI].G(T);
 268 }
 269
 270
 271 template<class ThermoType>
 272 Foam::scalar Foam::multiComponentMixture<ThermoType>::A
 273 (
 274     const label specieI,
 275     const scalar T
 276 ) const
 277 {
 278     return speciesData_[specieI].A(T);
 279 }
 280
 281
 282 template<class ThermoType>
 283 Foam::scalar Foam::multiComponentMixture<ThermoType>::mu
 284 (
 285     const label specieI,
 286     const scalar T
 287 ) const
 288 {
 289     return speciesData_[specieI].mu(T);
 290 }
 291
 292
 293 template<class ThermoType>
 294 Foam::scalar Foam::multiComponentMixture<ThermoType>::kappa
 295 (
 296     const label specieI,
 297     const scalar T
 298 ) const
 299 {
 300     return speciesData_[specieI].kappa(T);
 301 }
 302
 303
 304 template<class ThermoType>
 305 Foam::scalar Foam::multiComponentMixture<ThermoType>::alpha
 306 (
 307     const label specieI,
 308     const scalar T
 309 ) const
 310 {
 311     return speciesData_[specieI].alpha(T);
 312 }
 313
 314
 315 // ************************************************************************* //