Design Principle of Mixed-Wettability Surfaces for Inducing Droplet Directional Rebound

，，，，4*

1.School of Automotive and Traffic Engineering，Jiangsu University，Zhenjiang 212013，P.R.China；

2.R&D Center，Jiangsu Chaoli Electric Co.Ltd.，Zhenjiang 212321，P.R.China；

3.College of Aerospace Engineering，Nanjing University of Aeronautics and Astronautics，Nanjing 210016，P.R.China；

4.National Key Laboratory of Science and Technology on Aerodynamic Design and Research，Northwestern Polytechnical University，Xi’an 710072，P.R.China

Abstract: Directional rebound of droplet impacting on solid surfaces is of great significance in engineering applications such as anti-icing/fogging and self-cleaning.Mixed-wettability surfaces have been shown to be as an effective way for droplet manipulation.The impact of droplet on hydrophobic substrate decorated with a hydrophilic stripe is investigated numerically.The validated diffuse interface method is adopted for interface capture.First，the formation of satellite droplets during the impact process is explored，and the role of mixed-wettability surfaces in the droplet spreading，retraction and rebound stages is clarified by analyzing the vertical and lateral velocities of the droplet.After that，the effect of the stripe width on the droplet rebound form and the contact time is studied systematically.Special attention is paid to the mechanisms of dynamics and energy transfer during the evolution of liquid film and droplet bouncing.The obtained results can provide useful guidance for the design of mixed-wettability surfaces and further optimize the control of droplet directional rebound.

Key words：droplet impact；mixed-wettability surfaces；directional rebound；design principle

0 Introduction

Droplet impacting onto solid surfaces is a ubiquitous and important phenomenon in industrial processes，such as spray cooling，inkjet printing，selfcleaning and fluid transportation［1-6］.One of the key factors affecting the dynamic behavior of droplets is the surface wettability，and the droplets generally tend to rebound from the hydrophobic surfaces and adhere to the hydrophilic surfaces.Extensive works mainly focused on the contact timetcand spreading diameter of the droplet，which may produce a marked effect on the industrial process.In general，it is advantageous to reduce the contact time between the droplet and the surface，since it controls the exchange of mass，momentum and energy［7-9］.For example，rapid droplet detachment can effectively reduce the heat exchange and suppress the ice nucleation，thus preventing the aircraft and wire from icing at low temperature.Richard et al.［8］found thattconly scales with the inertia-capillarity time，whereρ1andR0are the density and the initial radius of the droplet，respectively，andσis the surface tension coefficient.The theoretical limit ofis approximately equal to 2.2 when the droplet undergoes the spread-retraction-rebound stage on flat superhydrophobic surfaces.Recently，the researchers reported that adding single or array structures on surfaces can reduce the contact time by 37%—80%，such as pancake bouncing［10］.

Although the contact time can be reduced by the chemical or physical treatments of the surface，the droplets usually rebound along the original route（in a perpendicular type）.This may limit its engineering applications，since the droplets will eventually fall back and adhere to the horizontal surface with the dissipation of kinetic energy.Therefore，the demand for droplet rebound with lateral motion（deflection）emerges，and many attempts have been made in this aspect.The intuitive idea is to exert a lateral force on the droplet during the rebound process，such as introducing the wind-driven，electric or magnetic fields［11-14］.However，these methods with external power input are sometimes unstable and hard to control.Recently，the use of engineering surfaces，such as roughness gradient，wettability gradient and bioinspired structures，to realize the spontaneous lateral motion of droplets has attracted wide attention［15-18］.

The development of mixed-wettability surfaces provides a fresh insight into the droplet manipulation，since the forces exerted on various wettability partitions are different.Farshchian et al.［19］reported a directional rebounding phenomenon when the droplet impacted on the borderline between the hydrophobic and hydrophilic portions of a heterogeneous surface.The experimental results showed that the lateral transportation distance can increase with the Weber number（We）in the range of 3.2 to 18.9.Chu et al.［20］further investigated the effect of impacting offset position on the deflection distance at higherWe.The maximum lateral transportation distance can be up to 13 times longer than the droplet diameter.In these studies，the nonorthogonal rebound of the droplet is realized by impacting on the junction region of a dual-textured substrate，and the lateral motion is mainly caused by asymmetric forces generated at the retraction stage.Based on this mechanism，the researchers tried to adopt the hydrophobic substrate decorated with hydrophilic patterns to control the droplets.Schutzius et al.［21］and Zhao et al.［22］designed a hydrophilic “arc” and “stripe” on a superhydrophobic surface，respectively.The impacting droplet could rebound toward the direction of the hydrophilic patterns instead of orthogonal to the surface.Then the simple practical applications were achieved by the well-controlled droplet motion，such as the on-demand chemical reaction system and droplet-driven micro-boat system.Furthermore，Li et al.［23］realized more various rebound behaviors of the droplets，such as gyration，rolling，and deflection，by using the heterogeneous surfaces with different high-adhesion（hydrophilic）patterns.In addition to experimental studies，Russo et al.［24］and Zhang et al.［25］recently numerically investigated the droplet impacting on a wettability-patterned surface.The numerical method coupled with the Kistler dynamic contact angle model and adaptive meshrefinement technique was fully validated，and the dynamic mechanisms of droplet rebound were revealed by analyzing the variation of the contact line and velocity profile.

In the abovementioned studies，various forms of hydrophilic patterns were designed，and the results showed that their effects on the dynamic behavior of droplets were nonlinear and complex.Taking the distanceabetween the impacting point and hydrophilic pattern as an example，the lateral velocityvxof the droplet continuously increases withafor the arc-shaped pattern，whilevxincreases first and then decreases with the increase ofafor the stripeshaped pattern.To clarify the relationship between hydrophilic patterns and the droplet lateral motion，the parameter of “effective retraction area”Sis proposed to quantify the geometric parameters and position of different patterns in our previous work［25］.Furthermore，we systematically investigated the effects of multiple parameters on the dynamic behaviors of the droplet impact，and some scaling laws for droplet rebound with the lateral motion were proposed［26］.Thus，the prediction and control of droplet directional rebound were preliminarily achieved.However，the design principle of mixed-wettability surfaces has never been established due to the following two problems.First，small satellite droplets may form during the impact process，which has occurred in both experimental observations and numerical simulations，but has not received much attention［21，24］.The energy transfer characteristics between the main and satellite droplets are not known yet，so the role of mixed-wettability surfaces in the formation of satellite droplets remains unclear.Second，although various forms of hydrophilic patterns have been studied in previous works，they were all carried out when the width of patterns was about 10% of the droplet diameter.It is not clear whether a better pattern width exists for droplet directional rebound.And the effect of the size（width）of the hydrophilic pattern on the adhesive force and droplet rebound modes（partial and complete rebound）is also unidentified.Therefore，the above two significant factors need to be further explored to establish the universal design principle.

In this paper，the impacting dynamics of the droplet on mix-wettability surfaces（hydrophobic substrate decorated with a hydrophilic stripe）are simulated by a reliable numerical method.The role of the hydrophilic stripe in the satellite droplet formation is clarified，and the dynamic and energy transfer mechanisms of the directional rebound of the main droplet are revealed.Then，the effect of the width of the hydrophilic stripe on adhesive force，liquid residual and rebound modes is systematically studied.On this basis，the comprehensive design principle of mix-wettability surfaces for inducing directional rebound is proposed，which provides valuable guidance for the droplet manipulation.

1 Problem Description and Methodology

The sketch of droplet impact is illustrated in Fig.1.The mixed-wettability surface is consisted of two parts：The superhydrophobic substrate（yellow region）and the hydrophilic stripe（red region）.The droplet with diameterD0is initially very close to the surface，and it impacts the surface with a predetermined impact velocityU0.The water and air are adopted as the working fluids.The size of the computational domain is taken to be 10D0×10D0×10D0so as to guarantee that the results are not affected by the boundaries，and the gravity accelerationgacts in the negativezdirection.In addition，the width of the hydrophilic stripe isw.In this paper，the impact velocity is in the range of ［1.25 m/s，2.08 m/s］and the droplet diameter is 2 mm.Since the droplet motion is mainly dominated by the surface tension and inertial forces，the dimensionless parameter ofWeis considered and defined asHere，ρlis the density of the water droplet.

Fig.1 Sketch of the droplet impacting on a mixed wettability surface

The droplet-solid interaction is considered as incompressible and isothermal multiphase flow problem in this paper.The evolution of water-air interface is captured by an improved diffuse interface（DI）method，which has been implemented in the OpenFOAM framework［27］

whereCis the order parameter in the range of ［0，1］anduthe fluid velocity.γis an adaptive parameter used to control the magnitude of the interfacecorrection term，which can improve the interface sharpness.In addition，the correction term is only taken into account in the diffuse layer，and it works normal to the interfaces，i.e.，in the directionMis the mobility andφthe chemical potential.Based on the free energy density model［28-29］，the chemical potential is defined as1)(2C-1)-ασεΔC.Here，εmeasures the interface thickness，andαis a positive constant，which is equal toThe equilibrium function of the one-dimensional interface profileCcan be written as

wherezdis the signed normal distance to the interface.

The flow field is obtained by solving the Navier-Stokes equations，which are written as

whereFσ=-C?φis the surface tension force.ρandμare the density and the dynamic viscosity，respectively，which depend on the order parameterC

where the subscripts “l(fā)” and “g” denote the water and air，respectively.

When the two-phase interface interacts with the solid surface，the unit normal vector of the interfacenshould be modified by the contact angleθd，and it is expressed as

wherenwandtware the unit normal and tangential vectors of the solid surface，respectively.In addition，when the contact line moves on the solid surface，the contact angle is no longer a constant due to the phenomenon of the contact angle hysteresis.In this paper，the variation of the contact angle in the impacting process is predicted by the Cox-Voinov model［32］，which is given in the form as

whereθappis the apparent contact angle；Ca=is the capillary number，hereuclis the contact line velocity；is the capillary length；χdenotes the grid size.Here，θappis obtained with Kistler’ law［33］，which is defined as

whereθadv，θstaandθrecare the advancing，static and receding contact angles，respectively.The Hoffman’s functionfHoffcan be written as

The Cox-Voinov dynamics contact angle model has extensively been used in the research of the droplet impact［26，32，34］.It can accurately predict the variation of the contact angle in the impacting process and obtain satisfactory numerical results.

The governing equations are discretized by the finite volume method.The detailed discretization procedure can be found in our previous work［26］.In addition，the adopted numerical model and grid independence have been well-validated by simulating the problems of droplet impacting on homogeneous and heterogeneous wettability surfaces［25］.Thus，the details are not presented in this paper.

2 Results and Discussion

In this section，the directional rebound of droplet impacting on mix-wettability surfaces is numerically investigated.The attention is focused on the formation of satellite droplets and the rebound mode of the main droplet.And the design principle of mixwettability surface is proposed for better droplet manipulation.It should be noted that the effects of multiple parameters on the droplet directional rebound have been systematically investigated in our previous works［25-26］，including the Weber number（droplet diameter and impact velocity），the arrangement of hydrophilic stripes（position and shape），and the values of contact angles.The obtained results show that the droplet deposits on the surface atWe≤5，and the droplet splashes atWe≥120.The droplet rebound usually occurs when the Weber number is in the range of 10 to 100.Therefore，the Weber numbers of 100 and 120 are selected to explore the droplet breakup in this paper，and a moderate value ofWe=43.3 is adopted to study the droplet rebound with lateral motion under various stripe widths.In addition，all cases are limited in the given contact angles.The advancing，static and receding contact angles on the hydrophobic part are set as 164°，162°，and 160°，respectively，and the corresponding contact angles on the hydrophilic part are 120°，80°，and 40°，respectively.The distanceabetween the impacting point and hydrophilic stripe is 1.5 mm.

2.1 Formation of satellite droplets

The formation of satellite droplets after impact is a well-known phenomenon.In order to explore the dynamic characteristics in this process，the droplet impacting on superhydrophobic and mixed-wettability surfaces are carried out，respectively.The evolutions of the liquid film for the cases ofWe=120（D0=2 mm andU0=2.08 m/s）andw=0.2 mm are plotted in Fig.2.After the droplet reaching its maximum spread（t=2 ms），it starts to retract under the action of superhydrophobic substrate.It can be seen that the liquid film on homogeneous surface moves uniformly and radially towards the impacting point，and no satellite droplets are formed in the retraction stage.On the mixed-wettability surface，the right liquid film can retract normally，but the left liquid film gets pinned by the hydrophilic stripe.This prevents the central liquid film from being fully filled，causing the liquid film to thin and form holes（t=5 ms）.Finally，small satellite droplets are generated with the breaking of the liquid film.

In addition，we can observe that there also are satellite droplets when the droplet rebound from the superhydrophobic surface（t=6.5 ms），as shown in Fig.2（a）.To further reveal the mechanisms of dynamics and energy transfer，the morphology of bouncing droplet is shown in Fig.3（a）.The result of the mixed-wettability surface is also presented for comparison.Here，the impact parameters are set as：We=100（D0=2 mm andU0=1.9 m/s）andw=0.4 mm.The iso-surface is colored by the vertical velocityUz.When the droplet rebounds from the superhydrophobic surface（t=8 ms），it is obvious that the gradient of the vertical velocity is larger than that from the mixed-wettability surface.As a result，the top liquid film is elongated and eventually separated from the main droplet.As shown in Fig.3（b），the liquid film retracts uniformly on the homogeneous surface，so the droplet rebounds vertically and the lateral velocity isUx≈0.However，the liquid film is pinned by the hydrophilic stripe of the mixed-wettability surface，and the rebound kinetic energy of the droplet is partially converted to the lateral direction.Thus，the bouncing droplet has the considerable lateral velocity.

Fig.2 Time evolution of the droplet impacting on superhydrophobic and mixed-wettability surfaces

Fig.3 Morphology of the droplet rebound from superhydrophobic and mixed-wettability surfaces

The formation of satellite droplets is related to two types of instability.One is the Kelvin-Helmholtz instability in the droplet spreading and retraction stages［35］.The liquid film is thin and further forms the hole，which will cause the droplet to break.For example，the formation of satellite droplets in Fig.2（b）occurs at the central liquid film，which is caused by the adhesion of the hydrophilic stripe in the retraction stage.Similarly，the splash may occur when the droplet impacts with a large Weber number，which is induced by the liquid film breaking at the edge of spreading droplet.The other is the Plateau-Rayleigh instability in the droplet rebound stage［36］.The rebound droplet is elongated under the action of the gradient of the vertical velocity，and excessive elongation will break the liquid film and generate the satellite droplets.Especially on the superhydrophobic surface，two or three satellite droplets may appear during the droplet rebounding process.For the mixed-wettability surface，we conclude that it can suppress the Plateau-Rayleigh instability in the droplet rebound stage，but promote the Kelvin-Helmholtz instability in the droplet retraction stage.No obvious effect is observed in the droplet spreading stage.

2.2 Effects of stripe width

In this subsection，the effect of the width of the hydrophilic stripe on the droplet rebound is studied.The value of Weber number is given as 43.3，and the stripe widthwvaries in the range of ［0.1 mm，1.6 mm］.

The evolutions of droplet rebounds from the mixed-wettability surfaces with different stripe widths are shown in Fig.4.In each case，the timetfrom left to right is 2，7，9，12，15，and 18 ms.It can be seen that the rebound form and contact time of the droplet vary greatly under different widths.For the case ofw=0.1 mm，the droplet leaves the surface immediately after retraction.The low adhesion exerted by the narrow stripe is not sufficient to hold the droplet，causing the droplet to rebound completely without the liquid residue.For the cases ofw=0.3 mm and 0.8 mm，a liquid bridge is formed between the rebound droplet and the residual liquid film.In the process of the lateral rebound，the liquid bridge is broken and the droplet leaves the surface.It should be noted that the “separate angle”of the droplet（marked with black lines in Figs.4（b，c）is different，which determines the contact time.With the further increase of the stripe width（w=1.2 mm and 1.6 mm），the droplet leaves the surface immediately after retraction，which is similar to the case ofw=0.1 mm.On the contrary，the high adhesion exerted by the wide stripe breaks the connection between the droplet and the liquid film，and a large amount of liquid remains on the hydrophilic stripe.

Fig.4 Evolutions of droplet rebounds from the mixed-wettability surfaces with different stripe widths

For the above three typical rebound forms of the droplet，the variation of the contact time with the stripe width is shown in Fig.5.The separation of droplet is mainly determined by the inertial force（rebound）and the surface tension force（adhesion）.In the phases of ①and ③，the inertial force and the surface tension force are dominate，respectively.Thus，the droplet can rebound from the surface rapidly，and the contact time is about 8 ms，which is similar to that of the droplet impacting on the superhydrophobic surface（w=0）.In the phase of ②，the contact time increases first and then decreases，which is related to the evolution of liquid bridge.For the small stripe width，the broken liquid bridge will leave the surface together with the rebound droplet.With the increase ofw，the adhesion of the hydrophilic stripe is enhanced，and the broken liquid bridge tends to remain on the surface.Here，the maximal contact time is obtained when the stripe width is 0.8 mm.

Fig.5 Variation of the contact time with the stripe width

2.3 Design principle

For the droplet directional rebound，three factors need to be considered in the design of the mixed-wettability surface，including the formation of satellite droplets，the lateral velocity of the rebound droplet and the liquid residue on the hydrophilic stripe.Among of them，the situations of satellite droplets and the liquid residue can be easily obtained from the above results.When the Weber number is in the range of ［10，100］，the mixed-wettability surface can suppress the formation of the satellite droplets in the rebound stage.However，it may induce the Kelvin-Helmholtz instability at the higher Weber number，thus promoting the formation of satellite droplets.In addition，the liquid residue is usually unfavorable in engineering applications，such as the anti-icing/fogging and self-cleaning.Fig.6 shows the percentage of residual liquid in the total mass of droplets.It is clear that the liquid residue increases with the stripe width.

Fig.6 Variation of the liquid residue with the stripe width

This subsection focuses on the effect of the stripe width on the performance of the droplet directional rebound.The kinetic energyEkof the droplet rebound from the mixed-wettability surface with different stripe widths are shown in Fig.7.The impacting droplet undergoes the spreading，retraction and rebound stages，in which part of the energy dissipates due to the surface friction and viscosity.It can be seen that the variation trend of rebounding kinetic energy is opposite to that of contact time（Fig.5）.This is mainly because the longer the droplet is pulled by the hydrophilic stripe，the more energy it loses during its movement.Therefore，the minimum kinetic energy is obtained atw=0.8 mm，which corresponds to the maximum contact time.

Fig.7 Variation of the rebounding kinetic energy Ek with the stripe width

Furthermore，the effect of the stripe width on the rebounding heightHand the lateral velocityvxof the droplet is shown in Fig.8.The lateral velocity can quantify the performance of droplet directional transportation，and the variation trend ofvxis similar to that of contact time（Fig.5）.This means that the energy continually transfer laterally before the droplet leaves the surface.The high lateral velocity can be obtained when the stripe width is in the range of ［0.3 mm，1.2 mm］.With the increase of the contact time，more energy is lost during the impact process and the more energy is transferred to the lateral direction.Thus，it is easily understood that the droplet rebounding heightHis lower at the medium stripe width.

Fig.8 Variations of the rebounding height H and the lateral velocity vx with the stripe width

In summary，two principles can be obtained to provide reference for the design of mixed-wettability surfaces.Firstly，the mixed-wettability surface is more suitable for the droplet directional rebound with Weber number between 5 and 100.In this range，the existence of the hydrophilic stripe can effectively suppress the formation of satellite droplets.Secondly，w=0.3 is a better choice for the purpose of obtaining the larger rebounding kinetic energy，the higher lateral velocity and the lesser liquid residue.

3 Conclusions

The dynamics of a droplet impacting on the mixed-wettability surface，which can induce the directional rebound，are systematically investigated.The numerical results show that the mixed-wettability surface can effectively suppress the formation of the satellite droplets in the rebound stage，but promote the Kelvin-Helmholtz instability in the retraction stage.Moreover，three typical droplet rebound forms are observed with different stripe widths.The larger lateral velocity and rebounding kinetic energy of the droplet，and the less liquid residue on the hydrophilic stripe can be obtained whenw=0.3 mm.These results may be useful for the design of the mixed-wettability surface，thus further promoting its application in engineering.