An official website of the United States government
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
This content will become publicly available on December 10, 2025
Unsteady evolution of slip and drag in surfactant-contaminated superhydrophobic channels
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
- Award ID(s):
- 2054894
- PAR ID:
- 10598757
- Publisher / Repository:
- Journal of Fluid Mechanics
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 1000
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
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
-
Drag for wall-bounded flows is directly related to the spatial flux of spanwise vorticity outward from the wall. In turbulent flows a key contribution to this wall-normal flux arises from nonlinear advection and stretching of vorticity, interpretable as a cascade. We study this process using numerical simulation data of turbulent channel flow at friction Reynolds number$$Re_\tau =1000$$. The net transfer from the wall of spanwise vorticity created by downstream pressure drop is due to two large opposing fluxes, one which is ‘down-gradient’ or outward from the wall, where most vorticity concentrates, and the other which is ‘up-gradient’ or toward the wall and acting against strong viscous diffusion in the near-wall region. We present evidence that the up-gradient/down-gradient transport occurs by a mechanism of correlated inflow/outflow and spanwise vortex stretching/contraction that was proposed by Lighthill. This mechanism is essentially Lagrangian, but we explicate its relation to the Eulerian anti-symmetric vorticity flux tensor. As evidence for the mechanism, we study (i) statistical correlations of the wall-normal velocity and of wall-normal flux of spanwise vorticity, (ii) vorticity flux cospectra identifying eddies involved in nonlinear vorticity transport in the two opposing directions and (iii) visualizations of coherent vortex structures which contribute to the transport. The ‘D-type’ vortices contributing to down-gradient transport in the log layer are found to be attached, hairpin-type vortices. However, the ‘U-type’ vortices contributing to up-gradient transport are detached, wall-parallel, pancake-shaped vortices with strong spanwise vorticity, as expected by Lighthill's mechanism. We discuss modifications to the attached eddy model and implications for turbulent drag reduction.more » « less
-
Recent simulations indicate that streamwise-preferential porous materials have the potential to reduce drag in wall-bounded turbulent flows (Gómez-de-Segura & García-Mayoral, J. Fluid Mech. , vol. 875, 2019, pp. 124–172). This paper extends the resolvent formulation to study the effect of such anisotropic permeable substrates on turbulent channel flow. Under the resolvent formulation, the Fourier-transformed Navier–Stokes equations are interpreted as a linear forcing–response system. The nonlinear terms are considered the endogenous forcing in the system that gives rise to a velocity and pressure response. A gain-based decomposition of the forcing–response transfer function – the resolvent operator – identifies response modes (resolvent modes) that are known to reproduce important structural and statistical features of wall-bounded turbulent flows. The effect of permeable substrates is introduced in this framework using the volume-averaged Navier–Stokes equations and a generalized form of Darcy's law. Substrates with high streamwise permeability and low spanwise permeability are found to suppress the forcing–response gain for the resolvent mode that serves as a surrogate for the energetic near-wall cycle. This reduction in mode gain is shown to be consistent with the drag reduction trends predicted by theory and observed in numerical simulations. Simulation results indicate that drag reduction is limited by the emergence of spanwise rollers resembling Kelvin–Helmholtz vortices beyond a threshold value of wall-normal permeability. The resolvent framework also predicts the conditions in which such energetic spanwise-coherent rollers emerge. These findings suggest that a limited set of resolvent modes can serve as the building blocks for computationally efficient models that enable the design and optimization of permeable substrates for passive turbulence control.more » « less
