Size effect on the coalescence-induced self-propelled droplet
Feng-Chao Wang, Fuqian Yang, and Ya-Pu Zhao
Citation: Appl. Phys. Lett. 98, 053112 (2011); doi: 10.1063/1.3553782
View online: http://dx.doi.org/10.1063/1.3553782
View Table of Contents: http://apl.aip.org/resource/1/APPLAB/v98/i5
Published by the American Institute of Physics.
Related Articles
A parameter-free, solid-angle based, nearest-neighbor algorithm
J. Chem. Phys. 136, 234107 (2012)
Surface charge density model for predicting the permittivity of liquid mixtures and composites materials
J. Appl. Phys. 111, 064101 (2012)
Energy dissipation in non-isothermal molecular dynamics simulations of confined liquids under shear
J. Chem. Phys. 135, 134708 (2011)
Extension of the Steele 10-4-3 potential for adsorption calculations in cylindrical, spherical, and other pore
geometries
J. Chem. Phys. 135, 084703 (2011)
Melting of iron at the Earth's core conditions by molecular dynamics simulation
AIP Advances 1, 032122 (2011)
Additional information on Appl. Phys. Lett.
Journal Homepage: http://apl.aip.org/
Journal Information: http://apl.aip.org/about/about_the_journal
Top downloads: http://apl.aip.org/features/most_downloaded
Information for Authors: http://apl.aip.org/authors
Downloaded 17 Sep 2012 to 166.111.54.144. Redistribution subject to AIP license or copyright; see http://apl.aip.org/about/rights_and_permissions
Size effect on the coalescence-induced self-propelled droplet
Feng-Chao Wang,1 Fuqian Yang,2,a� and Ya-Pu Zhao1,a�
1State Key Laboratory of Nonlinear Mechanics, Institute of Mechanics, Chinese Academy of Sciences,
Beijing 100190, People’s Republic of China
2Department of Chemical and Materials Engineering, University of Kentucky, Lexington,
Kentucky 40506, USA
�Received 8 November 2010; accepted 21 January 2011; published online 4 February 2011�
An analysis based on the energy conservation is presented for the self-propelled droplet during
coalescence of two droplets of the same size over a superhydrophobic rough surface. The
self-propelled behavior occurs only for the coalescence of droplets with a certain range of radius. An
analytical relation is established among the coalescence-induced velocity, surface energy, viscous
dissipation, and droplet size if gravity is negligible. The coalescence-induced velocity increases with
increasing droplet size to a maximum and then decreases with the size, which is in good accord with
the experimental observation reported in the literature. © 2011 American Institute of Physics.
�doi:10.1063/1.3553782�
Depending on the bonding energy between liquid mol-
ecules and surface atoms/molecules of a solid, a solid surface
can be wettable or nonwettable. Through chemical treatment,
such as the grafting of hydrogenated or fluorinated chlorosi-
lanes onto the Si–OH groups of a glass surface,1 one can
control the wettability of a solid surface. In general,
wettability determines the detachment and wetting of drop-
lets or bubbles on a solid surface, which is of practical im-
portance for many applications in self-cleaning, microfluid-
ics, lab-on-a-chip devices, and liquid-based cooling of
microelectronics.
Recently, there has been a great interest in developing
techniques to achieve the transition of the Wenzel state2 to
the Cassie–Baxter state3 for separating liquid from solid sur-
faces, which include vibration-induced dewetting4 and im-
pulse heating.5 All these works involved the use of external
energy, such as mechanical energy and electric energy, to
overcome the energy barrier between these two states. To
reduce the energy barrier and the wetting of liquid, hydro-
phobic or superhydrophobic surfaces are desired and differ-
ent methods have been developed to form synthetic superhy-
drophobic surfaces. The use of superhydrophobic surfaces in
the dropwise condensation has led to the observation of the
rapid motion of liquid droplet driven by the coalescence.6–9
Such a phenomenon is believed to be powered by the surface
energy released during the drop coalescence. Boreyko and
Chen6 did a simple scaling analysis by assuming that all the
released surface energy is converted to kinetic energy of the
merged droplet. However, a quantitative understanding of the
relationship between energy and kinetics has not been re-
sulted.
In this work, a quantitative approach is proposed,
whereby the kinetics of the coalescent droplet is governed by
an equilibrium balance of surface energy, kinetic energy, and
the viscous dissipation energy. The gravitational-potential
energy can be neglected in the analysis since the scale of the
droplets is much smaller than the capillary length
lc=��lv /�g �2.7 mm for water�, in which g is the gravita-
tional acceleration, � is the density of the liquid, and �lv is
the surface tension of the liquid-vapor interface. The interac-
tion energy between two droplets with radius of r1 and r2,
respectively, can be calculated as W=−AHr1r2 /6D�r1+r2�,10
where AH is the Hamaker constant and D is the separation
distance between the two droplets. The interaction energy is
about six to seven orders of magnitude less than the surface
energy and the viscous dissipation energy. Thus, the interac-
tion energy between two droplets is negligible in the analy-
sis.
The concept developed will be for the coalescence of
two quasistationary liquid droplets of the same size on a
superhydrophobic rough surface. The method can be applied
with some modifications to more complicated situations in-
volving the coalescence of multiple droplets. Consider two
droplets of the same radius r, which coalesce to form a new
droplet of radius R, as shown in Fig. 1. Under the assumption
of the absence of gravity �the typical Bond number Bo
=�gr2 /�lv�1�, the equilibrium shape of a droplet is a
spherical cap with an apparent contact angle �.
The surface energy of each stationary droplet on a hy-
drophobic rough surface is given by
ES = �lvAlv + �slAsl + �svAsv, �1�
where A is the interfacial area and the subscripts s, l, and v
denote the solid, liquid, and vapor, respectively. For a rough
a�Authors to whom correspondence should be addressed. Electronic ad-
dresses: fyang0@engr.uky.edu and yzhao@imech.ac.cn.
FIG. 1. �Color online� Schematic of the coalescence of two droplets to one
big droplet on a superhydrophobic rough surface. The transition of the Wen-
zel state to the Cassie–Baxter state is shown.
APPLIED PHYSICS LETTERS 98, 053112 �2011�
0003-6951/2011/98�5�/053112/3/$30.00 © 2011 American Institute of Physics98, 053112-1
Downloaded 17 Sep 2012 to 166.111.54.144. Redistribution subject to AIP license or copyright; see http://apl.aip.org/about/rights_and_permissions
Administrator
高亮
Administrator
高亮
Administrator
高亮
Administrator
高亮
surface, the surface energy can be written as11
ES = �lv�r2�2 − 2 cos � −��f�sin2 �� + �svAtot, �2�
in which Atot is the total area of the solid surface and
��f� = rf f cos �Y + f − 1, �3�
where rf is the roughness ratio of the wet surface area, f is
the fraction of the solid surface area wetted by the liquid, and
�Y is the Young contact angle as defined for an ideal surface.
Before coalescence, the droplets are in the Wenzel state,
which means f =1. A transition from the Wenzel to the
Cassie–Baxter state due to the coalescence can be observed.6
Thus, the surface energy of the droplet after coalescence has
a similar expression as Eq. �2�, while the parameters rf and f
should be changed. �W�f� and �C�f� would be used to dis-
tinguish the coefficients for the Wenzel and the Cassie–
Baxter state, respectively.
During the coalescence, there exists the energy dissipa-
tion due to the flow and mergence of the liquid droplets
against viscosity. Approximately, the viscous dissipation en-
ergy for each droplet can be estimated as12
Evis = �
0
��
�
d �dt �
� , �4�
in which
is the dissipation function
=
�
2 �vi�xj + �v j�xi
2
� 12� U
r
2, �5�
is the volume of each droplet, � is the viscosity of the
liquid, and the characteristic capillary time scale � is given
by
�� ��r3/�lv. �6�
As the coalescence starts, the capillary pressure inside the
droplet ��p=2�lv /r� will accelerate the droplet along the
horizontal direction.13 Thus, the average velocity U of each
droplet can be obtained as
U � � · �p · �r2
1
4��r3/3
=
3
2
��lv
�r
. �7�
Combining Eqs. �4�–�7� yields
Evis � 36����lvr3
�
. �8�
The total energy of the system after coalescence consists of
surface energy and kinetic energy of the coalescent droplet.
Then, balancing the change in the total energy of the system
due to the coalescence and introducing the mass conserva-
tion, i.e., R3=2r3, one obtains the kinetic energy of the coa-
lescent droplet as
EK = �lv�r2��2 − 22/3��2 − 2 cos �� + �22/3�C�f�
− 2�W�f��sin2 �� − 64����lvr3
�
, �9�
which depends on the size of the starting droplets and the
wetting condition of the rough surface. Some portions of the
surface energy released upon the coalescence will be used
to balance the viscous dissipation energy. To have the
coalescence-induced self-propelled droplet, it requires EK
0. It can be seen from Eq. �9� that the coalescence-induced
kinetic energy is the result of the competition between the
viscous energy and surface energy, which can be described
by the Ohnesorge number On=� /���lvr.14 For droplets of
small size �submicron�, the Ohnesorge number is larger than
the critical value 0.1,15 which indicates that the viscous dis-
sipation dominates during the coalescence and, consequently,
there is no energy available for the self-propelled behavior.
Interestingly, this corresponds to the “immobile coalescence”
reported by Boreyko and Chen.6
Once the self-propelled phenomenon emerges, the
coalescence-induced velocity V of the coalescent droplet can
be found from Eq. �9� as
V =��lv
�r
38C�f ,�� − 24 ����lvr�
1/2
, �10�
where C�f ,�� is the term in the curly brackets of Eq. �9�.
The water properties at 20 °C are used for the following
calculation, such as �=998.23 kg /m3, �lv=72.75 mN /m,
and �=1.0087 mPa s.16 Without loss of generality, the other
parameters are assumed to be �=170°, �Y =160°, rfW=1.5,
fW=1.0, rfC=1.0, and fC=0.5. Moreover, the choice of these
parameters has only a limited effect on the final results and
does not change the trend observed from the analysis. Using
these parameters, one can calculate the coalescence-induced
velocity of water droplets. Figure 2 shows the dependence of
the coalescence-induced velocity on the radius of water drop-
let. Obviously, the self-propelled behavior is not present for
droplets of small size. This is due to the larger Ohnesorge
number than the critical value 0.1. For water droplets of large
size, the coalescence-induced velocity first increases with the
droplet radius to a maximum velocity and then decreases
with the droplet radius. It should be noted that our analysis is
based on the assumption that the gravity can be neglected.
For droplets of larger size than the capillary length lc, the
variation of the gravitational-potential energy likely con-
sumes the surface energy released upon the coalescence of
droplets and then no coalescence-induced self-propelled be-
havior can be observed.
Our modeling analysis, with the viscous dissipation in-
volved, captures several features of the experimental results
FIG. 2. Dependence of the coalescence-induced velocity on the radius of the
droplet �the parameters used in the calculation are �=170°, �Y =160°, rfW
=1.5, fW=1.0, rfC=1.0, and fC=0.5; the properties of water at 20 °C are
used, such as �=998.23 kg /m3, �lv=72.75 mN /m, and �=1.0087 mPa s
�Ref. 16�.
053112-2 Wang, Yang, and Zhao Appl. Phys. Lett. 98, 053112 �2011�
Downloaded 17 Sep 2012 to 166.111.54.144. Redistribution subject to AIP license or copyright; see http://apl.aip.org/about/rights_and_permissions
reported by Boreyko and Chen,6 which includes: �i� the
coalescence-induced velocity increases with the droplet size
to a maximum and then decreases with the size, �ii� the
maximum velocity can be found for the coalescence of drop-
lets with the radius of about 50 �m, which is in accord with
the experimental observation, and �iii� although the theoret-
ical velocities are about two to three times larger than those
measured in Boreyko and Chen’s experiments, they are on
the same order of magnitude.
In summary, this work has attempted to emphasize the
importance of the viscous dissipation in examining the self-
propelled behavior during the coalescence of droplets over a
superhydrophobic rough surface. The theoretical analysis re-
veals that the self-propelled behavior can only occur for the
coalescence of droplets with the radius ranging from several
microns to a few millimeters, in which the surface energy
dominates compared with the viscous dissipation and the
gravitational-potential energy. In particular, the study has
shown that an equilibrium approach involving the conserva-
tion of energy can lead to a prediction of the coalescence-
induced velocity for the coalescence of two quasistationary
liquid droplets of the same size. The coalescence-induced
velocity is a function of the droplet size, which is consistent
with the experimental results. By modifying the surface
properties of liquid droplets through chemical treatment and
controlling the droplet size, one can use the coalescence-
induced self-propelled behavior for the applications in self-
cleaning, microfluidics, and lab-on-a-chip devices.
This work is jointly supported by the National Natural
Science Foundation of China �NSFC, Grant Nos. 60936001
and 11072244�, the National Basic Research Program of
China �973 Program, Grant No. 2007CB310500�, and Open-
ing Fund of State Key Laboratory of Nonlinear Mechanics.
1F. Brochard-Wyart, in Soft Matter Physics, edited by M. Daoud and C. E.
Williams �Springer, New York, 1999�.
2R. Wenzel, Ind. Eng. Chem. 28, 988 �1936�.
3A. Cassie and S. Baxter, Trans. Faraday Soc. 40, 546 �1944�.
4J. B. Boreyko and C. H. Chen, Phys. Rev. Lett. 103, 174502 �2009�.
5T. N. Krupenkin, J. A. Taylor, E. N. Wang, P. Kolodner, M. Hodes, and T.
R. Salamon, Langmuir 23, 9128 �2007�.
6J. B. Boreyko and C. H. Chen, Phys. Rev. Lett. 103, 184501 �2009�.
7C. H. Chen, Q. J. Cai, C. L. Tsai, C. L. Chen, G. Y. Xiong, Y. Yu, and Z.
F. Ren, Appl. Phys. Lett. 90, 173108 �2007�.
8C. Dorrer and J. Ruhe, Adv. Mater. 20, 159 �2008�.
9B. Mockenhaupt, H.-J. Ensikat, M. Spaeth, and W. Barthlott, Langmuir
24, 13591 �2008�.
10J. N. Israelachvili, Intermolecular and Surface Forces �Academic, San
Diego, 1992�.
11A. Marmur, Langmuir 19, 8343 �2003�.
12S. Chandra and C. T. Avedisian, Proc. R. Soc. London, Ser. A 432, 13
�1991�.
13S. T. Thoroddsen and K. Takehara, Phys. Fluids 12, 1265 �2000�.
14F. C. Wang, J. T. Feng, and Y. P. Zhao, J. Colloid Interface Sci. 326, 196
�2008�.
15R. Q. Duan, S. Koshizuka, and Y. Oka, Nucl. Eng. Des. 225, 37 �2003�.
16CRC Handbook of Physics and Chemistry, 63rd ed., edited by R. C. Weast
and M. J. Astle �CRC, Boca Raton, 1982�.
053112-3 Wang, Yang, and Zhao Appl. Phys. Lett. 98, 053112 �2011�
Downloaded 17 Sep 2012 to 166.111.54.144. Redistribution subject to AIP license or copyright; see http://apl.aip.org/about/rights_and_permissions
Administrator
高亮
Administrator
高亮
Administrator
高亮