Surface tension and wetting are dominating physical effects in microscale and nanoscale flows. We present an efficient and reliable model of surface tension and equilibrium contact angles in smoothed particle hydrodynamics for free-surface problems. We demonstrate its robustness and accuracy by simulating several three-dimensional free-surface flow problems driven by interfacial tension.

Impedance spectroscopy is widely adopted for probing the charge and charge mobility of soft ion-conducting media, such as synthetic membranes and biological tissue. The spectra exhibit a variety of distinctive signatures, but the physical basis of these is not well understood, e.g. models have not previously accounted for viscoelasticity, hydrodynamics or microstructural heterogeneity. This study explores a physically grounded continuum model that captures how these factors shape conductivity spectra. Nonlinear thermodynamics and linearised dynamics of a viscous electrolyte and compressible, elastic polymer network are coupled under the forcing of an oscillatory electric field. The model is solved in a one-dimensional spatially periodic unit cell, reporting conductivity and dielectric permittivity spectra, including Nyquist representations. Whereas rigid microstructures exhibit ion-diffusion-controlled relaxation, which manifests as a low-frequency dielectric ‘constant’, hydrodynamic and elastic forces contribute to a strongly diverging dielectric permittivity at low frequencies for heterogeneous anionic microstructures. The model also captures distinctive characteristics of experimentally reported impedance spectra for films bearing alternating layers of cationic and anionic charge, again highlighting the role of coupled hydrodynamic, elastic and electrical forces. Sufficiently thin and highly charged bilayers exhibit a notably low high-frequency conductivity. This is explained by strong low-frequency electrostatic polarisation and counter-ion release. The one-dimensional solutions computed herein provide a foundation for much more challenging computations in two and three dimensions.

Droplet spreading is ubiquitous and plays a significant role in liquid-based energy systems, thermal management devices and microfluidics. While the spreading of non-volatile droplets is quantitatively understood, the spreading and flow transition in volatile droplets remains elusive due to the complexity added by interfacial phase change and non-equilibrium thermal transport. Here we show, using both mathematical modelling and experiments, that the wetting dynamics of volatile droplets can be scaled by the spatial–temporal interplay between capillary, evaporation and thermal Marangoni effects. We elucidate and quantify these complex interactions using phase diagrams based on systematic theoretical and experimental investigations. A spreading law of evaporative droplets is derived by extending Tanner's law (valid for non-volatile liquids) to a full range of liquids with saturation vapour pressure spanning from $10^1$ to $10^4$ Pa and on substrates with thermal conductivity from $10^{-1}$ to $10^3\ {\rm W}\ {\rm m}^{-1}\ {\rm K}^{-1}$. In addition to its importance in fluid-based industries, the conclusions also enable a unifying explanation to a series of individual works including the criterion of flow reversal and the state of dynamic wetting, making it possible to control liquid transport in diverse application scenarios.

A theoretical analysis is presented of peristaltic pumping down a narrow conduit with permeable walls, motivated by the flushing action of lugworms and other marine organisms in sandy burrows. Flow in the conduit is dealt with using lubrication theory; the leakage into the surrounding medium is taken into account by exploiting slender-body theory to solve the associated Darcy problem. By adopting a model for the local force balance on the pumping surface, we bridge between the limits in which the pump operates with either fixed load or displacement. In the latter limit we characterize peristaltic waves with either fixed form or ones that partially collapse the conduit. We construct pump characteristics (the relation between the mean flux and net pressure drop) when the burrow wall is impermeable and pressures are fixed at each end, and compare the results with existing laboratory experiments performed on lugworms. We then consider how the peristaltic dynamics is changed when the wall is made permeable. Last, we consider pumping along an impermeable burrow into a leaky head shaft. The results reveal that the permeability of the conduit wall or end can greatly impact the direction and strength of the recirculating flow.

We consider the process of convective dissolution in a homogeneous and isotropic porous medium. The flow is unstable due to the presence of a solute that induces a density difference responsible for driving the flow. The mixing dynamics is thus driven by a Rayleigh–Taylor instability at the pore scale. We investigate the flow at the scale of the pores using Hele-Shaw type experiment with bead packs, two-dimensional direct numerical simulations and physical models. Experiments and simulations have been specifically designed to mimic the same flow conditions, namely matching porosities, high Schmidt numbers and linear dependency of fluid density with solute concentration. In addition, the solid obstacles of the medium are impermeable to fluid and solute. We characterise the evolution of the flow via the mixing length, which quantifies the extension of the mixing region and grows linearly in time. The flow structure, analysed via the centreline mean wavelength, is observed to grow in agreement with theoretical predictions. Finally, we analyse the dissolution dynamics of the system, quantified through the mean scalar dissipation, and three mixing regimes are observed. Initially, the evolution is controlled by diffusion, which produces solute mixing across the initial horizontal interface. Then, when the interfacial diffusive layer is sufficiently thick, it becomes unstable, forming finger-like structures and driving the system into a convection-dominated phase. Finally, when the fingers have grown sufficiently to touch the horizontal boundaries of the domain, the mixing reduces dramatically due to the absence of fresh unmixed fluid. With the aid of simple physical models, we explain the physics of the results obtained numerically and experimentally. The solute evolution presents a self-similar behaviour, and it is controlled by different length scales in each stage of the mixing process, namely the length scale of diffusion, the pore size and the domain height.

This paper reports on the efficacy of the Görtler number in scaling the laminar-turbulent boundary-layer transition on rotating cones facing axial inflow. Depending on the half-cone angle $\psi$ and axial flow strength, the competing centrifugal and cross-flow instabilities dominate the transition. Traditionally, the flow is evaluated by using two parameters: the local meridional Reynolds number $Re_l$ comparing the inertial versus viscous effects and the local rotational speed ratio $S$ accounting for the boundary-layer skew. We focus on the centrifugal effects, and evaluate the flow fields and reported transition points using Görtler number based on the azimuthal momentum thickness of the similarity solution and local cone radius. The results show that Görtler number alone dominates the late stages of transition (maximum amplification and turbulence onset phases) for a wide range of investigated $S$ and half-cone angle ($15^{\circ } \leq \psi \leq 50^{\circ }$), although the early stage (critical phase) seems to be not determined by the Görtler number alone on the broader cones ($\psi =30^{\circ }$ and $50^{\circ }$) where the primary cross-flow instability dominates the flow. Overall, this indicates that the centrifugal effects play an important role in the boundary-layer transition on rotating cones in axial inflow.

We investigate the three-dimensional (3-D) flow around and through a porous screen for various porosities at high Reynolds number $Re = {O}(10^4)$. Historically, the study of this problem has been focused on two-dimensional cases and for screens spanning completely or partially a channel. Since many recent problems have involved a porous object in a 3-D free flow, we present a 3-D model initially based on Koo & James (J. Fluid Mech., vol. 60, 1973, pp. 513–538) and Steiros & Hultmark (J. Fluid Mech., vol. 853, 2018 pp. 1–11) for screens of arbitrary shapes. In addition, we include an empirical viscous correction factor accounting for viscous effects in the vicinity of the screen. We characterize experimentally the aerodynamic drag coefficient for a porous square screen composed of fibres, immersed in a laminar air flow with various solidities and different angles of attack. We test various fibre diameters to explore the effect of the space between the pores on the drag force. Using PIV and hot wire probe measurements, we visualize the flow around and through the screen, and in particular measure the proportion of fluid that is deviated around the screen. The predictions from the model for drag coefficient, flow velocities and streamlines are in good agreement with our experimental results. In particular, we show that local viscous effects are important: at the same solidity and with the same air flow, the drag coefficient and the flow deviations strongly depend on the Reynolds number based on the fibre diameter. The model, taking into account 3-D effects and the shape of the porous screen, and including an empirical viscous correction factor that is valid for fibrous screens may have many applications including the prediction of water collection efficiency for fog harvesters.

The growth and characteristics of linear, oblique instabilities on a highly swept fin on a straight cone in Mach 6 flow are examined. Large streamwise pressure gradients cause doubly inflected cross-flow profiles and reversed flow near the wall, which necessitates using the harmonic linearized Navier–Stokes equations. The cross-flow instability is responsible for the most-amplified disturbances, however, not all disturbances show typical cross-flow characteristics. Distinct differences in perturbation structure are shown between small ($\sim$3–5 mm) and large ($\sim$10 mm) wavelength disturbances at the unit Reynolds number $Re' = 11 \times 10^6$ m$^{-1}$. As a result, amplification measurements based solely on wall quantities bias a most-amplified disturbance assessment towards larger wavelengths and lower frequencies than would otherwise be determined by an off-wall total-energy approach. A spatial-amplification energy-budget analysis demonstrates (i) that wall-normal Reynolds-flux terms dictate the local growth rate, despite other terms having a locally larger magnitude and (ii) that the Reynolds-stress terms are responsible for large-wavelength disturbances propagating closer to the wall compared with small-wavelength disturbances. Additionally, the effect of free-stream unit Reynolds number and small yaw angles on the perturbation amplification and energy budget is considered. At a higher Reynolds number ($Re' = 22 \times 10^6$ m$^{-1}$), the most-amplified wavelength shrinks. Perturbations do not behave self-similarly in the thinner boundary layer, and the shift in most-amplified wavelength is due to decreased dissipation relative to the lower-Reynolds-number case. Small yaw angles produce a streamwise shift in the boundary layer and disturbance amplification. The yaw results quantify a potential uncertainty source in experiments and flight.