1 Introduction
Superhydrophobic surfaces are typically no-slip surfaces endowed with additional microstructural features such as grooves, posts or holes and they can significantly reduce fluid drag in microchannels (Rothstein Reference Rothstein2010). In many cases, in the Cassie state the microstructural surface features are occupied by pockets of gas rather than being liquid-filled (as in the Wenzel state). The spanning gas–liquid interfaces, or menisci, that can protrude into or out of the fluid are then often close to being shear-free and allow slip. The drag reduction characteristics of such a surface have a strong dependence on the geometrical properties of these gas–liquid interfaces (Steinberger et al. Reference Steinberger, Cottin-Bizonne, Kleimann and Charlaix2007) and it is an important matter to be able to quantify them. For shear flows with shear rate $\dot{{\it\gamma}}$ over such a surface occupying the plane $y=0$ the velocity field far from the plane of the surface takes the form
where $\hat{\boldsymbol{x}}$ is the flow direction. The constant ${\it\lambda}$ is the effective slip length and is a measure of the frictional properties of the surface: it is the fictional distance below the surface at which the shear flow would extrapolate to zero.
In applications an important class of superhydrophobic surfaces comprises those with a periodic array of grooves aligned with the principal flow direction. Figure 1 shows a schematic diagram of such a surface in the Cassie state when a periodic array of parallel circular menisci of width $2c$ and period $2l$ , each having protrusion angle ${\it\theta}$ , span unidirectional grooves aligned along a $z$ -axis, say, with no-slip regions in between. Also shown is the period window $-l<x<l$ with the meniscus lying in the subinterval $-c<x<c$ . An important geometrical characteristic of such a ‘bubble mattress’ is the no-shear fraction $c/l$ . When no-shear conditions are imposed on the menisci a complete analytical characterization of the effective slip properties of such surfaces is known in the dilute limit of small no-shear fraction. Davis & Lauga (Reference Davis and Lauga2009) derived an explicit formula for transverse slip length over such a surface; Crowdy (Reference Crowdy2010) derived the corresponding result for longitudinal slip length. These results are valuable together because, for small no-shear fractions and zero capillary number, the slip length for a linear shear in a general direction over the surface is a linear combination of these results.
In other analytical work, Sbragaglia & Prosperetti (Reference Sbragaglia and Prosperetti2007) used a boundary perturbation analysis to investigate the effects of small interface curvature/protrusion on the effective slip length in a pressure-driven flow in channels. Their analysis is valid for any no-shear fraction, but it is limited to small protrusion angles of the menisci (close to flat). Numerical calculations of the slip lengths associated with such bubble mattresses have been performed by Teo & Khoo (Reference Teo and Khoo2010) and Ng & Wang (Reference Ng and Wang2011).
The present paper extends the author’s earlier result (Crowdy Reference Crowdy2010) beyond the dilute limit and produces formulae for the longitudinal slip length over these unidirectional surfaces that are accurate over a much larger range of no-shear fractions. Crowdy (Reference Crowdy2010) derived the longitudinal slip length ${\it\lambda}_{0}$ in the dilute limit as
The first of two new results of the present paper is to derive the more accurate formula:
where we use ${\it\lambda}_{1}$ to denote the improved result. We also write down the associated flow field. For no-shear fractions as high as $c/l=0.75$ (1.3) gives agreement with maximum relative errors (across protrusion angles) of 6–7 % compared to the numerical results of Teo & Khoo (Reference Teo and Khoo2010). Clearly (1.3) reduces to (1.2a−c ) to leading order in $c{\it\delta}$ (or, equivalently, to leading order in the no-shear fraction $c/l$ ). We have not presented (1.3) as an expansion in powers of $c/l$ because the functional form (1.3) arises naturally in our analysis and it turns out to give more accurate results than a truncated formal expansion.
The second result is to derive an even higher-order approximation, denoted by ${\it\lambda}_{2}$ , reported later in (4.13). It is not as simple to write down as (1.3), but it is nonetheless explicit and reduces the maximum relative errors in the slip length to as little as 1–2 % for no-shear fraction $c/l=0.75$ ; even for a no-shear fraction as high as $c/l=0.9$ , it is still accurate to within a maximum relative error of 8–9 % (see figure 3). Moreover, an approximation of the higher-order formula (4.13) (which ignores certain terms in (4.13) that are sixth order in the no-shear fraction) is the more concise, and only slightly less accurate, expression
with
This formula gives the required slip length with maximum relative error 2 % for no-shear fractions in the range $c/l\in [0,0.75]$ and 12 % for $c/l=0.9$ .
Formulae (1.3), (4.13) and (1.4) are the main results of this paper. They are valid without restriction on the protrusion angle ${\it\theta}$ . We show that they extend an analytical result of Sbragaglia & Prosperetti (Reference Sbragaglia and Prosperetti2007) beyond the limit of small protrusion angle. The formulae can, in principle, be derived using the usual formalisms of matched asymptotic expansions but we derive them using our own mathematical approach, which is of theoretical interest in itself.
2 The dilute limit
Crowdy (Reference Crowdy2010) has previously derived the formula (1.2a−c ) for the longitudinal slip length relevant in the dilute limit $c/l\rightarrow 0$ ; it is natural to think of fixing $c=1$ , say, and taking $l\rightarrow \infty$ . Those prior results can be summarized as follows. With $z=x+\text{i}y$ we consider a complex potential $h_{s}(z)$ , analytic in the fluid region, such that
The imposed far-field condition of simple shear requires that $w(x,y)\sim \dot{{\it\gamma}}y$ , hence
To satisfy the no-slip condition on the wall we need
while the no-stress condition on the meniscus requires that
which follows from the condition $\partial w/\partial n=0$ on use of the Cauchy–Riemann equations.
To solve this problem Crowdy (Reference Crowdy2010) employed conformal mapping techniques. With ${\it\beta}={\rm\pi}-{\it\theta}$ the conformal mapping transplanting the upper half of the unit disc in a parametric complex ${\it\zeta}$ -plane to the fluid region, and its inverse mapping function, are
Having derived these, Crowdy (Reference Crowdy2010) uses them to establish that
where
Actually, only ${\it\Delta}_{1}$ is needed to derive (1.2a−c ) based on a simple superposition argument.
Since his focus was to find the slip length Crowdy (Reference Crowdy2010) did not report the corresponding dilute-limit complex potential, but it will be useful for what follows. It turns out to be
We now explain this result because it helps to understand our derivations of the higher-order approximations. If $h(z)$ denotes the exact solution to the periodic problem then $h(z)$ must be analytic in the period window $D$ shown in figure 1 and satisfy the quasi-periodicity condition
in order that the velocity $w=\dot{{\it\gamma}}\text{Im}[h(z)]$ is periodic; the quasi-periodicity here is induced by the far-field $z$ behaviour. Similarly, $h_{0}(z)$ is required to be analytic in $D$ . This is true of $h_{0}(z)$ given in (2.8) since $h_{s}(z)$ is analytic there – it is analytic, by construction, everywhere in the upper half-plane and above the meniscus – while the apparent singularity of $h_{0}(z)$ at $z=0$ , which is inside the period window if ${\it\theta}<0$ , is in any case removable by the choice of ${\it\lambda}_{0}$ . Moreover, the cotangent function has a periodic array of singularities along the real axis but the two closest to the one at $z=0$ are in the two period windows neighbouring $D$ and not inside $D$ itself.
The first term added to $h_{s}(z)$ in (2.8) serves the purpose of removing the $O(1/z)$ behaviour of $h_{s}(z)$ as $z\rightarrow \infty$ so that, as $z\rightarrow \infty$ ,
where we have used the fact that $\cot ({\it\delta}z)\rightarrow -\text{i}$ as $y\rightarrow +\infty$ and the known far-field asymptotics (2.6) of $h_{s}(z)$ . This means that, on the edges of the period window where $|z|\geqslant l$ , the function $h_{0}(z)$ satisfies
which is a good approximation to the quasi-periodicity condition (2.9) if ${\it\delta}$ is small. It is easily checked that $h_{0}(z)$ is real when $z$ is real since this is true of $h_{s}(z)$ , implying that $h_{0}(z)$ satisfies the no-slip condition on the wall. Also, provided ${\it\delta}$ is small, then on the meniscus
Thus $h_{0}(z)$ also satisfies the meniscus boundary condition correct to $O(c^{2}{\it\delta}^{2})$ . In this way we have verified that $h_{0}(z)$ is the approximation to the required complex potential provided that ${\it\delta}$ is small.
3 An improved slip length formula
To produce an improved formula for the slip length, accurate at even larger no-shear fractions, it is natural to seek a higher-order approximation in ${\it\delta}c$ . Consider the modified complex potential
with ${\it\lambda}_{1}$ now chosen to satisfy
With this choice the $1/z$ term in the far-field expansion of the three non-cot terms
of $h_{1}(z)$ vanishes. This ensures that, as $z\rightarrow \infty$ ,
where we have again used the fact that $\cot ({\it\delta}z)\rightarrow -\text{i}$ as $y\rightarrow +\infty$ and the known far-field asymptotics (2.6) of $h_{s}(z)$ . Hence $h_{1}(z)$ has the correct far-field behaviour with slip length ${\it\lambda}_{1}$ . On the edges of the period window where $|z|\geqslant l$ , we still have
which is a good approximation to the required quasi-periodicity (2.9) if ${\it\delta}$ is small. The apparent singularity of the function (3.1) at $z=0$ is again removable while the nearest other singularities of the cotangent function are in the neighbouring period windows as before. It is therefore confirmed that $h_{1}(z)$ is analytic in $D$ . When $z$ is real, which is true on the no-slip wall, so is $h_{1}(z)$ , confirming that the no-slip condition is satisfied there. Finally, on expansion of the cotangent for small ${\it\delta}$ , notice that on the meniscus,
where we have used (2.4) twice as well as the Laurent expansion of $\cot ({\it\delta}z)$ about $z=0$ . Hence, (3.1) satisfies the boundary conditions on the meniscus, the no-slip wall and the edges of the period window correct to $O(c^{3}{\it\delta}^{3})$ . The quantity ${\it\lambda}_{1}$ as given by (3.2) is now the required slip length at this order of approximation and produces the result (1.3). Teo & Khoo (Reference Teo and Khoo2010) report the slip lengths ${\it\lambda}_{TK}$ , say, with the renormalizations
where we have taken $c=1$ . Figure 2 shows ${\it\lambda}_{TK}$ , as a function of protrusion angle ${\it\theta}$ , for $c/l=0.1$ , 0.25, 0.5, 0.75 and 0.9 together with corresponding data points from Teo & Khoo (Reference Teo and Khoo2010), which serve as our benchmark solution. The maximum relative error of this approximation for $c/l=0.75$ is between 6 and 7 % and for $c/l=0.9$ is around 25 %.
4 Higher-order analysis
The analysis of the previous section is readily extended to even higher order. We introduce $\tilde{h}_{s}(z)$ as the complex potential for a single protruding bubble with all the same boundary conditions as for $h_{s}(z)$ but now satisfying the modified far-field condition
This far-field flow is no longer a simple shear. By the same conformal mapping arguments (Crowdy Reference Crowdy2010) used to find $h_{s}(z)$ it can be shown that
where we have used (2.5) and
Now consider the new complex potential, constructed using $h_{s}(z)$ and $\tilde{h}_{s}(z)$ , given by
where the real parameters ${\it\lambda}_{2}$ and ${\it\mu}_{2}$ are chosen to satisfy
We claim that $h_{2}(z)$ is the required higher-order solution. To see this, the expansions
reveal that the third-order singularity of $h_{2}(z)$ at $z=0$ is removable. Furthermore,
where we have used the facts that $\cot ({\it\delta}z)\rightarrow -\text{i}$ and $\cot ({\it\delta}z)\text{cosec}^{2}({\it\delta}z)$ decays exponentially as $y\rightarrow +\infty$ , as well as the far-field forms (2.6) and (4.2) of $h_{s}(z)$ and $\tilde{h}_{s}(z)$ . Then on the edges of the period window where $|z|\geqslant l$ ,
so that $h_{2}(z)$ satisfies the required quasi-periodicity condition there (now to higher order than before). ${\it\lambda}_{2}$ is the required slip length. For $z$ on the no-slip wall $h_{2}(z)$ is real, implying that the no-slip condition is satisfied there. On expanding (4.4) for small ${\it\delta}$ , and on use of (4.5),
where we have used the Laurent expansions (4.6) and the fact that
Hence the meniscus boundary condition is satisfied correct to $O(c^{6}{\it\delta}^{6})$ . The solution of the $2\times 2$ system (4.5) gives the required slip length ${\it\lambda}_{2}$ . In matrix form it is
where
In view of (2.7) and (4.3) all matrix elements are known as explicit functions of $a$ and $c$ , consequently an explicit formula for the slip length ${\it\lambda}_{2}$ is
The normalized slip length (3.7) – but now with ${\it\lambda}_{1}$ replaced by ${\it\lambda}_{2}$ as given by (4.13) – is plotted in figure 3, as a function of ${\it\theta}$ for $c/l=0.1$ , 0.25, 0.5, 0.75 and 0.9. It gives markedly better agreement with the numerical solution at $c/l=0.9$ than that shown in figure 2. The maximum relative error of this approximation for $c/l=0.75$ is between 1 and 2 %, and for $c/l=0.9$ it is between 8 and 9 %.
Our analysis leads naturally to formula (4.13) but if we use the second equation in (4.5) to approximate
then, on substitution of this into the first equation in (4.5), we find the approximation
On eliminating $a$ in favour of ${\it\theta}$ , it can be shown that
Equation (4.15) is now equivalent to (1.4) reported earlier. The maximum relative error of this approximation is not quite as good as for formula (4.13), but is still impressive: at $c/l=0.75$ the maximum relative error is around 2 % (about the same as for (4.13)) and for $c/l=0.9$ it is around 12 % (slightly worse). The surprisingly wide range of no-shear fractions for which these formulae give good accuracy arguably obviates the need for higher-order approximations, but these are easily derivable in principle.
5 Verification against previous results
When formally expanded in powers of the no-shear fraction our formulae cross-check with other known results. First, Philip (Reference Philip1972) has shown that for a flat meniscus with ${\it\theta}=0$ , $c=1$ and at any no-shear fraction,
But for ${\it\theta}=0$ we find ${\it\alpha}(0)=1/2$ and ${\it\beta}(0)=1/8$ and it can be verified that (1.4) agrees with (5.1) to this order of expansion in ${\it\delta}$ .
By taking their microchannel width to infinity, Sbragaglia & Prosperetti (Reference Sbragaglia and Prosperetti2007) performed an analytical study for small protrusion angles ${\it\theta}$ in a semi-infinite domain. The expansion of (4.15) for small ${\it\delta}$ is
A Taylor expansion of ${\it\alpha}({\it\theta})$ and ${\it\beta}({\it\theta})$ to linear order in the protrusion angle ${\it\theta}$ gives
On substitution of (5.3) into (5.2), we find
The first term agrees with (5.1), valid for ${\it\theta}=0$ , once we set $c=1$ . On setting ${\it\delta}={\rm\pi}{\it\xi}/2c$ the normalized leading-order correction due to the meniscus curvature is
where, for small meniscus deflection downwards into the grooves, ${\it\theta}\approx \sin {\it\theta}=-c/R$ . We thus retrieve the result in (45)–(47) of Sbragaglia & Prosperetti (Reference Sbragaglia and Prosperetti2007), who derived it using quite different techniques. In summary, various series expansions of our formulae (1.3) and (4.13), or (1.4), give results that are consistent with other studies. Note, however, that we continue to present (1.3) and (1.4) in the unexpanded form arising naturally from our analysis.
6 A reciprocity result
In view of their convenient explicit forms we expect that the slip length formulae (1.3) and (1.4), and the associated complex potentials (3.1) and (4.4), will be useful in a variety of studies of superhydrophobic surfaces where, for example, additional physical effects are included (such as heat transfer, Marangoni or thermocapillary effects, or the influence of an enclosed gas phase). We end by showing how to combine the new formulae with a useful reciprocity result based on one of Green’s identities.
Let $w_{2}(z)=\text{Im}[\dot{{\it\gamma}}h_{2}(z)]$ , with associated slip length ${\it\lambda}_{2}$ . Suppose ${\hat{w}}$ is the solution of the same problem of shear flow over the bubble mattress with all the same boundary conditions imposed on $w_{2}(z)$ except on the meniscus, where we now require
where $\mathscr{A}(z,\overline{z})$ is some specified function (derived, say, from inclusion of additional physical effects). Let the slip length for this modified problem be $\hat{{\it\lambda}}$ . By Green’s identity, the harmonicity of $w_{2}$ and ${\hat{w}}$ , and the divergence theorem, we deduce that
where $\partial D$ is the boundary of $D$ . On the edges of the period window for which $|z|\geqslant l$ ,
where we have used (4.7). Hence the contributions to the right-hand side of (6.2) from the side edges are small (for small ${\it\delta}$ ) while the no-slip portions of the surface do not contribute at all. We then arrive at the approximation
where we have integrated around $\partial D$ in an anticlockwise direction and used the far-field conditions. By virtue of the results of this paper the right-hand side of (6.4) is an integral expression for the required slip length $\hat{{\it\lambda}}$ accurate to the same order in the no-shear fraction as the solution $h_{2}(z)$ and ${\it\lambda}_{2}$ used to obtain it.
7 Discussion
An asymptotic analysis similar to that presented here is, in principle, possible to generalize the dilute approximation to the transverse slip length over this class of surfaces found by Davis & Lauga (Reference Davis and Lauga2009); the biharmonic nature of the field equations there render the technical details more challenging. But complementary results in this direction would provide a fuller analytical description of the so-called slip tensor (Bazant & Vinogradova Reference Bazant and Vinogradova2008; Asmolov & Vinogradova Reference Asmolov and Vinogradova2012) for these surfaces, thereby generalizing the fairly complete analytical description in the dilute limit already available on combining the analytical formulae of Davis & Lauga (Reference Davis and Lauga2009) and Crowdy (Reference Crowdy2010).
The study here has imposed idealized shear-free boundary conditions at the liquid–gas interfaces and ignores additional dissipation associated with an enclosed gas. Incorporating dissipation in the gas subphase is a topic of much recent research: Schönecker & Hardt (Reference Schönecker and Hardt2013) and Schönecker, Baier & Hardt (Reference Schönecker, Baier and Hardt2014) have proposed a semianalytical method that approximates the liquid–gas interface as a constant shear-stress boundary (leading to a non-uniform local slip length); Nizkaya, Asmolov & Vinogradova (Reference Nizkaya, Asmolov and Vinogradova2014) have put forward a general ‘gas cushion model’ based on an operator method. Those authors have executed detailed studies of their models in the case of unidirectional (one-dimensional) surfaces with flat interfaces. Our results here – especially coupled with the reciprocity result of § 6 – might well be useful in extending those investigations to the important case of curved menisci.
Acknowledgements
This work is supported by EPSRC Fellowship EP/K019430/1, a Royal Society Wolfson Research Merit Award and EPSRC grant EP/K041134/1. The author is grateful to the authors of Teo & Khoo (Reference Teo and Khoo2010) for providing comparison data.