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
00036
00037 #ifndef swirlInjector_H
00038 #define swirlInjector_H
00039
00040 #include "injectorType.H"
00041 #include "vector.H"
00042
00043
00044
00045 namespace Foam
00046 {
00047
00048
00049
00050
00051
00052 class swirlInjector
00053 :
00054 public injectorType
00055 {
00056
00057 private:
00058
00059
00060
00061 typedef VectorSpace<Vector<scalar>, scalar, 2> pair;
00062
00063 dictionary propsDict_;
00064
00065 vector position_;
00066 vector direction_;
00067 scalar d_;
00068 scalar mass_;
00069 scalar injectionPressure_;
00070 scalar T_;
00071 label nParcels_;
00072 scalarField X_;
00073 List<pair> massFlowRateProfile_;
00074 List<pair> injectionPressureProfile_;
00075 List<pair> velocityProfile_;
00076 List<pair> CdProfile_;
00077 List<pair> TProfile_;
00078 scalar averageParcelMass_;
00079
00080 bool pressureIndependentVelocity_;
00081
00082
00083
00084 vector tangentialInjectionVector1_, tangentialInjectionVector2_;
00085
00086
00087
00088
00089
00090 swirlInjector(const swirlInjector&);
00091
00092
00093 void operator=(const swirlInjector&);
00094
00095
00096
00097 void setIntegralProfileValue();
00098
00099
00100
00101 void setTangentialVectors();
00102
00103
00104 scalar fractionOfInjection(const scalar time) const;
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117 public:
00118
00119
00120 TypeName("swirlInjector");
00121
00122
00123
00124
00125
00126 swirlInjector
00127 (
00128 const Time& t,
00129 const dictionary& dict
00130 );
00131
00132
00133
00134
00135 ~swirlInjector();
00136
00137
00138
00139
00140
00141 label nParcelsToInject
00142 (
00143 const scalar t0,
00144 const scalar t1
00145 ) const;
00146
00147
00148 const vector position() const;
00149
00150
00151 vector position
00152 (
00153 const scalar time,
00154 const bool twoD,
00155 const scalar angleOfWedge,
00156 const vector& axisOfSymmetry,
00157 const vector& axisOfWedge,
00158 const vector& axisOfWedgeNormal,
00159 Random& rndGen
00160 ) const;
00161
00162
00163 scalar d() const;
00164
00165
00166 const vector& direction() const;
00167
00168
00169 scalar mass
00170 (
00171 const scalar t0,
00172 const scalar t1,
00173 const bool twoD,
00174 const scalar angleOfWedge
00175 ) const;
00176
00177
00178 scalar mass() const;
00179
00180
00181 List<pair> massFlowRateProfile() const;
00182
00183
00184 scalar massFlowRate(const scalar time) const;
00185
00186
00187 List<pair> injectionPressureProfile() const;
00188
00189
00190 scalar injectionPressure(const scalar time) const;
00191
00192
00193 List<pair> velocityProfile() const;
00194
00195
00196 scalar velocity(const scalar time) const;
00197
00198
00199 List<pair> CdProfile() const;
00200
00201
00202 scalar Cd(const scalar time) const;
00203
00204
00205 const scalarField& X() const;
00206
00207
00208 List<pair> T() const;
00209
00210
00211 scalar T(const scalar time) const;
00212
00213
00214 scalar tsoi() const;
00215
00216
00217 scalar teoi() const;
00218
00219
00220 scalar injectedMass(const scalar t) const;
00221
00222 void correctProfiles
00223 (
00224 const liquidMixture& fuel,
00225 const scalar referencePressure
00226 );
00227
00228 bool pressureIndependentVelocity() const
00229 {
00230 return pressureIndependentVelocity_;
00231 }
00232 };
00233
00234
00235
00236
00237 }
00238
00239
00240
00241 #endif
00242
00243