An empirically validated analytical model of droplet dynamics in electrowetting on dielectric devices

Explicit analytical models that describe the capillary force on confined droplets actuated in electrowetting on dielectric devices and the reduction in that force by contact angle hysteresis as a function of the three-dimensional shape of the droplet interface are presented. These models are used to develop an analytical model for the transient position and velocity of the droplet. An order of magnitude analysis showed that droplet motion could be modeled using the driving capillary force opposed by contact angle hysteresis, wall shear, and contact line friction. Droplet dynamics were found to be a function of gap height, droplet radius, surface tension, fluid density, the initial and deformed contact angles, contact angle hysteresis, and friction coefficients pertaining to viscous wall friction and contact line friction. The first four parameters describe the device geometry and fluid properties; the remaining parameters were determined experimentally. Images of the droplet during motion were used to determine the evolution of the shape, position, and velocity of the droplet with time. Comparisons between the measured and predicted results show that the proposed model provides good accuracy over a range of practical voltages and droplet aspect ratios.     

Contact-angle hysteresis on super-hydrophobic surfaces

The relationship between perturbations to contact angles on a rough or textured surface and the super-hydrophobic enhancement of the equilibrium contact angle is discussed theoretically. Two models are considered. In the first (Wenzel) case, the super-hydrophobic surface has a very high contact angle and the droplet completely contacts the surface upon which it rests. In the second (Cassie-Baxter) case, the super-hydrophobic surface has a very high contact angle, but the droplet bridges across surface protrusions. The theoretical treatment emphasizes the concept of contact-angle amplification or attenuation and distinguishes between the increases in contact angles due to roughening or texturing surfaces and perturbations to the resulting contact angles. The theory is applied to predicting contact-angle hysteresis on rough surfaces from the hysteresis observable on smooth surfaces and is therefore relevant to predicting roll-off angles for droplets on tilted surfaces. The theory quantitatively predicts a "sticky" surface for Wenzel-type surfaces and a "slippy" surface for Cassie-Baxter-type surfaces.     

Anisotropic wetting characteristics on submicrometer-scale periodic grooved surface

Submicrometer-scale periodic structures consisting of parallel grooves were prepared on azobenzene-containing multiarm star polymer films by laser interference. The wetting characteristics on the patterned surfaces were studied by contact angle measurements. Macroscopic distortion of water drops was found on such small-scale surface structures, and the contact angles measured from the direction parallel to the grooves were larger than those measured from the perpendicular direction. A thermodynamic model was developed to calculate the change in the surface free energy as a function of the instantaneous contact angle when the three-phase contact line (TPCL) moves along the two orthogonal directions. It was found that the fluctuations, i.e., energy barriers, on the energy versus contact angle curves are crucial to the analysis of wetting anisotropy and contact angle hysteresis. The calculated advancing and receding contact angles from the energy versus contact angle curves were in good agreement with those measured experimentally. Furthermore, with the groove depth increasing, both the degree of wetting anisotropy and the contact angle hysteresis perpendicular to the grooves increased as a result of the increase in the energy barrier. The theoretical critical value of the groove depth, above which the anisotropic wetting appears, was determined to be 16 nm for the grooved surface with a wavelength of 396 nm. On the other hand, the effect of the groove wavelength on the contact angle hysteresis perpendicular to the grooves was also interpreted on the basis of the thermodynamic model. That is, with the wavelength decreasing, the contact angle hysteresis increased due to the increase in the number of energy barriers. These results may provide theoretical evidence for the design and application of anisotropic wetting surface.     

Macroscopic-wetting anisotropy on the line-patterned surface of fluoroalkylsilane monolayers

Micropatterned fluoroalkylsilane monolayer surfaces with liquidphilic/liquidphobic area (line width 1-20 microm) were prepared with few defects by vacuum ultraviolet (VUV) photolithography. The anisotropic wetting of a macroscopic droplet with a 0.5-5 mm diameter on the micropatterned surfaces was investigated. The strong anisotropy of the contact angle and the sliding angle and droplet distortion for fluoroalkylsilane/silanol patterned surfaces was attributed to the difference in the energy barrier of wetting between parallel and orthogonal lines. The wetting anisotropy decreased with decreases in the liquidphilic area. Fluoroalkylsilane/alkylsilane patterned surfaces with small differences in the surface free energies of the components showed anisotropic wetting only for the low-surface-tension liquids.     

Patterned nonadhesive surfaces: superhydrophobicity and wetting regime transitions

Nonadhesive and water-repellent surfaces are required for many tribological applications. We study mechanisms of wetting of patterned superhydrophobic Si surfaces, including the transition between various wetting regimes during microdroplet evaporation in environmental scanning electron microscopy (ESEM) and for contact angle and contact angle hysteresis measurements. Wetting involves interactions at different scale levels: macroscale (water droplet size), microscale (surface texture size), and nanoscale (molecular size). We propose a generalized formulation of the Wenzel and Cassie equations that is consistent with the broad range of experimental data. We show that the contact angle hysteresis involves two different mechanisms and how the transition from the metastable partially wetted (Cassie) state to the homogeneously wetted (Wenzel) state depends upon droplet size and surface pattern parameters.     

Super-hydrophobic, highly adhesive, polydimethylsiloxane (PDMS) surfaces

Super-hydrophobic surfaces have been fabricated by casting polydimethylsiloxane (PDMS) on a textured substrate of known surface topography, and were characterized using contact angle, atomic force microscopy, surface free energy calculations, and adhesion measurements. The resulting PDMS has a micro-textured surface with a static contact angle of 153.5° and a hysteresis of 27° when using de-ionized water. Unlike many super-hydrophobic materials, the textured PDMS is highly adhesive, allowing water drops as large as 25.0 μL to be inverted. This high adhesion, super-hydrophobic behavior is an illustration of the "petal effect". This rapid, reproducible technique has promising applications in transport and analysis of microvolume samples.     

Contact angle hysteresis of cylindrical drops on chemically heterogeneous striped surfaces

Contact angle hysteresis of a macroscopic droplet on a heterogeneous but flat substrate is studied using the interface displacement model. First, the apparent contact angle of a droplet on a heterogeneous surface under the condition of constant volume is considered. By assuming a cylindrical liquid-vapor surface (meniscus) and minimizing the total free energy, we derive an equation for the apparent contact angle, which is similar but different from the well-known Cassie's law. Next, using this modified Cassie's law as a guide to predict the behavior of a droplet on a heterogeneous striped surface, we examine several scenarios of contact angle hysteresis using a periodically striped surface model. By changing the volume of the droplet, we predict a sudden jump of the droplet edge, and a continuous change of the apparent contact angle at the edge of two stripes. Our results suggest that as drop volume is increased (advancing contact lines), the predominant drop configuration observed is the one whose contact angle is large; whereas, decreasing drop volume from a large value (receding contact lines) yields drop configuration that predominantly exhibit the smaller contact angle.     

Static and dynamic contact angles of water droplet on a solid surface using molecular dynamics simulation

The present study investigates the variation of static contact angle of a water droplet in equilibrium with a solid surface in the absence of a body force and the dynamic contact angles of water droplet moving on a solid surface for different characteristic energies using the molecular dynamics simulation. With increasing characteristic energy, the static contact angle in equilibrium with a solid surface in the absence of a body force decreases because the hydrophobic surface changes its characteristics to the hydrophilic surface. In order to consider the effect of moving water droplet on the dynamic contact angles, we apply the constant acceleration to an individual oxygen and hydrogen atom. In the presence of a body force, the water droplet changes its shape with larger advancing contact angle than the receding angle. The dynamic contact angles are compared with the static contact angle in order to see the effect of the presence of a body force.     

Condensation-induced wetting state and contact angle hysteresis on superhydrophobic lotus leaves

In this paper, we demonstrate how condensed moisture droplets wet classical superhydrophobic lotus leaf surfaces and analyze the mechanism that causes the increase of contact angle hysteresis. Superhydrophobic lotus leaves in nature show amazing self-cleaning property with high water contact angle (>150°) and low contact angle hysteresis (usually <10°), causing droplets to roll off at low inclination angles, in accordance with classical Cassie–Baxter wetting state. However, when superhydrophobic lotus leaves are wetted with condensation, the condensed water droplets are sticky and exhibit higher contact angle hysteresis (40–50°). Compared with a fully wetted sessile droplet (classical Wenzel state) on the lotus leaves, the condensed water droplet still has relatively large contact angle (>145°), suggesting that the wetting state deviates from a fully wetted Wenzel state. When the condensed water droplets are subjected to evaporation at room conditions, a thin water film is observed bridging over the micropillar structures of the lotus leaves. This causes the dew to stick to the surface. This result suggests that the condensed moisture does not uniformly wet the superhydrophobic lotus leaf surfaces. Instead, there occurs a mixed wetting state, between classical Cassie–Baxter and Wenzel states that causes a distinct increase of contact angle hysteresis. It is also observed that the mixed Cassie–Baxter/Wenzel state can be restored to the original Cassie–Baxter state by applying ultrasonic vibration which supplies energy to overcome the energy barrier for the wetting transition. In contrast, when the surface is fully wetted (classical Wenzel state), such restoration is not observed with ultrasonic vibration. The results reveal that although the superhydrophobic lotus leaves are susceptible to being wetted by condensing moisture, the configured wetting state is intermediate between the classical Cassie–Baxter and Wenzel states.     

Sliding of water droplets on hydrophobic surfaces with various hydrophilic region sizes

Four patterned surfaces with hydrophilic areas of different sizes were prepared using photolithography with a smooth octadecyltrimethoxysilane (ODS) hydrophobic coating. The hydrophilic area in the surfaces was aligned hexagonally with a constant area fraction. The sliding angle and contact angle hysteresis of the water droplets increased concomitantly with increasing pattern size. The increase of the contact line distortion between defects at the receding side plays an important role in this trend. The droplet sliding velocity also increased concomitantly with increasing pattern size. This trend was simulated by a simple flow model. The contribution of the interface between the ODS region and the hydrophilic area was deduced from this trend. This study demonstrated the different size dependency of the chemical surface defects for sliding behavior between the critical moment at which a droplet slides down and the period when a droplet is sliding.     

Theoretical consideration of wetting on a cylindrical pillar defect: pinning energy and penetrating phenomena

Wetting on a cylindrical pillar defect is discussed in terms of the free-energy difference ΔG. Wetting is divided into wetting on a flat surface, a pinning effect at the apex of the defect, and wetting on a pillar wall. First, we confirmed that ΔG between before and after ideal wetting on a flat surface can be derived as a function of the contact angle θ in which the free-energy minimum is obtained as the equilibrium contact angle θ(eq) described by Young's and Wenzel's laws. Second, the pinning effect at the apex in the cross section of the pillar defect is discussed in ΔG, where the pinning effect is shown to originate from the energy barrier by an increase in the air-liquid interfacial area of a pinned droplet induced by deformation. Next, the ΔG profiles of wetting on the pillar wall are drawn based on the theory of Carroll (Carroll, B. J. J. Colloid Interface Sci.1976, 57, 488-495) to better understand the ΔG profile during penetration. Differences in the manner of wetting between the wetting state on a flat surface and the pillar wall are reflected in ΔG. Finally, penetration of a droplet into a pillar defect is comprehensively discussed on the basis of wetting on a flat surface and a pillar wall. If we consider a simple manner of penetration, another type of energy barrier resulting from an anomalous deformation of the air-liquid interface of the penetrating droplet can be theoretically suggested. Consequently, two types of energy barrier are found. These energy barriers should play a significant role in the hysteresis of wetting, the liquid-repellent Cassie-Baxter state (CB), and the CB-Wenzel wetting transition on a microtextured surface.     

Contact angle hysteresis on regular pillar-like hydrophobic surfaces

A series of pillar-like patterned silicon wafers with different pillar sizes and spacing are fabricated by photolithography and further modified by a self-assembled fluorosilanated monolayer. The dynamic contact angles of water on these surfaces are carefully measured and found to be consistent with the theoretical predictions of the Cassie model and the Wenzel model. When a water drop is at the Wenzel state, its contact angle hysteresis increases along with an increase in the surface roughness. While the surface roughness is further raised beyond its transition roughness (from the Wenzel state to the Cassie state), the contact angle hysteresis (or receding contact angle) discontinuously drops (or jumps) to a lower (or higher) value. When a water drop is at the Cassie state, its contact angle hysteresis strongly depends on the solid fraction and has nothing to do with the surface roughness. Even for a superhydrophobic surface, the contact angle hysteresis may still exhibit a value as high as 41 degrees for the solid fraction of 0.563.     

Model description of contact angles in electrowetting on dielectric layers

The electrowetting on dielectric (EWOD) technique has considerable potential for microfluidic and biomedical applications. The Lippmann-Young model based on the force balance concept has long been used to predict the contact angles of droplets under electrowetting. However, recent experimental evidence has indicated that this model fails to provide accurate predictions of the lower contact angles associated with saturation conditions at higher electric potentials. Hence, the study simulates the internal flow in an actuated droplet and treats it as stagnation-point flow. This kinetic energy is then taken into consideration while calculating the contact angles using an energy balance model. The energy of an actuated droplet is contributed by the combination of the side surface tension energy, the base tension energy, the dielectric energy, and the kinetic energy when deriving the energy balance model. Consequently, the new energy balance model modifies the Lippmann-Young equation, thereby providing enhanced reasonable predictions of the droplet contact angle across the higher electric potential where the contact angles are close to the saturated condition.     

低能量链端对聚苯乙烯薄膜铺展和润湿行为的影响

通过低能量功能端基的表面富集作用,研究了聚苯乙烯（ＰＳ）薄膜在聚甲基丙烯酸甲酯（ＰＭＭＡ）表面上的铺展和润湿动力学．用光学显微镜跟踪了ＰＳ薄膜的润湿行为,并对高分子熔体膜中非连续部分尺寸的增大速率进行了测定．分别用ＸＰＳ和ＡＦＭ对ＰＳ薄膜的表面组成和ＰＳ液滴的平衡接触角进行了测定．发现具有低表面能的氟碳端基在薄膜表面富集使ＰＳ薄膜的表面张力下降,并使ＰＳ液滴在ＰＭＭＡ表面上的平衡接触角减小,从而使高分子熔体膜中非连续部分尺寸的增长速率下降,得到了与液液界面铺展和润湿理论一致的实验结果．     

Modeling contact angle hysteresis of a liquid droplet sitting on a cosine wave-like pattern surface

A liquid droplet sitting on a hydrophobic surface with a cosine wave-like square-array pattern in the Wenzel state is simulated by using the Surface Evolver to determine the contact angle. For a fixed drop volume, multiple metastable states are obtained at two different surface roughnesses. Unusual and non-circular shape of the three-phase contact line of a liquid droplet sitting on the model surface is observed due to corrugation and distortion of the contact line by structure of the roughness. The contact angle varies along the contact line for each metastable state. The maximum and minimum contact angles among the multiple metastable states at a fixed viewing angle correspond to the advancing and the receding contact angles, respectively. It is interesting to observe that the advancing/receding contact angles (and contact angle hysteresis) are a function of viewing angle. In addition, the receding (or advancing) contact angles at different viewing angles are determined at different metastable states. The contact angle of minimum energy among the multiple metastable states is defined as the most stable (equilibrium) contact angle. The Wenzel model is not able to describe the contact angle along the three-phase contact line. The contact angle hysteresis at different drop volumes is determined. The number of the metastable states increases with increasing drop volume. Drop volume effect on the contact angles is also discussed.     

Contact angle hysteresis at the nanoscale: a molecular dynamics simulation study

In this paper, the contact angle hysteresis (CAH) of nanodroplets on both rigid and flexible substrates with different wettabilities was studied using molecular dynamics (MD) simulations. The critical shear stress (CSS) that determines the motion of the contact line (CL) was investigated. A theoretical correlation between CAH and CSS was proposed. Both CAH and CSS reflect the energy dissipation at the CL of the droplet in response to the exerted force. MD results of CAH are qualitatively consistent with the theoretical model. Simulation results also show that, for the same liquid–solid interactions, CAH on the flexible substrate is larger than that on the rigid substrate. These findings aim to enhance our understanding of the mechanism of the CAH at the nanoscale.     

Electrowetting-based control of static droplet states on rough surfaces

Electrowetting (EW) is a powerful tool to control fluid motion at the microscale and has promising applications in the field of microfluidics. The present work analyzes the influence of an electrowetting voltage in determining and altering the state of a static droplet resting on a rough surface. An energy-minimization-based modeling approach is used to analyze the influence of interfacial energies, surface roughness parameters, and electric fields in determining the apparent contact angle of a droplet in the Cassie and Wenzel states under the influence of an EW voltage. The energy-minimization-based approach is also used to analyze the Cassie-Wenzel transition under the influence of an EW voltage and estimate the energy barrier to transition. The results obtained show that EW is a powerful tool to alter the relative stabilities of the Cassie and Wenzel states and enable dynamic control of droplet morphology on rough surfaces. The versatility and generalized nature of the present modeling approach is highlighted by application to the prediction of the contact angle of a droplet on an electrowetted rough surface consisting of a dielectric layer of nonuniform thickness.     

Frictional behaviors of three kinds of nanotextured surfaces

Nanodot‐textured surface, nanorod‐textured surface and nanocomposite‐textured surface were prepared by the hydrothermal technique and successive ion layer absorption and reaction technique. The adhesion and friction properties of the three kinds of nanotextured surfaces were investigated using an atomic force microscope colloidal probe. Experimental results revealed that the nanorod‐textured surface and nanocomposite‐textured surface can significantly reduce adhesive and friction forces compared with a nanodot‐textured surface. The main reason for this phenomenon was that the nanorod and nanocomposite textures can reduce contact area between the sample surface and the colloidal probe. The effects of surface root mean square roughness, applied load and sliding velocity on the adhesion and friction behaviors of the nanotextured surfaces were investigated. The adhesive and friction forces of the nanotextured surfaces decreased with the increasing surface root mean square roughness. Compared with the nanocomposite‐textured surface, the nanorod‐textured surface has better adhesion and frictional performance. Copyright © 2016 John Wiley & Sons, Ltd.     

Shape-gradient composite surfaces: water droplets move uphill

The approach of water droplets self-running horizontally and uphill without any other forces was proposed by patterning the shape-gradient hydrophilic material (i.e., mica) to the hydrophobic matrix (i.e., wax or low-density polyethylene (LDPE)). The shape-gradient composite surface is the best one to drive water droplet self-running both at the high velocity and the maximal distance among four different geometrical mica/wax composite surfaces. The driving force for the water droplets self-running includes: (1) the great difference in wettability of surface materials, (2) the low contact angle hysteresis of surface materials, and (3) the space limitation of the shape-gradient transportation area. Furthermore, the average velocity and the maximal distance of the self-running were mainly determined by the gradient angle (alpha), the droplet volume, and the difference of the contact angle hysteresis. Theoretical analysis is in agreement with the experimental results.     

Adhesion and debonding of soft elastic films: crack patterns, metastable pathways, and forces

We study the phenomenon of debonding in a thin soft elastic film sandwiched between two rigid plates as one of the plates is brought into intimate contact and then pulled away from contact proximity by application of a normal force. Nonlinear simulations based on minimization of total energy (composed of stabilizing elastic strain energy and destabilizing adhesive interaction energy) are employed to address the problems of contact hysteresis, cavitation, crack morphology, variation of contact area, snap-off distance, pull-off force, work done, and energy loss. Below a critical distance (d(c)) upon approach, simulations show the formation of columnar structures and nonrandom, regularly arranged nanocavities at the soft interface at a length scale of approximately 3h (h being the thickness of the film). The persistence of such instability upon withdrawal (distance >d(c)) indicates a contact hysteresis, which is caused by an energy barrier that separates the metastable states of the patterned configuration and the global minimum state of the flat film. The energy and the pull-off force are found to be nonequilibrium and nonunique properties depending on the initial contact, defects, noise, etc. Three broad pathways of debonding leading to adhesive failure of the interface, depending on the stiffness of the film, step size of withdrawal, and the imposed noise, are identified: a catastrophic column collapse mode, a peeling mode involving a continuous decrease in the contact area, and a column splitting mode. The first two modes are caused by a very high stress concentration near the cavity edges. These metastable patterned configurations engender pull-off forces that are orders of magnitude smaller than that required to separate two flat surfaces from contact.