Behavior of a Droplet Impacting on a Horizontal Plate

Size: px

Start display at page:

Download "Behavior of a Droplet Impacting on a Horizontal Plate"

Dorcas Bridges
5 years ago
Views:

1 Journal of JSEM, Vol., Special Issue (0) SS5-SS30 Copyright C 0 JSEM Behavior of a Droplet Impacting on a Horizontal Plate Yusuke SAKAI, Yoshiaki UEDA, Toshio ISHII and Manabu IGUCHI Graduate School of Engineering, Hokkaido University, Sapporo , Japan JFE Steel Corporation, Kanagawa, Japan (Received 7 December 00; received in revised form 7 April 0; accepted 5 June 0) Abstract The aim of this study is to investigate an unsteady behavior of a droplet impacting on a horizontal plate with poor wettability. Visualization techniques employed in this study include Computational Fluid Dynamics (CFD) and experimental one using a high-speed camera. The contact angle of the plate is varied by coating a water repellent material on its surface. This paper presents snapshots of bouncing motion of a droplet, the deformation rate at the impact on the plate, the coefficient of restitution, and flow in the droplet. Key words Droplet, Wettability, Deformation Rate, Coefficient of Restitution, Surface Tension, Experimental Visualization, CFD. Introduction Information on a single droplet impacting on a solid surface is of practical importance in many engineering fields. Thus, many researchers have investigated the dynamic behavior of a droplet impacting on a solid surface [ 3]. The results of these studies are applicable to fabrication of semiconductor, cleaning of ink-jet printers, development of steelmaking processes and so on [4 6]. In the steelmaking processes, for example, the refractory is usually chosen to be poorly wetted by molten steel to prevent metallurgical reactions with the steel [7]. The wettability of the refractory therefore plays an important role for the refining of molten steel. A recent review on the wettability and surface tension physics is found in the textbook of de Gennes, Brochard- Wyart & Quéré [8]. Kato et al. investigated the critical inclination angle of a plate for sliding down of a liquid drop [9 ]. Sonoyama and Iguchi calculated the advancing and receding contact angles of a bubble or a droplet using Finite Element Method (FEM) []. Kagawa and Iguchi investigated a detachment shape of a single silicone oil or water droplet from the exit of a wetted or a poorly wetted single-hole nozzle [3]. Akao et al., Wachters et al. and Qiao et al. investigated the behaviors of a droplet impacting on a heated substrate surface [4 6]. Many researchers have investigated the maximum spread of the droplet impacting on a rigid surface [7 9]. Gunjal et al. researched the influence of the diameter of a droplet, impact velocity, surface tension, and contact angle on the dynamic behavior of the droplet by experiment and numerical simulation using the VOF method [0]. Mao et al. compared the experimentally measured maximum spread of a droplet impacting on a flat surface of c = 30 ~ 0 deg. with the numerically calculated one [], where c is the equilibrium contact angle. Previous researchers focused on the dynamic behavior of a droplet impacting on a rigid surface. Investigations on the bounce-off height after the impact and coefficient of restitution are very limited. For example, one of research groups in France dealt with dynamic behavior of a droplet on a super-hydrophobic inclined plate of a contact angle of approximately 70 deg. [], the detailed mechanism is still open as far as the authors know. This study investigates the dynamic behavior of a water droplet impacting and bouncing on a hydrophobic horizontal plate and then compares the experimental results obtained in this study with the existing equation proposed by Mao et al. Furthermore, a particular attention is paid to the influence of the number of impacts on the bounce-off height and the coefficient of restitution. Empirical equations are proposed for the coefficient of restitution as a function of the Weber number.. Experimental Procedure The experimental apparatus and physical setting is shown in Fig.. The flat plate made of transparent acrylic resin was 00 mm long and 00 mm wide. The original contact angle, c, of de-ionized water on the acrylic plate in the air was 70 deg. This experiment used a water repellent to vary the contact angle c to 44 deg. A single water droplet was generated at the exit of a single-hole nozzle. The inner nozzle diameter, d ni, was 0.0,.0, and.7 mm and the corresponding outer nozzle diameter was 0.6,.0, and 3. mm. The volume of the droplet, V D, was calculated by measuring the weight of 0 droplets. The volume of the droplet tank cock Fig. Experimental apparatus and physical setting h 0 high-speed camera -SS5-

2 Y. SAKAI, Y. UEDA, T. ISHII and M. IGUCHI symmetric face droplet hydrophobic plate 6.5d D tension. The wettability at the meniscus was taken into account using the geometrical condition n = n w cos c + t w sin c () on the contact line where n w and t w are the normal and tangential vectors on a wall. We employed the following Laplace-Young equation. p = ij (/R + /R ) (3) z y x 5.0d D 3.5 d D grids points Fig. Computational grids and initial setting (V D = 33. mm 3, h 0 = 5.4 mm) droplet was found to be 4.3, 33., and 45.4 mm 3. The droplet was released from a predetermined height, h 0. The Weber number, We, was changed from 0.00 to 36.. The Weber number is defined as follows: We = L d D v D / where L is the liquid density, d D is the equivalent volume diameter, v D is the velocity of the droplet at the impact, and is the surface tension. The Weber number means the ratio of the inertial force to the force due to the surface tension. 3. Numerical Procedure In all numerical simulations, FLUENT TM ver.6..6 was employed on.66 GHz Intel (R) Core (TM) Quad CPU Q9400 processor with 3.5 GB RAM. The construction of the three-dimensional computational grids was made by GAMBIT. The computational grid was made of cubic 3,30 elements (.5d D < x <.5d D and 0 < y < 6.5 d D and 0 < z < 3.5d D, d D : diameter of an initial droplet), as shown in Fig.. To save CPU time, the droplet and flow field were assumed to be symmetric with respect to the symmetry face shown in Fig.. The contact angle of de-ionized water on the plate was 44 deg. The spherical water droplet was set at the certain height from the horizontal plate, which was identical with the above-mentioned experimental condition. In the FLUENT code, we employed a segregated implicit solver and first-order upwind interpolation scheme. A small time-step of t = (s) was adapted to achieve a convergence in every iteration. A free surface behavior (i.e. gas-liquid interface behavior) was tracked by the VOF model. The contribution of a volume force from a wall adhesion was added to the momentum equation as the source term, i.e., F vol = ij ( i grad i ) / [(/) ( i + j )] with the curvature i = div n and the surface normal n = grad i. Here, i is the density of the liquid i, j is the density of the gas j, is the volume-averaged density, i is the volume fraction of the liquid i, and ij is the surface where p is the pressure difference between the two fluids and, /R and /R are the surface curvatures as measured by two radii in orthogonal directions (see FLUENT 6. User s Guide [3] for more details). Laplace Young equation is available only in the case of the static condition. However, if the time-step t is chosen to be much shorter than the whole deformation period of the droplet, a quasistatic condition would be realized during the time-step t because the droplet would deform a little. This is the reason why Eq. (3) was adopted in this study. In our previous study, we affirmed this numerical procedure was in agreement with experimental results of dynamic behavior of a bouncing droplet on an inclined plate [4]. 4. Results and Discussion 4. Comparison of the measured diameter of a falling droplet from the nozzle exit and calculated value from an empirical equation Table shows the relationship between the nozzle diameter and the diameter of a single droplet falling in the atmosphere. The following empirical equation proposed by Ohtake et al. [5] could predict the diameter, d D, of a single droplet falling in the atmosphere under the critical flow rate. d D =.6 d n [ / ( gd n ) ] 0.35 (4) where ( = ij ) is the surface tension, is the density difference between liquid and gas, and g is the acceleration due to gravity. When the wettability of the nozzle is good, the outer diameter, d no, should be substituted into d n in Eq. Table Comparison of the measured diameter of falling droplet from the exit nozzle and empirical equation (4) outer nozzle diameter d no = 0.6 mm d no =.0 mm d no = 3. mm empirical equation d D =.7 d D = 3.99 d D = 4.59 experiment d D =.0 d D = 3.99 d D = SS6-

3 Journal 5th of JSEM, Vol., Special Issue (0)Science and Technology in Experimental Mechanics, 4-7 November, 00, Kyoto, Japan International Symposium on Advanced (4), while when the wettability of the nozzle is poor, the inner diameter, dni, should be into dn. In this study, the wettability of the nozzle was good, so we adopted dno as dn. As shown in Table, all experimental results for dd were in good agreement with Eq. (4). 4. Behavior of a droplet impacting on a horizontal plate Details of the behavior of a droplet impacting on the hydrophobic plate can be observed by a high-speed camera. The photographs displayed in Fig. 3 were taken at 500 frames per second. The volume of the droplet was 33. mm3, and it was set at a height of 9.7 mm from the horizontal plate. The Weber number was 7.3. At first, the droplet impacted on the hydrophobic plate (Fig. 3(a)), and then deformed like a pan-cake (Fig. 3(c)). After the impact, the droplet bounced off the hydrophobic plate several times (Fig. 3(i)). On the other hand, a droplet was observed to pinch off when the Weber number was 6.7, as shown in Fig. 4(m). A similar phenomenon was reported by Akao et al [4]. They investigated the behavior of a droplet impacting on a hot plate. They concluded that disintegration of the droplet was caused by the instability of the liquid column. The deformation rate, ed, of the droplet was defined as follows: ed = dmax / dd We We )e d ]e d ( [ ( cos T c ) Re (6) Figure 7 shows the comparison between the measured values of ed and the calculated ones. The deformation rate can be predicted by Eq. (6) for ed approximately greater than.6. Figure 8 shows the relationship between the Weber number and the deformation rate ratio, edexp / edcal, where the subscripts EXP and CAL denote the measured and calculated values, respectively. This figure demonstrates that Eq. (6) is valid for the Weber number greater than about 0, although Eq. (6) was originally derived for Tc d 0 deg. (a) (b) (c) (d) (e) (f) (g) (h) (i) (j) (k) (l) (m) (n) (o) (5) where dmax denotes the maximum width of the droplet during the impact motion and dd is the droplet diameter. Figure 5 shows a comparison of ed with the number of impacts, N. When the Weber number was large (We = 3.5), ed decreased with an increase in the number of impacts. The values of ed, however, were kept constant when the Weber number was small (We =.6). The reason is that the droplet hardly deforms like the feature in Fig. 6. Mao et al. proposed the following empirical equation for the deformation rate, ed. The range of the impact (a) velocity is from 000 to 6000 mm/s. (b) (c) (d) (e) Fig. 4 Behavior of a droplet pinching off after the impact on the hydrophobic plate. The volume of the droplet is 33. mm3, the Weber number is 6.7, and the time interval of each photograph is 4 ms ed.5 (f) (g) (h) (i) (j) (k) (l) (m) (n) (o).5 0 Fig. 3 Photographs of a droplet impacting on the hydrophobic plate. The volume of the droplet is 33. mm3, the Weber number is 7.3, and the time interval of each photograph is ms 3 4 N Fig. 5 Deformation rate of a droplet impacting on the hydrophobic plate. The Weber number is.6 (circle),.6 (triangle), and 3.5 (square), respectively -SS7-

Y. SAKAI, Y. UEDA, T. ISHII and M. IGUCHI (a) (b) (c) (d) (e) (f) (g) (h) (i) (j) (k) (l) (m) (n) (o) e 0.8 0.6 h 0 = 4.8mm 0.4 h 0 = 8.0mm h 0 = 9.8mm 0. h 0 = 3.0mm h 0 = 5.6mm 0 0 3 4 N Fig.

4 Y. SAKAI, Y. UEDA, T. ISHII and M. IGUCHI (a) (b) (c) (d) (e) (f) (g) (h) (i) (j) (k) (l) (m) (n) (o) e h 0 = 4.8mm 0.4 h 0 = 8.0mm h 0 = 9.8mm 0. h 0 = 3.0mm h 0 = 5.6mm N Fig. 6 Photographs of a droplet impacting on the hydrophobic plate taken by high-speed camera. The volume of the droplet is 4.3 mm 3, the Weber number is 0.63, and the time interval is ms Calculated e d Measured e d Fig. 7 Comparison between the measured deformation rate, e d, and the calculated one proposed by Mao et al e dexp / e dcal We Fig. 8 Relationship between the Weber number and the deformation rate ratio The coefficient of restitution (COR), e, of a droplet was defined as e = (h n+ / h n ) / (7) where h n is the n-th peak height from the center of the droplet to the plate. Figure 9 shows the relationship between the number of impacts, N, and COR, e. The volume of the Fig. 9 The coefficient of restitutions of a droplet of 33.mm 3 at the impact on the hydrophobic plate e % 0. e e e 3 Richard et al. (d D = 0.8mm) Richard et al. (d D =.0mm) We +0% Fig. 0 Relationship between the Weber number and the coefficient of restitution. Circles indicate the value of e, triangles the value of e, and squares the value of e 3. The values of e for We.6 is fitted as e = We 0. droplet was 33. mm 3. As shown in Fig. 9, e became greater with a decrease in h 0. Here, the subscript denotes the first impact. The reason is that the droplet deforms more significantly at the impact for higher falling velocity (see Fig. 4 (a) (g)) and dissipates energy due to the intense surface deformation. The value of e 3 was not dependent on the number of impacts. This is because energy loss at the third impact is very small, as the droplet hardly deforms (Fig. 6). Figure 0 shows the relationship between the Weber number, We, and COR, e, obtained in the experiments. The solid circle and solid triangle are the results obtained by Richard et al. [8]. A water droplet, whose diameter was 0.8mm or.0mm, impacted on a super-hydrophobic surface of c = 70 deg. The plane was slightly tilted at an angle of about deg. They mentioned that when the Weber number was smaller than around 0.4, the error increased. The reason is that it becomes hard to determine the heights from the plate from which the drops are detaching. The present data are in good agreement with those of Richard et al. except the Weber number smaller than 0.4. When the Weber number was smaller than unity, the force due to the surface tension was dominant to the behavior of the droplet. In this region (We <.6), the droplet hardly loses -SS8-

5 Journal of JSEM, Vol., Special Issue (0) its energy due to the viscosity of the droplet; the droplet slightly deforms, as shown in Fig. 6. As seen in Fig. 0, when the Weber number was smaller than about unity, e was slightly smaller than unity. When the Weber number was larger than about unity, the experimental results of the first impact decreased sharply. The following empirical equations therefore were derived. e = 0.95 (We <.6). (8) e = We 0. (We.6). (9) The scattering of the data around Eq. (9) was 0%. The dynamic viscosity of water ranged from kg / (m s) to kg / (m s). 4.3 Internal flow of a droplet impacting on a hydrophobic plate In this subsection, we focus on the internal flow of a droplet impacting on a hydrophobic plate. Figures (a) ~ (e) show the dynamic behavior of small particles pre-mixed in the droplet. The volume of the droplet was 33.mm 3, and the Weber number was.. In the experiment, tracer particles of a density s of 040 kg/m 3 were mixed in the droplet. The mass percentage of the particles was 0. %. We confirmed that a small amount of the particles does not affect the dynamic behavior of the droplet. The behavior of the particles in the droplet can not be observed clearly in the experiment. This is because the free surface of the droplet is distorted due to the impact. In the numerical simulation, the particles were successively put in the droplet to observe the internal flow of the droplet impacting on the hydrophobic plate. The physical properties of the particles are identical with those used in the above-mentioned experiment. Figures (f) ~ (j) show the results obtained by the numerical simulation. The profile of the droplet at the symmetric face is drawn to understand an instantaneous situation of the impact motion of the droplet. At first, the droplet impacts on the hydrophobic plate (Fig. (f)), and then it deforms like a disk, as shown in Fig. (h) of the computational results. After deforming, the droplet is ringshaped pan-cake due to gravity, as seen in Fig. (i). The droplet returns to nearly spherical shape due to the surface tension (Fig. (j)) and then it bounces off the hydrophobic plate. 5. Conclusions This study has investigated the dynamic behavior of a water droplet impacting on and bouncing off a horizontal hydrophobic plate. The visualization techniques employed were the high-speed camera and the computation with the aid of the FLUENT numerical software. Main findings are summarized as follows: () The empirical equation proposed by Mao et al. was applicable to predicting the deformation rate, e d, for the Weber number greater than about 0. () When the Weber number was large, the deformation rate, e d, decreased with an increase in the number of impacts, N. The reason is that the droplet deforms and then dissipates energy due to the intense surface deformation. The values of e d, however, was kept constant when the Weber number was smaller than about unity. This is because the droplet hardly deforms. (3) The relationship between the Weber number, We, and the coefficient of restitution, e, was described by the following equations. The dynamic viscosity of water ranged from kg/(m s) to kg / (m s). e = 0.95 (We <.6). (8) e = We 0. (We.6). (9) (a) (b) (c) (d) (e) (f) (g) (h) (i) (j) Fig. Comparison of the particle trajectory in the droplet between (top) experimental results and (bottom) CFD ones. The volume of the droplet is 33. mm 3 and the Weber number is.. In the CFD, the photographs show the profile of the droplet at the symmetric face. The line in the droplet consists of the particles ( s = 040 kg/m 3 ). The particles did not affect the behavior of the droplet -SS9-

6 Y. SAKAI, Y. UEDA, T. ISHII and M. IGUCHI References [] Benintendi, S. W. and Smith, M. K.: The Spreading of a Non-Isothermal liquid Droplet, Phys. Fluids, (999), [] Xiying, L., Xuehu, M and Zhong, L.: Dynamic Behavior of the Water Droplet Impact on a Textured Hydrophobic/Superhydrophobic Surface: The Effect of the Remaining Liquid Film Arising on the Pillars' Tops on the Contact Time, Langmuir, 6(00), [3] Okumura, K., Chevy, F., Richard, D., Quéré, D., and Clanet, C.: Water Spring: A Model for Bouncing Drops, Europhys. Lett., 6 (003), [4] Chen, R. H. and Wang, H.W.: Effects of Tangential Speed on Low-Normal-Speed Liquid Drop Impact on a Non-Wettable Solid Surface, Exp. Fluids., 39(005), [5] Heinlein, J. and Fritsching, U.: Droplet Clustering in Sprays, Exp. Fluids, 40(006), [6] Meacham, J. M., Varady, M. J., Degertekin, F. L. and Fedorov, A. G.: Droplet Formation and Ejection from a Micromachined Ultrasonic Droplet Generator: Visualization and Scaling, Phys. Fluids, 7(005), 00605()-00605(8). [7] Ichikawa, K and Lin, W.: Interaction of Molten Steel with Refractories Alumina Clogging in Submerged Nozzles (in Japanese), CAMP-ISIJ, 3 (000), 70. [8] De Gennes, P. G., Brochard-Wyart, F. and Quéré, D.: Gouttes, Bulles, Perles et Ondes (translated by Okuyama, K.), Yoshioka Shoten, Kyoto, (003). [9] Kato, K., Azuma, T., Higashine, M. and Miyamoto, Y.: On the Sliding Down of Liquid Drops on Inclined Plates (st report, Critical Inclination Angle of Plates) (in Japanese), Trans. JSME, (B) 7-77(006), [0] Kato, K., Higashine, M., Nakamoto, N. and Azuma, T.: On the Sliding Down of Liquid Drops on Inclined Plates (nd Report, Sliding Angles on SAMs Surfaces) (in Japanese), Trans. JSME, (B) 73-73(007), [] Kato, K., Nishida, T. and Higashine, M.: On the Sliding of Liquid Droplets on Inclined Plates (Critical Volume of Liquid Droplets) (in Japanese), Trans. JSME, (B) (009), [] Sonoyama, N. and Iguchi, M.: Bubble Formation and Detachment on Nonwetted Surfaces, Metallurgical and Materials Transactions, 33B(00), [3] Kagawa, T. and Iguchi, M.: Model Experiment on the Production of Silicon Droplet, Materials Transactions, 45-8(004), [4] Akao, F., Araki, K., Mori, S. and Moriyama, A.: Deformation Behaviors of a Liquid Droplet Impinging onto Hot Metal Surface, Trans. ISIJ, 0(980), [5] Wachters, L. H. J. and Westerling, N. A. J.: The Heat Transfer from a Hot Wall to Impinging Water Drops in the Spheroidal State, Chem. Eng. Sci., (966), [6] Qiao, Y. M. and Chandra, S.: Boiling of Droplets on a Hot Surface in Low Gravity, Int. J. Heat Mass Trans., 39(996), [7] Park, H., Carr, W.W., Zhu, J. and Morris, J.: Single Drop Impaction on a Solid Surface, AIChE Journal, 49(003), [8] Li, X., Ma, X. and Lan, Z.: Behavioral Patterns of Drop Impingement onto Rigid Substrates with a Wide Range of Wettability and Different Surface Temperatures, AIChE Journal, 55(009), [9] Megaridis, C. M, Boomsma, K. and Bayer, I. S.: Partial Rebound of Molten-Metal Droplets Impacting on Solid Substrates, AIChE Journal, 50(004), [0] Gunjal, P. R., Ranade, V. V. and Chaudhari, R. V.: Dynamics of Drop Impact on Solid Surface: Experiments and VOF simulations, AIChE Journal, 5(005), [] Mao, T, Kuhn, D. C. S and Tran, H.: Spread and Rebound of Liquid Droplets upon Impact on Flat Surfaces, AIChE Journal, 43(997), [] Richard, D. and Quéré, D.: Bouncing Water Drops, Europhys. Lett., 50-6(000), [3] Fluent Inc.: FLUENT 6. User s Guide, Lebanon, NH (005). [4] Sakai, Y., Ueda, Y., Ishii, T. and Iguchi, M.: Behavior of a Droplet on an Inclined Plate under Various Wettability Conditions, J. JSEM., 0(00), 6-3. [5] Ohtake, T. and Fujita, S.: Kagaku-kikai, 3(949), (in Japanese). -SS30-

Maximal deformation of an impacting drop

J. Fluid Mech. (24), vol. 57, pp. 99 28. c 24 Cambridge University Press DOI:.7/S222494 Printed in the United Kingdom 99 Maximal deformation of an impacting drop By CHRISTOPHE CLANET, CÉDRIC BÉGUIN, DENIS