00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035 #ifndef veryInhomogeneousMixture_H
00036 #define veryInhomogeneousMixture_H
00037
00038 #include "combustionMixture.H"
00039
00040
00041
00042 namespace Foam
00043 {
00044
00045
00046
00047
00048
00049 template<class ThermoType>
00050 class veryInhomogeneousMixture
00051 :
00052 public combustionMixture
00053 {
00054
00055
00056 static const int nSpecies_ = 3;
00057 static const char* specieNames_[3];
00058
00059 dimensionedScalar stoicRatio_;
00060
00061 ThermoType fuel_;
00062 ThermoType oxidant_;
00063 ThermoType products_;
00064
00065 mutable ThermoType mixture_;
00066
00067
00068 volScalarField& ft_;
00069
00070
00071 volScalarField& fu_;
00072
00073
00074 volScalarField& b_;
00075
00076
00077 veryInhomogeneousMixture(const veryInhomogeneousMixture<ThermoType>&);
00078
00079
00080 public:
00081
00082
00083 typedef ThermoType thermoType;
00084
00085
00086
00087
00088
00089 veryInhomogeneousMixture(const dictionary&, const fvMesh&);
00090
00091
00092
00093
00094 const dimensionedScalar& stoicRatio() const
00095 {
00096 return stoicRatio_;
00097 }
00098
00099 const ThermoType& mixture(const scalar, const scalar) const;
00100
00101 const ThermoType& cellMixture(const label celli) const
00102 {
00103 return mixture(ft_[celli], fu_[celli]);
00104 }
00105
00106 const ThermoType& patchFaceMixture
00107 (
00108 const label patchi,
00109 const label facei
00110 ) const
00111 {
00112 return mixture
00113 (
00114 ft_.boundaryField()[patchi][facei],
00115 fu_.boundaryField()[patchi][facei]
00116 );
00117 }
00118
00119 const ThermoType& cellReactants(const label celli) const
00120 {
00121 return mixture(ft_[celli], ft_[celli]);
00122 }
00123
00124 const ThermoType& patchFaceReactants
00125 (
00126 const label patchi,
00127 const label facei
00128 ) const
00129 {
00130 return mixture
00131 (
00132 ft_.boundaryField()[patchi][facei],
00133 ft_.boundaryField()[patchi][facei]
00134 );
00135 }
00136
00137 const ThermoType& cellProducts(const label celli) const
00138 {
00139 scalar ft = ft_[celli];
00140 return mixture(ft, fres(ft, stoicRatio().value()));
00141 }
00142
00143 const ThermoType& patchFaceProducts
00144 (
00145 const label patchi,
00146 const label facei
00147 ) const
00148 {
00149 scalar ft = ft_.boundaryField()[patchi][facei];
00150 return mixture(ft, fres(ft, stoicRatio().value()));
00151 }
00152
00153
00154 void read(const dictionary&);
00155 };
00156
00157
00158
00159
00160 }
00161
00162
00163
00164 #ifdef NoRepository
00165 # include "veryInhomogeneousMixture.C"
00166 #endif
00167
00168
00169
00170 #endif
00171
00172