EXP 1-1 Introduction: Difference between revisions
Jan.Jedelsky (talk | contribs) |
Jan.Jedelsky (talk | contribs) |
||
Line 72: | Line 72: | ||
</math>|{{EquationRef|5}}}} | </math>|{{EquationRef|5}}}} | ||
<br/><br/> | <br/><br/> | ||
semi-empirical '''Equation (3)''' in <ref name="Rezaei26"> S. Rezaei, F. Vashahi, G. Ryu, and J. Lee, Fuel 258, 116094 (2019) </ref>, or analytically using linear stability analysis (LISA) <ref name="Arun27"> G. Arun Vijay and N. Shenbaga Vinayaga Moorthi, Journal of Propulsion and Power 32 (2), 448 (2015) </ref>. Countless correlations have been developed to describe droplet size as a function of atomizer operating parameters. The one by Wang and Lefebvre <ref name="Lefebvre28"> X. F. Wang and A. H. Lefebvre, Journal of Propulsion and Power 3 (1), 11 (1987) </ref> calculates <math>D_{32}</math> considering relevant physical phenomena during atomization:<br/><br/> | semi-empirical '''Equation (3)''' in <ref name="Rezaei26"> S. Rezaei, F. Vashahi, G. Ryu, and J. Lee, Fuel 258, 116094 (2019) </ref>, or analytically using linear stability analysis (LISA) <ref name="Arun27"> G. Arun Vijay and N. Shenbaga Vinayaga Moorthi, Journal of Propulsion and Power 32 (2), 448 (2015) </ref>. Countless correlations have been developed for swirl atomizers spraying into the steady environment to describe droplet size as a function of atomizer operating parameters. The one by Wang and Lefebvre <ref name="Lefebvre28"> X. F. Wang and A. H. Lefebvre, Journal of Propulsion and Power 3 (1), 11 (1987) </ref> calculates <math>D_{32}</math> considering relevant physical phenomena during atomization:<br/><br/> | ||
{{NumBlk|:|<math> | {{NumBlk|:|<math> | ||
D_{32}=4.52\left(\frac{\sigma u_{l}^2}{\rho_g p_{i n}^2}\right)^{0.25}\left(\tau_0 \cos (S C A)\right)^{0.25}+0.39\left(\frac{\sigma \varrho_l}{\rho_g p_{i n}}\right)^{0.25}\left(\tau_0 \cos (S C A)\right)^{0.75} | D_{32}=4.52\left(\frac{\sigma u_{l}^2}{\rho_g p_{i n}^2}\right)^{0.25}\left(\tau_0 \cos (S C A)\right)^{0.25}+0.39\left(\frac{\sigma \varrho_l}{\rho_g p_{i n}}\right)^{0.25}\left(\tau_0 \cos (S C A)\right)^{0.75} |
Revision as of 11:20, 24 July 2023
Introduction
The subject of the case is a PSA spray exposed to cross-flowing air. A small low-pressure atomizer was used for the study. This atomizer was developed for spraying aviation fuel Jet A-1 (kerosene) into the combustion chamber of a small gas turbine (GT) engine. The here documented operation conditions of the atomizer and the flow velocity correspond to the engine's low-power or steady-flight conditions. The airflow is forced perpendicularly to the main spraying axis, which is considered a cross-flow case. The flow is homogeneous, isothermal and with low turbulence intensity, . Similar atomizers of this type and size used together with the operating pressure and cross-flow air velocity conditions cover many industrial spray applications ranging from small GT combustors to chemical spray reactors. The conditions are also relevant for agriculture and domestic sprayers. The processed results of the present case were published in [1], with work carried out in the frame of projects №. GA18-15839S and GA 22-17806S funded by Czech Science Foundation and project “Computer Simulations for Effective Low-Emission Energy Engineering” No. CZ.02.1.01/0.0/0.0/16_026/0008392 funded by Operational Programme Research, Development and Education, Priority axis 1: Strengthening capacity for high-quality research. The present case is one of several cases measured and studied in [1]. The data are relevant to CFD engineers and scientists. They can differentiate the crucial phenomena to be considered in their numerical simulations of that disperse two-phase flow case. The modellers can highlight the important features of the complex two-phase flows and use the data for validation purposes. The case is equally interesting to engineers dealing with the processes where the gas–liquid energy transfer and droplet transport are important. The case data can be used for further processing to obtain new findings of the problem, derive empirical models and serve as benchmark data.
Main characteristics of the flow and spray
The PSA sprays water (which represents low viscosity liquid) into cross-flowing air with low turbulence. There are several forces relevant to the case. Cohesive and consolidating forces acting on the liquid film are the surface tension force and the viscosity force . These are counteracted by disruptive compressive and momentum forces . Apart from these also the gravity force acts. We can neglect the other forces possibly acting on the droplets and other liquid structures, such as stochastic force that accounts for Brownian collisions of the droplet with surrounding fluid molecules, or Basset force. The case can be decomposed into several consequent stages with different relevant phenomena, due to the physical actions of these forces, as shown in Figure 2a:
- Liquid flow inside the atomizer, its discharge,
- Sheet formation and the primary break-up of the liquid sheet,
- Liquid secondary break-up and spray formation,
- Interaction of the sprayed liquid with surrounding air: gas–liquid mixing, droplet collisions, droplet clustering, and droplet repositioning.
From a thermodynamic point of view, the case is isothermal and isobaric, except for possible evaporation, which can modify the droplet size [2] and can introduce thermal effects, such as the exchange of heat between the discharged liquid and the surrounding air, which are otherwise unimportant. For the purpose of numerical simulations, the case features a two-way to four-way coupling between the gas and liquid phases depending on the position in the spray [3].
Underlying flow physics which characterise this case
The four stages of this case are explained in the consequent subsections.
Liquid flow inside the atomizer and its discharge
The formation of the liquid film and the resulting spray depend mainly on the internal flow, the geometry of the outlet and the interaction with the surrounding environment being additional factors. The liquid is forced into rotational motion inside the swirl chamber due to the tangentially oriented inlet ports, see Figure 7. The swirling flow reduces the pressure near the atomizer axis, and when it drops below the air pressure at the exit, an air core establishes along the main atomizer axis. The shape and stability of the air core inside the nozzle directly affect the formation, geometrical characteristics and stability of the liquid sheet that emerges from the nozzle [4]. Therefore, the flow field inside the swirl chamber is key for understanding these processes.
The liquid, pumped under pressure through the tangentially oriented inlet ports, creates a swirling flow inside the swirl chamber, see Figure 7. Its main purpose is to form a thin liquid film at the nozzle exit. The swirl momentum is determined by the swirl number, , which represents the ratio of the momentum from the swirl component of the velocity to the axial component
-
(1)
where is the cross-section of the inlet ports, is the fluid flow rate, and are the radius of the swirl chamber and the exit orifice respectively (as shown in Figure 7). Thus, the swirl number depends on the nozzle geometry and can be expressed by the nozzle dimension constant .
The nature of the flow depends on Reynolds number, (ratio of momentum and viscous forces), and insignificantly on Weber number, , (ratio of momentum and surface tension forces), see Equations (11), (12), where is considered as the velocity of the fluid entering the chamber and is its diameter. For larger nozzles and enlarged models, the Froude number:
-
(2)
which describes the effect of gravitational forces on the flow, is also important. Here is the air core radius, and is the gravitational acceleration.
Earlier theoretical works assumed the internal flow as a non-viscous free vortex [5], and with the consideration of Bernoulli's equation for an ideal fluid neglecting the potential term and the radial velocity component, the resulting equation reads:
-
(3)
where is the radius of the inlet ports and is the air core radius in the position of the exit orifice. For the solution, the principle of maximum flow is assumed, i.e., is adjusted so that the flow rate is always maximum: . The non-viscous description of the flow was revised in [5], [6] and recognised as suitable only for understanding the nature of the flow or preliminary nozzle design. Discrepancies of experiments with the non-viscous theory led to corrections of this model [7] [8] [9] [10] [11] and more complex analytical approaches [12] [13] [14], which, however, do not reach the accuracy of CFD models.
The liquid swirling inside the nozzle discharges from the exit orifice at a high velocity into the surrounding air (Figure 2a). The annular liquid structure formed at the orifice features a relatively low discharge coefficient, , which according Rizk and Lefebvre [7] is: for this case. This value agrees well with the experimental data in Table 2. The efficiency of the conversion of inlet potential energy into kinetic energy at the nozzle exit is . This so-called nozzle efficiency was studied and estimated by several authors. Horvay with Leuckel [15] found , Yule with Chinn [16] reported , and we [3] for similarly sized atomizer found .
Formation and primary break-up of the liquid sheet
The conical liquid sheet spreads in the axial and radial direction, attenuates downstream the nozzle and undergoes a dynamic liquid–gas interaction. That depends on the airflow conditions, which can be distinguished into simple categories of still, co-, counter- or crossflowing air.
A high-velocity shear between the discharged liquid and the surrounding air produces Kelvin–Helmholtz-type instabilities on the sheet. These add to turbulent perturbations induced by the swirling motion inside the chamber and deform the sheet. The disrupting gas forces and the consolidating surface tension forces of the liquid film are compared using the gas Weber number, , where the indices and stand for air and liquid, respectively, and is the sheet thickness. A critical Weber number [17] distinguishes domination of long-wave or short-wave growth on the sheet; long waves prevail when and short waves in the opposite case. The actual (see Table 2), compared with , shows that long-wave growth appears at lower air velocity, and the transition to the short-wave happens at higher air velocity. The sheet thickness reduces to its critical value, and the surface tension forces perforate the perturbed sheet. The sheet then disrupts or tears into fragments at the break-up distance.
The internal flow is complex, and the internal disturbances can turbulise the emerging liquid sheet and these disturbances may, depending on their frequency and intensity, reduce the break-up length. This scenario is supported by Sharief et al. [18] and Yule and Chinn [16], in contradiction to the numerical findings of Deng et al. [19]. The primary break-up features a contraction and ordering of detached sheet fragments into irregularly shaped filaments. These, due to the capillary instability [20], break down into single droplets that form a hollow-cone spray. The relative importance of internal viscous and surface tension forces during the sheet disintegration is indicated by the ratio of to of the liquid phase at the discharge orifice after Yule and Dunkley [21] : . The originally two-dimensional sheet breaks down, and its oscillations and mixing with air result in a radial redistribution of the liquid fragments and droplets according to their size classes. The spray acquires a Gaussian velocity profile normal to the sheet surface [22]. The moving liquid film, fragments, and droplets experience mechanical interactions with the air through viscous drag. The droplets, moving with low , typically below 100, decelerate according to Stokes' law as , and establish a positive size–velocity correlation which contribute to droplet collisions in the dense spray region [23]. The gas–liquid interaction is described in detail in [3].
Spray formation
The produced droplets cover a wide size range and form a single or double-peak size distribution (in Figure 2b), depending on the position in the spray. The droplets with 20 μm decelerate fast to the airflow velocity, medium-size droplets (20 μm 50 μm) feature a positive size–velocity correlation, and the largest droplets up to 100 µm keep the original velocity of the discharged liquid. The smallest droplets follow the local air velocity closely, and so the velocity of droplets sized below 5 µm can serve as the air velocity estimate.
The size distribution of the droplets can be represented simply at each position by a suitable mean droplet diameter, for which a general expression is
-
(4)
where is the diameter of individual droplet and is the total number of droplets at the position. The most often used mean droplet diameters are: , which is the arithmetic mean diameter (here and , this diameter is used for comparison of disperse systems), (surface mean diameter for vaporization studies) and is called the Sauter mean diameter (or volume/surface mean diameter) and this one is used for mass and heat transfer evaluations. The spray itself, if sprayed into still air can be considered roughly axially symmetrical with large size and velocity variability in the radial direction. The radial profiles of the mean liquid velocity are self-similar along the axial locations with a peak close to the sheet position. The sprayed mass is mostly distributed along the sheet trajectory, and it forms a hollow-cone spray. The inner region contains only small droplets that are driven there due to the air drag. The main semi-conical spray region behind the disintegrated liquid sheet contains larger high-energetic droplets with high penetration ability. The outer spray periphery covers a small portion of droplets with velocity decreasing with radial distance.
The liquid sheet, its fractions and larger droplets in the near-nozzle area follow the trajectory given by the discharge conditions, while the ensuing flow and motion of small droplets in the far field are more influenced by the interaction with the surrounding gas.
One of the PSA's main parameters is the break-up distance, , which determines the volume of the ligaments and the size of the resulting droplets. It can be determined from the empirical Equation (5) [24]:
-
(5)
semi-empirical Equation (3) in [25], or analytically using linear stability analysis (LISA) [26]. Countless correlations have been developed for swirl atomizers spraying into the steady environment to describe droplet size as a function of atomizer operating parameters. The one by Wang and Lefebvre [27] calculates considering relevant physical phenomena during atomization:
-
(6)
Available correlations are not fully reliable or universal, so further experiments are required.
Interaction of the sprayed liquid with the surrounding air
All the above descriptions and most published studies have considered atomizer spraying in the absence of ambient flow. Though liquid spraying into still surrounding air is the most frequently investigated configuration, many applied atomizers work in a flowing environment. Cross-flow spray configuration is relevant, e.g. for Venturi scrubbers [28] where typical velocities range from 10 m/s [29] or 20 m/s [30] up to 30 m/s and even 80 m/s [31] [32] [33].
Droplet size at higher ambient flow velocities is affected by the secondary aerodynamic break-up. If we consider a maximum droplet size in the spray of 100 µm, then after [34]
-
(7)
a relative velocity between the droplet and surrounding gas, , of about 78 m/s is required for the aerodynamic break-up of a droplet. Thus, the ambient flow velocities tested here do not allow for the secondary break-up to apply. Their resulting size depends on the primary break-up and the design of the atomizer itself. The combination of the forces , , and acting on the liquid sheet results in its disintegration [35]. The presence of transverse flow disturbs the flow field around the atomizer and can change compressive force and momentum force which are derived in detail in [35].
-
(8)
-
(9)
Here is the wave number, is the velocity of the liquid sheet, and is the velocity of the air on the front and rear side of the sheet respectively. With the presence of cross-flow, the air velocity around the liquid sheet changes, which then affects the acting on the liquid sheet, see Figure 3. The cross-flow represents an "additional resistance" of the surrounding environment that the liquid sheet must overcome. With this also, the rate of change of momentum of the liquid film and the ratios in Equation (9) change. The increased and in the cross-flow reduce the break-up distance, , and leads to the formation of larger droplets, as observed in [36].
The studies dealing with sprays in cross-flow express the effect of the ambient flow on the spray by the ratio of the liquid and air momentum (), the aerodynamic Weber number () and the relative Weber number () using Equations ((10) – (12)).
The in Equation (10) denotes the velocity of the liquid sheet at the discharge point and it is calculated after [37]. The Weber number is defined in two ways here. The first definition, used in [37], incorporates the cross-flow velocity () and the diameter of the discharge orifice (). The second one, , contains the relative velocity of the liquid sheet () to the cross-flow velocity (), which is denoted . The determination of is shown in Figure 4 and its calculation given by Equation (13), where is the spray cone angle.
-
(10)
-
(11)
-
(12)
-
(13)
![]() |
Figure 4: Graphical representation of in cross-flow, the right image is a magnified view of the velocity vectors |
Main quantities of interest
The PDA measurements produce data allowing the calculation of droplet size and velocity statistics and, to some extent estimating the local airflow velocity. These data give detailed information on the velocity field of the sprayed liquid and surrounding air. The HSV provides photogrammetric information on the discharged liquid. The data can be used for the estimation of relevant dimensionless criteria that characterise the individual processes involved in the studied case, as summarised in Table 1. The table also contains information on experimental and simulation techniques and approaches used and applicable to study these processes by different researchers.
Process |
Output, parameters questioned |
Relevant criteria |
Approaches | |||
---|---|---|---|---|---|---|
Experiment |
Simulation | |||||
1 |
Internal flow |
Velocity field, air core properties |
, , |
LDA,HSV |
Laminar, URANS, LES | |
2 |
Discharge and liquid film formation |
, , velocity, stability, liquid film thickness |
, , , |
HSV | ||
3 |
Break-up into smaller structures (primary) |
break-up character, |
, , , , |
LIF |
HSV |
(LES) DNS |
4 |
Subsequent disintegration into droplets (secondary) |
Droplet size, velocity, concentration |
, , , |
PDA, HSV | ||
5 |
Interaction of droplets with the surrounding environment and with each other |
Character of interaction, energy transfer, droplet collision, evaporation |
, , , , |
URANS, LES, Stat |
Laser Doppler anemometry, phase Doppler anemometry, high-speed visualisation, laser-induced fluorescence, collision (for definition see [2]), concentration of droplets in spray, statistical approaches, Bond number (also called Eötvös number) , can only be significant at very low discharge velocity, is the difference in density between the liquid and the gas, Ohnesorge number is the characteristic dimension, , and are defined similarly to , and in Equations (11) (12) (14) with appropriate characteristic lengths and velocities and the index and denote the gas and liquid respectively. Froude number according to Equation (2). Stokes number is the difference between the gas and droplet velocity, is characteristic distance.
References
- ↑ 1.0 1.1 O. Cejpek, M. Maly, J. Slama, M. M. Avulapati, and J. Jedelsky, Continuum Mechanics and Thermodynamics 34 (6), 1497 (2022)
- ↑ 2.0 2.1 J. Jedelský, M. Malý, S. K. Vankeswaram, M. Zaremba, R. Kardos, D. Csemány, A. Červenec, and V. Józsa, (http://dx.doi.org/10.2139/ssrn.4385285)
- ↑ 3.0 3.1 3.2 J. Jedelsky, M. Maly, N. Pinto del Corral, G. Wigley, L. Janackova, and M. Jicha, International Journal of Heat and Mass Transfer 121, 788 (2018)
- ↑ M. Maly, O. Cejpek, M. Sapik, V. Ondracek, G. Wigley, and J. Jedelsky, Experimental Thermal and Fluid Science 120, 110210 (2021)
- ↑ 5.0 5.1 A. J. Yule and J. Chinn, presented at the International Conference on Liquid Atomization and Sprays, ICLASS-94, Rouen, France, 1994 (unpublished).
- ↑ J. J. Chinn, Atomization and Sprays 19 (3) (2009); J. J. Chinn, Atomization and Sprays 19 (3) (2009)
- ↑ 7.0 7.1 N. K. Rizk and A. H. Lefebvre, Journal of Propulsion and Power 1 (3), 193 (1985)
- ↑ A. Jones, presented at the Proceedings of the Second International Conference on Liquid Atomization and Spray Systems, 1982
- ↑ J. Ballester and C. Dopazo, Atomization and sprays 4 (3) (1994)
- ↑ M. Benjamin, A. Mansour, U. Samant, S. Jha, Y. Liao, T. Harris, and S. Jeng, presented at the ASME 1998 International Gas Turbine and Aeroengine Congress and Exhibition, 1998
- ↑ A. Sakman, M. Jog, S. Jeng, and M. Benjamin, AIAA journal 38 (7), 1214 (2000)
- ↑ L. Craig, N. Barlow, S. Patel, B. Kanya, and S. P. Lin, Atomization and Sprays 19 (12), 1113 (2009)
- ↑ M. Malý, L. Janáčková, J. Jedelský, M. Jícha, R. Lenhard, and K. Kaduchová, presented at the AIP Conference Proceedings, 2016
- ↑ G. Amini, International Journal of Multiphase Flow 79, 225 (2016)
- ↑ M. Horvay and W. Leuckel, German chemical engineering 9 (5), 276 (1986)
- ↑ 16.0 16.1 A. Yule and J. Chinn, Atomization and Sprays 10 (2), 121 (2000)
- ↑ P. K. Senecal and D. P. Schmidt, IRutland, C.JReitz, R.DCorradini, M.L, International Journal of Multiphase Flow 25 (6–7), 1073 (1999)
- ↑ R. Sharief, J. Jeong, and D. James, Atomization and Sprays 10 (6) (2000)
- ↑ H.-Y. Deng, F. Feng, and X.-S. Wu, Atomization and Sprays 26 (4) (2016)
- ↑ E. Villermaux, New Journal of Physics 6 (1), 125 (2004)
- ↑ A. J. Yule and J. J. Dunkley, Atomization of Melts: For Powder Production and Spray Deposition. (Oxford University Press, USA, 1984)
- ↑ H. Liu, Science and engineering of droplets: fundamentals and applications. (Noyes Publications; William Andrew Pub., Park Ridge, N.J (1999)
- ↑ J. L. Santolaya, J. A. García, E. Calvo, and L. M. Cerecedo, International Journal of Multiphase Flow 56 (0), 160 (2013)
- ↑ T. Arai and H. Hashimoto, Transactions of the Japan Society of Mechanical Engineers. B 51, 3336 (1985)
- ↑ S. Rezaei, F. Vashahi, G. Ryu, and J. Lee, Fuel 258, 116094 (2019)
- ↑ G. Arun Vijay and N. Shenbaga Vinayaga Moorthi, Journal of Propulsion and Power 32 (2), 448 (2015)
- ↑ X. F. Wang and A. H. Lefebvre, Journal of Propulsion and Power 3 (1), 11 (1987)
- ↑ C. M, H. R, G. S, and J. R. Coury, Brazilian Journal of Chemical Engineering 21 (2004)
- ↑ T. Mi and X. M. Yu, Chemical Engineering and Processing: Process Intensification 62, 159 (2012)
- ↑ N. Abbaspour, M. Haghshenasfard, M. R. Talaei, and H. Amini, Journal of Molecular Liquids 303, 112689 (2020)
- ↑ D. Breitenmoser, P. Papadopoulos, T. Lind, and H.-M. Prasser, International Journal of Multiphase Flow 142, 103694 (2021)
- ↑ M. Ali, Y. Qi, and K. Mehboob, Research Journal of Applied Sciences, Engineering and Technology 4 (2012)
- ↑ J. A. Gonçalves, M. A. Costa, M. L. Aguiar, and J. R. Coury, J Hazard Mater 116 (1-2), 147 (2004)
- ↑ W. R. Lane, Industrial & Engineering Chemistry 43 (6), 1312 (1951)
- ↑ 35.0 35.1 N. Dombrowski and W. R. Johns, Chemical Engineering Science 18, 203 (1963)
- ↑ S. Lee, W. Kim, and W. Yoon, Journal of Mechanical Science and Technology 24 (2), 559 (2010)
- ↑ 37.0 37.1 R. Surya Prakash, H. Gadgil, and B. N. Raghunandan, International Journal of Multiphase Flow 66, 79 (2014)
Contributed by: Ondrej Cejpek, Milan Maly, Ondrej Hajek, Jan Jedelsky — Brno University of Technology
© copyright ERCOFTAC 2024