FreeFOAM The Cross-Platform CFD Toolkit
veryInhomogeneousMixture.H
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) 1991-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 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 
39 
40 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
41 
42 namespace Foam
43 {
44 
45 /*---------------------------------------------------------------------------*\
46  Class veryInhomogeneousMixture Declaration
47 \*---------------------------------------------------------------------------*/
48 
49 template<class ThermoType>
51 :
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)
78 
79 
80 public:
81 
82  //- The type of thermodynamics this mixture is instantiated for
84 
85 
86  // Constructors
87 
88  //- Construct from dictionary and mesh
90 
91 
92  //- Destructor
94  {}
95 
96 
97  // Member functions
98 
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 
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 
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 
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 
162 
163 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
164 
165 } // End namespace Foam
166 
167 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
168 
169 #ifdef NoRepository
171 #endif
172 
173 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
174 
175 #endif
176 
177 // ************************ vim: set sw=4 sts=4 et: ************************ //