Anisotropic flow in striped 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.

Directional Droplet Transport Mediated by Circular Groove Arrays. Part I: Experimental Findings

Directional transport of liquid droplets is crucial for various applications including water harvesting, anti-icing, and condensation heat transfer. Here, bouncing of water droplets with patterned superhydrophobic surfaces composed of circular equidistant grooves was studied. The directional transport of droplets toward the pole of the grooves was observed. The impact of the Weber number, initial polar distance r, and geometrical parameters of the surface on the directional droplet bouncing was experimentally explored. The nature of bouncing was switched when the Weber numbers exceeded We ≅ 20-25. The rebouncing height was slightly dependent on the initial polar coordinate of the impact point for a fixed We, whereas it grew for We > 20. The weak dependence of the droplet spreading time on the Weber number was close to the dependence predicted by the Hertz bouncing, thus evidencing the negligible influence of viscosity in the process. Change in the scaling exponent describing the dependence of the normalized spreading time on the Weber number was registered for We ≅ 25. The universal dependence of the offset distance ΔL of the droplets on the Weber number ΔL/D₀ ∼ We^1.5 was established. The normalized offset distance decreased with the normalized initial polar distance as ΔL/D₀ ∼ (r/S)^-1, where D₀ and S are the droplet diameter and groove width, respectively. This research may yield more insights into droplet bouncing on patterned surfaces and offer more options in directed droplet transportation.

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.

Hydrodynamic lift forces on solutes in a tilted nanopillar array: A computer simulation study

We study solute transport in a microfluidic channel, where the walls hold an array of tilted rigid nanopillars. By solving numerically the flow equations in the channel, we show that a combination of hydrodynamic effects with excluded volume interactions between the solute particles and the pillars leads to a hydrodynamic lift effect, which varies with the particle size, and depends in a strongly nonlinear fashion on the flow rate. We show that the lift force can be sufficiently strong to drive the solute accumulation or removal from the pillar region and can be switched to the opposite direction by variation of the shear rate or driving pressure. We also demonstrate that the nanopillar array can be used to selectively attract particles of certain size and enhance solute trapping at the surface.

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.

Electrohydrodynamics near hydrophobic surfaces

We show that an electro-osmotic flow near the slippery hydrophobic surface depends strongly on the mobility of surface charges, which are balanced by counterions of the electrostatic diffuse layer. For a hydrophobic surface with immobile charges, the fluid transport is considerably amplified by the existence of a hydrodynamic slippage. In contrast, near the hydrophobic surface with mobile adsorbed charges, it is also controlled by an additional electric force, which increases the shear stress at the slipping interface. To account for this, we formulate electrohydrodynamic boundary conditions at the slipping interface, which should be applied to quantify electro-osmotic flows instead of hydrodynamic boundary conditions. Our theoretical predictions are fully supported by dissipative particle dynamics simulations with explicit charges. These results lead to a new interpretation of zeta potential of hydrophobic surfaces.

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.

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.

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.