The ability of streamwise-travelling waves of spanwise velocity to reduce the turbulent skin-friction drag is assessed in the compressible regime. Direct numerical simulations are carried out to compare drag reduction in subsonic, transonic and supersonic channel flows. Compressibility improves the benefits of the travelling waves, in a way that depends on the control parameters: drag reduction becomes larger than the incompressible one for small frequencies and wavenumbers. However, the improvement depends on the specific procedure employed for comparison. When the Mach number is varied and, at the same time, wall friction is changed by the control, the bulk temperature in the flow can either evolve freely in time until the aerodynamic heating balances the heat flux at the walls, or be constrained such that a fixed percentage of kinetic energy is transformed into thermal energy. Physical arguments suggest that, in the present context, the latter approach should be preferred. This provides a test condition in which the wall-normal temperature profile more realistically mimics that in an external flow, and also leads to a much better scaling of the results, over both the Mach number and the control parameters. Under this comparison, drag reduction is only marginally improved by compressibility.

We use video footage of a water-tunnel experiment to construct a 2-D reduced-order model of the flapping dynamics of an inverted flag in uniform flow. The model is obtained as the reduced dynamics on a 2-D attracting spectral submanifold (SSM) that emanates from the two slowest modes of the unstable fixed point of the flag. Beyond an unstable fixed point and a limit cycle expected from observations, our SSM-reduced model also confirms the existence of two unstable fixed points for the flag, which were found by previous studies. Importantly, the model correctly reconstructs the dynamics from a small number of general trajectories and no further information on the system. In the chaotic flapping regime, we construct a 4-D SSM-reduced model that captures the system's chaotic attractor.

The concept of added mass is generalized to stratified fluids, accounting for the presence of internal waves. Once the added mass of a moving body is known, so is the hydrodynamic force exerted on it by the fluid, and the energy imparted by it to the fluid. As a function of frequency, added mass is complex. Its real part is associated with inertia and its imaginary part, only present in the frequency range of propagating waves, with wave damping. Owing to causality, these two parts satisfy Kramers–Kronig relations. The added masses of an elliptic cylinder of horizontal axis, typical of two-dimensional bodies, and a spheroid of vertical axis, typical of three-dimensional bodies, are deduced from their dipole strengths, themselves deduced from their representations as single layers. The wave power is shown to be a maximum, for fixed oscillation amplitude, at approximately $0.8$ times the buoyancy frequency. In the temporal domain, added mass appears as a new memory force taking the form of a convolution integral. The kernel of this integral combines algebraically decaying oscillations at the buoyancy frequency on the one hand; and an exponentially damped oscillation for the horizontal motion of the spheroid, implying short-term memory, an aperiodic algebraic decay for its vertical motion, implying long-term memory, and a constant for the motion of the cylinder, implying everlasting memory, on the other hand. A limitation of the study is its restriction to translational motion.

The passage of single air bubbles through the horizontal interface between viscoplastic and Newtonian fluids, considering various combinations of densities and viscosities for the fluid layers, is studied computationally. The primary focus is on the quantity of liquid transferred from the lower layer (viscoplastic fluid) to the upper layer (Newtonian fluid) as a result of the bubble's ascent, a factor with significant implications for the turbidity of methane-emitting lakes and water bodies. The entrainment characteristics are observed to vary considerably with the bubble shape, within the lower layer and as the bubble approaches the interface. The results show that at Bond number $(Bo)>1$ and moderate Archimedes ($Ar$), prolate-shaped bubbles crossing the interface undergo elongation in the direction of their poles. This elongation is further accentuated when the viscosity of upper layer is less than the plastic viscosity of the lower layer. The bubble is found to break up when leaving the lower layer, of a critical capillary number, $Ca_c \approx 5$. The results show a significant reduction in the volume of entrainment compared with the Newtonian counterpart. This suggests disturbances caused by the rising bubble at the interface dissipate over a smaller region. Four distinct entrainment regimes are identified, mainly indicating the height to which the entrained fluid can be transported away from the interface. In contrast to Newtonian fluids, the volume of entrainment increases by decreasing the viscosity of the upper layer. Interestingly, the heavy viscoplastic fluid that has been dragged up into the light Newtonian fluid does not recede with time.

New experiments in highly turbulent, steady, subcritical and uniform water open-channel flows have been carried out to measure the mean turbulent kinetic energy (TKE) budget of sediment-laden boundary layer flows with two sizes (dp = 3 mm and 1 mm) of Plexiglas particles $({\rm relative\ density}\ = 1.192)$. The experiments covered energetic sediment transport conditions (Shields number of $0.35 < \theta < 1.2$) ranging from non-capacity to full-capacity flows in bedload-to-suspension-dominated transport modes (suspension number of $0.5 < {w_s}/{u_\ast } < 1.3$ where ${w_s}$ is the settling velocity and $u_*$ is the friction velocity) and for weakly to highly inertial, finite size turbulence-particle conditions (Stokes number of 0.1 < St < 3.5 and dp/η > 10 where $\eta$ is the Kolmogorov length scale). It was shown that the effects of sediments on the TKE budget are very pronounced in all large particle experiments for which a bedload layer of several grain diameter thickness is developed above the channel bed. When compared with the corresponding reference clear-water flows, the TKE shear-production rate for the 3 mm particle flows is strongly reduced in the wall region corresponding to the bedload layer. This turbulence damping is seen to increase with sediment load until full capacity for flows with constant Shields value, as well as with Shields number value. Inside this damped TKE shear-production zone, a distinct peak of maximal turbulence production appears to coincide with the upper edge of the bedload layer delimited by a sharp gradient in mean sediment concentration. This vertically upshifted peak of TKE production is accompanied by an enhanced net downward oriented TKE flux when compared with the reference clear-water flows. The downward diffused TKE is found to act in the bedload layer as a local energy source in reasonable balance with the sediment transport term. The mechanism behind this downward TKE transport was further analysed on the basis of coherent flow structure dynamics controlled by ejection- and sweep-type events. The agreement between the height of downward directed mean TKE flux and the height below which sweep-type events dominate the Reynolds shear-stress contribution over ejections, revealed the leading role played by sweeps in mean TKE transport. This agreement holds for all reference clear-water flows supporting the well-known wall-roughness-induced dominance of the sweep contribution in turbulent, rough clear-water boundary layer flows. Furthermore, for all 3 mm particle flows, the two referred to transition levels were significantly and similarly upshifted to the upper edge of the bedload layer. Only for these sediment-laden flows, the bedload layer thickness is seen to exceed the wall-roughness sublayer of the reference clear-water flows. This supports a strong analogy between wall-roughness effects in clear-water flows and bedload layer effects in sediment-laden flows, on the mean TKE budget induced by a similarly modified coherent flow structure dynamics. The bedload layer-controlled wall roughness is finally confirmed by the good prediction of the wall-roughness parameter ks of the logarithmic velocity distribution. An empirical formulation fitting the presented measurements is presented, valid over the range of Shields number values covered herein.

The geometrical properties of streamlines, such as the curvatures, directions and positions, are studied in steady inviscid compressible flow fields via differential geometry theories and conservation laws. The influences of the streamline geometries on the flow speeds and pressures are also identified and discussed. By transforming the streamlines to fill the domain and satisfy the boundary conditions, a unified geometry-based solver, the streamline transformation method, is proposed for both subsonic and supersonic regions. The governing equations and boundary conditions along streamlines and shock waves are also derived. This method is verified by numerical results of three typical flow fields, including the subsonic channel flow, the supersonic downstream of attached shock waves and especially the subsonic/supersonic downstream of detached bow shock waves. Both two-dimensional planar and axisymmetric flow fields are considered. Compared with the results from computational fluid dynamics, good agreements are achieved by this method, while fewer computational resources, by an order of magnitude, are consumed. Features of these flow fields are also analysed from a geometrical perspective, such as flow speeds and pressures deviated by the wall curvatures, and three-dimensional effects in the after-shock flow fields. For a hyperbolic-shaped bow shock wave, the stand-off distances and the transitions from subsonic to supersonic regions are also discussed. As indicated by the accuracy, efficiency and applicability in a wide range of flow speeds, the streamline transformation method would be a potential candidate for the theoretical analysis and inverse design of high-speed flow fields, especially where the subsonic regions exist downstream of strong shock waves.

In the interaction of water waves with marine structures, the interplay between wave diffraction and drag-induced dissipation is seldom, if ever, considered. In particular, linear hydrodynamic models, and extensions thereof through the addition of a quadratic force term, do not represent the change in amplitude of the waves diffracted and radiated to the far field, which should result from local energy dissipation in the vicinity of the structure. In this work, a series of wave flume experiments is carried out, whereby waves of increasing amplitude impinge upon a vertical barrier, extending partway through the flume width. As the wave amplitude increases, the effect of drag – which is known to increase quadratically with the flow velocity – is enhanced, thus allowing the examination of the far-field effect of drag-induced dissipation, in terms of wave reflection and transmission. A potential flow model is proposed, with a simple quadratic pressure drop condition through a virtual porous surface, located on the sides of the barrier (where dissipation occurs). Experimental results confirm that drag-induced dissipation has a marked effect on the diffracted flow, i.e. on wave reflection and transmission, which is appropriately captured in the proposed model. Conversely, when diffraction becomes dominant as the barrier width becomes comparable to the incoming wavelength, the diffracted flow must be accounted for in predicting drag-induced forces and dissipation.

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.