Tensorial slip of superhydrophobic channels

Modeling of Shear Flows over Superhydrophobic Surfaces: From Newtonian to Non-Newtonian Fluids

The design and use of superhydrophobic surfaces have gained special attentions due to their superior performances and advantages in many flow systems, e.g., in achieving specific goals including drag reduction and flow/droplet handling and manipulation. In this work, we conduct a brief review of shear flows over superhydrophobic surfaces, covering the classic and recent studies/trends for both Newtonian and non-Newtonian fluids. The aim is to mainly review the relevant mathematical and numerical modeling approaches developed during the past 20 years. Considering the wide ranges of applications of superhydrophobic surfaces in Newtonian fluid flows, we attempt to show how the developed studies for the Newtonian shear flows over superhydrophobic surfaces have been evolved, through highlighting the major breakthroughs. Despite the fact that, in many practical applications, flows over superhydrophobic surfaces may show complex non-Newtonian rheology, interactions between the non-Newtonian rheology and superhydrophobicity have not yet been well understood. Therefore, in this Review, we also highlight emerging recent studies addressing the shear flows of shear-thinning and yield stress fluids in superhydrophobic channels. We focus on reviewing the models developed to handle the intricate interaction between the formed liquid/air interface on superhydrophobic surfaces and the overlying flow. Such an intricate interaction will be more complex when the overlying flow shows nonlinear non-Newtonian rheology. We conclude that, although our understanding on the Newtonian shear flows over superhydrophobic surfaces has been well expanded via analyzing various aspects of such flows, the non-Newtonian counterpart is in its early stages. This could be associated with either the early applications mainly concerning Newtonian fluids or new complexities added to an already complex problem by the nonlinear non-Newtonian rheology. Finally, we discuss the possible directions for development of models that can address complex non-Newtonian shear flows over superhydrophobic surfaces.

Flow Structure and Mixing Efficiency of Viscous Fluids in Microchannel with a Striped Superhydrophobic Wall

The peculiarities of a Newtonian fluid flow structure in microchannels with a striped superhydrophobic lower wall texture are studied by means of numerical modeling. In the Cassie-Baxter state, an oblique orientation of such a texture induces helicoidal streamlines with micro spirals. Such a flow structure favors the enhancement of fluid mixing efficiency, which can be quantified using the total root-mean-square deviation of streamlines from the microchannel axis. This characteristic was shown to be a nonmonotonic function of the striped texture tilt angle and to depend strongly on microchannel thickness. The mechanisms of micro and macro helicoidal flow structure formation are investigated, and the mixing quality of miscible fluids is estimated for various Peclet numbers and texture tilt angles. It was found that the striped superhydrophobic wall leads to a notable enhancement in the microchannel mixing efficiency at sufficiently large Peclet numbers.

Apparent slip of shear thinning fluid in a microchannel with a superhydrophobic wall

The peculiarities of simple shear flow of shear thinning fluids over a superhydrophobic wall consisting of a set of parallel gas-filled grooves and solid stripes (domains with slip and stick boundary conditions) are studied numerically. The Carreau-Yasuda model is used to provide further insight into the problem of the slip behavior of non-Newtonian fluids having a decreasing viscosity with a shear rate increase. This feature is demonstrated to cause a nonlinear velocity profile leading to the apparent slip. The corresponding transverse and longitudinal apparent slip lengths of a striped texture are found to be noticeably larger than the respective effective slip lengths of Newtonian liquids in microchannels of various thicknesses and surface fractions of the slip domains. The viscosity distribution of the shear thinning fluid over the superhydrophobic wall is carefully investigated to describe the mechanism of the apparent slip. Nonmonotonic behavior of the apparent slip length as a function of the applied shear rate is revealed. This important property of shear thinning fluids is considered to be sensitive to the steepness of the viscosity flow curve, thus providing a way to decrease considerably the flow resistance in microchannels.

Advective superdiffusion in superhydrophobic microchannels

We consider pressure-driven flows in wide microchannels, and discuss how a transverse shear, generated by misaligned superhydrophobic walls, impacts cross-sectional spreading of Brownian particles. We show that such a transverse shear can induce an advective superdiffusion, which strongly enhances dispersion of particles compared to a normal diffusion, and that maximal cross-sectional spreading corresponds to a crossover between its subballistic and superballistic regimes. This allows us to argue that an advective superdiffusion can be used for boosting dispersion of particles at smaller Péclet numbers compared to known concepts of passive microfluidic mixing. This implies that our superdiffusion scenario allows one efficient mixing of much smaller particles or using much thinner microchannels than methods, which are currently being exploited.

Study of the Relationship between Boundary Slip and Nanobubbles on a Smooth Hydrophobic Surface

Surface nanobubbles, which are nanoscopic or microscopic gaseous domains forming at the solid/liquid interface, have a strong impact on the interface by changing the two-phase contact to a three-phase contact. Therefore, they are believed to affect the boundary condition and liquid flow. However, there are still disputes in the theoretical studies as to whether the nanobubbles can increase the slip length effectively. Furthermore, there are still no direct experimental studies to support either side. Therefore, an intensive study on the effective slip length for flows over bare surfaces with nanobubbles is essential for establishing the relation between nanobubbles and slip length. Here, we study the effect of nanobubbles on the slippage experimentally and theoretically. Our experimental results reveal an increase from 8 to 512 nm in slip length by increasing the surface coverage of nanobubbles from 1.7 to 50.8% and by decreasing the contact angle of nanobubbles from 42.8 to 16.6°. This is in good agreement with theoretical results. Our results indicate that nanobubbles could always act as a lubricant and significantly increase the slip length. The surface coverage, height, and contact angle are key factors for nanobubbles to reduce wall friction.

Probing effective slippage on superhydrophobic stripes by atomic force microscopy

While the effective slippage of water past superhydrophobic surfaces has been studied over a decade, theoretical predictions have never been properly confirmed by experiments. Here we measure a drag force on a sphere approaching a plane decorated by superhydrophobic grooves and compare the results with the predictions of semi-analytical theory developed here, which employs the gas cushion model to calculate the local slip length at the gas sectors. We demonstrate that at intermediate and large (compared to a texture period) separations the half-sum of longitudinal and transverse effective slip lengths can be deduced from the force-distance curve by using the known analytical theory of hydrodynamic interaction of a sphere with a homogeneous slipping plane. This half-sum is shown to depend on the fraction of gas sectors and its value is in excellent agreement with theoretical predictions. At small distances the half-sum of effective longitudinal and transverse slip lengths becomes separation-dependent, and is in quantitative agreement with the predictions of our semi-analytical theory.

Principles of transverse flow fractionation of microparticles in superhydrophobic channels

We propose a concept of fractionation of micron-sized particles in a microfluidic device with a bottom wall decorated by superhydrophobic stripes. The stripes are oriented at an angle α to the direction of a driving force, G, which generally includes an applied pressure gradient and gravity. Separation relies on the initial sedimentation of particles under gravity in the main forward flow, and their subsequent lateral deflection near a superhydrophobic wall due to generation of a secondary flow transverse to G. We provide some theoretical arguments allowing us to quantify the transverse displacement of particles in the microfluidic channel, and confirm the validity of theoretical predictions in test experiments with monodisperse fractions of microparticles. Our results can guide the design of superhydrophobic microfluidic devices for efficient sorting of microparticles with a relatively small difference in size and density.

Flows and mixing in channels with misaligned superhydrophobic walls

Aligned superhydrophobic surfaces with the same texture orientation reduce drag in the channel and generate secondary flows transverse to the direction of the applied pressure gradient. Here we show that a transverse shear can be easily generated by using superhydrophobic channels with misaligned textured surfaces. We propose a general theoretical approach to quantify this transverse flow by introducing the concept of an effective shear tensor. To illustrate its use, we present approximate theoretical solutions and Dissipative Particle Dynamics simulations for striped superhydrophobic channels. Our results demonstrate that the transverse shear leads to complex flow patterns, which provide a new mechanism of a passive vertical mixing at the scale of a texture period. Depending on the value of Reynolds number two different scenarios occur. At relatively low Reynolds number the flow represents a transverse shear superimposed with two corotating vortices. For larger Reynolds number these vortices become isolated, by suppressing fluid transport in the transverse direction.

Green-Kubo relation for friction at liquid-solid interfaces

We have developed a Green-Kubo relation that enables accurate calculations of friction at solid-liquid interfaces directly from equilibrium molecular dynamics (MD) simulations and that provides a pathway to bypass the time-scale limitations of typical nonequilibrium MD simulations. The theory has been validated for a number of different interfaces and it is demonstrated that the liquid-solid slip is an intrinsic property of an interface. Because of the high numerical efficiency of our method, it can be used in the design of interfaces for applications in aqueous environments, such as nano- and microfluidics.

Electrokinetics over charge-modulated surfaces in the presence of patterned wettability: role of the anisotropic streaming potential

In the present study, we focus on evaluating the induced streaming electric field along the orthogonal directions in a narrow fluidic confinement in the presence of patterned surface wettability and modulated surface charges. We attempt to assess the implications of such modulations on the related important quantities and pinpoint the regimes of improved induced streaming potential field and the resulting anisotropy in the induced potential. Our results reveal that for certain combinations of the parameters characterizing the modulated slip, a significant amount of augmentation in the streaming electric field might be obtained, whereas in other cases the effects may lead to adverse consequences. We further demonstrate that the presence of anisotropic modulations on the channel walls give rise to considerable off-diagonal effects, which makes the streaming potential "disoriented" with the applied pressure gradient, when the same is not applied along one of the orthogonal directions. Our analysis also shows that one can remove such "mis-orientations" by finely tuning several relevant flow and geometric parameters, which may bear immense scientific and technological consequences towards an improved design of miniaturized energy conversion devices.

Effective slip-length tensor for a flow over weakly slipping stripes

We discuss the flow past a flat heterogeneous solid surface decorated by slipping stripes. The spatially varying slip length, b(y), is assumed to be small compared to the scale of the heterogeneities, L, but finite. For such weakly slipping surfaces, earlier analyses have predicted that the effective slip length is simply given by the surface-averaged slip length, which implies that the effective slip-length tensor becomes isotropic. Here we show that a different scenario is expected if the local slip length has steplike jumps at the edges of slipping heterogeneities. In this case, the next-to-leading term in an expansion of the effective slip-length tensor in powers of max[b(y)/L] becomes comparable to the leading-order term, but anisotropic, even at very small b(y)/L. This leads to an anisotropy of the effective slip and to its significant reduction compared to the surface-averaged value. The asymptotic formulas are tested by numerical solutions and are in agreement with results of dissipative particle dynamics simulations.

Interfacial instability of compressible slip flows in a microchannel

In microfluidics the instability of an interface between two fluids may be favorable or unfavorable. Such an instability for incompressible flows has been extensively studied in the literature; in this paper we extend it to compressible flows which slip on the solid boundary. A generalized Darcy's law taking into consideration the compressibility and velocity slip is obtained. Using a conformal-mapping method we derive an interface equation on a fixed domain, the solutions of which determine the moving interface. We also examine the linear stability of the base-state flow by perturbing the corresponding conformal map and show that the velocity slip has a stabilizing effect.

Control of slippage with tunable bubble mattresses

Tailoring the hydrodynamic boundary condition is essential for both applied and fundamental aspects of drag reduction. Hydrodynamic friction on superhydrophobic substrates providing gas-liquid interfaces can potentially be optimized by controlling the interface geometry. Therefore, establishing stable and optimal interfaces is crucial but rather challenging. Here we present unique superhydrophobic microfluidic devices that allow the presence of stable and controllable microbubbles at the boundary of microchannels. We experimentally and numerically examine the effect of microbubble geometry on the slippage at high resolution. The effective slip length is obtained for a wide range of protrusion angles, θ, of the microbubbles into the flow, using a microparticle image velocimetry technique. Our numerical results reveal a maximum effective slip length, corresponding to a 23% drag reduction at an optimal θ ≈ 10°. In agreement with the simulation results, our measurements correspond to up to 21% drag reduction when θ is in the range of -2° to 12°. The experimental and numerical results reveal a decrease in slip length with increasing protrusion angles when >/~ 10°. Such microfluidic devices with tunable slippage are essential for the amplified interfacial transport of fluids and particles.

Tensorial slip theory for gas flows and comparison with molecular dynamics simulations using an anisotropic gas-wall collision mechanism

In this paper we examine the anisotropic slip theory for gas flows based on tangential accommodation coefficients and compare it with molecular dynamics (MD) results. A special gas-wall boundary condition is employed within MD simulations to mimic the anisotropic gas-wall collision mechanism. Results from MD simulations with different surface orientations show good agreement with the slip quantification proposed in this work.

Flow past superhydrophobic surfaces with cosine variation in local slip length

Anisotropic superhydrophobic surfaces have the potential to greatly reduce drag and enhance mixing phenomena in microfluidic devices. Recent work has focused mostly on cases of superhydrophobic stripes. Here, we analyze a relevant situation of cosine variation of the local slip length. We derive approximate formulas for maximal (longitudinal) and minimal (transverse) directional effective slip lengths that are in good agreement with the exact numerical solution and lattice-Boltzmann simulations. Compared to the case of superhydrophobic stripes, the cosine texture can provide a very large effective slip. However, the difference between eigenvalues of the slip-length tensor is smaller, indicating that the flow is less anisotropic.