src/thermophysicalModels/basicSolidThermo/solidMixtureThermo/mixtures/multiComponentSolidMixture/multiComponentSolidMixture.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2010-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 #include "multiComponentSolidMixture.H"
  27
  28 // * * * * * * * * * * * * * Private Member Functions  * * * * * * * * * * * //
  29
  30 template<class ThermoSolidType>
  31 void Foam::multiComponentSolidMixture<ThermoSolidType>::correctMassFractions()
  32 {
  33     volScalarField Yt("Yt", Y_[0]);
  34
  35     for (label n=1; n<Y_.size(); n++)
  36     {
  37         Yt += Y_[n];
  38     }
  39
  40     forAll(Y_, n)
  41     {
  42         Y_[n] /= Yt;
  43     }
  44
  45
  46 }
  47
  48
  49 template<class ThermoSolidType>
  50 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::X
  51 (
  52     label iComp, label celli, scalar T
  53 ) const
  54 {
  55     scalar rhoInv = 0.0;
  56     forAll(solidData_, i)
  57     {
  58         rhoInv += Y_[i][celli]/solidData_[i].rho(T);
  59     }
  60
  61     scalar X = Y_[iComp][celli]/solidData_[iComp].rho(T);
  62
  63     return (X/rhoInv);
  64 }
  65
  66
  67 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
  68
  69 template<class ThermoSolidType>
  70 Foam::multiComponentSolidMixture<ThermoSolidType>::multiComponentSolidMixture
  71 (
  72     const dictionary& thermoSolidDict,
  73     const fvMesh& mesh
  74 )
  75 :
  76     basicSolidMixture
  77     (
  78         thermoSolidDict.lookup("solidComponents"),
  79         mesh
  80     ),
  81     solidData_(components_.size())
  82 {
  83
  84     forAll(components_, i)
  85     {
  86         solidData_.set
  87         (
  88             i,
  89             new ThermoSolidType
  90             (
  91                 thermoSolidDict.subDict(components_[i] + "Coeffs")
  92             )
  93         );
  94     }
  95     correctMassFractions();
  96 }
  97
  98
  99 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 100
 101 template<class ThermoSolidType>
 102 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::rho
 103 (
 104     scalar T, label celli
 105 ) const
 106 {
 107     scalar tmp = 0.0;
 108     forAll(solidData_, i)
 109     {
 110         tmp += solidData_[i].rho(T)*X(i, celli, T);
 111     }
 112     return tmp;
 113 }
 114
 115
 116 template<class ThermoSolidType>
 117 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::hf
 118 (
 119     scalar, label celli
 120 ) const
 121 {
 122     scalar tmp = 0.0;
 123     forAll(solidData_, i)
 124     {
 125         tmp += solidData_[i].hf()*Y_[i][celli];
 126     }
 127     return tmp;
 128 }
 129
 130
 131 template<class ThermoSolidType>
 132 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::hs
 133 (
 134     scalar T, label celli
 135 ) const
 136 {
 137     scalar tmp = 0.0;
 138     forAll(solidData_, i)
 139     {
 140         tmp += solidData_[i].hs(T)*Y_[i][celli];
 141     }
 142     return tmp;
 143 }
 144
 145
 146 template<class ThermoSolidType>
 147 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::h
 148 (
 149     scalar T, label celli
 150 ) const
 151 {
 152     scalar tmp = 0.0;
 153     forAll(solidData_, i)
 154     {
 155         tmp += solidData_[i].h(T)*Y_[i][celli];
 156     }
 157     return tmp;
 158 }
 159
 160
 161 template<class ThermoSolidType>
 162 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::kappa
 163 (
 164     scalar T, label celli
 165 ) const
 166 {
 167     scalar tmp = 0.0;
 168     forAll(solidData_, i)
 169     {
 170         tmp += solidData_[i].kappa(T)*X(i, celli, T);
 171     }
 172     return tmp;
 173 }
 174
 175
 176 template<class ThermoSolidType>
 177 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::sigmaS
 178 (
 179     scalar T, label celli
 180 ) const
 181 {
 182     scalar tmp = 0.0;
 183     forAll(solidData_, i)
 184     {
 185         tmp += solidData_[i].sigmaS(T)*X(i, celli, T);
 186     }
 187     return tmp;
 188 }
 189
 190
 191 template<class ThermoSolidType>
 192 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::K
 193 (
 194     scalar T, label celli
 195 ) const
 196 {
 197     scalar tmp = 0.0;
 198     forAll(solidData_, i)
 199     {
 200         tmp += solidData_[i].K(T)*X(i, celli, T);
 201     }
 202     return tmp;
 203 }
 204
 205
 206 template<class ThermoSolidType>
 207 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::emissivity
 208 (
 209     scalar T, label celli
 210 ) const
 211 {
 212     scalar tmp = 0.0;
 213     forAll(solidData_, i)
 214     {
 215         tmp += solidData_[i].emissivity(T)*Y_[i][celli];
 216     }
 217     return tmp;
 218 }
 219
 220
 221 template<class ThermoSolidType>
 222 Foam::scalar Foam::multiComponentSolidMixture<ThermoSolidType>::Cp
 223 (
 224     scalar T, label celli
 225 ) const
 226 {
 227     scalar tmp = 0.0;
 228     forAll(solidData_, i)
 229     {
 230         tmp += solidData_[i].Cp(T)*Y_[i][celli];
 231     }
 232     return tmp;
 233 }
 234
 235
 236 template<class ThermoSolidType>
 237 void Foam::multiComponentSolidMixture<ThermoSolidType>::read
 238 (
 239     const dictionary& thermoDict
 240 )
 241 {
 242     forAll(components_, i)
 243     {
 244         solidData_[i] =
 245             ThermoSolidType(thermoDict.subDict(components_[i] + "Coeffs"));
 246     }
 247 }
 248
 249
 250 // ************************************************************************* //