src/thermophysicalModels/reactionThermo/mixtures/veryInhomogeneousMixture/veryInhomogeneousMixture.H

   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 Class
  25     Foam::veryInhomogeneousMixture
  26
  27 Description
  28     Foam::veryInhomogeneousMixture
  29
  30 SourceFiles
  31     veryInhomogeneousMixture.C
  32
  33 \*---------------------------------------------------------------------------*/
  34
  35 #ifndef veryInhomogeneousMixture_H
  36 #define veryInhomogeneousMixture_H
  37
  38 #include "basicMultiComponentMixture.H"
  39
  40 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  41
  42 namespace Foam
  43 {
  44
  45 /*---------------------------------------------------------------------------*\
  46                   Class veryInhomogeneousMixture Declaration
  47 \*---------------------------------------------------------------------------*/
  48
  49 template<class ThermoType>
  50 class veryInhomogeneousMixture
  51 :
  52     public basicMultiComponentMixture
  53 {
  54     // Private data
  55
  56         static const int nSpecies_ = 3;
  57         static const char* specieNames_[3];
  58
  59         dimensionedScalar stoicRatio_;
  60
  61         ThermoType fuel_;
  62         ThermoType oxidant_;
  63         ThermoType products_;
  64
  65         mutable ThermoType mixture_;
  66
  67         //- Mixture fraction
  68         volScalarField& ft_;
  69
  70         //- Fuel mass fraction
  71         volScalarField& fu_;
  72
  73         //- Regress variable
  74         volScalarField& b_;
  75
  76         //- Construct as copy (not implemented)
  77         veryInhomogeneousMixture(const veryInhomogeneousMixture<ThermoType>&);
  78
  79
  80 public:
  81
  82     //- The type of thermodynamics this mixture is instantiated for
  83     typedef ThermoType thermoType;
  84
  85
  86     // Constructors
  87
  88         //- Construct from dictionary and mesh
  89         veryInhomogeneousMixture(const dictionary&, const fvMesh&);
  90
  91
  92     //- Destructor
  93     virtual ~veryInhomogeneousMixture()
  94     {}
  95
  96
  97     // Member functions
  98
  99         const dimensionedScalar& stoicRatio() const
 100         {
 101             return stoicRatio_;
 102         }
 103
 104         const ThermoType& mixture(const scalar, const scalar) const;
 105
 106         const ThermoType& cellMixture(const label celli) const
 107         {
 108             return mixture(ft_[celli], fu_[celli]);
 109         }
 110
 111         const ThermoType& patchFaceMixture
 112         (
 113             const label patchi,
 114             const label facei
 115         ) const
 116         {
 117             return mixture
 118             (
 119                 ft_.boundaryField()[patchi][facei],
 120                 fu_.boundaryField()[patchi][facei]
 121             );
 122         }
 123
 124         const ThermoType& cellReactants(const label celli) const
 125         {
 126             return mixture(ft_[celli], ft_[celli]);
 127         }
 128
 129         const ThermoType& patchFaceReactants
 130         (
 131             const label patchi,
 132             const label facei
 133         ) const
 134         {
 135             return mixture
 136             (
 137                 ft_.boundaryField()[patchi][facei],
 138                 ft_.boundaryField()[patchi][facei]
 139             );
 140         }
 141
 142         const ThermoType& cellProducts(const label celli) const
 143         {
 144             scalar ft = ft_[celli];
 145             return mixture(ft, fres(ft, stoicRatio().value()));
 146         }
 147
 148         const ThermoType& patchFaceProducts
 149         (
 150             const label patchi,
 151             const label facei
 152         ) const
 153         {
 154             scalar ft = ft_.boundaryField()[patchi][facei];
 155             return mixture(ft, fres(ft, stoicRatio().value()));
 156         }
 157
 158         //- Read dictionary
 159         void read(const dictionary&);
 160
 161         //- Return thermo based on index
 162         const ThermoType& getLocalThermo(const label specieI) const;
 163
 164
 165         // Per specie properties
 166
 167             //- Number of moles []
 168             virtual scalar nMoles(const label specieI) const;
 169
 170             //- Molecular weight [kg/kmol]
 171             virtual scalar W(const label specieI) const;
 172
 173
 174         // Per specie thermo properties
 175
 176             //- Heat capacity at constant pressure [J/(kg K)]
 177             virtual scalar Cp(const label specieI, const scalar T) const;
 178
 179             //- Heat capacity at constant volume [J/(kg K)]
 180             virtual scalar Cv(const label specieI, const scalar T) const;
 181
 182             //- Enthalpy [J/kg]
 183             virtual scalar H(const label specieI, const scalar T) const;
 184
 185             //- Sensible enthalpy [J/kg]
 186             virtual scalar Hs(const label specieI, const scalar T) const;
 187
 188             //- Chemical enthalpy [J/kg]
 189             virtual scalar Hc(const label specieI) const;
 190
 191             //- Entropy [J/(kg K)]
 192             virtual scalar S(const label specieI, const scalar T) const;
 193
 194             //- Internal energy [J/kg]
 195             virtual scalar E(const label specieI, const scalar T) const;
 196
 197             //- Gibbs free energy [J/kg]
 198             virtual scalar G(const label specieI, const scalar T) const;
 199
 200             //- Helmholtz free energy [J/kg]
 201             virtual scalar A(const label specieI, const scalar T) const;
 202
 203
 204         // Per specie transport properties
 205
 206             //- Dynamic viscosity [kg/m/s]
 207             virtual scalar mu(const label specieI, const scalar T) const;
 208
 209             //- Thermal conductivity [W/m/K]
 210             virtual scalar kappa(const label specieI, const scalar T) const;
 211
 212             //- Thermal diffusivity [kg/m/s]
 213             virtual scalar alpha(const label specieI, const scalar T) const;
 214 };
 215
 216
 217 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 218
 219 } // End namespace Foam
 220
 221 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
 222
 223 #ifdef NoRepository
 224 #   include "veryInhomogeneousMixture.C"
 225 #endif
 226
 227 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 228
 229 #endif
 230
 231 // ************************************************************************* //