![Arxiv:1907.02962V1 [Physics.Flu-Dyn] 3 Jul 2019 Cr a Critical Weber Number](http://data.docslib.org/img/bfc21c48633fcfac39c343dd582555e2-1.webp)
Deformation and breakup of droplets in an oblique continuous air stream
Surendra Kumar Soni†, Pavan Kumar Kirar††, Pankaj Kolhe†∗, and Kirti Chandra Sahu††∗ †Department of Mechanical and Aerospace Engineering, Indian Institute of Technology Hyderabad, Sangareddy 502 285, Telangana, India ††Department of Chemical Engineering, Indian Institute of Technology Hyderabad, Sangareddy 502 285, Telangana, India
We experimentally investigate the deformation and breakup of droplets interacting with an oblique continuous air stream. A high-speed imaging system is employed to record the trajectories and topological changes of the droplets of different liquids. The droplet size, the orientation of the air nozzle to the horizontal and fluid properties (surface tension and viscosity) are varied to study different breakup modes. We found that droplet possessing initial momentum prior to entering the continuous air stream exhibits a variation in the required Weber number for the vibrational to the bag breakup transition with a change in the angle of the air stream. The critical Weber numbers (W ecr) for the bag-type breakup are obtained as a function of the E¨otv¨osnumber (Eo), angle of inclination of the air stream (α) and the Ohnesorge number (Oh). It is found that although the droplet follows a rectilinear motion initially that transforms to a curvilinear motion at later times when the droplet undergoes topological changes. The apparent acceleration of the droplet and its size influence the critical Weber number for the bag breakup mode. The departure from the cross- flow arrangement shows a sharp decrease in the critical Weber number for the bag breakup which asymptotically reaches to a value associated with the in-line (opposed) flow configuration for the droplet breakup.
I. INTRODUCTION
Fragmentation of liquid droplets due to the flow of high-speed air stream is observed in many industrial applications [1] and natural phenomena [2]. In atomization, initially fuel jet forms ligaments and breaks into droplets (primary atomization). These droplets become unstable when the aerodynamic force overcomes the influence of the surface tension force and the viscous force, and subsequently break into smaller satellite droplets (secondary atomization). The surface area to volume ratio significantly increases in the secondary atomization process, which in turn increases its efficiency (see e.g. Refs. [3, 4]). Thus fragmentation of a liquid droplet into tiny satellite droplets in in-line and cross-flow configurations has been a subject of research for the past several decades (see for instance [2, 5–10]). Previous experiments suggest that there exists different types of mode when a liquid droplet interacts with high speed air stream, namely, the vibration, bag, bag-stamen, dual-bag, multi-mode, shear and catastrophic breakup modes [7, 11–13] owing to different mechanisms and flow characteristics. The Weber number that measures of the relative importance of the inertia over the surface tension force is mainly used to demonstrate different breakup 2 modes. The Weber number is given by W e ≡ ρaU D/σ. The other dimensionless numbers which have been used to study droplet breakup phenomenon are the Reynolds number, Re ≡ ρlUD/µl (inertia force/viscous force), the E¨otv¨osnumber ≡ ( − ) 2 (gravitational force/surface tension force) and the Ohnesorge number, ≡ √ √Eo ρl ρa gD /σ Oh W e/Re = µl/ ρlσD, wherein, σ is the interfacial tension, g is the acceleration due to gravity, D is the diameter of the droplet, ρa and ρl are the density of the air and liquid, respectively, and U is the average velocity of the air stream. At low Weber numbers, a droplet undergoes shape oscillations at a certain frequency (vibrational regime). As the Weber number is increased slowly by increasing the aerodynamic force, and keeping the surface tension force constant, the droplet exhibits a transition from the vibrational mode to the bag breakup mode. The value of the Weber number at which this transition occurs (i.e. the Weber number at which the droplet just starts to form a bag) is termed as the critical Weber number (W e ). Taylor [5] was the first to predict that a liquid droplet undergoes breakup above arXiv:1907.02962v1 [physics.flu-dyn] 3 Jul 2019 cr a critical Weber number. Since then many researchers [6–10] have reported different values of the critical Weber number in different flow configurations. In cross-flow configuration, Table I presents the range of the Weber numbers associated with the bag breakup region obtained by previous researchers. The discrepancies in obtaining the critical Weber number can be attributed to the definition of the breakup regime and the different experimental facilities (shock tubes, wind tunnels, convergent nozzles, etc.) used in their investigations. Moreover, the fragmentation process, due to its inherent nature, is highly sensitive to the flow conditions. Among all the fragmentation regimes described above,
∗Electronic address: [email protected];[email protected] 2
Table I: The Weber number range for the bag type breakup. References Weber number range Pilch and Erdman [8] 12 < W e ≤ 50 Guildenbecher et al. [11] 11 < W e < 35 Krzeczkowski [18] 10 < W e < 18 Jain et al. [6] 12 < W e < 24 Hsiang and Faeth [19] W e ≤ 11 ± 2 Wierzba [17] 13.7 < W e < 14.07 Kulkarni and Sojka [14] 12 < W e < 16 Wang et al. [20] 10 < W e < 35 the bag breakup regime has the widest range of applications since it occurs at relatively low Weber numbers (see for instance, Refs. [2, 6, 7, 11, 12, 14, 15]). The bag breakup occurs at the beginning of the secondary atomization process and shows similar characteristics with that of the bag-stamen, the multi-bag and the shear stripping breakup modes [6]. Hanson et al. [16] and Wierzba [17] demonstrated that the critical Weber number decreases as the droplet size increases in cross-flow configuration. Their results are summarised in Table II. Three types approaches have been used to study the droplet breakup phenomena: (i) the shock tube approach [18, 19, 21, 22] (ii) the continuous air jet method [6, 12, 14, 15, 23] and (iii) the droplet tower approach [2]. In the present study, we implement a continuous air jet method. The main components of this method are an air nozzle that can produce a top-hat velocity profile and a compressed air measuring and controlling device. The results obtained using the first and the second approaches can be made similar by ensuring that the droplet interacts in the continuous air stream. This is possible when the time taken by the droplet to cross the shear boundary layer is much smaller than the resident time of the droplet when the droplet deforms and breaks [11].
Table II: The critical Weber number for different droplet sizes reported by Hanson et al. [16] and Wierzba [17].
References Fluids D (mm) W ecr 0.426 6.62 0.285 7.18 Silicone oil (10 cSt) 0.180 7.93 0.117 7.94 0.531 10.4 0.273 11.8a Silicone oil (50 cSt) 0.210 13.5 0.126 15.4 0.540 13.1 Hanson et al. [16] 0.338 14.3 0.239 15.5 0.25 16.2 Silicone oil (100 cSt) 0.185 22.6 0.150 23.8 0.391 4.79 0.180 6.37 Water 0.0945 7.14 0.471 6.76 Methyl alcohol 0.186 7.45 Wierzba [17] Water 2.22-3.9 13.7 to 14.07
As the above review indicates all the previous studies considered small droplets either in the cross-flow (e.g. [16, 17] or in the in-line (e.g. [2, 24]) configurations. However, in many situations, such as swirling spray applications and during rain, droplets may interact with the continuous air stream at an angle. In such situations, apart from the aerodynamic, the viscous and the surface tension forces, the gravitational force also influences the droplet breakup p morphology when the droplet size is comparable or larger than the capillary length scale, given by lc = σ/(ρl − ρa)g [25]. When D ≥ lc, the surface tension force will not be able to hold the drop in a spherical shape and the droplet will deform due to the effect of gravity, and thereby changing the critical force requirement for different types of breakup. To the best of our knowledge, such a study has not been conducted yet. In the present work, we experimentally investigate the deformation and breakup of a droplet interacting with an air stream flowing at an angle, α to the horizontal direction. A high-speed imaging system is employed to record the 3 trajectory and topological change of the droplet and the effect of obliquity of the air stream on the droplet breakup process has been investigated. The droplet size, the orientation of the air nozzle, fluid properties (surface tension and viscosity) are varied to study different breakup modes. The critical Weber number has been obtained for different values of the E¨otv¨osnumber (Eo), the angle of inclination of the air stream (α) and the Ohnesorge number (Oh). It is found that, although the droplet exhibits rectilinear motion during its entry into the oblique air stream and prior to the formation of the bag breakup mode, a curvilinear motion is depicted by the droplet at later times, which indicates the variation in forces acting on the droplet. The apparent acceleration experienced by the droplet and its size are linked with the variation in the critical Weber number requirement for the bag breakup mode. The departure from the cross-flow arrangement shows a sharp decrease in the critical Weber number for the bag breakup with asymptotic behaviour at higher oblique angles. At high obliquity, the critical Weber number approaches a value of 6 as reported in the literature (see for instance, Ref. [2]) for the inline (opposed) flow configuration for the droplet breakup. It is also noted that in the present study, the effort has been made to give a detailed description of the experimental procedure, the cautions taken during the experiments, the uncertainty analysis that confirms the repeatability of the experimental results. The rest of the paper is organized as follows. A brief description of the mathematical analysis is given in Section II, the experimental procedure is discussed in Section III, which is followed by the discussion of our experimental results in Section IV. Finally, concluding remarks are given in Section V.
II. MATHEMATICAL FORMULATION
A schematic diagram showing a liquid droplet freely falling under the action of gravity and subjected to an air stream at an angle α to the horizontal is shown in Fig. 1. In the crossed flow configuration, α = 0, such that the air stream flows in the horizontal (X) direction and the droplet falls in the vertical (Y ) direction. For any π 3 oblique angle, α of the air stream, the various forces acting on the droplet are the gravity force (Fg = 6 D ρlg), the π 3 π 3 dud 1 2 buoyant force (FB = 6 D ρag), the inertia force (Fi = 6 D ρa dt ) and the drag force (FD = 2 CdρaU A), such that Fi = Fg − Fb + FD. Here ud is the velocity of the droplet, A is the projected area of the deformed droplet and Cd is drag coefficient and t is time. Buoyancy force
F
F
Air stream F
Y F
↵
Figure 1: Schematic diagram showing a liquid droplet freely falling under the action of gravity and subjected to an air stream at an angle α to the horizontal.
In the absence of the air stream, balancing the gravity force, Fg with the surface tension force, Fs = πDσ, one can get the following expression for maximum droplet size, Dc for which the surface tension dominates the other forces acting on the droplet while maintaining its sphericity. s 6σ Dc = . (1) (ρl − ρa)g
In case of an oblique air stream, the droplet follows a curvilinear motion and experiences an apparent weight due 4 to the rate of change of its acceleration. In this case, the resultant acceleration acting on the droplet (ba) is given by F F F sinα F cosα = i = + D = + D + D (2) ba π 3 g π 3 g π 3 π 3 6 D ρl 6 D ρl 6 D ρl 6 D ρl | {z } | {z } Vertical component Horizontal component FD sinα The net resultant acceleration have two components: the vertical, g + π 3 and the horizontal component, D ρl 6 FD cosα ◦ π 3 . Thus, for the cross-flow condition (α = 0 ), the the vertical and the horizontal components of the D ρl 6 FD ◦ acceleration are g and π 3 , respectively. For the in-line (counter) flow configuration (α = 90 ), the component D ρl 6 FD of acceleration in the vertical direction is g + π 3 and the droplet does not encounter any acceleration in the 6 D ρl horizontal direction. Using the apparent weight of the droplet, the modified critical droplet diameter to maintained its sphericity can be obtained as s 6σ Dcc = . (3) (ρl − ρa)ba
Thus for a droplet with D > Dcc, the surface tension is weaker than the external force and the droplet deforms from its spherical shape. Now we like to turn the attention to the time scales used in the previous studies. Different researchers have used different time scales to non-dimensionalise their results [18, 26]. Most of the researchers non-dimensionalised the breakup time with the transport time (a ratio of droplet diameter to its velocity) as suggested by Nicholls and Ranger [27]. Such nondimensionalisation provides an opposite trend as compared to the actual observed breakup time in case of oblique flow configurations. Considering the droplet breakup process being dictated by forces acting on the droplet, we propose and use the following time scale, ts ≡ ρaU/(ρl − ρa)g to nondimensionalise time, in the present study.
III. EXPERIMENTAL SET-UP AND PROCEDURE
5 1. High speed camera 2. Air nozzle 3. Diffuser sheet 4. Light source 5. Dispensing needle H Y X Z 1 4 α 2 3
Figure 2: Schematic of the experimental set-up. It consists of a high-speed camera, a nozzle, a diffuser sheet, a light source and a dispensing needle. The angle of inclination of the air nozzle to the horizontal is denoted by α. H (≈ 40 mm) is the height of the needle tip from the centre of the nozzle. Gravity, g acts in the negative Y direction. The orientation of the nozzle and the deflection of the nozzle is shown in X − Y plane. 5
Figure 2 depicts the schematic of the experimental set-up used in the present investigation. A metallic circular nozzle has been used to create a continuous and uniform air jet flow. This nozzle is hinged in such a way that the orientation of the nozzle can be fixed to any angle, α to the horizontal. A set of honeycomb pallet is placed near the upstream end of the nozzle to straighten the flow and to minimise disturbances. The internal section of the nozzle is designed to obtain a top-hat velocity profile at the nozzle exit, which has been verified with the aid of a pitot tube and differential head of the manometer measurements. The internal and external diameters of the pitot tube are 1.5 mm and 4 mm, respectively. The horizontal distance between the nozzle and the pitot tube is kept within 10 mm (< 0.5 diameter of the nozzle). For every experiment, it is ensured that the droplet enters in the continuous velocity profile of air or the potential jet core in the immediate downstream region of the nozzle exit before undergoing deformation and breakup. The air flow rate has been measured and controlled using an accurate and reliable ALICAT digital mass flow controller (range 0 to 2000 liter per minute) which is installed before the air nozzle. For generating the same size of liquid droplets, a regular infusion pump (SP-810 from Medical Point) is used. The drop size is varied in two ways: (i) by changing the mass flow rate setting in the infusion pump and (ii) by changing the spout having different exit diameters. The droplet generator and the circular nozzle are arranged in the cross-flow configuration (α = 0◦). The liquid droplet drips from the tip of the spout under the action of gravity and enters the air stream. A high-speed camera (Phantom VEO 710) with a 50 mm Nikon lens is used to capture the droplet deformation and breakup phenomena occurring near the nozzle. The images are recorded at 15000 frames per second (fps) with an exposure time of 10 µs for all the experiments conducted in the present study. The camera is connected to a computer and is operated via Phantom Camera Control (PCC) 3.1 software. The camera has a built-in memory card that allows for direct recording of the images. A diffused back-lit illumination (Light Emitting Diode (LED) based lighting, model 900445, 12000 lm, Visual Instrumentation Corporation) shadowgraph technique has been used to illuminate the background for droplet identification. In this technique, a light source is placed behind a diffusion screen to diffuse the backlit into uniform light. The recorded images are analysed using the Matlab software. The geometric size of the droplet has been calculated by image processing. The droplet deforms while suspended in syringe spout and no longer remains perfectly spherical, hence the equivalent droplet initial diameter has been calculated as [28]
1 2 3 D = (Da × Db) , (4) where Da and Db are droplet diameter corresponding to the major axis and minor axis of the slightly deformed droplet. The droplet velocity is found to be insignificant as compared to air velocity that allows us to use the absolute velocity of the air stream instead of the relative velocity in our calculations. Note that free jet is employed to investigate the droplet deformation and breakup. The uniform flow requirement for droplet breakup study is ensured through droplet entering the potential core region immediate downstream of the nozzle exit. For oblique air stream angles, if sufficient care is not taken, the droplet deformation and breakup process will occur in the shear layer instead of the jet potential core region. This will result in a higher critical Weber number as compared to the actual value. Thus, in the present study, sufficient care has been taken to ensure that droplet enters the potential core region of the air stream. In our experiments, the uniform size of droplets is created with the help of a syringe pump. For calculating the dimensionless numbers, the standard properties of the fluids are taken at 30◦C. To see the effect of surface tension and viscosity on the droplet breakup phenomenon, four different fluids, namely, deionised water, 50 %, 70 % and 80 % aqueous mixtures of glycerin (by volume) have been considered. The viscosities of the fluids are measured using a rotational rheometer (Model: RheolabQC). The density and surface tension of the fluids are taken from Ref. [29]. The density (kg/m3), the surface tension (N/m) and the dynamic viscosity (mPa·s) of the working fluids are given in Table III.
Table III: Physical properties of the working fluids. Working fluid Dynamic viscosity Surface tension Density Droplet diameter 3 µl (mPa·s) σ (mN/m) ρl (kg/m ) D (mm) Water 0.78 72 998 2.0-6.3 Glycerine (50%) 4.85 67.5 1127 2.5 -5.75 Glycerine (70%) 17.48 65.8 1182 2.8-5.75 Glycerine (80%) 43.08 64.8 1209 3-5.6
For current experiments, the nozzle orientation is varied from α = 0◦ to α = 60◦ with an interval of 10◦ in order to investigate the droplet deformation and breakup. 6
The uncertainty in the calculation of the field variables, like droplet diameter (D), the air phase velocity and other dimensionless numbers are listed in Table IV. The maximum uncertainty propagation rule has been applied to calculate the uncertainty. According to this rule, if p, q, r, ... are the field variables having uncertainties Up, Uq, Ur, ... in the measurements, the uncertainty in Q = Q(p, q, r, ...) can be calculated as follows s δQ 2 δQ 2 δQ 2 δQ = · U + · U + · U + ...... (5) δp p δq q δr r
Table IV: Uncertainty in the calculations of different variables in the present study. Variable Maximum uncertainty α (degree) 5 × 10−1 Air velocity (m/s) 6 × 10−3 Weber number (W e) 3.3 × 10−1 E¨otv¨osnumber (Eo) 8 × 10−2 Nozzle diameter (m) 5 × 10−4 Droplet diameter (m) 10−4
IV. RESULTS AND DISCUSSION
We begin the presentation of our results by demonstrating the topological changes of a freely falling water droplet under the action of gravity when subjected to an air stream flowing at an angle α = 60◦ to the horizontal. The dynamics is compared again the cross-flow configuration (α = 0◦), which has been previously studied by several researchers (see for instance, Refs. [6, 11, 12, 14, 15]). Figs. 3 and 4 present the temporal variations of a droplet undergoing the vibrational, the transitional and the bag breakup for α = 0◦ and 60◦, respectively. By varying the droplet size and air-stream velocity, W e and Eo are varied to demonstrate the dynamics for α = 0◦ and 60◦. The dimensionless time, τ (≡ (ρl − ρa)gt/Uρa) is used to show the temporal evolutions of the droplet, such that τ = 0 represents the dimensionless time when the droplet just enters the uniform air-stream.
Vibrational breakup τ = 0 τ = 8.55 τ = 19.99 τ = 25.06 τ = 30.76
Transition τ = 0 τ = 11.05 τ = 20.29 τ = 27.17 τ = 31.95
Bag breakup τ = 0 τ = 11.24 τ = 16.64 τ = 18.58 τ = 19.01
Figure 3: The breakup dynamics of a water droplet in the cross-flow condition (α = 0◦) for Eo ≈ 2.4 demonstrating the vibrational (W e = 7.87), the transitional (W e = 8.27) and the bag breakup (W e = 9.13) modes. The dimensionless times, τ are shown at the top each panel, such that τ = 0 represents the time when the droplet just enters the uniform air-stream.
For α = 0◦ and Eo ≈ 2.4, it can be seen in the first row of Fig. 3 that (for W e = 7.87) the aerodynamic force is not enough to overcome the effect of the surface tension and viscous force, and the droplet undergoes shape oscillations 7
(the vibrational regime). Due to the aerodynamic force, the amplitude of oscillations continue to increase and the droplet decomposes into comparable size droplets (not shown). In this case, it is not necessary that the disturbance created by the aerodynamic force is sufficient enough for the fragmentation process. The deformation time scale of the droplet undergoing the vibrational mode is also high as compared to the other modes. For a higher value of the Weber number (W e = 8.27), the droplet exhibits the transitional mode (second row of Fig. 3) which is an intermediate stage between the vibrational mode and the bag breakup mode. As we increase the value of Weber number further (W e = 9.13) keeping the rest of the parameters fixed, the droplet exhibits the bag breakup mode as seen in the third row of Fig. 3. The minimum value of the Weber number at which the droplet undergoes bag breakup phenomenon is known as the critical Weber number, W ecr. It can be concluded that for a fixed set of other parameters, the bag breakup mode can be achieved by increasing the air stream velocity only. The temporal evolutions for the vibrational, the transitional and the bag breakup modes for α = 60◦ and Eo ≈ 1.8 are shown in Fig. 4 for W e = 5.13, 5.73 and 6.35, respectively. It can be seen that the phenomena are qualitatively similar to those observed in the cross-flow configuration. In all these cases, the droplet enters inside the air nozzle at τ ≈ 10, 10.38 and 11.56 and comes out at τ ≈ 42.72, 31.13 and 31.07 with oblate, disk and bag-type shapes for the vibrational, the transitional and the bag breakup modes, respectively. This confirms that even in the oblique configuration the droplet was interacting with the continuous air stream in the potential core region.
Vibrational breakup τ = 0 9.99 26.91 42.72 51.04
Transition τ = 0 10.38 20.75 31.13 36.13
Bag breakup τ = 0 11.56 21.24 31.07 32.26
Figure 4: The breakup dynamics of a water droplet in an oblique configuration (α = 60◦) for Eo ≈ 1.8 demonstrating the vibrational (W e = 5.13), the transitional (W e = 5.73) and the bag breakup (W e = 6.35) modes. The dimensionless times, τ are shown at the top each panel, such that τ = 0 represents the time when the droplet just enters the uniform air-stream. The inset in each panel of the last two columns represent the zoomed view of the droplet.
Next, we investigate the effect of the air stream directionality on the droplet deformation and breakup phenomena. The drop size, height from the nozzle to tip, air flow velocity are kept constant and only the angle of the air stream, α is varied. Figs. 5(a), (b) and (c) present the temporal evolutions of the droplet along with their trajectories for α = 0◦, 30◦ and 60◦, respectively at Eo = 1.25. Note that the Weber number can be calculated based on the relative velocity of the air stream with respect to the droplet or only the average velocity of the air stream from the nozzle. The parameters used to generate Fig. 5 are given in Table V. It can be seen that the droplet velocity, ud (calculated by the image processing method) is much smaller than the air stream velocity in all the cases. Two sizes of water droplet with D = 3.03 mm and 4.23 mm are considered for which the values of the E¨otv¨osnumber are 1.25 and 2.44, respectively. It can be observed that the Weber number calculated based on the air stream velocity, U is constant for each droplet size, but the Weber number calculated based on the relative velocity between the air stream and the droplet decreases slightly with the increase in the oblique angle of the air stream, α. Thus, the Weber number, W e calculated based on the average air stream velocity is used to present the rest of the presented in this study. Note that in Fig. 5, only α is varied while the rest of the parameters remain fixed. It can be seen in Fig. 5(a) for α = 0◦ that the droplet undergoes the vibrational mode. As we increase the angle of the air stream, the droplet exhibits a transitional mode for α = 30◦ (see Fig. 5(b)) and the bag breakup mode for α = 60◦ (see Fig. 5(c)). At τ=0, the droplet is almost spherical for all the cases. In the vibrational mode, the non-uniform pressure distribution around the interface of the spherical drop initiates an internal flow within the drop, which in turn deforms the droplet to an oblate shape at τ ≈ 6.44 and a disk type shape at τ ≈ 14.79. In this case, the 8 aerodynamic force is not enough to disintegrate the droplet into multiple pieces. Due to the competition between the aerodynamic force and the surface tension force, the droplet undergoes shape oscillations with a certain frequency. In some situations, these oscillations may lead to a breakup of the droplet into large ligament and satellite droplets of comparable sizes.
Table V: The parameters considered to generate Fig. 5. Here, W er corresponds the Weber number defined based on the relative velocity between the air stream and the droplet.
α D (mm) U (m/s) ud (m/s) Eo W e W er 0◦ 3.03 12.33 0.3 1.25 8.0 8.40 4.23 9.82 0.3 2.44 7.1 7.54 30◦ 3.03 12.33 0.3 1.25 8.0 8.34 4.23 9.82 0.3 2.44 7.1 7.48 60◦ 3.03 12.33 0.3 1.25 8.0 8.20 4.23 9.82 0.3 2.44 7.1 7.32
For the same operating condition, as the air stream angle is increased to α=30◦ the droplet undergoes a transitional region (Fig. 5(b)). In this case, as the droplet enters the air stream, the unequal distribution of the pressure changes the droplet shape from the spherical to a disk-like shape at τ ≈ 10.69. At τ ≈ 15.85, a bag-like structure appears in the air flow direction, but due to insufficient aerodynamic force, the droplet bag re-tracked back to a disk-like shape. Subsequent at later times τ > 20.53, the droplet breaks into similar size droplets (not shown). The resultant satellite droplets exhibit vibrational modes. At α = 60◦, the droplet clearly shows the bag-type breakup (Fig. 5(c)). In this case, when the spherical drop enters the region of the steady air stream, it is influenced by the distribution of air pressure around the drop. At an equilibrium condition, the internal and the surface tension forces are just balanced with the external aerodynamic force. However, when a steady air stream flows around a drop, the air velocity distribution and the air pressure distribution at any point on the drop surface are not uniform (see Fig. 6). The air velocity is maximum at point ii and equals to zero at point i (the stagnation point) in Fig. 6. Thus, by Bernoulli’s equation, the air pressure becomes highest at point i, and lowest at point ii. In this situation, the external aerodynamic pressure causes the drop to deform to an oblate ellipsoid shape in the direction normal to the air flow. As the velocity at the equator of drop increases, the Bernoulli pressure difference increases accordingly, which flattens the drop more. Finally, the drop becomes a pancake-like shape (not shown as the droplet is inside the air nozzle at this instant). Due to the contraction of the drop, a wake region [15] is developed on the back side of the droplet, which facilities the bag formation at τ ≈ 19. The bag carries more mass at the periphery due to the pressure distribution along the interface, that suppresses the instability while at the centre of the bag grows parallel to the air stream direction. The hollow portion of the bag becomes thin, like a membrane, and finally ruptures after τ ≈ 23.62 for this set of parameters. Due to this rupture, many satellite droplets are generated. Then the toroidal ring expands in the radial direction and breaks into relatively bigger child droplets in the downstream direction (not shown). Comparison of the trajectories of the droplet in Figs. 5(a), (b) and (c) also reveals that as the oblique angle of the air stream, α increases, the droplet follow a curvilinear motion with a sharp turn and the distance travelled by droplet from the nozzle decreases. The change in projected area of the droplet changes the drag coefficient. As the tangent to the path of the droplet gives its velocity. The acceleration, thus the inertia force acting on the droplet, also changes in an oblique configuration. Thus, using such oblique configurations, one can design compact combustors with a small primary reaction zone. In Fig. 5(d), the trajectories of the droplets exhibiting the vibrational, the transitional and the bag breakup modes for Eo = 1.25 (D = 3.03 mm, solid lines) and 2.44 (D = 4.23 mm, dashed lines) are shown. For both values of Eo, at α = 0◦, the droplet enters the aerodynamic field, deforms into a disk-like shape and dragged away in the direction of the air flow. The shape of the droplet for α = 0◦ is shown at τ=25.59. At this time, the droplet is in the shear layer region outside the potential core region of the air stream. For higher values of α, the droplet decelerates while migrating in the downward direction. When the droplet reaches the equilibrium stage, i.e., when the forces acting on the droplet in downward direction are balanced by the forces in the upward direction, the droplet comes to a stationary position for a very short time and starts to move in the opposite direction. At this stage, the droplet deforms to a shape that will facilitate it to migrate in the air flow direction. The droplet remains in contact with the strong aerodynamic force for a longer time as compared to the low α case (this residence time is proportional to the value of α). This increases the probability to break the droplet even with smaller air velocity. Therefore, based on the above discussion, it can be concluded that the critical Weber number is sensitive to the air flow direction, α. Increasing α decreases the critical Weber number for the bag breakup. This observation can play a significant role in the atomization process where the co and counter swirling secondary air streams are used to promote the 9 secondary atomization process. The critical Weber number also sensitive to the droplet size. It is observed that for water droplets of diameters 3.03 (Eo = 1.25) and 4.23 mm (Eo = 2.44), the values of W ecr are 8 and 7.1, respectively. (a) (b)
1 1 1 τ = 0 2 2 3 τ = 0 4 τ = 20.53 5 τ = 6.44 1 2 3 τ = 25.59 5 τ = 4.47 τ = 15.85 4 τ = 14.79 τ = 10.69 3 τ = 22.07 4 2 5 5 3 4
(c) (d)
τ = 23.62 1 5 τ = 0 1 훼 = 60° 5 휏 = 23.61
τ = 21.33 4 1 2 훼 = 30° 1 휏 = 20.53 4 2 2 τ = 17.76 3 3
3 3 τ = 5.00 2 훼 = 0° 휏 = 26.61
Figure 5: Trajectory and the shape of the droplet at different times for Eo = 1.25 that corresponds to D = 3.03 mm for different values of α. (a) α = 0◦, (b) α = 30◦ (c) α = 60◦. The panel (d) combines the trajectories for Eo = 1.25 (D = 3.03 mm) and Eo = 2.44 (D = 4.23 mm) at different oblique angles of the air stream. The numbers in each panel show the instants at which the droplet shapes are presented.
Figure 6: Pressure distribution on the droplet subjected to air stream. 10
(a) (b) 40 40 Eo We Eo We 1.10 7.09 1.03 9.66 1.60 7.56 1.83 10.18 20 2.23 7.87 20 2.29 10.23
Y (mm) 0 Y (mm) 0
-20 -20 0 5 10 15 20 25 30 35 0 5 10 15 20 25 30 35 X (mm) X (mm) (c) (d) 40 40 Eo We Eo We 0.77 5.96 0.85 6.94 1.59 5.70 1.71 6.70 20 2.21 5.59 20 2.33 6.10
Y (mm) 0 Y (mm) 0
-20 -20 0 5 10 15 20 25 30 35 0 5 10 15 20 25 30 35 X (mm) X (mm) (e) (f) 40 40 Eo We 1.06 6.14 1.65 5.13 2.32 4.97 20 Predicted 20
Eo We
Y (mm) 0 Y (mm) 0 1.12 6.77 1.82 6.35 2.54 6.23 Predicted -20 -20 0 5 10 15 20 25 30 35 0 5 10 15 20 25 30 35 X (mm) X (mm)
Figure 7: Trajectories of a water droplet (Oh ≈ 0.0015) exhibiting (a,c,e) the vibrational and (b,d,f) the bag breakup modes for (a,b) α = 0◦, (c,d) α = 30◦ and (e,f) α = 60◦.
Next we investigate the path followed by the droplet subjected to an oblique air flow for different sets of Eo and W e as it plays important role in combustors designing and optimisation. Figs. 7(a,b), (c,d) and (e,f) correspond to α = 0◦ (cross-flow configuration), α = 30◦ and α = 60◦, respectively. Figs. 7 (a,c,e) and (b,d,f) show the trajectories for the droplet exhibiting the vibrational mode and the bag breakup, respectively. Here, Y = 0 represents the centre line of the nozzle. The trajectories are plotted before fragmentation of the droplet. It can be seen that droplet travel a short distance from the nozzle and undergoes a bigger turning in the bag breakup cases as compared to the vibrational cases. It is observed that for all α values considered, when the droplet velocity vector (tangent to the trajectory) 11 aligns closer to the air flow direction, the drag force exerted on the droplet is more, which increases the stretching of the droplet, and the surface tension force can no longer keep the droplet intact and it disintegrates into smaller droplets. Increasing the droplet size (increasing Eo) increases the curvature of the droplet trajectory. It is to be noted here that the value of W e considered is close to the critical Weber number for each case. However, the fact that the larger droplets, although smaller than the critical diameter based on the capillary length scale, break up even when the trajectory of the droplet is not aligned with the air stream direction. This indicates that the gravity force, aiding into droplet distortion, leads to increase the drag to overcome the cohesive force that keeps the droplet intact.
(a) α = 0◦ (b) α = 30◦ (c) α = 60◦
D (mm) Wecr (i) (ii) D (mm) We (i) (ii) cr D (mm) Wecr (i) (ii)
2.62 14.56 2.29 10.19 2.60 6.93
τ = 6.77 τ = 8.05 τ = 15.38
2.67 10.69 2.70 7.73 2.79 6.80
τ = 8.48 τ = 10.59 τ = 18.71
3.67 9.57 3.54 6.89 3.58 6.27
τ = 15.23 τ = 18.35 τ = 25.18
4.07 9.13 4.09 6.82 4.32 6.23
τ = 17.80 τ = 24.30 τ = 28.77
5.81 8.74 6.03 6.46 6.08 5.61
τ = 0 τ = 31.57 τ = 0 τ = 38.62 τ = 0 τ = 57.55
Figure 8: Breakup of a water droplet (Oh ≈ 0.0015) for (a) α = 0◦, (b) α = 30◦ and (c) α = 60◦. For each angle, the left side images are corresponding to drop shape just before reaching the air stream, and the right images are the same droplet corresponding to τ (written below the images) just before the onset of the bag breakup.
Figs. 8(a), (b) and (c) show instantaneous snapshots of the droplet of different sizes just before entering the air stream (panel (i); τ = 0) and just before the onset of bag breakup (panel (ii); the value of τ is written below the droplets) for α = 0◦, 30◦ and 60◦, respectively. The value of τ mentioned in panel (ii) for each α essentially represents the time required to complete the bag formation for different sizes of the droplet at its critical Weber number, W ecr. It can be observed that for each value of α, increasing the droplet size increases the value of the dimensionless breakup time. Increasing Eo (increasing droplet size) decreases the value of W ecr, and increases the breakup time because of the decrease in the critical air velocity requirement for the larger droplets. It can also be seen that with the increase in the obliquity angle of the air nozzle, the dimensionless breakup time, τ shows a significant increase. This can be attributed to the droplet deceleration prior to aligning closely with the air jet direction. Interestingly, the droplet breakup occurs much closer to the air nozzle for high obliquity angle (α = 60◦) with significantly larger residence time. This is an important requirement from the design perspective of a swirl stabilised combustors as the flame would anchor much closer to the swirler warranting accelerating swirl flow (or converging) to increase flame standoff distance. Such a practice have been adopted in previous studies (see Refs. [30, 31]). For α = 0◦ (cross-flow configuration), the path traveled by the droplet is large as compared to the path traveled by the droplet at any other obliquity angle conditions. Note that the droplet traveled more distance from the nozzle in α = 0◦ case than other cases with α values, but still the non-dimensional time is lower, which can be attributed to the acceleration of the droplet in the cross-flow condition. This is unlike the oblique cases where the deceleration of droplet is not observed as droplet motion is not totally against gravity. In order to demonstrate that small droplets require less time to disintegrate as compared to bigger droplets at their respective values of the critical Weber number, in Fig. 9, we present the shapes of different size droplets at two 12 dimensional times, t = 0 and t = 23.67 ms in the cross-flow configuration (α = 0). It can be seen that the droplets are almost spherical when they enter the air stream. We observed that the small droplets (with D = 2.62 mm and 2.76 mm) at their critical Weber numbers completely disintegrate, the intermediate size droplets (D = 3.37 mm and 4.10 mm) form bags and the bigger droplet (D = 5.91 mm) appears as a disk at t = 23.67 ms. Although the bigger droplet breaks at smaller Weber number, the time required for the breakup increases with the droplet size. As in this case, the velocity requirement is lower, the rate of energy transfer to the droplet from the air stream is reduced resulting into larger residence time to overcome the effect due to the surface tension and the viscous forces.
(i) D (mm) (ii)
2.62
2.76
3.37
4.10
5.91
t = 0 ms t = 23.67 ms Figure 9: Breakup of a water droplet (Oh ≈ 0.0015) at α = 0 (cross-flow configuration). The left side images are corresponding to drop shape just before entering the air stream and the right set of images are for the same droplet t = 23.67 ms.
Then we investigate the effect of angle of the air stream, α on the critical Weber number W ecr. In Fig. 10, we present the region maps showing the vibrational and the bag breakup modes in the W e − Eo plane for different values of α. Here, the Eo is changed by changing the droplet size in our experiments. In each panel, the red dashed line with filled square symbols depicts the boundary between the vibrational (open circle symbols) and the bag breakup ◦ (open square symbols) regimes. For α = 0 (cross-flow configuration), at Eo = 1.1, the value of the W ecr ≈ 14, which conforms the finding of previous studies [16, 17]. Increasing the value of Eo decreases the value of W ecr initially and then decreases at a slower rate. For high values of Eo, D ≥ lc. Thus the gravity force starts to dominate the surface tension force, and the droplet deforms from its initially spherical shape. This, in turn, increases the drag force exerted on droplet due to the increase in its surface area, which helps the droplet to deform more [26]. Thus, the critical Weber number decreases with the increase in Eo (or size of the droplet). Note that in the cross-flow condition (α = 0◦), air flows in the horizontal direction and the component of net acceleration in the vertical direction remains as g. Hence the reported E¨otv¨osnumber is calculated based on the gravitational acceleration component only. For the cases with the oblique air flow, the droplet undergoes a curvilinear motion with a rapid deceleration followed by acceleration with a small radius of curvature in its trajectory (see Fig. 5). The curvilinear trajectory makes the internal pressure distribution to fluctuate as the droplet is continuously changing its direction. The pressure variation in a curvilinear motion is normal to the instantaneous velocity vector, which increases the stretching of the droplet, and hence the frontal area and the drag force also increase. Comparison of the results presented in Fig. 10(a)-(f) reveals that increasing the value of α decreases the value of W ecr. For all values of α considered, W ecr drops sharply as we increase Eo initially and then decreases slightly or remains constant. This behaviour can also be seen in Fig. 11, where the variations of W ecr versus α are plotted for different values of Eo. It can be seen that W ecr decreases to an asymptotic value of 6 (approximately). Interestingly, Villermaux and Bossa [2] also theoretically found that 13
◦ ◦ W ecr ≈ 6 for α = 90 (in case of the opposed flow configuration). It can be seen that Fig. 5 that at α = 60 the ◦ droplet almost aligned with the direction of air stream. This is the reason why W ecr for α = 60 is approximately equal to 6, like in the case of the opposed flow configuration considered by Villermaux and Bossa [2].
(a) α = 0◦ (b) α = 10◦ 16 16
Bag breakup 12 12 Bag breakup
We 8 We 8
4 4 Vibration Vibration
0 0 1 2 3 4 5 1 2 3 4 5 Eo Eo
(c) α = 20◦ (d) α = 30◦ 16 16
Bag breakup Bag breakup 12 12
We 8 We 8
4 4 Vibration Vibration
0 0 1 2 3 4 5 1 2 3 4 5 Eo Eo (e) α = 40◦ (f) α = 60◦ 16 16 Bag breakup Bag breakup Vibration 12 12 Bag breakup Transition line
We 8 We 8
4 4 Vibration Vibration
0 0 1 2 3 4 5 1 2 3 4 5 Eo Eo
Figure 10: Region maps showing the vibrational and the bag breakup modes in W e − Eo plane for a water droplet (Oh ≈ 0.0015) at different angles of inclination of the air stream, α. The red dashed line in each panel shows the transition between the vibrational and bag breakups. 14
Eo 1.1 12 1.4 2.4 4.8 We 9 cr
6
3 0 10 20 30 40 50 60 α
Figure 11: Variation of the transitional Weber number, W ecr with α (in degree) for water droplet (Oh ≈ 0.0015) and different values of E¨otv¨osnumber. The uncertainty in the measurement in angles is ±0.5◦.
(a) (b) 20 α α o 9 0 0o o o 30 4 30 o o 15 60 60 /10 6 cr Re cr Uc (m/s) 3
10 We
0 2 3 4 5 6 7 1 2 3 4 5 6 D (mm) Eo (c) (d) 20 α 12 0 ° 30 °
4 9 60 ° 15 /10
cr 6 Re cr
Uc (m/s) α
10 We 3 0° 30° 60° 0 2 3 4 5 6 7 1 2 3 4 5 6 D (mm) Eo
Figure 12: (a,c) The critical dimensional air velocity, Uc (m/s) versus D (mm), and (b,d) W ecrRecr versus Eo for different values of α. Panels (a,b) and (c,d) correspond to the droplets of water (Oh ≈ 0.0015) and 80% glycerol + 20% water solution (Oh ≈ 0.078), respectively.
In Figs. 12(a) and (c), the critical air stream velocity, Uc that is just sufficient for the bag breakup is plotted against the droplet diameter, D for three different values of α with droplets of water and 80% glycerol+ 20 %water solution, 15 respectively. It can be seen that increasing the size of the droplet decreases the value of the velocity of the air stream required for the bag breakup. This behaviour can be attributed to the large frontal cross-section area and the large drag coefficient due to the deformation associated with bigger droplets. Additionally, the wake region [15] developed behind the deformed droplet facilitates the skewed internal pressure distribution leading to excessive stretching. A similar finding was also reported in Refs. [9, 17]. Close inspection of Fig. 12(a) and (c) also reveals that increasing the viscosity of the droplet the critical velocity requirement increases slightly for the same size of the droplet, which is also observed by previous researchers [8, 11, 16]. This is due to the suppression of the surface instability as the viscosity is increased. It can be seen that increasing α for the same droplet size decreases the critical velocity requirements for the droplet to undergo the bag breakup. This can be attributed to the effective energy dissipation to overcome the surface tension and the viscous effects through the larger turning in the droplet trajectory. Note that the rate of change of the angular momentum provides a measure of the energy transfer from air stream to the droplet to overcome the viscous and surface energy dissipation. Thus decreasing the radius of curvature in the curvilinear trajectory of the droplet increases the change in its angular momentum leading to a more effective energy transfer in accordance with the conservation of energy principle. It is obvious that as the droplet size increases, the critical Weber number (W ecr) decreases but the corresponding critical Reynolds number (Recr) of the air stream increases even though the air velocity is decreased (except when the stretching is excess). However, in Fig. 12(b) and (d), it can be seen that the product of W ecr and Recr does not vary much with the increase in the droplet size. This is observed for all angles of the air stream. The product of W ecr and Recr signifies the inertial force requirement to overcome the viscous and the surface tension forces, such that a low value of W ecr × Recr represents a situation when a less inertia force can lead to the droplet breakup. Increasing α decreases W ecr × Recr as the curvature in the droplet trajectory increases with increasing α leading to more effective energy transfer from the air stream to the droplet to overcome viscous and surface tension forces. Hence, it can be concluded that the opposed flow configuration is an effective way to promote the secondary atomization at a low cost of input energy [30, 31]).
(a) α = 0◦ (b) α = 30◦ 16 16
14 14
We We cr12 cr12 Eo Eo 10 1.45 10 1.45 2.78 2.78 5.36 5.36 8 8 0 0.025 0.05 0.075 0.1 0 0.025 0.05 0.075 0.1 Oh Oh (c) α = 60◦ 16 Eo 1.45 14 2.78 5.36 We cr12
10
8 0 0.025 0.05 0.075 0.1 Oh
◦ Figure 13: The critical Weber number, W ecr versus Ohnesorge number, Oh for different values of Eo: (a) α = 0 , (b) α = 30◦ and (c) α = 60◦. 16
Two dimensionless numbers, namely the Weber number, W e and the Ohnesorge number, Oh are found to be effective in characterising the breakup regimes and transition between the vibrational and the bag breakup regimes. For relatively low Ohnesorge numbers (Oh < 0.1), the Weber number has a dominant influence on droplet breakup, and it determines the breakup regimes as discussed above for Oh ≈ 0.0015 (water droplets). In Fig. 13(a-c), we present the variations of W ecr versus Oh for different values of Eo and α. The value of Oh is varied by changing the concentration of the glycerol-water solution. As expected. it can be seen that W ecr increases with increasing Oh for all the cases. Aalburg et al. [32] also found that the effect of surface tension, which is responsible to oppose the fragmentation, decreases with increasing Oh. Finally, we investigate the variations of W ecr with α for different values of Eo and for different glycerol-water solutions (for different values of Oh). It can be seen that W ecr decreases with the increase in α irrespective of the liquids considered. However, as observed for a low value of Oh (see Fig. 11), W ecr did not approach the asymptotic value of 6 observed in the opposed flow configuration [2]. The Ohnesorge number, Oh has been accounted in the m correlation to predict the modified critical Weber number, W ecr of a viscous droplet as follows [8, 11, 14]
m c2 W ecr = W ecr(1 + c1Oh ), (6) where the values of c1 and c2 are 0.776 and 0.338 respectively, which are obtained using a regression analysis. In the present study, the mean absolute percentage error between the predicted and experimental value is found to be < 6%. (a) Oh ≈ 0.009 (b) Oh ≈ 0.032 16 16 Eo Eo 1.45 1.45 14 14 2.78 2.78 5.36 5.36 12 12 We We cr cr 10 10
8 8
6 6 0 10 20 30 40 50 60 0 10 20 30 40 50 60 α α (c) Oh ≈ 0.078 16
14
12 We cr 10 Eo 1.45 8 2.78 5.36 6 0 10 20 30 40 50 60 α
Figure 14: Variations of the critical Weber number, W ecr versus α (in degree) for (a) 50% glycerol + 50% water solution (Oh ≈ 0.009), (b) 70% glycerol + 30 %water solution (Oh ≈ 0.032) and (c) 80% glycerol + 20% water solution (Oh ≈ 0.078).
V. CONCLUDING REMARKS
We have conducted experiments to investigate the effect of the air stream obliquity (α), the droplet diameter (D) and the fluid properties of the liquid on the deformation and the breakup phenomena of droplets experiencing an 17 oblique continuous air stream. A high-speed imaging system is used to record the trajectories and topological changes of the droplets. Four different liquids, namely, deionised water, 50%, 70% and 80% aqueous mixtures of glycerol (by volume) are considered in our experiments. The angle of the continuous air stream is varied from α = 0◦ to α = 60◦ ◦ with an interval of 10 . The values of the critical Weber number (W ecr) for the bag-type breakup are obtained as a function of the E¨otv¨osnumber (Eo), angle of inclination of the air stream (α) and the Ohnesorge number (Oh). It is found that, in case of the oblique air stream, the droplet entering the air stream follows a rectilinear motion that transforms into a curvilinear motion when the droplet undergoes topological changes. As we gradually change the cross-flow configuration to the in-line configuration, a sharp decrease in the critical Weber number is observed, which asymptotically reaches to the same value as that observed in an in-line (opposed) flow configuration. The droplet breakup time also increases with the increase in the oblique angle of the continuous air stream. This behaviour can be attributed to the deceleration of the droplet as it follows a curvilinear trajectory and shifts its direction to the nearly in-line (co-flow) arrangement. Thus increasing the residence time in a short path length leads to a more effective energy transfer from the air stream to the droplets, which dissipates to overcome the viscous and the surface tension forces, to form smaller satellite droplets. As expected, increasing the droplet diameter (or increasing Eo, which signifies a decrease in the relative influence of the surface tension force over the gravitational force) decreases the value of W ecr. Increasing the viscosity of the droplet (which is equivalent to increasing the value of Oh) requires higher W e for the droplet to undergo bag-type breakup.
Acknowledgments
Authors would like to acknowledge the financial support from the Department of Science and Technology, India (SERB grant No. ECR/2015/000365). SKS and PKK also thank the Ministry of Human Resource Development, India for providing research fellowships.
[1] T. Sikroria, A. Kushari, S. Syed, and J. A. Lovett, “Experimental investigation of liquid jet breakup in a cross flow of a swirling air stream,” J. Eng. Gas Turb. Power 136, 061501 (2014). [2] E. Villermaux and B. Bossa, “Single-drop fragmentation determines size distribution of raindrops,” Nature Physics 5, 697 (2009). [3] C. M. Varga, J. C. Lasheras, and E. J. Hopfinger, “Initial breakup of a small-diameter liquid jet by a high-speed gas stream,” J. Fluid Mech. 497, 405 (2003). [4] A. H. Lefebvre and V. G. McDonell, Atomization and sprays (CRC press, ADDRESS, 2017). [5] G. I. Taylor, “The shape and acceleration of a drop in a high speed air stream,” The scientific papers of GI Taylor 3, 457 (1963). [6] M. Jain, R. S. Prakash, G. Tomar, and R. V. Ravikrishna, “Secondary breakup of a drop at moderate Weber numbers,” Proc. R. Soc. Lond. A 471, 20140930 (2015). [7] Z. Dai and G. M. Faeth, “Temporal properties of secondary drop breakup in the multimode breakup regime,” Int. J. Multiphase Flow 27, 217 (2001). [8] M. Pilch and C. A. Erdman, “Use of breakup time data and velocity history data to predict the maximum size of stable fragments for acceleration-induced breakup of a liquid drop,” Int. J. Multiphase Flow 13, 741 (1987). [9] W.-H. Chou and G. M. Faeth, “Temporal properties of secondary drop breakup in the bag breakup regime,” Int. J. Multiphase Flow 24, 889 (1998). [10] S. A. Krzeczkowski, “Measurement of liquid droplet disintegration mechanisms,” Int. J. Multiphase Flow 6, 227 (1980). [11] D. R. Guildenbecher, C. L´opez-Rivera, and P. E. Sojka, “Secondary atomization,” Exp. Fluids 46, 371 (2009). [12] X.-K. Cao et al., “A new breakup regime of liquid drops identified in a continuous and uniform air jet flow,” Phys. Fluids 19, 057103 (2007). [13] S. R. and G. Tomar, “Secondary atomization,” J. Indian I. Sci. 99, 77 (2019). [14] V. Kulkarni and P. E. Sojka, “Bag breakup of low viscosity drops in the presence of a continuous air jet,” Phys. Fluids 26, 072103 (2014). [15] A. K. Flock et al., “Experimental statistics of droplet trajectory and air flow during aerodynamic fragmentation of liquid drops,” Int. J. Multiphase Flow 47, 37 (2012). [16] A. R. Hanson, E. G. Domich, and H. S. Adams, “Shock tube investigation of the breakup of drops by air blasts,” Phys. Fluids 6, 1070 (1963). [17] A. Wierzba, “Deformation and breakup of liquid drops in a gas stream at nearly critical Weber numbers,” Exp. Fluids 9, 59 (1990). [18] S. A. Krzeczkowski, “Measurement of liquid droplet disintegration mechanisms,” Int. J. Multiphase Flow 6, 227 (1980). [19] L. P. Hsiang and G. M. Faeth, “Drop properties after secondary breakup,” Int. J. Multiphase Flow 19, 721 (1993). 18
[20] C. Wang, S. Chang, H. Wu, and J. Xu, “Modeling of drop breakup in the bag breakup regime,” Appl. Phys. Lett. 104, 154107 (2014). [21] Z. Dai and G. M. Faeth, “Temporal properties of secondary drop breakup in the multimode breakup regime,” Int. J. Multiphase Flow 27, 217 (2001). [22] W. H. Chou and G. M. Faeth, “Temporal properties of secondary drop breakup in the bag breakup regime,” Int. J. Multiphase Flow 24, 889 (1998). [23] H. Zhao et al., “Influence of surfactant on the drop bag breakup in a continuous air jet stream,” Phys. Fluids 28, 054102 (2016). [24] T. Inamura, H. Yanaoka, and T. Kawada, “Visualization of airflow around a single droplet deformed in an airstream,” Atomization and Sprays 19, 667 (2009). [25] M. Brenner and A. Bertozzi, “Spreading of droplets on a solid surface,” Physical review letters 71, 593 (1993). [26] M. Jalaal and K. Mehravaran, “Fragmentation of falling liquid droplets in bag breakup mode,” Int. J. Multiphase Flow 47, 115 (2012). [27] J. A. Nicholls and A. A. Ranger, “Aerodynamic shattering of liquid drops.,” AIAA J. 7, 285 (1969). [28] R. Rioboo, M. Marengo, and C. Tropea, “Time evolution of liquid drop impact onto solid, dry surfaces,” Exp. Fluids 33, 112 (2002). [29] “http://www.aciscience.org/docs/physical-properties-of-glycerine-and-its-solutions.pdf,” . [30] K. Merkle, H. Haessler, H. B¨uchner, and N. Zarzalis, “Effect of co-and counter-swirl on the isothermal flow-and mixture- field of an airblast atomizer nozzle,” International Journal of Heat and Fluid Flow 24, 529 (2003). [31] A. K. Gupta, M. J. Lewis, and M. Daurer, “Swirl effects on combustion characteristics of premixed flames,” J. Eng. Gas Turbines Power 123, 619 (2001). [32] C. Aalburg, B. v. Leer, and G. Faeth, “Deformation and drag properties of round drops subjected to shock-wave distur- bances,” AIAA journal 41, 2371 (2003).