skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on May 10, 2026

Title: Drag reduction in surfactant-contaminated superhydrophobic channels at high Péclet numbers
Motivated by microfluidic applications, we investigate drag reduction in laminar pressure-driven flows in channels with streamwise-periodic superhydrophobic surfaces (SHSs) contaminated with soluble surfactant. We develop a model in the long-wave and weak-diffusion limit, where the streamwise SHS period is large compared with the channel height and the Péclet number is large. Using asymptotic and numerical techniques, we determine the influence of surfactant on drag reduction in terms of the relative strength of advection, diffusion, Marangoni effects and bulk–surface exchange. In scenarios with strong exchange, the drag reduction exhibits a complex dependence on the thickness of the bulk-concentration boundary layer and surfactant strength. Strong Marangoni effects immobilise the interface through a linear surfactant distribution, whereas weak Marangoni effects yield a quasi-stagnant cap. The quasi-stagnant cap has an intricate structure with an upstream slip region, followed by intermediate inner regions and a quasi-stagnant region that is mediated by weak bulk diffusion. The quasi-stagnant region differs from the immobile region of a classical stagnant cap, observed for instance in surfactant-laden air bubbles in water, by displaying weak slip. As exchange weakens, the bulk and interface decouple: the surfactant distribution is linear when the surfactant is strong, whilst it forms a classical stagnant cap when the surfactant is weak. The asymptotic solutions offer closed-form predictions of drag reduction across much of the parameter space, providing practical utility and enhancing understanding of surfactant dynamics in flows over SHSs.  more » « less
Award ID(s):
2054894
PAR ID:
10598756
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Journal of Fluid Mechanics
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
1010
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Although superhydrophobic surfaces (SHSs) show promise for drag reduction applications, their performance can be compromised by traces of surfactant, which generate Marangoni stresses that increase drag. This question is addressed for soluble surfactant in a three-dimensional laminar channel flow, with periodic SHSs made of long finite-length longitudinal grooves located on both walls. We assume that bulk diffusion is sufficiently strong for cross-channel concentration gradients to be small. Exploiting long-wave theory and accounting for the difference between the rapid transverse and slower longitudinal Marangoni flows, we derive a one-dimensional model for surfactant transport from the full three-dimensional transport equations. Our one-dimensional model allows us to predict the drag reduction and surfactant distribution across the parameter space. The system exhibits multiple regimes, involving competition between Marangoni effects, bulk and interfacial diffusion, bulk and interfacial advection, shear dispersion and surfactant exchange between the bulk and the interface. We map out asymptotic regions in the high-dimensional parameter space, and derive explicit closed-form approximations of the drag reduction, without any fitting or empirical parameters. The physics underpinning the drag reduction effect and the negative effect of surfactant is discussed through analysis of the velocity field and surfactant concentrations, which show both uniform and non-uniform stress distributions. Our theoretical predictions of the drag reduction compare well with results from the literature solving numerically the full three-dimensional transport problem. Our atlas of maps provides a comprehensive analytical guide for designing surfactant-contaminated channels with SHSs, to maximise the drag reduction in applications. 
    more » « less
  2. Recognising that surfactants can impede the drag reduction resulting from superhydrophobic surfaces (SHS), we investigate the impact of spatio–temporal fluctuations in surfactant concentration on the drag-reduction properties of SHS. We model the unsteady transport of soluble surfactant in a channel flow bounded by two SHS. The flow is laminar, pressure driven and the SHS are periodic in the streamwise and spanwise directions. We assume that the channel length is much longer than the streamwise period, the streamwise period is much longer than the channel height and spanwise period, and bulk diffusion is sufficiently strong for cross-channel concentration gradients to be small. By combining long-wave and homogenisation theories, we derive an unsteady advection–diffusion equation for surfactant-flux transport over the length of the channel, which is coupled to a quasi-steady advection–diffusion equation for surfactant transport over individual plastrons. As diffusion over the length of the channel is typically small, the surfactant flux is governed by a nonlinear wave equation. In the fundamental case of the transport of a bolus of surfactant, we predict its propagation speed and describe its nonlinear evolution via interaction with the SHS. The propagation speed can fall below the average streamwise velocity as the surfactant adsorbs and rigidifies the plastrons. Smaller concentrations of surfactant are advected faster than larger ones, so that wave-steepening effects can lead to shock formation in the surfactant-flux distribution. Our asymptotic results reveal how unsteady surfactant transport can affect the spatio–temporal evolution of the slip velocity, drag reduction and effective slip length in SHS channels. 
    more » « less
  3. Recent experimental and computational investigations have shown that trace amounts of surfactants, unavoidable in practice, can critically impair the drag reduction of superhydrophobic surfaces (SHSs), by inducing Marangoni stresses at the air–liquid interface. However, predictive models for realistic SHS geometries do not yet exist, which has limited the understanding and mitigation of these adverse surfactant effects. To address this issue, we derive a model for laminar, three-dimensional flow over SHS gratings as a function of geometry and soluble surfactant properties, which together encompass 10 dimensionless groups. We establish that the grating lengthgis the key geometric parameter and predict that the ratio between actual and surfactant-free slip increases withg2. Guided by our model, we perform synergistic numerical simulations and microfluidic experiments, finding good agreement with the theory as we vary surfactant type and SHS geometry. Our model also enables the estimation, based on velocity measurements, of a priori unknown properties of surfactants inherently present in microfluidic systems. For SHSs, we show that surfactant effects can be predicted by a single parameter, representing the ratio between the grating length and the interface length scale beyond which the flow mobilizes the air–water interface. This mobilization length is more sensitive to the surfactant chemistry than to its concentration, such that even trace-level contaminants may significantly increase drag if they are highly surface active. These findings advance the fundamental understanding of realistic interfacial flows and provide practical strategies to maximize superhydrophobic drag reduction. 
    more » « less
  4. In this work, we experimentally investigated the impact of surface roughness on drag reduction as well as the plastron stability of superhydrophobic surfaces (SHSs) in turbulent flows. A series of SHSs were fabricated by spraying hydrophobic nanoparticles on sandpapers. By changing the grit size of sandpapers from 240 to 1500, the root mean square roughness height (krms) of the SHSs varied from 4 to 14 μm. The experiments were performed in a turbulent channel flow facility, where the mean flow speed (Um) varied from 0.5 to 4.4 m/s, and the Reynolds number (Rem) based on Um and channel height changed from 3400 to 26 400. The drag reduction by SHSs was measured based on pressure drops in the fully developed flow region. The plastron status and gas fraction (φg) were simultaneously monitored by reflected-light microscopy. Our results showed a strong correlation between drag reduction and krms+ = krms/δv, where δv is the viscous length scale. For krms+ < 1, drag reduction was independent of krms+. A maximum 47% drag reduction was observed. For 1 < krms+ < 2, less drag reduction was observed due to the roughness effect. And for krms+ > 2, the SHSs caused an increase in drag. Furthermore, we found that surface roughness influenced the trend of plastron depletion in turbulent flows. As increasing Rem, φg reduced gradually for SHSs with large krms, but reduced rapidly and maintained as a constant for SHSs with small krms. Finally, we found that as increasing Rem, the slip length of SHS reduced, although φg was nearly a constant. 
    more » « less
  5. Marangoni flow is the motion induced by a surface tension gradient along a fluid–fluid interface. In this study, we report a Marangoni flow generated when a bath of surfactant contacts a pre-wetted film of deionized water on a vertical substrate. The thickness profile of the pre-wetted film is set by gravitational drainage and so varies with the drainage time. The surface tension is lower in the bath due to the surfactant, and thus a liquid film climbs upwards along the vertical substrate due to the surface tension difference. Particle tracking velocimetry is performed to measure the dynamics in the film, where the mean fluid velocity reverses direction as the draining film encounters the front of the climbing film. The effect of the surfactant concentration and the pre-wetted film thickness on the film climbing is then studied. High-speed interferometry is used to measure the front position of the climbing film and the film thickness profile. As a result, higher surfactant concentration induces a faster and thicker climbing film. Also, for high surfactant concentrations, where Marangoni driving dominates, increasing the film thickness increases the rise speed of the climbing front, since viscous resistance is less important. In contrast, for low surfactant concentrations, where Marangoni driving balances gravitational drainage, increasing the film thickness decreases the rise speed of the climbing front while enhancing gravitational drainage. We rationalize these observations by utilizing a dimensionless parameter that compares the magnitudes of the Marangoni stress and gravitational drainage. A model is established to analyse the climbing front, either in the Marangoni-driving-dominated region or in the Marangoni-balanced drainage region. Our work highlights the effects of the gravitational drainage on the Marangoni flow, both by setting the thickness of a pre-wetted film and by resisting the film climbing. 
    more » « less