Research Article  Open Access
Simulation of Effective Slip and Drag in PressureDriven Flow on Superhydrophobic Surfaces
Abstract
The flow on superhydrophobic surfaces was investigated using finite element modeling (FEM). Surfaces with different textures like grooves, square pillars, and cylinders immersed in liquid forming Cassie state were modeled. Nonslip boundary condition was assumed at solidliquid interface while slip boundary condition was supposed at gasliquid interface. It was found that the flow rate can be affected by the shape of the texture, the fraction of the gasliquid area, the height of the channel, and the driving pressure gradient. By extracting the effective boundary slip from the flow rate based on a model, it was found that the shape of the textures and the fraction of the gasliquid area affect the effective slip significantly while the height of the channel and the driving pressure gradient have no obvious effect on effective slip.
1. Introduction
Development in microfluidics has stimulated interest in manipulating flows in micro/nanoscale channels [1, 2]. Drag reduction at the solidliquid interface is an important challenge in the research and application of microfluidics. Boundary slip can reduce the drag at the interface [3–7]. Boundary slip is quantified by slip length , which is the distance beyond the solidliquid interface at which the velocity profile extrapolates to zero [7–9]. Researches have shown that flat, smooth, and chemically homogenous hydrophilic surfaces exhibit no or less slip [10, 11]. However, if the surface is hydrophobic, there can be slip with the magnitude up to tens of nanometers [6, 12–15]. In the meanwhile, superhydrophobic surfaces, on which the contact angle exceeds 150° and the contact angle hysteresis is less than 5°, have been proved to be able to significantly amplify the magnitude of slip. These surfaces have composite surface structures, which can trap gas at the grooves or pits [16, 17] and form a Cassie state. The slip length on superhydrophobic surface or other complex heterogeneous surfaces is represented by effective slip . The is obtained by averaging the flow over the characteristic length scale of the configuration (e.g., a channel) [4], which is usually defined as where means the average value in the supposed plane on surfaces [18]. Experiments of flow on superhydrophobic surfaces have indicated that the can reach to the order of several or tens of microns [19–24].
As most of the “super” properties of superhydrophobic surfaces are due to the gas entrapment, the Cassie state is preferred rather than the Wenzel state. Numerous research efforts have been done during the last few decades in order to stabilize the Cassie state [21, 25, 26]. Since that it is hard to measure the actual fraction of solid in contact with liquid for arbitrary rough or fractional surface, surfaces with well controlled periodic textures have been widely used [27].
In recent years, the researches focused on quantitative understanding of the effective slippage on the superhydrophobic surfaces. Attention has been paid to surfaces with directional patterns, such as the arrays of parallel superhydrophobic grooves, square pillars, and cylinders that generate anisotropic effective slip in the Cassie regime experimentally and theoretically [17, 28–31]. This paper extended the work to investigate the effective slip as well as drag in pressuredriven flow on superhydrophobic surfaces decorated with arrays of grooves, square pillars, and cylinders, respectively, through simulation. The flow has been studied with the help of FEM technique and the results as well as discussion are given in the following sections.
2. Models and Calculation
2.1. Models
Surfaces in Cassie state with textures in micro scale usually lead to the socalled superhydrophobicity. Three different textures, that is, grooves, square pillars, and cylinders, were chosen to decorate the surface.
In most of the microfluidic applications, the flow is driven through pressure. Therefore, pressuredriven fluid that flows through a straight channel was modeled, in which the upper plate is an ideal flat one, while the lower one is the superhydrophobic surfaces with special textures (Figure 1). Fluid flows through the channel, fully wetting the upper plate but partially wetting the lower plate, which form the Cassie state. To simplify the model, the meniscus curvature at the gasliquid interface was ignored, as considered in most previous publications [6, 13, 32–34]. The length and width of the channel, the period length , and the height of the textures are in constant number as shown in Table 1. The size of the textures , the width of the interval (), and the height of the channel are main variables. Using these parameters, the fraction of gasliquid area can be calculated.
 
, , and in “range (, , )” represent the initial value, step length, and final value. 
Figure 2 shows that the fluid flows from one side to the opposite side. According to the flow direction, there can be five typical flow regimes over these textures: (a) longitudinal direction of grooves, (b) transverse direction of the grooves, (c) parallel (or perpendicular) to the array of the square pillars, (d) 45° to the array of the square pillars, and (e) parallel (or perpendicular) to the array of the cylinders.
(a)
(b)
(c)
(d)
(e)
2.2. Fluid Properties
Incompressible water was chosen as the fluid, of which the flow follows the NavierStokes equations and the mass continuity equation:where is the density of the fluid, is the of velocity vector, is the driven pressure, is the momentum vector, is the dynamic viscosity of the fluid, and is the external force. The flow field remains laminar due to microscale and the low driven pressure. The external reference pressure is 1 bar ( Pa), and room temperature is 293.15 K, the density of water g/cm^{3}, and dynamic viscosity Pas [32].
2.3. Boundary Condition
In the socalled “gas cushion model,” the interfacial region is described as a lubricating “gas film” of thickness and viscosity , while the bulk value is . The slip can be calculated [33]When the thickness of gas film ranges between 0.1 and 10 µm, the slip at the gasliquid interface is about 5–500 µm (dynamics viscosity of gas Pas at 293.15 K [33]). This model is justified for a continuous “gas film” on a homogeneous surface [32] or surface with shallow textures [34]. Note that it has been verified that the slip length at the smooth solidliquid interface is no more than tens of nanometers. Therefore, the slip at the upper flat plate (Figure 1) can be ignored. Accordingly, this simulation has set the boundary conditions to nonslip at the upper plate as well as at the top of the textures, as described by (4). It has been assumed that the slip only occurs at the space between the textures, expressed by (5). The slip condition depends on the contact state at the superhydrophobic surface. It is confined by the geometry of gasliquid interface rather than the depth of the textures. Equations (4) and (5) take the formwhere is the normal vector and . These boundary conditions assume that the gas layer has no friction force to the liquid flow but confine the fluid velocity direction parallel to the wall.
Fluids at the inlet and the outlet are all of normal flow without backflow, followingwhere is the pressure at the outside of the fluid field near inlet or outlet, at the inlet, and at the outlet; is the tangential vector to the wall. In many previous publications about pressuredriven microchannel flow [15, 35, 36], the pressure gradient is assumed . With reference to this, in this model, is set a variable parameter ranging from 6 to 600 Pa with the step length of 6 Pa, while ; in which case, ranges from to with the step length . Moreover, the lateral sides of the fluid flow field are defined by symmetric plane.
The flow was simulated in Comsol Multiphysics. The effect of the shape of textures, fraction of gas, height of channel, and driving pressure gradient on effective slip as well as drag was studied. Five different cases, as shown in Figure 1, were simulated. The results are given in the corresponding sections. Moreover, it has set that the flow driven by pressure in microchannel with superhydrophobic surface still follows NavierStokes equation, but with an effective slip on the superhydrophobic surface, described asThe flow rate in the channel is calculated byCombining (7) and (8), it can get relationship between effective slip and flow rate :Therefore, the effective slip can be obtained through above equation as the flow rate is obtained through the simulation.
3. Results and Discussions
3.1. Effect of the Surface Textures
The interface at flat surface or the top of the textures refers to a solidliquid contact while the interface of the gap and the liquid refers to a gasliquid contact. A nonslip condition was supposed at the solidliquid interface while a slip condition was assumed at gasliquid interface. As an example, the distributions of velocity and pressure of the fluid at the interface of superhydrophobic surfaces with , µm, and are shown in Figure 3. The velocity at the gasliquid interface is higher than that at the solidliquid interface; however, the shape of the textures still affects the distribution and magnitude of velocity of the fluid. If the textures have small angle of inclination to the direction of flow, there will be lower drag; thus the flow can be enhanced. It can be noticed that the liquid flows smoothly on the longitudinally grooved surface. However, the flow on surfaces textured with square pillars or cylinders changes with the angle of inclination. The drag will be larger with larger angle of inclination.
(a)
(b)
(c)
(d)
(e)
Figure 4 shows the flow distribution in the direction perpendicular to the flow. The velocity near the walls is lower while it is higher at the center of the channel. Noticeably, at the contact line of superhydrophobic surface and liquid, the velocity is zero at the solidliquid interface while it is large at the gasliquid interface where it shows a finite slip.
3.2. Effect of the GasLiquid Fraction
Different sized textures were modeled on the surfaces shown in Table 1. The effect of the gasliquid fraction on the flow rate and effective slip are shown in Figures 5 and 6, respectively. is used to show the change of the flow rate where is the flow rate in microchannels with textured surfaces and is the flow rate in microchannels with smooth surfaces in the same condition. Due to the constraint of the shape, the fraction of gasliquid area of surfaces decorated with square pillars with the angle of inclination 45° cannot be lower than 50%, while that surface textured with cylinders should be larger than about 21%. Also, surfaces decorated with certain patterns, on which the fraction of gasliquid area can be up to 99%, have been realized and applied in experiment [20].
As shown in Figures 5 and 6, the lines of flow rate and effective slip in the same regime show a similar trend and increase exponentially with the increase of . When is lower than about 84%, the flow rate and effective slip on the surfaces with longitudinal grooves are larger than the other surfaces. When increases over 90%, the enhancement of the flow in the other surfaces is more significant than the channel with transverse grooves decorated surface. The obstruction of the flow by the transverse grooves may be the main reason of the weak enhancement of flow rate. Moreover, the magnitude of the effective slip ranges in nanometers when the fraction is lower than 50%, and it reaches to about 3–57 µm when the fraction is larger than 90%, which is in agreement with the previous experimental as well as theoretical results [19–24]. Maali et al. [19] got 291 nm effective slip length on cylinders decorated surface with a = 0.63 measured by AFM, while it shows about 700 nm effective slip in this simulation in the same state. Choi et al. [20] studied the effect of gas fraction of textured hydrophobic surfaces on slip using rheometer system, and they obtained about 18 µm, 25 µm, 50 µm, 80 µm, and 100 µm slip length at the gas fraction of 50%, 85%, 95%, 98%, and 99% on surface textured with cylinders, respectively. At the same fraction of gasliquid area, this simulation shows 350 nm, 3 µm, 8 µm, 17 µm, and 29 µm slip length, respectively. Although the numbers of effective slips have deviations, the trends of the boundary slip with the fraction are in agreement. The inaccuracy of experiments [12, 37], the existence of meniscus at the gasliquid interface [38–40], and the simplifications in this simulation might be the reason of the deviations on the numbers.
3.3. Effect of the Height of the Channel and Driving Pressure
The effect of the height of the channel on the flow rate was also studied. The gasliquid fraction was fixed at 0.5 while the driving pressure gradient . ranged from 1 to 100 µm with the step length 1 µm. The results are shown in Figures 7 and 8 and Table 2.

As shown in Figure 7, the flow rate in channels characterized with textures can reach twice the flow rate on smooth surface when is as small as 1 µm, and the enhancement decreases rapidly with the increase of and the flow rates on all textured surfaces are nearly the same as the flow rate on smooth surface when reaches 100 µm. It strongly suggests that textures on surfaces can remarkably boost the flow in channels with height in small scales. Among the five cases, surface decorated with longitudinal grooves has the largest flow rate, which shows a better lubrication.
In Figure 8, the effective slip shows a clear distinction: a sharp increase when increases from 1 µm to about 5 µm; then the effective slip remains approximately constant as increase to 100 µm but with visible fluctuations. It is noticeable that when is small, the upper wall with nonslip will have a strong effect on the flow near the lower superhydrophobic surfaces. With the increase of the height, this effect will decay until it is close to zero. When is larger than 20 µm, fluctuates around a certain number (details shown in Table 2). The meshing in this FEM simulation should be the reason, with an increasing ; the simulation needs to reconstruct the mesh in every step, which leads to the deviations. Therefore, it can be believed that the height of the channel has no obvious effect on the effective slip at the interface of superhydrophobic surface and liquid. The results show that channel characterized with longitudinal grooves has the largest effective slip while channel with square pillars has the smallest effective slip. This can be explained by the obstruction effect of the textures as mentioned in effect of surface textures.
Fixing the fraction of gasliquid area at 0.5 and the height of the channel at 20 µm, the effect of driving pressure gradient on the flow rate and the effective slip were studied based on the models shown in Figure 2. ranges from to with the step length . The results are shown in Figures 9 and 10.
Figure 9 shows that the superhydrophobic surfaces significantly promote the flow, and the flow rate is larger in the channel characterized with longitudinal grooves which suggests lower drag and larger effective slip at the same driving pressure gradient, as shown in Figure 10. It is obvious that the effective slip on all kinds of surfaces descends slightly with the increase of the driving pressure gradient. The data essentially reveal the degree of lubrication of the textures. Order of the intercepts of the five cases also can be explained by the obstruction effect of their shapes.
4. Conclusion
The flow rate in the channel with different textures, fraction of gasliquid area, height of the channel, and driving pressure gradient are analyzed. Unlike most previous research, the independent effect of the four factors was studied, respectively. Textures with small angle of inclination to the direction of flow can reduce the drag and enhance the flow. The flow rate will increase more sharply with larger fraction of gasliquid which is in agreement with previous studies [19–24].
The results also show that the effective slip is affected by the surface textures and the gasliquid fraction. The height of the channel and the change in the driving pressure gradient have little or no effect on the effective slip. The flow analysis showed that the flow on the surface with longitudinal grooves and higher fraction of gasliquid exhibit less drag.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This research was supported by National Natural Science Foundation of China (Grants nos. 51475118 and 51505108).
References
 H. A. Stone, A. D. Stroock, and A. Ajdari, “Engineering flows in small devices: microfluidics toward a labonachip,” Annual Review of Fluid Mechanics, vol. 36, pp. 381–411, 2004. View at: Publisher Site  Google Scholar
 T. M. Squires and S. R. Quake, “Microfluidics: fluid physics at the nanoliter scale,” Reviews of Modern Physics, vol. 77, no. 3, pp. 977–1026, 2005. View at: Publisher Site  Google Scholar
 D. Jing and B. Bhushan, “The coupling of surface charge and boundary slip at the solid–liquid interface and their combined effect on fluid drag: a review,” Journal of Colloid and Interface Science, vol. 454, pp. 152–179, 2015. View at: Publisher Site  Google Scholar
 O. I. Vinogradova and A. V. Belyaev, “Wetting, roughness and flow boundary conditions,” Journal of Physics: Condensed Matter, vol. 23, no. 18, Article ID 184104, 16 pages, 2011. View at: Publisher Site  Google Scholar
 G. McHale, M. I. Newton, and N. J. Shirtcliffe, “Immersed superhydrophobic surfaces: gas exchange, slip and drag reduction properties,” Soft Matter, vol. 6, no. 4, pp. 714–719, 2010. View at: Publisher Site  Google Scholar
 O. I. Vinogradova and G. E. Yakubov, “Dynamic effects on force measurements. II. Lubrication and the atomic force microscope,” Langmuir, vol. 19, no. 4, pp. 1227–1234, 2003. View at: Publisher Site  Google Scholar
 O. I. Vinogradova, “Slippage of water over hydrophobic surfaces,” International Journal of Mineral Processing, vol. 56, no. 1, pp. 31–60, 1999. View at: Publisher Site  Google Scholar
 L. Bocquet and J.L. Barrat, “Flow boundary conditions from nano to microscales,” Soft Matter, vol. 3, no. 6, pp. 685–693, 2007. View at: Publisher Site  Google Scholar
 E. Lauga, M. P. Brenner, and H. A. Stone, “Microfluidics: the noslip boundary condition,” in Handbook of Experimental Fluid Dynamics, vol. 19, pp. 1219–1240, Springer, 2007. View at: Google Scholar
 C. Sendner, D. Horinek, L. Bocquet, and R. R. Netz, “Interfacial water at hydrophobic and hydrophilic surfaces: slip, viscosity, and diffusion,” Langmuir, vol. 25, no. 18, pp. 10768–10781, 2009. View at: Publisher Site  Google Scholar
 P. A. Thompson and S. M. Troian, “A general boundary condition for liquid flow at solid surfaces,” Nature, vol. 389, no. 6649, pp. 360–362, 1997. View at: Publisher Site  Google Scholar
 C. CottinBizonne, B. Cross, A. Steinberger, and E. Charlaix, “Boundary slip on smooth hydrophobic surfaces: intrinsic effects and possible artifacts,” Physical Review Letters, vol. 94, no. 5, Article ID 056102, 4 pages, 2005. View at: Publisher Site  Google Scholar
 C. Ybert, C. Barentin, C. CottinBizonne, P. Joseph, and L. Bocquet, “Achieving large slip with superhydrophobic surfaces, scaling laws for generic geometries,” Physics of Fluids, vol. 19, no. 12, Article ID 123601, 10 pages, 2007. View at: Publisher Site  Google Scholar
 O. I. Vinogradova, K. Koynov, A. Best, and F. Feuillebois, “Direct measurements of hydrophobic slippage using doublefocus fluorescence crosscorrelation,” Physical Review Letters, vol. 102, no. 11, Article ID 118302, 5 pages, 2009. View at: Publisher Site  Google Scholar
 Y. Pan, B. Bhushan, and X. Zhao, “The study of surface wetting, nanobubbles and boundary slip with an applied voltage: a review,” Beilstein Journal of Nanotechnology, vol. 5, no. 1, pp. 1042–1065, 2014. View at: Publisher Site  Google Scholar
 O. I. Vinogradova, N. F. Bunkin, N. V. Churaev, O. A. Kiseleva, A. V. Lobeyev, and B. W. Ninham, “Submicrocavity structure of water between hydrophobic and hydrophilic walls as revealed by optical cavitation,” Journal of Colloid and Interface Science, vol. 173, no. 2, pp. 443–447, 1995. View at: Publisher Site  Google Scholar
 C. CottinBizonne, J.L. Barrat, L. Bocquet, and E. Charlaix, “Lowfriction flows of liquid at nanopatterned interfaces,” Nature Materials, vol. 2, no. 4, pp. 237–240, 2003. View at: Google Scholar
 S. Schmieschek, A. V. Belyaev, J. Harting, and O. I. Vinogradova, “Tensorial slip of superhydrophobic channels,” Physical Review E, vol. 85, no. 1, Article ID 016324, 11 pages, 2012. View at: Publisher Site  Google Scholar
 A. Maali, Y. Pan, B. Bhushan, and E. Charlaix, “Hydrodynamic dragforce measurement and slip length on microstructured surfaces,” Physical Review E, vol. 85, no. 6, Article ID 066310, 2012. View at: Publisher Site  Google Scholar
 C.H. Choi, C.J. Kim, and C. Lee, “Structured surfaces for a giant liquid slip,” Physical Review Letters, vol. 101, no. 6, Article ID 064501, 2008. View at: Publisher Site  Google Scholar
 C.H. Choi, U. Ulmanella, J. Kim, C.M. Ho, and C.J. Kim, “Effective slip and friction reduction in nanograted superhydrophobic microchannels,” Physics of Fluids, vol. 18, Article ID 087105, 8 pages, 2006. View at: Publisher Site  Google Scholar
 P. Joseph, C. CottinBizonne, J.M. Benoît et al., “Slippage of water past superhydrophobic carbon nanotube forests in microchannels,” Physical Review Letters, vol. 97, no. 15, Article ID 156104, 4 pages, 2006. View at: Publisher Site  Google Scholar
 J. Ou and J. P. Rothstein, “Direct velocity measurements of the flow past dragreducing ultrahydrophobic surfaces,” Physics of Fluids, vol. 17, no. 10, Article ID 103606, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. R. Philip, “Flows satisfying mixed noslip and noshear conditions,” Zeitschrift für Angewandte Mathematik und Physik, vol. 23, no. 3, pp. 353–372, 1972. View at: Google Scholar
 D. Quéré, “Wetting and roughness,” Annual Review of Materials Research, vol. 38, pp. 71–99, 2008. View at: Publisher Site  Google Scholar
 J. Bico, C. Marzolin, and D. Quéré, “Pearl drops,” Europhysics Letters, vol. 47, no. 2, pp. 220–226, 1999. View at: Publisher Site  Google Scholar
 A. L. Dubov, J. Teisseire, and E. Barthel, “Elastic instability and contact angles on hydrophobic surfaces with periodic textures,” EPL, vol. 97, no. 2, Article ID 26003, 2012. View at: Publisher Site  Google Scholar
 C. Lee, C. H. Choi, and C. J. Kim, “Structured surfaces for a giant liquid slip,” Physical Review Letters, vol. 101, no. 6, Article ID 064501, 4 pages, 2008. View at: Publisher Site  Google Scholar
 M. Z. Bazant and O. I. Vinogradova, “Tensorial hydrodynamic slip,” Journal of Fluid Mechanics, vol. 613, pp. 125–134, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 F. Feuillebois, M. Z. Bazant, and O. I. Vinogradova, “Effective slip over superhydrophobic surfaces in thin channels,” Physical Review Letters, vol. 102, Article ID 026001, 4 pages, 2009. View at: Google Scholar
 J. P. Rothstein, “Slip on superhydrophobic surfaces,” Annual Review of Fluid Mechanics, vol. 42, pp. 89–109, 2010. View at: Publisher Site  Google Scholar
 W. M. Haynes, Handbook of Chemistry and Physics, Academic Press, New York, NY, USA, 91st edition, 2010.
 O. I. Vinogradova, “Drainage of a thin liquid film confined between hydrophobic surfaces,” Langmuir, vol. 11, no. 6, pp. 2213–2220, 1995. View at: Publisher Site  Google Scholar
 T. V. Nizkaya, E. S. Asmolov, and O. I. Vinogradova, “Flow in channels with superhydrophobic trapezoidal textures,” Soft Matter, vol. 9, no. 48, pp. 11671–11679, 2013. View at: Publisher Site  Google Scholar
 D. Jing and B. Bhushan, “Electroviscous effect on fluid drag in a microchannel with large zeta potential,” Beilstein Journal of Nanotechnology, vol. 6, no. 1, pp. 2207–2216, 2015. View at: Publisher Site  Google Scholar
 H. Ban, B. C. Lin, and Z. R. Song, “Effect of electrical double layer on electric conductivity and pressure drop in a pressuredriven microchannel flow,” Biomicrofluidics, vol. 4, Article ID 014104, 14 pages, 2010. View at: Publisher Site  Google Scholar
 K. Ahmad, X. Zhao, Y. Pan, W. Wang, and Y. Huang, “Atomic force microscopy measurement of slip on smooth hydrophobic surfaces and possible artifacts,” The Journal of Physical Chemistry C, vol. 119, no. 22, pp. 12531–12537, 2015. View at: Publisher Site  Google Scholar
 A. Gaddam, M. Garg, A. Agrawal, and S. S. Joshi, “Modeling of liquidgas meniscus for textured surfaces: effects of curvature and local slip length,” Journal of Micromechanics and Microengineering, vol. 25, no. 12, Article ID 125002, 13 pages, 2015. View at: Publisher Site  Google Scholar
 A. M. J. Davis and E. Lauga, “Geometric transition in friction for flow over a bubble mattress,” Physics of Fluids, vol. 21, no. 1, Article ID 011701, 2009. View at: Publisher Site  Google Scholar
 A. Steinberger, C. CottinBizonne, P. Kleimann, and E. Charlaix, “High friction on a bubble mattress,” Nature Materials, vol. 6, no. 9, pp. 665–668, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2016 Yuanding Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.