Newer
Older
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software; you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by the
Free Software Foundation; either version 2 of the License, or (at your
option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM; if not, write to the Free Software Foundation,
Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
Class
Foam::spray
Description
A spray is a cloud of parcels
\*---------------------------------------------------------------------------*/
#ifndef spray_H
#define spray_H
#include "parcel.H"
#include "injector.H"
#include "IOPtrList.H"
#include "interpolation.H"
#include "liquid.H"
#include "autoPtr.H"
#include "liquidMixture.H"
#include "Random.H"
#include "thermoPhysicsTypes.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
namespace Foam
{
class atomizationModel;
class breakupModel;
class collisionModel;
class dispersionModel;
class dragModel;
class evaporationModel;
class injectorModel;
class heatTransferModel;
class wallModel;
Andrew Heather
committed
class basicMultiComponentMixture;
/*---------------------------------------------------------------------------*\
Class spray Declaration
\*---------------------------------------------------------------------------*/
class spray
:
public Cloud<parcel>
{
// Private data
// References to the database and meshes
const Time& runTime_;
scalar time0_;
const fvMesh& mesh_;
//- Random number generator
Random rndGen_;
//- Acceleration due to gravity
const vector& g_;
// References to the physical fields
const volVectorField& U_;
const volScalarField& rho_;
const volScalarField& p_;
const volScalarField& T_;
//- The spray properties
IOdictionary sprayProperties_;
scalar ambientPressure_;
scalar ambientTemperature_;
IOPtrList<injector> injectors_;
// References to the spray sub-models
autoPtr<atomizationModel> atomization_;
autoPtr<dragModel> drag_;
autoPtr<evaporationModel> evaporation_;
autoPtr<heatTransferModel> heatTransfer_;
autoPtr<wallModel> wall_;
autoPtr<breakupModel> breakupModel_;
autoPtr<collisionModel> collisionModel_;
autoPtr<dispersionModel> dispersionModel_;
autoPtr<liquidMixture> fuels_;
autoPtr<injectorModel> injectorModel_;
//- Minimum number of lagrangian subcycles
const label subCycles_;
// Composition properties
const PtrList<gasThermoPhysics>& gasProperties_;
Andrew Heather
committed
const basicMultiComponentMixture& composition_;
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
List<label> liquidToGasIndex_;
List<label> gasToLiquidIndex_;
List<bool> isLiquidFuel_;
// Necessary 2D-information
bool twoD_;
vector axisOfSymmetry_;
vector axisOfWedge_;
vector axisOfWedgeNormal_;
scalar angleOfWedge_;
// Interpolation
dictionary interpolationSchemes_;
autoPtr<interpolation<vector> > UInterpolator_;
autoPtr<interpolation<scalar> > rhoInterpolator_;
autoPtr<interpolation<scalar> > pInterpolator_;
autoPtr<interpolation<scalar> > TInterpolator_;
// Spray Source Terms
//- Momentum
vectorField sms_;
//- Enthalpy
scalarField shs_;
//- Mass
PtrList<scalarField> srhos_;
//- The total mass of the injected liquid
scalar totalInjectedLiquidMass_;
//- The (total added) injected kinetic energy of the liquid
scalar injectedLiquidKE_;
// Private Member Functions
//- Disallow default bitwise copy construct
spray(const spray&);
//- Disallow default bitwise assignment
void operator=(const spray&);
public:
// Constructors
//- Construct from components
spray
(
const volVectorField& U,
const volScalarField& rho,
const volScalarField& p,
const volScalarField& T,
Andrew Heather
committed
const basicMultiComponentMixture& composition,
const PtrList<gasThermoPhysics>& gasProperties,
const dictionary& thermophysicalProperties,
const dimensionedVector& g
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
);
// Destructor
~spray();
// Member Functions
// Spray tracking and evolution functions
//- Evolve the spray (move, inject and breakup)
void evolve();
//- Move the spray parcels
void move();
//- Inject more parcels
void inject();
//- Primary breakup droplets
void atomizationLoop();
//- Secondary breakup droplets
void breakupLoop();
// Access
inline const Time& runTime() const;
inline const fvMesh& mesh() const;
inline const volVectorField& U() const;
inline const volScalarField& rho() const;
inline const volScalarField& p() const;
inline const volScalarField& T() const;
inline PtrList<injector>& injectors();
inline const PtrList<injector>& injectors() const;
inline const atomizationModel& atomization() const;
inline const breakupModel& breakup() const;
inline const collisionModel& collisions() const;
inline const dispersionModel& dispersion() const;
inline const dragModel& drag() const;
inline const evaporationModel& evaporation() const;
inline const heatTransferModel& heatTransfer() const;
inline const injectorModel& injection() const;
inline const wallModel& wall() const;
inline tmp<volVectorField> momentumSource() const;
inline tmp<volScalarField> evaporationSource(const label i) const;
inline tmp<volScalarField> heatTransferSource() const;
inline Random& rndGen();
inline label subCycles() const;
inline const vector& g() const;
inline const liquidMixture& fuels() const;
inline const PtrList<gasThermoPhysics>& gasProperties() const;
Andrew Heather
committed
inline const basicMultiComponentMixture& composition() const;
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
inline const List<label>& liquidToGasIndex() const;
inline const List<label>& gasToLiquidIndex() const;
inline const List<bool>& isLiquidFuel() const;
inline const bool& twoD() const;
inline const vector& axisOfSymmetry() const;
inline const vector& axisOfWedge() const;
inline const vector& axisOfWedgeNormal() const;
inline const scalar& angleOfWedge() const;
inline const interpolation<vector>& UInterpolator() const;
inline const interpolation<scalar>& rhoInterpolator() const;
inline const interpolation<scalar>& pInterpolator() const;
inline const interpolation<scalar>& TInterpolator() const;
inline vectorField& sms();
inline const vectorField& sms() const;
inline scalarField& shs();
inline const scalarField& shs() const;
inline PtrList<scalarField>& srhos();
inline const PtrList<scalarField>& srhos() const;
inline const scalar& ambientPressure() const;
inline const scalar& ambientTemperature() const;
// Check
//- Returns the liquid mass that has been injected
scalar injectedMass(const scalar t) const;
//- Returns the liquid mass that will be injected by the injectors
scalar totalMassToInject() const;
//- Returns the injected enthalpy
scalar injectedEnthalpy(const scalar t) const;
//- Returns current total liquid mass in the domain
scalar liquidMass() const;
//- Returns the enthalpy of all the liquid in the domain
// Hdrop = Hgas - Hlat
scalar liquidEnthalpy() const;
//- Returns the enthalpy (total) of all the liquid in the domain
// Hdrop = Hgas - Hlat + (P-Psat)/rhoDrop;
scalar liquidTotalEnthalpy() const;
//- Returns the kinetic energy of the liquid phase
scalar liquidKineticEnergy() const;
//- Returns the injected kinetic energy of the liquid phase
scalar injectedLiquidKineticEnergy() const;
//- Returns the droplet penetration for 'prc' percent of the
// liquid from nozzle 'nozzlei'
scalar liquidPenetration
(
const label nozzlei,
const scalar prc
) const;
//- Returns the droplet penetration for 'prc' percent of the
// liquid from nozzle 0
scalar liquidPenetration(const scalar prc) const;
//- Return Sauter Mean Diameter
scalar smd() const;
//- Return Maximum Diameter
scalar maxD() const;
//- Return Ambient Pressure
void calculateAmbientPressure();
//- Return Ambient Temperature
void calculateAmbientTemperature();
// I/O
//- Write fields
};
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
} // End namespace Foam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#include "sprayI.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#endif
// ************************************************************************* //