Hostname: page-component-78c5997874-v9fdk Total loading time: 0 Render date: 2024-11-05T04:20:39.718Z Has data issue: false hasContentIssue false

Analytical formulae for longitudinal slip lengths over unidirectional superhydrophobic surfaces with curved menisci

Published online by Cambridge University Press:  24 February 2016

Darren G. Crowdy*
Affiliation:
Department of Mathematics, Imperial College London, 180 Queen’s Gate, London SW7 2AZ, UK
*
Email address for correspondence: [email protected]

Abstract

This paper reports new analytical formulae for the longitudinal slip lengths for simple shear over a superhydrophobic surface, or bubble mattress, comprising a periodic array of unidirectional circular menisci, or bubbles, protruding into, or out of, the fluid. The accuracy of the formulae is tested against results from full numerical simulations; they are found to give small relative errors even at large no-shear fractions. In the dilute limit the formulae reduce to an earlier result by Crowdy (Phys. Fluids, vol. 22, 2010, 121703). They also extend analytical results of Sbragaglia & Prosperetti (Phys. Fluids, vol. 19, 2007, 043603) beyond the limit of a small protrusion angle.

Type
Rapids
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
© 2016 Cambridge University Press

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

(1.1) $$\begin{eqnarray}\boldsymbol{u}=\dot{{\it\gamma}}(y+{\it\lambda})\hat{\boldsymbol{x}},\end{eqnarray}$$

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.

Figure 1. Longitudinal shear flow over a surface comprising a unidirectional periodic array, with period $2l$ , of menisci of width $2c$ protruding at angle ${\it\theta}$ with intermediate no-slip zones.

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

(1.2ac ) $$\begin{eqnarray}\frac{{\it\lambda}_{0}}{c}=c{\it\delta}{\it\alpha}({\it\theta}),\quad {\it\alpha}({\it\theta})=\frac{3{\rm\pi}^{2}-4{\rm\pi}{\it\theta}+2{\it\theta}^{2}}{6({\rm\pi}-{\it\theta})^{2}},\quad {\it\delta}\equiv \frac{{\rm\pi}}{2l}.\end{eqnarray}$$

The first of two new results of the present paper is to derive the more accurate formula:

(1.3) $$\begin{eqnarray}\frac{{\it\lambda}_{1}}{c}=\frac{c{\it\delta}{\it\alpha}({\it\theta})}{\displaystyle \left[1-\frac{(c{\it\delta})^{2}{\it\alpha}({\it\theta})}{3}\right]},\end{eqnarray}$$

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.2ac ) 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

(1.4) $$\begin{eqnarray}\frac{{\it\lambda}_{2}}{c}\approx \frac{\displaystyle c{\it\delta}{\it\alpha}({\it\theta})\left[1+\frac{(c{\it\delta})^{4}{\it\beta}({\it\theta})}{15}\right]}{\displaystyle \left[1-\frac{(c{\it\delta})^{2}{\it\alpha}({\it\theta})}{3}-\frac{(c{\it\delta})^{4}{\it\beta}({\it\theta})}{15}\right]},\end{eqnarray}$$

with

(1.5) $$\begin{eqnarray}{\it\beta}({\it\theta})=\frac{1}{360({\rm\pi}-{\it\theta})^{4}}(32{\it\theta}^{4}-128{\rm\pi}{\it\theta}^{3}+212{\rm\pi}^{2}{\it\theta}^{2}-168{\rm\pi}^{3}{\it\theta}+45{\rm\pi}^{4}).\end{eqnarray}$$

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.2ac ) 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

(2.1) $$\begin{eqnarray}w(x,y)=\dot{{\it\gamma}}\,\text{Im}[h_{s}(z)].\end{eqnarray}$$

The imposed far-field condition of simple shear requires that $w(x,y)\sim \dot{{\it\gamma}}y$ , hence

(2.2) $$\begin{eqnarray}h_{s}(z)\sim z+O(1/z),\quad \text{as}~z\rightarrow \infty .\end{eqnarray}$$

To satisfy the no-slip condition on the wall we need

(2.3) $$\begin{eqnarray}\text{Im}[h_{s}(z)]=0,\quad \text{on}~y=0,x\notin [-c,c],\end{eqnarray}$$

while the no-stress condition on the meniscus requires that

(2.4) $$\begin{eqnarray}\text{Re}[h_{s}(z)]=0,\quad \text{on the meniscus},\end{eqnarray}$$

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

(2.5a,b ) $$\begin{eqnarray}z({\it\zeta})=c\left[\frac{(1-{\it\zeta})^{2{\it\beta}/{\rm\pi}}+(1+{\it\zeta})^{2{\it\beta}/{\rm\pi}}}{(1-{\it\zeta})^{2{\it\beta}/{\rm\pi}}-(1+{\it\zeta})^{2{\it\beta}/{\rm\pi}}}\right],\quad {\it\zeta}(z)=\frac{(z/c-1)^{{\rm\pi}/2{\it\beta}}-(z/c+1)^{{\rm\pi}/2{\it\beta}}}{(z/c-1)^{{\rm\pi}/2{\it\beta}}+(z/c+1)^{{\rm\pi}/2{\it\beta}}}.\end{eqnarray}$$

Having derived these, Crowdy (Reference Crowdy2010) uses them to establish that

(2.6) $$\begin{eqnarray}h_{s}(z)=-\frac{2{\rm\pi}c}{{\it\beta}}\left[\frac{((z/c)^{2}-1)^{{\rm\pi}/2{\it\beta}}}{(z/c-1)^{{\rm\pi}/{\it\beta}}-(z/c+1)^{{\rm\pi}/{\it\beta}}}\right]\sim z+\frac{{\it\Delta}_{1}}{z}+\frac{{\it\Delta}_{3}}{z^{3}}+\cdots ,\quad \text{as}~|z|\rightarrow \infty ,\end{eqnarray}$$

where

(2.7ac ) $$\begin{eqnarray}a=\frac{{\rm\pi}}{2({\rm\pi}-{\it\theta})},\quad {\it\Delta}_{1}=-\frac{c^{2}}{3}(2a^{2}+1),\quad {\it\Delta}_{3}=-\frac{2c^{4}}{45}(-7a^{4}+5a^{2}+2).\end{eqnarray}$$

Actually, only ${\it\Delta}_{1}$ is needed to derive (1.2ac ) 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

(2.8a,b ) $$\begin{eqnarray}h_{0}(z)=h_{s}(z)-\frac{{\it\Delta}_{1}}{z}-{\it\lambda}_{0}\cot ({\it\delta}z),\quad -\frac{{\it\lambda}_{0}}{{\it\delta}}={\it\Delta}_{1}.\end{eqnarray}$$

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

(2.9) $$\begin{eqnarray}h(z+2l)=h(z)+2l\end{eqnarray}$$

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$ ,

(2.10) $$\begin{eqnarray}h_{0}(z)\sim z+\text{i}{\it\lambda}_{0}+O(1/z^{3}),\end{eqnarray}$$

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

(2.11) $$\begin{eqnarray}h_{0}(z+2l)=h_{0}(z)+2l+O({\it\delta}^{3}),\end{eqnarray}$$

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

(2.12) $$\begin{eqnarray}\displaystyle \text{Re}[h_{0}(z)] & = & \displaystyle \text{Re}\left[h_{s}(z)-\frac{{\it\Delta}_{1}}{z}-{\it\lambda}_{0}\cot ({\it\delta}z)\right]\nonumber\\ \displaystyle & {\sim} & \displaystyle \text{Re}\left[-\frac{{\it\Delta}_{1}}{z}-{\it\lambda}_{0}\left\{\frac{1}{{\it\delta}z}-\frac{{\it\delta}z}{3}+\cdots \right\}\right]=\text{Re}\left[\frac{{\it\lambda}_{0}{\it\delta}z}{3}+\cdots \right]\nonumber\\ \displaystyle & = & \displaystyle O({\it\lambda}_{0}c{\it\delta})=O(c^{2}{\it\delta}^{2}).\end{eqnarray}$$

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

(3.1) $$\begin{eqnarray}h_{1}(z)=h_{s}(z)+\frac{{\it\lambda}_{1}}{{\it\delta}}\frac{1}{z}-{\it\lambda}_{1}\cot ({\it\delta}z)+\frac{{\it\lambda}_{1}{\it\delta}}{3}[h_{s}(z)-z],\end{eqnarray}$$

with ${\it\lambda}_{1}$ now chosen to satisfy

(3.2) $$\begin{eqnarray}{\it\Delta}_{1}+\frac{{\it\lambda}_{1}}{{\it\delta}}+\frac{{\it\lambda}_{1}{\it\delta}{\it\Delta}_{1}}{3}=0.\end{eqnarray}$$

With this choice the $1/z$ term in the far-field expansion of the three non-cot terms

(3.3ac ) $$\begin{eqnarray}h_{s}(z),\quad +\frac{{\it\lambda}_{1}}{{\it\delta}}\frac{1}{z},\quad +\frac{{\it\lambda}_{1}{\it\delta}}{3}[h_{s}(z)-z]\end{eqnarray}$$

of $h_{1}(z)$ vanishes. This ensures that, as $z\rightarrow \infty$ ,

(3.4) $$\begin{eqnarray}h_{1}(z)\sim z+\text{i}{\it\lambda}_{1}+O(1/z^{3}),\end{eqnarray}$$

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

(3.5) $$\begin{eqnarray}h_{1}(z+2l)=h_{1}(z)+2l+O({\it\delta}^{3}),\end{eqnarray}$$

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,

(3.6) $$\begin{eqnarray}\displaystyle \text{Re}[h_{1}(z)] & = & \displaystyle \text{Re}\left[h_{s}(z)+\frac{{\it\lambda}_{1}}{{\it\delta}}\frac{1}{z}-{\it\lambda}_{1}\cot ({\it\delta}z)+\frac{{\it\lambda}_{1}{\it\delta}}{3}[h_{s}(z)-z]\right]\nonumber\\ \displaystyle & {\sim} & \displaystyle \text{Re}\left[+\frac{{\it\lambda}_{1}}{{\it\delta}}\frac{1}{z}-{\it\lambda}_{1}\left\{\frac{1}{{\it\delta}z}-\frac{{\it\delta}z}{3}+O(c^{3}{\it\delta}^{3})\right\}-\frac{{\it\lambda}_{1}{\it\delta}z}{3}\right]\nonumber\\ \displaystyle & = & \displaystyle O({\it\lambda}_{1}c^{3}{\it\delta}^{3})=O(c^{4}{\it\delta}^{4}),\end{eqnarray}$$

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

(3.7) $$\begin{eqnarray}{\it\lambda}_{TK}=\frac{{\it\lambda}_{1}}{2/l}=\frac{{\it\lambda}_{1}}{4{\it\delta}/{\rm\pi}},\end{eqnarray}$$

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 %.

Figure 2. Normalized slip length (3.7) as a function of protrusion angle ${\it\theta}$ for $c/l=0.1$ (a), 0.25 (b), 0.5 (c), 0.75 (d) and 0.9 (e). Crosses show numerical data from Teo & Khoo (Reference Teo and Khoo2010).

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

(4.1) $$\begin{eqnarray}\tilde{h}_{s}(z)\sim z^{3}+O(1/z),\quad \text{as}~z\rightarrow \infty .\end{eqnarray}$$

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

(4.2) $$\begin{eqnarray}\tilde{h}_{s}(z)=-a^{3}c^{3}\left[\frac{1}{{\it\zeta}^{3}}-{\it\zeta}^{3}\right]+c^{2}(1-a^{2})\left[\frac{1}{{\it\zeta}}-{\it\zeta}\right]\sim z^{3}+\frac{\tilde{{\it\Delta}}_{1}}{z}+\frac{\tilde{{\it\Delta}}_{3}}{z^{3}}+\cdots ,\quad \text{as}~|z|\rightarrow \infty ,\end{eqnarray}$$

where we have used (2.5) and

(4.3) $$\begin{eqnarray}\left.\begin{array}{@{}c@{}}\displaystyle \tilde{{\it\Delta}}_{1}=-\frac{2c^{4}}{15}(-7a^{4}+5a^{2}+2),\\ \displaystyle \tilde{{\it\Delta}}_{3}=-\frac{c^{6}}{3780}(4855a^{6}+945a^{5}-5817a^{4}+6615a^{3}-6090a^{2}+3780a-508).\end{array}\right\}\end{eqnarray}$$

Now consider the new complex potential, constructed using $h_{s}(z)$ and $\tilde{h}_{s}(z)$ , given by

(4.4) $$\begin{eqnarray}\displaystyle h_{2}(z) & = & \displaystyle h_{s}(z)+\frac{{\it\lambda}_{2}}{{\it\delta}}\frac{1}{z}+\frac{{\it\mu}_{2}}{{\it\delta}^{3}}\frac{1}{z^{3}}-{\it\lambda}_{2}\cot ({\it\delta}z)-{\it\mu}_{2}\cot ({\it\delta}z)\text{cosec}^{2}({\it\delta}z)\nonumber\\ \displaystyle & & \displaystyle +\left[\frac{1}{3}{\it\lambda}_{2}{\it\delta}+\frac{1}{15}{\it\mu}_{2}{\it\delta}\right](h_{s}(z)-z)+\left[\frac{1}{45}{\it\lambda}_{2}{\it\delta}^{3}+\frac{20}{945}{\it\mu}_{2}{\it\delta}^{3}\right](\tilde{h}_{s}(z)-z^{3}),\qquad\end{eqnarray}$$

where the real parameters ${\it\lambda}_{2}$ and ${\it\mu}_{2}$ are chosen to satisfy

(4.5) $$\begin{eqnarray}\left.\begin{array}{@{}c@{}}\displaystyle {\it\Delta}_{1}+\frac{{\it\lambda}_{2}}{{\it\delta}}+\left[\frac{1}{3}{\it\lambda}_{2}{\it\delta}+\frac{1}{15}{\it\mu}_{2}{\it\delta}\right]{\it\Delta}_{1}+\left[\frac{1}{45}{\it\lambda}_{2}{\it\delta}^{3}+\frac{20}{945}{\it\mu}_{2}{\it\delta}^{3}\right]\tilde{{\it\Delta}}_{1}=0,\\ \displaystyle {\it\Delta}_{3}+\frac{{\it\mu}_{2}}{{\it\delta}^{3}}+\left[\frac{1}{3}{\it\lambda}_{2}{\it\delta}+\frac{1}{15}{\it\mu}_{2}{\it\delta}\right]{\it\Delta}_{3}+\left[\frac{1}{45}{\it\lambda}_{2}{\it\delta}^{3}+\frac{20}{945}{\it\mu}_{2}{\it\delta}^{3}\right]\tilde{{\it\Delta}}_{3}=0.\end{array}\right\}\end{eqnarray}$$

We claim that $h_{2}(z)$ is the required higher-order solution. To see this, the expansions

(4.6a,b ) $$\begin{eqnarray}\cot z=\frac{1}{z}-\frac{z}{3}-\frac{z^{3}}{45}-\frac{2}{945}z^{5}+\cdots ,\quad \cot (z)\text{cosec}^{2}(z)=\frac{1}{z^{3}}-\frac{z}{15}-\frac{20}{945}z^{3}+\cdots ,\end{eqnarray}$$

reveal that the third-order singularity of $h_{2}(z)$ at $z=0$ is removable. Furthermore,

(4.7) $$\begin{eqnarray}h_{2}(z)\rightarrow z+\text{i}{\it\lambda}_{2}+O(1/z^{5}),\end{eqnarray}$$

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$ ,

(4.8) $$\begin{eqnarray}h_{2}(z+2l)=h_{2}(z)+2l+O({\it\delta}^{5})\end{eqnarray}$$

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),

(4.9) $$\begin{eqnarray}\displaystyle \text{Re}[h_{2}(z)] & = & \displaystyle \text{Re}\left[\frac{{\it\lambda}_{2}}{{\it\delta}}\frac{1}{z}+\frac{{\it\mu}_{2}}{{\it\delta}^{3}}\frac{1}{z^{3}}-{\it\lambda}_{2}\left[\frac{1}{{\it\delta}z}-\frac{{\it\delta}z}{3}-\frac{{\it\delta}^{3}z^{3}}{45}+O({\it\delta}^{5}c^{5})\right]\right.\nonumber\\ \displaystyle & & \displaystyle -\,{\it\mu}_{2}\left[\frac{1}{{\it\delta}^{3}z^{3}}-\frac{{\it\delta}z}{15}-\frac{20}{945}{\it\delta}^{3}z^{3}+O({\it\delta}^{5}c^{5})\right]-\left[\frac{1}{3}{\it\lambda}_{2}{\it\delta}+\frac{1}{15}{\it\mu}_{2}{\it\delta}\right]z\nonumber\\ \displaystyle & & \displaystyle -\left.\left[\frac{1}{45}{\it\lambda}_{2}{\it\delta}^{3}+\frac{20}{945}{\it\mu}_{2}{\it\delta}^{3}\right]z^{3}\right]=O({\it\lambda}_{1}c^{5}{\it\delta}^{5}),\qquad\end{eqnarray}$$

where we have used the Laurent expansions (4.6) and the fact that

(4.10) $$\begin{eqnarray}\text{Re}[h_{s}(z)]=\text{Re}[\tilde{h}_{s}(z)]=0,\quad \text{on the meniscus}.\end{eqnarray}$$

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

(4.11) $$\begin{eqnarray}\left(\begin{array}{@{}cc@{}}a_{11} & a_{12}\\ a_{21} & a_{22}\end{array}\right)\left(\begin{array}{@{}c@{}}{\it\lambda}_{2}\\ {\it\mu}_{2}\end{array}\right)=\left(\begin{array}{@{}c@{}}-{\it\Delta}_{1}{\it\delta}\\ -{\it\Delta}_{3}{\it\delta}^{3}\end{array}\right),\end{eqnarray}$$

where

(4.12) $$\begin{eqnarray}\left.\begin{array}{@{}c@{}}\displaystyle a_{11}=1+\frac{{\it\Delta}_{1}{\it\delta}^{2}}{3}+\frac{\tilde{{\it\Delta}}_{1}{\it\delta}^{4}}{45},\quad a_{12}=\frac{{\it\Delta}_{1}{\it\delta}^{2}}{15}+\frac{20\tilde{{\it\Delta}}_{1}{\it\delta}^{4}}{945},\\ \displaystyle a_{21}=\frac{{\it\Delta}_{3}{\it\delta}^{4}}{3}+\frac{\tilde{{\it\Delta}}_{3}{\it\delta}^{6}}{45},\quad a_{22}=1+\frac{{\it\Delta}_{3}{\it\delta}^{4}}{15}+\frac{20\tilde{{\it\Delta}}_{3}{\it\delta}^{6}}{945}.\end{array}\right\}\end{eqnarray}$$

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

(4.13) $$\begin{eqnarray}{\it\lambda}_{2}=\frac{{\it\Delta}_{3}{\it\delta}^{3}a_{12}-{\it\Delta}_{1}{\it\delta}a_{22}}{a_{11}a_{22}-a_{12}a_{21}}.\end{eqnarray}$$

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 %.

Figure 3. Normalized slip length, now based on the approximation ${\it\lambda}_{2}$ in (4.13), as a function of ${\it\theta}$ for $c/l=0.1$ (a), 0.25 (b), 0.5 (c), 0.75 (d) and 0.9 (e). Crosses show numerical data from Teo & Khoo (Reference Teo and Khoo2010).

Our analysis leads naturally to formula (4.13) but if we use the second equation in (4.5) to approximate

(4.14) $$\begin{eqnarray}{\it\mu}_{2}\approx -{\it\Delta}_{3}{\it\delta}^{3}\end{eqnarray}$$

then, on substitution of this into the first equation in (4.5), we find the approximation

(4.15) $$\begin{eqnarray}{\it\lambda}_{2}\approx -\frac{\displaystyle {\it\delta}{\it\Delta}_{1}\left(1-\frac{{\it\delta}^{4}{\it\Delta}_{3}}{15}\right)}{\displaystyle 1+\frac{{\it\Delta}_{1}{\it\delta}^{2}}{3}+\frac{{\it\delta}^{4}{\it\Delta}_{3}}{15}}.\end{eqnarray}$$

On eliminating $a$ in favour of ${\it\theta}$ , it can be shown that

(4.16) $$\begin{eqnarray}\left.\begin{array}{@{}c@{}}\displaystyle {\it\Delta}_{1}=-\frac{c^{2}}{6({\rm\pi}-{\it\theta})^{2}}(2{\it\theta}^{2}-4{\rm\pi}{\it\theta}+3{\rm\pi}^{2}),\\ \displaystyle {\it\Delta}_{3}=-\frac{c^{4}}{360({\rm\pi}-{\it\theta})^{4}}(32{\it\theta}^{4}-128{\rm\pi}{\it\theta}^{3}+212{\rm\pi}^{2}{\it\theta}^{2}-168{\rm\pi}^{3}{\it\theta}+45{\rm\pi}^{4}).\end{array}\right\}\end{eqnarray}$$

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,

(5.1) $$\begin{eqnarray}{\it\lambda}=\frac{2l}{{\rm\pi}}\log \sec \left[\frac{{\rm\pi}}{2l}\right]=-\frac{1}{{\it\delta}}\log \cos {\it\delta}\approx {\it\delta}\left[\frac{1}{2}+\frac{{\it\delta}^{2}}{12}+\frac{{\it\delta}^{4}}{45}+\cdots \right],\quad {\it\delta}\equiv \frac{{\rm\pi}}{2l}.\end{eqnarray}$$

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

(5.2) $$\begin{eqnarray}\frac{{\it\lambda}_{2}}{c}={\it\delta}c{\it\alpha}({\it\theta})\left[1+\frac{c^{2}{\it\delta}^{2}{\it\alpha}({\it\theta})}{3}+{\it\delta}^{4}c^{4}\left(\frac{2{\it\beta}({\it\theta})}{15}+\frac{{\it\alpha}({\it\theta})^{2}}{9}\right)+\cdots \right].\end{eqnarray}$$

A Taylor expansion of ${\it\alpha}({\it\theta})$ and ${\it\beta}({\it\theta})$ to linear order in the protrusion angle ${\it\theta}$ gives

(5.3a,b ) $$\begin{eqnarray}{\it\alpha}({\it\theta})=\frac{1}{2}+\frac{{\it\theta}}{3{\rm\pi}}+\cdots ,\quad {\it\beta}({\it\theta})=\frac{1}{8}+\frac{{\it\theta}}{30{\rm\pi}}+\cdots .\end{eqnarray}$$

On substitution of (5.3) into (5.2), we find

(5.4) $$\begin{eqnarray}\frac{{\it\lambda}_{2}}{c}={\it\delta}c\left[\frac{1}{2}+\frac{{\it\delta}^{2}c^{2}}{12}+\frac{{\it\delta}^{4}c^{4}}{45}+\cdots \right]+{\it\delta}c\left(\frac{{\it\theta}}{{\rm\pi}}\right)\left[\frac{1}{3}+\frac{c^{2}{\it\delta}^{2}}{9}+\frac{24}{675}c^{4}{\it\delta}^{4}+\cdots \right]+\cdots .\end{eqnarray}$$

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

(5.5) $$\begin{eqnarray}\frac{c{\it\theta}{\it\xi}}{2}\left[\frac{1}{3}+\frac{{\it\xi}^{2}{\rm\pi}^{2}}{36}+\frac{{\it\xi}^{4}{\rm\pi}^{4}}{450}\right]=-\frac{c^{2}}{2R}\left[\frac{{\it\xi}}{3}+\frac{{\it\xi}^{3}{\rm\pi}^{2}}{36}+\frac{{\it\xi}^{4}{\rm\pi}^{4}c^{4}}{450}\right],\end{eqnarray}$$

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

(6.1) $$\begin{eqnarray}\frac{\partial {\hat{w}}}{\partial n}=\mathscr{A}(z,\overline{z}),\end{eqnarray}$$

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

(6.2) $$\begin{eqnarray}0=\int _{D}[w_{2}{\rm\nabla}^{2}{\hat{w}}-{\hat{w}}{\rm\nabla}^{2}w_{2}]\,\text{d}A=\oint _{\partial D}w_{2}\frac{\partial {\hat{w}}}{\partial n}\,\text{d}s-\oint _{\partial D}{\hat{w}}\frac{\partial w_{2}}{\partial n}\,\text{d}s,\end{eqnarray}$$

where $\partial D$ is the boundary of $D$ . On the edges of the period window for which $|z|\geqslant l$ ,

(6.3a,b ) $$\begin{eqnarray}w_{2}=\dot{{\it\gamma}}(y+{\it\lambda}_{2})+O({\it\delta}^{5}),\quad \frac{\partial w_{2}}{\partial n}\,\text{d}s=\pm \frac{\partial w_{2}}{\partial x}\,\text{d}y,\end{eqnarray}$$

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

(6.4) $$\begin{eqnarray}\hat{{\it\lambda}}\approx {\it\lambda}_{2}-\frac{1}{2l}\int _{\mathit{meniscus}}\mathscr{A}(z,\overline{z})\,\text{Im}[h_{2}(z)]\,\text{d}s,\end{eqnarray}$$

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.

References

Asmolov, E. S. & Vinogradova, O. I. 2012 Effective slip boundary conditions for arbitrary one-dimensional surfaces. J. Fluid Mech. 706, 108117.Google Scholar
Bazant, M. Z. & Vinogradova, O. I. 2008 Tensorial hydrodynamic slip. J. Fluid Mech. 613, 125134.Google Scholar
Crowdy, D. G. 2010 Slip length for longitudinal shear flow over a dilute periodic mattress of protruding bubbles. Phys. Fluids 22, 121703.Google Scholar
Davis, A. M. J. & Lauga, E. 2009 Geometric transition in friction for flow over a bubble mattress. Phys. Fluids 21, 011701.Google Scholar
Ng, C.-O. & Wang, C. Y. 2011 Effective slip for Stokes flow over a surface patterned with two- or three-dimensional protrusions. Fluid Dyn. Res. 43, 065504.Google Scholar
Nizkaya, T. V., Asmolov, E. S. & Vinogradova, O. I. 2014 Gas cushion model and hydrodynamic boundary conditions for superhydrophobic textures. Phys. Rev. E 90, 043017.Google Scholar
Philip, J. R. 1972 Flows satisfying mixed no-slip and no-shear conditions. Z. Angew. Math. Phys. 23, 353372.Google Scholar
Rothstein, J. P. 2010 Slip on superhydrophobic surfaces. Annu. Rev. Fluid Mech. 42, 89109.Google Scholar
Sbragaglia, M. & Prosperetti, A. 2007 A note on the effective slip properties for microchannel flows with ultrahydrophobic surfaces. Phys. Fluids 19, 043603.CrossRefGoogle Scholar
Schönecker, C., Baier, T. & Hardt, S. 2014 Influence of the enclosed fluid on the flow over a microstructured surface in the Cassie state. J. Fluid Mech. 740, 168195.Google Scholar
Schönecker, C. & Hardt, S. 2013 Longitudinal and transverse flow over a cavity containing a second immiscible fluid. J. Fluid Mech. 717, 376394.Google Scholar
Steinberger, A., Cottin-Bizonne, C., Kleimann, P. & Charlaix, E. 2007 High friction on a bubble mattress. Nat. Mater. 6, 665668.Google Scholar
Teo, C. J. & Khoo, B. C. 2010 Flow past superhydrophobic surfaces containing longitudinal grooves: effects of interface curvature. Microfluid. Nanofluid. 9, 499511.Google Scholar
Figure 0

Figure 1. Longitudinal shear flow over a surface comprising a unidirectional periodic array, with period $2l$, of menisci of width $2c$ protruding at angle ${\it\theta}$ with intermediate no-slip zones.

Figure 1

Figure 2. Normalized slip length (3.7) as a function of protrusion angle ${\it\theta}$ for $c/l=0.1$ (a), 0.25 (b), 0.5 (c), 0.75 (d) and 0.9 (e). Crosses show numerical data from Teo & Khoo (2010).

Figure 2

Figure 3. Normalized slip length, now based on the approximation ${\it\lambda}_{2}$ in (4.13), as a function of ${\it\theta}$ for $c/l=0.1$ (a), 0.25 (b), 0.5 (c), 0.75 (d) and 0.9 (e). Crosses show numerical data from Teo & Khoo (2010).