egrMixture.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011-2015 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 "egrMixture.H"
27 #include "fvMesh.H"
28 
29 // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
30 
31 template<class ThermoType>
32 const char* Foam::egrMixture<ThermoType>::specieNames_[3] = {"ft", "b", "egr"};
33 
34 
35 // * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
36 
37 template<class ThermoType>
39 (
40  const dictionary& thermoDict,
41  const fvMesh& mesh,
42  const word& phaseName
43 )
44 :
46  (
47  thermoDict,
48  speciesTable(nSpecies_, specieNames_),
49  mesh,
50  phaseName
51  ),
52 
53  stoicRatio_(thermoDict.lookup("stoichiometricAirFuelMassRatio")),
54 
55  fuel_(thermoDict.subDict("fuel")),
56  oxidant_(thermoDict.subDict("oxidant")),
57  products_(thermoDict.subDict("burntProducts")),
58 
59  mixture_("mixture", fuel_),
60 
61  ft_(Y("ft")),
62  b_(Y("b")),
63  egr_(Y("egr"))
64 {}
65 
66 
67 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
68 
69 template<class ThermoType>
71 (
72  const scalar ft,
73  const scalar b,
74  const scalar egr
75 ) const
76 {
77  if (ft < 0.0001)
78  {
79  return oxidant_;
80  }
81  else
82  {
83 
84  scalar fu = b*ft + (1.0 - b)*fres(ft, stoicRatio().value());
85  scalar ox = 1 - ft - (ft - fu)*stoicRatio().value();
86 
87  fu *= (1.0 - egr);
88  ox *= (1.0 - egr);
89 
90  scalar pr = 1 - fu - ox;
91 
92  mixture_ = fu/fuel_.W()*fuel_;
93  mixture_ += ox/oxidant_.W()*oxidant_;
94  mixture_ += pr/products_.W()*products_;
95 
96  return mixture_;
97  }
98 }
99 
100 
101 template<class ThermoType>
103 {
104  stoicRatio_ = thermoDict.lookup("stoichiometricAirFuelMassRatio");
105 
106  fuel_ = ThermoType(thermoDict.subDict("fuel"));
107  oxidant_ = ThermoType(thermoDict.subDict("oxidant"));
108  products_ = ThermoType(thermoDict.subDict("burntProducts"));
109 }
110 
111 
112 template<class ThermoType>
114 (
115  const label speciei
116 ) const
117 {
118  if (speciei == 0)
119  {
120  return fuel_;
121  }
122  else if (speciei == 1)
123  {
124  return oxidant_;
125  }
126  else if (speciei == 2)
127  {
128  return products_;
129  }
130  else
131  {
133  << "Unknown specie index " << speciei << ". Valid indices are 0..2"
134  << abort(FatalError);
135 
136  return fuel_;
137  }
138 }
139 
140 
141 // ************************************************************************* //
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
Foam::egrMixture.
Definition: egrMixture.H:49
error FatalError
A list of keyword definitions, which are a keyword followed by any number of values (e...
Definition: dictionary.H:137
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
Definition: error.H:319
const dictionary & subDict(const word &) const
Find and return a sub-dictionary.
Definition: dictionary.C:633
hashedWordList speciesTable
A table of species as a hashedWordList.
Definition: speciesTable.H:41
void read(const dictionary &)
Read dictionary.
Definition: egrMixture.C:102
const dimensionedScalar b
Wien displacement law constant: default SI units: [m.K].
Definition: createFields.H:27
A class for handling words, derived from string.
Definition: word.H:59
Specialization of the basicSpecieMixture for combustion.
errorManip< error > abort(error &err)
Definition: errorManip.H:131
const ThermoType & getLocalThermo(const label speciei) const
Return thermo based on index.
Definition: egrMixture.C:114
const ThermoType & mixture(const scalar, const scalar, const scalar) const
Definition: egrMixture.C:71
PtrList< volScalarField > & Y
Mesh data needed to do the Finite Volume discretisation.
Definition: fvMesh.H:78
ITstream & lookup(const word &, bool recursive=false, bool patternMatch=true) const
Find and return an entry data stream.
Definition: dictionary.C:451