ILASS-Americas 29th Annual Conference on Liquid Atomization and Spray Systems, Atlanta, GA, May 2017 Size-velocity pdfs for Drop Fragments Formed via Bag Breakup C.M.L. White,* G. Sondgeroth,+ W. Shang,+ L. Yao,+ J. Chen+, and Paul E. Sojka+ *School of Aeronautics and Astronautics +School of Mechanical Engineering all Maurice J. Zucrow Laboratories Purdue University West Lafayette, IN 47907 USA and D.R. Guildenbecher Sandia National Laboratories Albuquerque, NM Abstract Digital-inline holography (DIH) was used to measure size probability distribution functions, pdf(d), for fragments formed via bag breakup. A MatLab script was used to reconstruct drop fragmentation dynamics, from which fragment sizes, as well as velocities (and accelerations), can be extracted. Results, which demonstrate bimodal fragment size distributions, are reported in terms of Weber and Ohnesorge numbers for the ranges 13<We<30 and 0.002<Oh<0.45. Physical explanations for the presence of the two-and three-peaked size distribution are presented. The data will be useful to those modeling sprays in gas turbine engines, pharmaceutical tablet coaters, and spray dryers.
ILASS-Americas 29th Annual Conference on Liquid Atomization and Spray Systems, Atlanta, GA, May 2017 Introduction and Literature Review Drop breakup occurs in many multi-phase flows. Analytical or numerical models of that breakup can be used to optimize procedures such as agricultural spray drift, pharmaceutical tablet coating, and engine efficiency and emissions. The advantage of modeling drop break up lies in prediction of the relative quantity of different sized drops so their effects on a larger system can be predicted and system performance optimized. The majority of drop breakup research has focused on Newtonian drops. Nicholls and Ranger (1969), Krzeczkowski (1980), Pilch and Erdman (1987), Wierzba and Takayama (1988), Hsiang and Faeth (1992, 1995), and Liu and Reitz (1997) are among those who identified the five modes of breakup for low Ohnesrge numbers (Oh < 0.1). These modes include vibration (We < 11), bag (11 < We < 35), multi-mode (35 < We < 80), sheet thinning (80 < We < 350), and catastrophic (We > 350). These findings are summarized in the review of Guildenbecher et al. (2009). In order to predict drop breakup, the aerodynamic fragmentation is often modelled by solving gas-phase dynamics over a fixed grid. This grid aids in tracking drop dynamics, which, when coupled with models for aerodynamic drag, drop deformation, and drop fragmentation, results in predictions for breakup. Earlier models were criticized because they were based on boundary layer effects. Current versions employ Kelvin-Helmholtz and Rayleigh-Taylor instabilities (Theofanous, 2011). Deciding which instability is dominant during different time intervals within the breakup event can t be completed because of insufficient experimental data for this dynamic process. Such information must have sufficient temporal and spatial resolution, as well as dynamic range, to capture both large- (initial drop, rim or core) and small-scale (fragments) structures. Previous efforts on bag regime breakup observed the formations of complex structures, but did not address them quantitatively. This was due to a lack of spatiotemporal resolving capacity of the available measuring techniques, which, as a result, greatly hindered drop breakup modelling. Consequentially, a comprehensive and fundamental characterization of Newtonian drop fragmentation and experimental datasets to support model development remain absent. The majority of experimental data acquired for the modelling of liquid drop breakup were collected by means of 2D high speed imaging and/or PDA. However, these measuring techniques are incapable of recording the 3D spatial distribution of particle sizes and velocities. Furthermore, 2D imaging techniques may result in erroneous detection if the object in question is out of focus. Digital in-line holography (DIH) overcomes these limitations, as it provides 3D, time-resolved information. The measuring technique has been applied to the characterization of multiphase droplets (Tian et al., 2010; Gopalan and Katz, 2010; Lebrun et. al. 2011; Gao et al., 2013) Guildenbecher et al., 2016). While in use, DIH has demonstrated the ability to measure the sizes and 3D locations of the droplets without knowing the refractive index of the liquid. For this study, DIH was used to measure the size pdf for fragments formed during the bag, rim, and stamen stages of breakup. These measurements are significant because they have the potential to promote the development of new techniques for the optimal design of atomization systems that use droplet breakup. For example, this research may lead to new designs of fuel injectors in gas turbine and other engines, by creating optimal conditions for control of fuel droplet sizes, which may eventually lead to more efficient energy production and lower pollution levels. It may also contribute to the design of agricultural spray atomizers that optimize the consumption of herbicides and pesticides while also reducing drift to surrounding areas. This research is therefore expected to not only improve the fundamental understanding of drop breakup, but also have a broad and positive impact on society. Experimental Apparatus The drop generator and air nozzle used to study fragmentation are described by Guildenbecher and Sojka (2011). The air nozzle produces a nearly uniform jet velocity profile while the drop generator operates in the periodic dripping regime and is positioned to release drops of 3 mm diameter above the centerline of the air jet at approximately three second intervals. In order to mimic plug flow air jet initial conditions, a drop initial height of ~0.5 m is used. This provides a vertical drop velocity ~3 m/s. During this study DIH was used to investigate fragmentation of DI-water drops as a function of We. Weber number was controlled by varying the air jet velocity. The DIH experimental apparatus follows that of Guildenbecher et al. (2016) and consists of a Big Sky Laser Evergreen (m/n 00248300-2) frequency doubled Nd:YAG laser that produces a double-pulse (each 90 nsec) train at 2 Hz. The Nd:YAG beam first enters a spatial filter (25 μm pinhole) and lens (125 mm focal length) that expands and collimates it to a diameter of 5.0 cm. The 5.0 cm beam then illuminates the breakup field before being recorded by a Phantom VR2512 camera with pixel size 25x25 μm. Typical raw holograms, shadowgraphs, and reconstructed holograms are presented below. Data extraction from these images follows and is discussed in the next section.
Results and Discussions Data were collected at three different times during the breakup events. The first time was after the bag had fragmented and before the rim did. The second time was after the rim had broken up, while the third time was after the stamen broke up. Figure 2 shows typical bag breakup data; We=15. pdf(d) has a single maximum and goes to zero at a diameter of about 400 m. Similar behavior was observed at We 30. This suggests a single mechanism for fragment formation, a conclusion supported by shadowgraphs which show bag breakup resulting in a plethora of small drops. mechanism for stamen breakup, a conclusion again supported by shadowgraphs. Figure 3. Fragment size pdf(d) for rim breakup at We=13. Figure 2. Fragment size pdf(d) for bag breakup at We=15. Figure 3 shows typical data after rim breakup has occurred. The majority of cases for 13 We 30 showed a second pdf(d) peak, which occurred at a larger diameter than for the bag fragments (between 500 and 700 m). The maximum diameter for pdf(d) is about 800 m in all cases. This suggests a two separate mechanisms for bag and rim fragment formation, a conclusion again supported by shadowgraphs where rim fragments were observed to be much larger than those from bag breakup. Figure 4 shows typical pdf(d) after stamen breakup. This only occurred for We 25. Here, the bag and rim pdf(d) are essentially the same as for the cases discussed above. However, a third pdf(d) peak emerges at about 900 m. This new peak suggests yet another physical Figure 4. Fragment size pdf(d) for stamen breakup at We=25. Summary and Conclusions Drop breakup was investigated at 13 We 30 using DIH. Bag breakup yielded a mono-modal pdf whose peak was near 100 m and whose maximum was approximately 400 m. Rim breakup usually resulted in a second pdf(d) peak between 500 and 700 m; the rim pdf(d) dropped to zero at about 800 m. Shadowgraph data support the conclusion that separate physical mechanisms are responsible for bag and rim breakup.
Bag-and-stamen breakup was observed for We 25. In these cases there was a third pdf(d) peak, due to stamen fragmentation, which usually occurred at about 800 m. The maximum non-zero pdf(d) value rose to about 1 mm. Acknowledgements This work was supported by the Purdue University Provost Office, the Purdue University School of Mechanical engineering. References Gao, J., D.R. Guildenbecher, P.L. Reu, V. Kulkarni, P.E. Sojka, and J. Chen, Quantitative 3D diagnostics of multiphase drop fragmentation vis digital inline holography, Optics Letters 38 (2013): 1893 1895. Gopalan, B. and Katz, J., Turbulent shearing of crude oil mixtures with dispersants generates long microthreads and microdroplets, Phys Rev Lett, 104, 054501 (2010). Guildenbecher, D.R., Lopez-Rivera, C, and Sojka, P.E., Secondary atomization, Experiments Fluids, 46, 371 (2009). Guildenbecher, D.R., and Sojka, P.E., Experimental investigation of aerodynamic fragmentation of liquid drops modified by electrostatic surface charge, Atomization Sprays, 21(2), 139-147 (2011). Guildenbecher, D.R., M.A. Cooper, and P.E. Sojka, High-speed (20 khz) digital in-line holography for transient particle tracking and sizing in multiphase flows, Applied Optics 55(11) 2892-2903 (2016). Hsiang, L.P., and Faeth, G.M., Near-limit drop deformation and secondary breakup, Int l J Multiphase Flow, 18(5), 635-652 (1992). Hsiang, L.P., and Faeth, G.M., Drop deformation and breakup due to shock wave and steady disturbances, Int l J Multiphase Flow, 21(4), 545-560 (1995). Krzeczkowski, S.A., Measurement of liquid droplet disintegration mechanisms, Int l J Multiphase Flow, 6(3), 227-239 (1980). Lebrun, D., Allano, D., Mees, L., Walle, F., Corbin, F., Boucheron, R., and Frechou, D., Size measurements of bubbles in a cavitation tunnel by digital in-line holography, Applied Optics, 50(34), H1-H9 (2011). Liu, Z., and Reitz, R.D., An analysis of the distortion and breakup mechanisms of high speed liquid drops, Int l J Multiphase Flow, 23(4), 631-650 (1997). Pilch, M., and Erdman, C.A., 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 l J Multiphase Flow, 13(6), 741-757 (1987)). Nicholls, J.A. and Ranger, A.A., Aerodynamic shattering of liquid drops, AIAA J 7(2), 285-290 (1969). Theofanous, T.G., Aerobreakup of Newtonian and viscoelastic liquids, Ann Rev Fluid Mechanics, 43, 661-690 (2011). Tian, L., Loomis, N., Dominguez-Caballero, J.A., and Barbastahis, G., Quantitative measurement of size and three-dimensional position of fast-moving bubbles in air-water mixture flows using digital holography, Applied Optics, 49(9), 1549-1554 (2010). Wierzba, A., and Takayama, K., Experimental investigation of the aerodynamic breakup of liquid drops, AIAA J, 26(11), 1329-1335 (1988).