# Effective slippage on superhydrophobic trapezoidal grooves

## Abstract

We study the effective slippage on superhydrophobic grooves with trapezoidal cross-sections of various geometries (including the limiting cases of triangles and rectangular stripes), by using two complementary approaches. First, dissipative particle dynamics (DPD) simulations of a flow past such surfaces have been performed to validate an expression [E. S. Asmolov and O. I. Vinogradova, J. Fluid Mech. 706, 108 (2012)] that relates the eigenvalues of the effective slip-length tensor for one-dimensional textures. Second, we propose theoretical estimates for the effective slip length and calculate it numerically by solving the Stokes equation based on a collocation method. The comparison between the two approaches shows that they are in excellent agreement. Our results demonstrate that the effective slippage depends strongly on the area-averaged slip, the amplitude of the roughness, and on the fraction of solid in contact with the liquid. To interpret these results, we analyze flow singularities near slipping heterogeneities, and demonstrate that they inhibit the effective slip and enhance the anisotropy of the flow. Finally, we propose some guidelines to design optimal one-dimensional superhydrophobic surfaces, motivated by potential applications in microfluidics.

## I Introduction

The design and fabrication of micro- and nanotextured surfaces has received much attention during the past decade. In case of a hydrophobic texture a modified surface profile can lead to a very large contact angle, which induces exceptional wetting properties Quere (2005). A remarkable mobility of liquids on such superhydrophobic surfaces is observed in the Cassie state, where the textures are filled with gas, which renders them “self-cleaning” and causes droplets to roll (rather than slide) under gravity Richard and Quere (1999) and rebound (rather than spread) upon impact Richard and Quere (2000); Tsai et al. (2010). Furthermore, patterned superhydrophobic materials are important in context of fluid dynamics and their superlubricating properties Vinogradova and Dubov (2012); McHale et al. (2010). In particular, superhydrophobic heterogeneous surfaces in the Cassie state exhibit very low friction, and this drag reduction is associated with the large slippage of liquids Bocquet and Barrat (2007); Rothstein (2010); Vinogradova and Belyaev (2011).

It is very difficult to quantify the flow past heterogeneous surfaces. However, analytical results can often be obtained using an effective slip boundary condition, , at the imaginary smooth homogeneous, but generally anisotropic surface Vinogradova and Belyaev (2011); Kamrin et al. (2010). For anisotropic textures, the effective slip generally depends on the direction of the flow and is a tensor, represented by a symmetric, positive definite matrix Bazant and Vinogradova (2008)

(1) |

diagonalized by a rotation

(2) |

Equation (1) allows us to calculate an effective slip in any direction given by an angle , provided the two eigenvalues of the slip-length tensor, () and (), are known. The concept of an effective slip length tensor is general and can be applied for arbitrary channel thickness Schmieschek et al. (2012). It is a global characteristic of a channel Vinogradova and Belyaev (2011), and the eigenvalues usually depend not only on the parameters of the heterogeneous surfaces, but also on the channel thickness. However, for a thick channel (compared to a texture period, ) they become a characteristic of a heterogeneous interface solely.

An important type of such surfaces are highly anisotropic textures, where a (scalar) partial local slip varies in only one direction. Such surfaces can be created experimentally by making textured surfaces with one dimensional surface profiles, and in the Cassie state, the local slip length will always reflect the relief of the texture. This is the consequence of the “gas cushion model”, which takes into account that the dissipation at the gas/liquid interface is dominated by the shearing of a continuous gas layer Vinogradova (1995); Andrienko et al. (2003)

(3) |

where and are dynamic viscosities of a gas and a liquid, and is the thickness of the gas layer. Eq. (3) represents an upper bound for a local slip at the gas area, which is attained in the limit of a small fraction of the solid phase. At a relatively large solid fraction Ybert et al. (2007) or if the flux in the “pockets” of superhydrophobic surfaces is equal to zero Maynes et al. (2007); Nizkaya et al. (2013), Eq. (3) could overestimate the local slip. However, even in all situations the local slip length profile, , necessarily follows the relief of the texture provided it is shallow enough, Nizkaya et al. (2013). Since in case of water Eq. (3) leads to , it is expected that the “gas cushion model” could be safely applied up to or so.

The flow properties on highly anisotropic surfaces can be highly peculiar. In particular, off-diagonal effects may be obtained where a driving force (such as pressure gradient, shear rate, or electric field) in one direction induces a flow in a different direction, with a measurable and useful perpendicular component.This could be exploited for implementing transverse pumps, mixers, flow detectors, and more Stroock et al. (2002); Ajdari (2001). The main task is to find the connection between the eigenvalues of the slip-length tensor and the parameters of such one-dimensional surface textures. Up to now analytical results for the effective slip length have only been obtained for the simplest geometry of rectangular grooves Lauga and Stone (2003); Belyaev and Vinogradova (2010a); Ng et al. (2010); Zhou et al. (2012); Schmieschek et al. (2012); Asmolov et al. (2013a). We are unaware of any prior work that has obtained analytical results for other types of one-dimensional textures.

However, some recently derived relations for the case of an anisotropic isolated surface (corresponding to the thick channel limit) allow one to simplify the analysis of more complex geometries, and to rationalize the limiting behavior of the effective slip length in certain regimes. We mention first that the transverse component of the slip-length tensor was predicted to be exactly half the longitudinal component that one would obtain if the local slip were multiplied with a factor two, Asmolov and Vinogradova (2012)

(4) |

This relation can greatly simplify the analysis since it indicates that the flow along any direction of the one-dimensional surface can be easily determined, once the longitudinal component of the effective slip tensor is found. Equation (4) has been recently verified for a texture with sinusoidal superhydrophobic grooves by using a Lattice Boltzmann (LB) method Asmolov et al. (2013b), and for weakly slipping stripes by using dissipative particle dynamics (DPD) simulations Asmolov et al. (2013a).

Another useful relation addresses the particular case where the local slip is large, . In that case, the effective slip length tensor becomes isotropic and is given by Asmolov et al. (2013b); Hendy and Lund (2007)

(5) |

where is the mean inverse local slip. Finally, for weakly slipping patterns with , the flow again becomes isotropic, and its value is equal to the surface average of the local slip length Belyaev and Vinogradova (2010a); Kamrin et al. (2010):

(6) |

Note that Eq. (6) is expected to be accurate for a continuous local slip profile only. If the local slip length exhibits step-like jumps at the edges of heterogeneities, the expansions leading to Eq. (6) become questionable. Higher order contributions to the effective slip become comparable to the first-order term given by Eq. (6) and can no longer be ignored Asmolov et al. (2013a).

Computer simulations and numerical modeling can shed light on flow phenomena past one-dimensional surfaces. In efforts to better understand the connection between the parameters of the texture and the effective slip lengths, several groups have performed continuum and Molecular Dynamics simulations of flow past anisotropic surfaces. Most of these studies focused on calculating eigenvalues of the effective slip-length tensor for superhydrophobic stripes Cottin-Bizonne et al. (2004); Priezjev et al. (2005); Priezjev (2011). To study various aspects of a flow past striped walls, recent works employed Dissipative particle dynamics (DPD) Zhou et al. (2012); Asmolov et al. (2013a); Tretyakov and Müller (2013) and Lattice Boltzmann (LB) Schmieschek et al. (2012) methods. In the case of stripes with a piecewise constant local slip profile (i.e., the profile features steplike jumps at the boundaries of the stripes), the singularity in the slip profile induces singularities both in the pressure and velocity gradient, which generates an additional mechanism for dissipation even at small local slip Asmolov and Vinogradova (2012); Asmolov et al. (2013a). Anisotropic one-dimensional texture with continuously varying patterns can produce a very large effective slip. This has been confirmed in recent LB simulations studies of sinusoidal textures Asmolov et al. (2013b). One-dimensional continuous surfaces may also exhibit singularities (for example, the first derivative may be discontinuous). However, we are unaware of any previous work that has addressed the question of effective slip lengths past such textures.

In this paper we study the friction properties on trapezoidal superhydrophobic grooves [see Fig. 1(a)]. This texture is more general than rectangular grooves since it provides an additional geometrical structure parameter, i.e., the trapezoid base width. Also, the trapezoid texture can be extended to its extremes, such as a triangular patterns, where the liquid/solid contact area vanishes, and rectangular reliefs. Superhydrophobic trapezoidal textures were shown to give large contact angles, and to provide stable Cassie states Li et al. (2010). An obvious advantage of the trapezoidal surface relief is that the manufactured texture becomes mechanically more stable against bending compared to rectangular grooves. Furthermore, it can be naturally formed from dilute rectangular grooves as a result of elasto-capillarity Tanaka et al. (1993). This type of surfaces has already been intensively used in slippage experiments Choi et al. (2006). However, the impact of trapezoidal texture reliefs on transport and hydrodynamic slippage remains largely unknown.

Our study is based on the use of DPD simulations, which we have put forward in our previous papers Zhou et al. (2012); Asmolov et al. (2013a), on the numerical solution of the Stokes equations (a collocation method), and on theoretical analysis. In particular, we validate Eq. (4) for trapezoidal textures with various parameters, we apply Eq. (5) to evaluate an effective slip for strongly slipping patterns, we discuss the singularities in velocities caused by inflection points in the local slip profiles, which in particular lead to deviations from Eq. (6), and we finally suggest the optimal design of trapezoidal patterns.

## Ii Model and Governing Equations

We consider creeping flow along a plane anisotropic superhydrophobic wall, using a Cartesian coordinate system (Fig. 1). As in most previous publications Bocquet and Barrat (2007); Belyaev and Vinogradova (2010b); Priezjev et al. (2005); Lauga and Stone (2003); Ybert et al. (2007); Asmolov et al. (2011); Feuillebois et al. (2009); Schmieschek et al. (2012), we model the superhydrophobic plate as a flat heterogeneous interface. In such a description, we neglect an additional mechanism for a dissipation connected with the meniscus curvature Sbragaglia and Prosperetti (2007); Hyväluoma and Harting (2008); Ybert et al. (2007). Note however that such an ideal situation has been achieved in many recent experiments Steinberger et al. (2007); Karatay et al. (2013); Haase et al. (2013). The origin of coordinates is placed at the flat interface, characterized by a slip length , spatially varying in one direction, and the texture varies over a period .

The local slip length is taken to be a piecewise linear function. It includes two regions with a constant slip: one where the slip length assumes its minimum value, , which has area fraction of , and another with area fraction where the slip length is maximal. These regions are connected by two transition regions with equal surface fractions , where the local slip length varies linearly from to . If we rewrite the coordinate in the form , where is an integer number and , the local slip-length profile can be expressed in terms of :

(7) |

For , we get a triangular profile

(8) |

Note that regular trapezoidal and triangular profiles are continuous functions (albeit with a discontinuous first derivative). However, in the limit , the profile reduces to the discontinuous rectangular slip length profile (alternating stripe profile) studied in earlier work

(9) |

Our results apply to a single surface in a thick channel, where effective hydrodynamic slip is determined by flow at the scale of roughness, so that the velocity profile sufficiently far above the surface may be considered as a linear shear flow Bocquet and Barrat (2007); Vinogradova and Belyaev (2011). However it does not apply to a thin Feuillebois et al. (2009) or an arbitrary Schmieschek et al. (2012) channel situation, where the effective slip scales with the channel width. Note that the flow in thin channels with superhydrophobic trapezoidal textures has been considered in recent work Nizkaya et al. (2013).

Dimensionless variables are defined by using as a reference length scale, the inverse shear rate sufficiently far from the surface as a reference time scale, and setting the fluid kinematic viscosity, , to unity (this sets the scale of mass). At small Reynolds number, the dimensionless velocity and pressure fields, and , satisfy Stokes equations

(10) | |||

Due to Eq. (4), it is sufficient to consider the flow parallel to the pattern (-axis). In this case, the velocity has only one component in the -direction , and we seek a solution for the velocity profile of the form

(11) |

where is the undisturbed linear shear flow. The perturbation of the flow, , which is caused by the presence of the texture and decays far from the interface. It follows from Stokes’ equations that . Since the pressure is constant far from the superhydrophobic surface this implies that it also remains constant in the entire flow region. The velocity perturbation satisfies a dimensionless Laplace equation

(12) |

The boundary conditions at the wall and at infinity are defined in the usual way

(13) | |||||

(14) |

The general solution to Eq. (12), decaying at infinity, can be represented in terms of a cosine Fourier series as

(15) |

where are coefficients to be found by applying the appropriate boundary condition (13). It can be written in terms of the Fourier coefficients as

(16) |

Equation (16) was used to determine the longitudinal effective slip numerically. This was done using a method similar to Asmolov et al. (2013a). We truncate the sum in Eq. (16) at some cut-off number (usually ) and evaluate it in the points where are numbers varying from to . The problem is then reduced to a linear system of equations, where ; , and The system is solved using the IMSL routine LSARG. The transverse effective slip was then calculated using Eq. (4).

## Iii Dissipative Particle Dynamics Simulation

Mesoscale fluid simulations were performed using the dissipative particle dynamics (DPD) method Hoogerbrugge and Koelman (1992); Español and Warren (1995); Groot and Warren (1997). More specifically, we use a version of DPD without conservative interactions Soddemann et al. (2003) and combine that with a tunable-slip method Smiatek et al. (2008) to model patterned surfaces. The tunable-slip method allows one to implement arbitrary hydrodynamic boundary condition in coarse-grained simulations. In the following, we give a brief description of the model and introduce the implementation of planar surfaces with arbitrary one-dimensional slip profile.

The basic DPD equations involve pair interactions between fluid particles. The pair forces are equal in magnitude but opposite in direction; thus the momentum is conserved for the whole system. This leads to the correct long-time hydrodynamic behavior (i.e. Navier-Stokes). The pair force consists of a dissipative contribution and a random component

(17) |

The dissipative part is proportional to the relative velocity between two particles ,

(18) |

with a local distance-dependent friction coefficient . The weight function is a monotonically decreasing function of and vanishes at a given cutoff, and is a parameter that controls the strength of the dissipation. The random force has the form

(19) |

where is a symmetric random number with a zero mean and a unit variance. The amplitude of the stochastic contribution is related to the dissipative contribution by a fluctuation-dissipation theorem in order to ensure that the model reproduces the correct equilibrium distribution.

The tunable-slip method Smiatek et al. (2008) introduces the wall interaction in a similar fashion. The force exerted by the wall on particle is given by

(20) |

The first term is a repulsive interaction that prevents the fluid particles from penetrating the wall. It can be written in term of the gradient of a Weeks-Chandler-Andersen potential Weeks et al. (1971):

(21) |

where is the distance between the fluid particle and the wall, assuming the wall lies in plane, and and set the energy and length scales. In the following, physical quantities will be reported in a model unit system of (length), (mass), and (energy). The dissipative contribution is similar to Eq. (18), with the velocity replaced by the relative particle velocity with respect to the wall. For immobile walls, the dissipative part can be written as

(22) |

The parameter characterizes the strength of the wall friction and can be used to tune the value of the slip length. For example, results in a perfectly slippery wall, whereas a positive value of leads to a finite slip length. The random force has the form

(23) |

where each component of is a random variable. Again, the magnitude of the random force is chosen such that the fluctuation-dissipation theorem is satisfied.

fluid density | |
---|---|

friction coefficient for DPD interaction | |

cutoff for DPD interaction | 1.0 |

no-slip boundary condition () | |

cutoff for wall interaction | 2.0 |

dynamic viscosity | |

hydrodynamic boundary from the simulation wall | |

simulation box | |

pattern period, |

For homogeneous surfaces, a useful analytical expression can be derived which establishes a relation between the slip length and the simulation parameters (the wall friction parameter and the cutoff ) to a good approximation Smiatek et al. (2008). However, the accuracy of the analytic result is not satisfactory at small slip lengths. For the purpose of the present work, we need an accurate expression for the relation between the slip length and the tuning parameter . Therefore, we have simulated plane Poiseuille and Couette flows with different to obtain the slip length and the position of the hydrodynamic boundary. The simulation results were fitted to a truncated Laurent series up to second order in and first order in . The simulation data and the fitted curve are shown in Fig. 2. The fit yields the expression

(24) |

This formula relates the wall friction coefficient , which is set in the simulation, to the slip length . Also shown in Fig. 2 is the analytic result in dashed line, which agrees with the simulation data quite well at large slip length.

One important observation is that both the slip length and the friction coefficient are local quantities. Thus, it is reasonable to assume that the relation (24) also holds for inhomogeneous surfaces, as long as the slip length varies on length scales much larger than the range of the weight functions, and . For a patterned surface with much larger than , the local friction coefficient can be calculated from the slip profile by inverting Eq. (24). The resulting profile is then used in simulations to set up the patterned surfaces.

The simulations are carried out using the open source simulation package ESPResSo Limbach et al. (2006). Modifications have been made to incorporate the patterned surfaces. Table 1 summarizes the simulation parameters. In simulations, the patterned surface has a period of . The simulation box is a rectangular cuboid of size , periodic boundary conditions are assumed in the plane, and the two symmetric patterned surfaces are located at and . The separation between the two surfaces is about twice the pattern periodicity. We have checked that increasing the channel thickness does not change the results.

The simulation starts with randomly distributed particles inside the channel, and flow is induced by applying a body force to all particles. Body forces used in this simulation were chosen small enough that the flow velocity near the wall was less than , in order to ensure that the flow is strictly laminar and effects of finite Reynolds numbers do not affect the simulation results. The magnitude of the slip velocity also affects the spatial resolution in modeling the variation of the slip length. Since the trajectories of the fluid particles are integrated at discrete time steps, the friction experienced by the particle between time steps is determined from the early instant. Therefore, we used a small time step of . Due to the large system size, the relaxation time required for reaching a steady state was very large as well (over time steps). The flow velocity is small in comparison with the thermal fluctuation; thus many time steps are necessary to gather sufficient statistics. In this work, velocity profiles were averaged over time steps. A fit to the plane Poiseuille flow then gives the effective slip length Zhou et al. (2012). The error bars are obtained by six independent runs with different initialization.

Finally, we note that in case of a triangular local slip, only a very small region near the no-slip point contributes to the shear stress at large . This presents a complication for simulations since fluid particles with a finite translational velocity travel a small distance over the no-slip region before reaching the zero averaged velocity. This distance is related to the initial velocity of the fluid particles when they enter the no-slip region, which is in general proportional to and the pressure gradient (the body force in simulation). To obtain accurate simulation data at large , one has to reduce the body force significantly, and the necessary simulation times increase prohibitively. Therefore, only simulation results for and smaller are presented for this texture.

## Iv Eigenvalues of the slip-length tensor

Here we present the results for the eigenvalues of the slip-length tensor obtained by simulations of different local slip profiles. We first verify Eq. (4), derived earlier by some of us Asmolov and Vinogradova (2012) to relate the longitudinal and transverse effective slip lengths. We use local slip profiles with . In that special case, the profile can be characterized by three parameters: the area fraction of the no-slip region, the surface-averaged slip length , and the heterogeneity factor, , which is a measure for the relative variation of the slip length and defined as

(25) |

A homogeneous surface then corresponds to , and low values of corresponds to small . Patterns involving no-slip regions, (i.e. ) have . Such a no-slip assumption can be justified provided , which is typical for many situations. Indeed, for most smooth hydrophobic surfaces, the slip length is of the order of a few tens of nm Vinogradova and Yakubov (2003); Vinogradova et al. (2009); Cottin-Bizonne et al. (2005); Joly et al. (2006), which is very small compared to the slip length values of up to tens or even hundreds of m slip lengths that may be obtained on the gas areas.

### iv.1 Fixed average slip

Let us first investigate the effect of on the effective slip at fixed . Figure 3(a) shows the simulation results for eigenvalues of the slip length tensor as a function of obtained for a trapezoidal local slip (). Also included are numerical results (solid curves) calculated with Eqs. (4) and (16). The agreement between the simulation and numerical results is quite good for all . This demonstrates both the accuracy of our simulations, and the validity of Eq. (4). The data presented in Fig. 3(a) show that the longitudinal slip is larger than the transverse slip, i.e., the flow is anisotropic, and the anisotropy increases with . The effective slip assumes its maximum value at , when the surface is homogeneous, and then decreases monotonously. This implies that the effective slip is found to be always smaller than the surface-averaged value, and the difference from the average value increases with the heterogeneity. A similar behavior was reported in prior work for other types of patterns Vinogradova (1999); Asmolov et al. (2013b). In the case of large local slip, the effective slip length can be obtained from Eq. (5), which gives for a trapezoidal slip profile, Eq. (7):

(26) |

Figure 3(a) includes the theoretical curve expected for this (isotropic) case (dash-dotted line). We note that the slip length is always smaller than the average slip, . Perhaps the most interesting and important result here is that Eq. (26) remains surprisingly accurate even in the case of moderate, but not vanishing, local slip, provided is small enough. We remind the reader that Eq. (26) may not be used when , i.e., for textures with no-slip regions, , where it would simply predicts the effective slip length to be zero.

Qualitatively similar results were obtained for rectangular stripes and triangular local slip profiles. Therefore we do not show these results here. We only recall that in the case of stripes, the asymptotic result for large slip gives (see Asmolov et al. (2013b); Cottin-Bizonne et al. (2004) for the original derivations),

(27) |

and for triangular profiles, one can derive

(28) |

To examine the difference between the different slip profiles, the results for a longitudinal effective slip for a texture with a trapezoidal local slip from Fig. 3(a) are reproduced in Fig. 3(b). Also included in Fig. 3(b) are numerical and simulation data for rectangular () and triangular () profiles with the same average slip. The data show that the triangular profile gives a larger, and the rectangular one a smaller effective slip than that resulting from the trapezoidal local slip. We note that is very different for these slip-length profiles. This suggests that one of the key parameters determining the effective slip is the area fraction of solid in contact with the liquid. If it is very small, the effective slip length becomes large. Another important factor is the presence of flow singularities near slipping heterogeneities as we will discuss below.

### iv.2 Fixed heterogeneity factor

We now fix and consider flows that are maximally anisotropic for a given pattern, which corresponds to (no-slip areas) and . It is instructive to focus first on the role of . Fig. 4 shows the eigenvalues of the slip-length tensor as a function of as obtained from simulations (symbols). As before we include numerical results (solid curves) calculated with Eqs. (4) and (16). We remind that the “gas cushion model”, Eq. (3), becomes very approximate at , and the local slip profile could not necessarily follow the relief of the texture. It is, however, instructive to include large values of into our consideration below since this will provide us with some guidance.

The eigenvalues of the slip-length tensor for a superhydrophobic trapezoidal texture are shown in Fig. 4(a). The agreement between the theoretical and simulation data is very good for all , again confirming the accuracy of the simulations and the validity of Eq. (4). Qualitatively, the results are similar to those obtained earlier for a striped surfaces Belyaev and Vinogradova (2010a): At very small , the effective slip is small and appears to be isotropic, in accordance to predictions of Eq. (6). It becomes anisotropic and increases for larger , and saturates at . In case of stripes with alternating regions of no-slip () and partial-slip () areas the situation can be described by simple analytic formulae Belyaev and Vinogradova (2010a)

(29) | |||||

(30) |

In the limit of , these expressions turn into those derived earlier for perfect slip stripes Philip (1972); Lauga and Stone (2003)

(31) |

For comparison, we plot the above expressions (using the value of the trapezoidal texture) in Fig. 4(a) as dash-dotted [Eqs. (29,30)] and dashed [Eq. (31] lines. One can see that with these parameters, the results for stripes practically coincides with those for trapezoidal slip at very small and very large . At intermediate local slip there is a small discrepancy, suggesting that the equations for partial slip stripes overestimate the effective slip for trapezoidal texture. We will return to a discussion of these results later in the text.

Figure 4(b) plots the simulation results for a triangular local slip profile, which is particularly interesting since . As explained above, simulation data could only be obtained for small values of in this case. In contrast, our numerical approach allowed us to calculate the effective slip length over a much larger interval of . These numerically computed slip lengths are also included in Fig. 4(b). In the regime where simulation data are available, the agreement with the simulation and numerical results is excellent, therefore validating Eq.(4). It is interesting to note that according to the numerical calculations, the effective slip lengths grow weakly with at large and likely do not saturate, in contrast to rectangular and trapezoidal textures.

Finally we compare the effective slip of surfaces with the same and , but different . Fig. 5 shows the corresponding curves for a maximal heterogeneity factor, (), and a relatively small . The results for larger are qualitatively similar and not shown here. Only the values of the effective slip lengths become larger.) We note that by varying from 0 to 0.5, we change the local slip profile from triangular () to rectangular (stripes, ). Regular trapezoidal textures correspond to intermediate values of . Fig. 5 also shows the analytical prediction of Eq. (6) for weakly slipping surfaces, i.e., an effective slip which is isotropic and independent on . The effective slip according to simulations and numerical results is always smaller than the first-order theoretical prediction, and the deviations from the area-average slip model increase with . Furthermore, the data show that the longitudinal and transverse slip lengths are different from each other, i.e., the flow is anisotropic, even though it should be nearly isotropic according to the asymptotic theory, Eq. (6). For striped patterns, similar deviations from the asymptotic theory were shown to result from flow perturbations in the vicinity of jumps in the slip length profile Asmolov et al. (2013a). Here, we find that the friction and anisotropy of the flow past triangular and trapezoidal, weakly slipping surfaces is increased, compared to the expectation from the asymptotic theory, Eq. (6). This somewhat counter-intuitive result suggests that similar, although perhaps weaker, viscous dissipation and singularities of the shear stress can also appear if singularities in the slip length profile only appear in its first derivative. Below we will focus on this phenomenon.

## V Flow singularities near slipping heterogeneities

Motivated by the findings from the previous section, we now discuss the flow near slipping heterogeneity in more detail. We remind the reader that the shear stress is known to be singular near sharp corners of rectangular hydrophilic grooves, being proportional to for longitudinal configurations, and to for transverse configurations Wang (2003) ( is the distance from the corner). In the case of a striped superhydrophobic surface, the shear stress behaves as Sbragaglia and Prosperetti (2007); Asmolov and Vinogradova (2012); Asmolov et al. (2013a), i.e., the singularity is stronger and creates a source of additional viscous dissipation that reduces effective slip Asmolov et al. (2013a). We are now in a position to prove that similar singularities also arise at the border between no-slip and linear local slip regions.

As above, we develop our theory for a longitudinal configuration only. The solution of Laplace equation, Eq. (12), can be constructed in the vicinity of the edge of slipping regions by using polar coordinates Wang (2003); Asmolov et al. (2013a). It is convenient to chose different coordinates for trapezoidal and triangular textures, and to consider these two cases separately. For a trapezoidal slip profile, we use the same coordinates as Wang (2003); Asmolov et al. (2013a) (illustrated in Fig. 6(a)). In the case of triangular textures, for symmetry reasons it is more convenient to use the coordinates shown in Fig. 6(b).

We begin by studying the trapezoidal slip profile. For the trapezoid, the origin of coordinates is at , with and , , so that the n o-slip and the slip regions correspond to and . The slip condition (13) becomes

(32) |

where is defined as

(33) |

and can vary from to , depending on the local slip, , and on the value of . In the case of stripes, one has , resulting in at finite .

The solution of Eq. (12) that satisfies the no-slip boundary condition at is given by

(34) |

and the components of the velocity gradient are

(35) | |||||

(36) |

The velocity at the edge is finite provided but its gradient becomes singular when . In this case the two terms in the left-hand side of Eq. (32) are of the same order, , but the right-hand side is smaller and can be neglected, since as . As a result Eq. (32) can be rewritten, taking into account (34-36) as

or equivalently as,

(37) |

The theoretical prediction of Eq. (37) for the behavior of as a function of is shown in Fig. 7 (solid curve). Note that decays monotonically from down to . In other words, for any finite , we have and the shear stress becomes singular. As shown in Fig. 7 these theoretical predictions are in excellent quantitative agreement with the simulation results, where we have measured the slip velocity along the surface and then obtained by using Eq. (34). Altogether the theoretical and numerical results confirm that the singular is weaker for trapezoidal profiles than for stripes, such that the viscous dissipation is smaller. This explains qualitatively the results presented in Fig. 5, where the slip length deviates more strongly from the area-averaged slip length on striped than on trapezoidal textures.

In the theoretical analysis of the triangular slip profile () it is convenient to shift the origin of the angular variable compared to the trapezoidal case (see Fig. 6), so that

The even solution of the Laplace equation (12) is

(38) |

and the components of the velocity gradient are then

(39) | |||||

(40) |

Applying the condition , we can neglect the right-hand side to obtain

(43) |

The numerical solution of this equation is included in Fig. 7 (dashed curve), and supported by simulation data. We see that again decays monotonically from , but at large the solution of Eq. (43) is

i.e., different from expected for a trapezoidal profile, where would tend to 0.5. This implies that for given , the triangular slip profile always gives a stronger singularity compared to a trapezoid slip profile. At small , the singularity is weaker than for a striped surface.

A striking conclusion from our analysis is that non-smooth surface textures with continuous local slip can generate stronger singularities in flow and shear stress than surfaces with a discontinuous, piece-wise constant slip length profile. Indeed, we have for (or ), i.e., the singularity for a triangular texture becomes stronger than for stripes and should significantly reduce the longitudinal effective slip. This is likely a reason why the effective slip in Fig. 4(b) grows only weakly at very large . Coming back to Fig. 5, we remark that the results for triangular profile there correspond to and , i.e., the singularity is weaker in this case than for stripes with . The singularities for trapezoidal profiles may be even weaker: The exponents vary from to . Note however that since the no-slip area remains the key factor determining the effective slip, always decreases with .

To calculate the exponents for a transverse configuration, one can use the relation between the velocity fields in eigendirections Asmolov and Vinogradova (2012):

(44) | |||

(45) |

where is the velocity field for the longitudinal pattern with a twice larger local slip length. It follows from Eqs. (44) and (45), that the velocity gradients and the pressure in the transverse configuration for both trapezoidal and triangular textures have stronger singularities at the edges of slipping regions than those for the longitudinal configuration:

where . This leads to larger viscous dissipation. As a result the anisotropy of the slip-length tensor is finite even at small and is approximately the same for triangular, trapezoidal and stripe profiles ( see Fig. 5). For example, the eigenvalues shown there for a triangular texture correspond to and . Therefore, in the transverse configuration even for a shallow triangular texture, , the singularity becomes stronger than for stripes (.

## Vi Optimization of the effective slip

Based on the results obtained in the previous sections, we now discuss how to optimize the effective slippage, and what maximum slip lengths may actually be expected, taking into account the actual limitation in surface engineering.

Very large effective slip lengths were recently found for sinusoidal profiles with a no-slip point Asmolov et al. (2013b). If we request that the local slip pattern includes a no-slip area, we also find here that the largest values of are obtained for triangular profiles, where the no-slip area reduces to just one point. However, real textures always have finite contact area with liquid. The effective slip length of a trapezoidal texture with large is close to that of stripes [see Fig. 4(a)]. In this final section, we therefore focus on the case of a trapezoidal profile with small , where both and the slope of the inclined regions is large. Such a texture should lead to large . To find one-dimensional textures which optimize the effective forward slip we consider profiles of the following form

(46) |

The bottom base of the trapezoid for such textures is wide, , [see Fig. 8(a)]. For comparison, we also consider stripes with the same maximal local slip length, :

(47) |

Both profiles are characterized by only one parameter .

The simulation and numerical results for and are shown in Fig. 8(b). The effective slip lengths for the trapezoid and stripe profiles are nearly equal for . For example, for the stripe profile with is only slightly larger than the corresponding value for the trapezoid profile with the same and wide bottom base [see Fig. 8(a)]. Note that such narrow trapezoids with wide bottom base allowed researchers to achieve a very large effective slip in recent experiments Choi et al. (2006). The same value of can be achieved with a triangular profile with . Thus the effect of the inclined regions is relatively small for large slopes, and the main contribution to the effective slip length stems from the no-slip region and the singularity of the shear stress.

## Vii Concluding Remarks

In conclusion, we have presented simulation and numerical results and some asymptotic law analysis that allowed us to assess the frictional properties of superhydrophobic surfaces as a function of the texture geometric parameters. We have focused on one-dimensional superhydrophobic surfaces with trapezoidal grooves and systematically investigated the effect of several parameters, such as the maximum local slip and the heterogeneity factor. The slip properties on a trapezoidal surface were compared with two limiting cases: striped and triangular slip length profiles. Our simulation and numerical results suggested that the effective slip length is mainly determined by (i) the area fraction of the no-slip region and (ii) the singularity of the shear stress at the no-slip edge. The first factor is obvious; a smaller fraction of no-slip area results in less friction of the surface, which in turn leads to large slip length. The second effect is less intuitive, especially for a continuous slip profile. Here we demonstrated that the singularity develops for the trapezoidal and triangular local slip profiles, and in the latter case can be even stronger than for discontinuous slip of striped textures. The singularity is stronger for a transverse configuration in comparison to a longitudinal configuration, and as a result, it always contributes to the flow anisotropy. Our results can guide the design of superhydrophobic surfaces for micro/nanofluidics.

###### Acknowledgements.

This research was supported by the RAS through its priority program “Assembly and Investigation of Macromolecular Structures of New Generations”, and by the DFG through SFB TR6 and SFB 985. The simulations were carried out using computational resources at the John von Neumann Institute for Computing (NIC Jülich), the High Performance Computing Center Stuttgart (HLRS) and Mainz University (MOGON).### References

- D. Quere, Rep. Prog. Phys. 68, 2495 (2005).
- D. Richard and D. Quere, Europhys. Lett. 48, 286 (1999).
- D. Richard and D. Quere, Europhys. Lett. 50, 769 (2000).
- P. Tsai, R. C. A. van der Veen, M. van de Raa, and D. Lohse, Langmuir 26, 16090 (2010).
- O. I. Vinogradova and A. L. Dubov, Mendeleev Commun. 19, 229 (2012).
- G. McHale, M. I. Newton, and N. J. Schirtcliffe, Soft Matter 6, 714 (2010).
- L. Bocquet and J. L. Barrat, Soft Matter 3, 685 (2007).
- J. P. Rothstein, Annu. Rev. Fluid Mech. 42, 89 (2010).
- O. I. Vinogradova and A. V. Belyaev, J. Phys.: Cond. Matter 23, 184104 (2011).
- K. Kamrin, M. Bazant, and H. A. Stone, J. Fluid Mech. 658, 409 (2010).
- M. Z. Bazant and O. I. Vinogradova, J. Fluid Mech. 613, 125 (2008).
- S. Schmieschek, A. V. Belyaev, J. Harting, and O. I. Vinogradova, Phys. Rev. E 85, 016324 (2012).
- O. I. Vinogradova, Langmuir 11, 2213 (1995).
- D. Andrienko, B. Dünweg, and O. I. Vinogradova, J. Chem. Phys. 119, 13106 (2003).
- C. Ybert, C. Barentin, C. Cottin-Bizonne, P. Joseph, and L. Bocquet, Phys. Fluids 19, 123601 (2007).
- D. Maynes, K. Jeffs, B. Woolford, and B. W. Webb, Phys. Fluids 19, 093603 (2007).
- T. V. Nizkaya, E. S. Asmolov, and O. I. Vinogradova, Soft Matter (2013), doi:10.1039/C3SM51850G.
- A. D. Stroock, S. K. Dertinger, G. M. Whitesides, and A. Ajdari, Anal. Chem. 74, 5306 (2002).
- A. Ajdari, Phys. Rev. E 65, 016301 (2001).
- E. Lauga and H. A. Stone, J. Fluid Mech. 489, 55 (2003).
- A. V. Belyaev and O. I. Vinogradova, J. Fluid Mech. 652, 489 (2010a).
- C. O. Ng, H. C. W. Chu, and C. Y. Wang, Phys. Fluids 22, 102002 (2010).
- J. Zhou, A. V. Belyaev, F. Schmid, and O. I. Vinogradova, J. Chem. Phys. 136, 194706 (2012).
- E. S. Asmolov, J. Zhou, F. Schmid, and O. I. Vinogradova, Phys. Rev. E 88, 023004 (2013a).
- E. S. Asmolov and O. I. Vinogradova, J. Fluid Mech. 706, 108 (2012).
- E. S. Asmolov, S. Schmieschek, J. Harting, and O. I. Vinogradova, Phys. Rev. E 87, 023005 (2013b).
- S. C. Hendy and N. J. Lund, Phys. Rev. E 76, 066313 (2007).
- C. Cottin-Bizonne, C. Barentin, E. Charlaix, L. Bocquet, and J.-L. Barrat, Eur. Phys. J. E 15, 427 (2004).
- N. V. Priezjev, A. A. Darhuber, and S. M. Troian, Phys. Rev. E 71, 041608 (2005).
- N. V. Priezjev, J. Chem. Phys. 135, 204704 (2011).
- N. Tretyakov and M. Müller, Soft Matter 9, 3613 (2013).
- W. Li, X. S. Cui, and G. P. Fang, Langmuir 26, 3194 (2010).
- T. Tanaka, M. Morigami, and N. Atoda, Jpn. J. Appl. Phys. 32, 6069 (1993).
- C.-H. Choi, U. Ulmanella, J. Kim, C.-M. Ho, and C.-J. Kim, Phys. Fluids 18, 087105 (2006).
- A. V. Belyaev and O. I. Vinogradova, Soft Matter 6, 4563 (2010b).
- E. S. Asmolov, A. V. Belyaev, and O. I. Vinogradova, Phys. Rev. E 84, 026330 (2011).
- F. Feuillebois, M. Z. Bazant, and O. I. Vinogradova, Phys. Rev. Lett. 102, 026001 (2009).
- M. Sbragaglia and A. Prosperetti, Phys. Fluids 19, 043603 (2007).
- J. Hyväluoma and J. Harting, Phys. Rev. Lett. 100, 246001 (2008).
- A. Steinberger, C. Cottin-Bizonne, P. Kleimann, and E. Charlaix, Nature Materials 6, 665 (2007).
- E. Karatay, A. S. Haase, C. W. Visser, C. Sun, D. Lohse, P. A. Tsai, and R. G. H. Lammertink, PNAS 110, 8422 (2013).
- A. S. Haase, E. Karatay, P. A. Tsai, and R. G. H. Lammertink, Soft Matter 9, 8949 (2013).
- P. J. Hoogerbrugge and J. M. V. A. Koelman, Europhys. Lett. 19, 155 (1992).
- P. Español and P. Warren, Europhys. Lett. 30, 191 (1995).
- R. D. Groot and P. B. Warren, J. Chem. Phys. 107, 4423 (1997).
- T. Soddemann, B. Dünweg, and K. Kremer, Phys. Rev. E 68, 046702 (2003).
- J. Smiatek, M. Allen, and F. Schmid, Eur. Phys. J. E 26, 115 (2008).
- J. D. Weeks, D. Chandler, and H. C. Andersen, J. Chem. Phys. 54, 5237 (1971).
- H. Limbach, A. Arnold, B. Mann, and C. Holm, Comput. Phys. Commun. 174, 704 (2006).
- O. I. Vinogradova and G. E. Yakubov, Langmuir 19, 1227 (2003).
- O. I. Vinogradova, K. Koynov, A. Best, and F. Feuillebois, Phys. Rev. Lett. 102, 118302 (2009).
- C. Cottin-Bizonne, B. Cross, A. Steinberger, and E. Charlaix, Phys. Rev. Lett. 94, 056102 (2005).
- L. Joly, C. Ybert, and L. Bocquet, Phys. Rev. Lett. 96, 046101 (2006).
- O. I. Vinogradova, Int. J. Miner. Proc. 56, 31 (1999).
- J. R. Philip, J. Appl. Math. Phys. 23, 353 (1972).
- C. Y. Wang, Phys. Fluids 15, 1114 (2003).