One-degree-of-freedom flow-induced vibration (FIV) and energy harvesting through FIV of an elastically mounted circular cylinder with mechanically coupled rotation were investigated numerically for low Reynolds number 100, mass ratio 8 and a wide range of reduced velocities. The aims of this study are to investigate the effect of the flow direction angle $\beta$ on the vibration and energy harvesting through FIV. Two types of lock-in are found: vortex-induced vibration (VIV) and galloping. The response amplitude increases with the increase of $\beta$ in both regimes. Both VIV response and galloping regimes are found for $\beta$ = 45° to $\beta$ = 90°. For $\beta$ = −90° to $\beta$ = 0°, only VIV response regimes are found. The fluid force and fluid torque play different roles in exciting/damping the vibration. In the high-amplitude gallop regime, the fluid force excites the vibration, and the torque damps the vibration. Energy harvesting at flow direction angle 90° is investigated as this flow direction has the maximum galloping amplitude. The energy harvesting is achieved by a linear electric damping coefficient in the numerical model. The maximum harvestable power in the galloping regime is significantly greater than that in the VIV regime, and it increases with the increase of the reduced velocity. When the reduced velocity is 20, the harvested power is over 20 times that in the VIV regime, and can further increase if reduced velocity further increases. The maximum efficiency over all simulated parameters is 0.424, occurring when the reduced velocity is 20, and electric damping factor is 0.04.

An arbitrary Lagrangian–Eulerian finite element method and numerical implementation for curved and deforming lipid membranes is presented here. The membrane surface is endowed with a mesh whose in-plane motion need not depend on the in-plane flow of lipids. Instead, in-plane mesh dynamics can be specified arbitrarily. A new class of mesh motions is introduced, where the mesh velocity satisfies the dynamical equations of a user-specified two-dimensional material. A Lagrange multiplier constrains the out-of-plane membrane and mesh velocities to be equal, such that the mesh and material always overlap. An associated numerical inf–sup instability ensues, and is removed by adapting established techniques in the finite element analysis of fluids. In our implementation, the aforementioned Lagrange multiplier is projected onto a discontinuous space of piecewise linear functions. The new mesh motion is compared to established Lagrangian and Eulerian formulations by investigating a pre-eminent numerical benchmark of biological significance: the pulling of a membrane tether from a flat patch and its subsequent lateral translation.

A new statistical definition for the mean turbulent boundary layer (TBL) thickness is introduced, based on identification of the wall-normal location where the streamwise velocity skewness changes sign, from negative to positive, in the outermost region of the boundary layer. Importantly, this definition is independent of arbitrary thresholds, and broadly applicable, including to past single-point measurements. Furthermore, this definition is motivated by the phenomenology of streamwise velocity fluctuations near the turbulent/non-turbulent interface (TNTI), whose local characteristics are shown to be universal for TBLs under low free-stream turbulence conditions (i.e. with or without pressure gradients, surface roughness, etc.) through large-scale experiments, simulations and coherent structure-based modelling. The new approach yields a TBL thickness that is consistent with previous definitions, such as those based on Reynolds shear stress or ‘composite’ mean velocity profiles, and which can be used practically, e.g. to calculate integral thicknesses. Two methods are proposed for estimating the TBL thickness using this definition: one based on simple linear interpolation and the other on fitting a generalised Fourier model to the outer skewness profile. The robustness and limitations of these methods are demonstrated through analysis of several published experimental and numerical datasets, which cover a range of canonical and non-canonical TBLs. These datasets also vary in key characteristics such as wall-normal resolution and measurement noise, particularly in the critical TNTI region.

Floating particles deform the liquid–gas interface, which may lead to capillary repulsion or attraction and aggregation of nearby particles (e.g. the Cheerios effect). Previous studies employed the superposition of capillary multipoles to model interfacial deformation for circular or ellipsoidal particles. However, the induced interfacial deformation depends on the shape of the particle and becomes more complex as the geometric complexity of the particle increases. This study presents a generalised solution for the liquid–gas interface near complex anisotropic particles using the domain perturbations approach. This method enables a closed-form solution for interfacial deformation near particles with an anisotropic shape, as well as the varying height of the pinned liquid–gas contact line. We verified the model via experiments performed with fixed particles held at the water level with shapes such as a circle, hexagon and square, which have either flat or sinusoidal pinned contact lines. Although in this study we concentrate on the equilibrium configuration of the liquid–gas interface in the vicinity of particles placed at fixed positions, our methodology paves the way to explore the interactions among multiple floating anisotropic particles and, thus, the role of particle geometry in self-assembly processes of floating particles.

This research examines in detail the complex nonlinear forces generated when steep waves interact with vertical cylindrical structures, such as those typically used as offshore wind turbine foundations. These interactions, particularly the nonlinear wave forces associated with the secondary load cycle, present unanswered questions about how they are triggered. Our experimental campaigns underscore the occurrence of the secondary load cycle. We also investigate how the vertical distributions of the scattering force, pressure field and wave field affect the nonlinear wave forces associated with the secondary load cycle phenomena. A phase-based harmonic separation method isolates harmonic components of the scattering force’s vertical distribution, pressure field and wave field. This approach facilitates the clear separation of individual harmonics by controlling the phase of incident waves, which offers new insights into the mechanisms of the secondary load cycle. Our findings highlight the importance of complex nonlinear wave–structure interactions in this context. In certain wave regimes, nonlinear forces are locally larger than the linear forces, highlighting the need to consider the secondary load cycle in structural design. In addition, a novel discovery emerges from our comparative analysis, whereby very high-frequency (over the fifth in harmonic and order) oscillations, strongly correlated to wave steepness, have the potential to play a role in structural fatigue. This new in-depth analysis provides a unique insight regarding the complex interplay between severe waves and typical cylindrical offshore structures, adding to our understanding of the secondary load cycle for applications related to offshore wind turbine foundations.

Flows enabled by phoretic mechanisms are of significant interest in several biological and biomedical processes, such as bacterial motion and targeted drug delivery. Here, we develop a homogenisation-based macroscopic boundary condition that describes the effective flow across a diffusio-phoretic microstructured membrane, where the interaction between the membrane walls and the solute particles is modelled via a potential approach. We consider two cases where potential variations occur (i) at the pore scale and (ii) only in the close vicinity of the boundary, allowing for a simplified version of the macroscopic flow description, in the latter case. Chemical interactions at the microscale are rigorously upscaled to macroscopic phoretic solvent velocity and solute flux contributions, and added to the classical permeability and diffusivity properties of the membrane. These properties stem from the solution of Stokes advection–diffusion problems at the microscale, some of them forced by an interaction potential term. Eventually, we show an application of the macroscopic model to develop minimal phoretic pumps, showcasing its suitability for efficient design and optimisation procedures.

Surface quasi-geostrophic (SQG) theory describes the two-dimensional active transport of a scalar field, such as temperature, which – when properly rescaled – shares the same physical dimension of length/time as the advecting velocity field. This duality has motivated analogies with fully developed three-dimensional turbulence. In particular, the Kraichnan – Leith – Batchelor similarity theory predicts a Kolmogorov-type inertial range scaling for both scalar and velocity fields, and the presence of intermittency through multifractal scaling was pointed out by Sukhatme & Pierrehumbert (2002 Chaos 12, 439–450), in unforced settings. In this work, we refine the discussion of these statistical analogies, using numerical simulations with up to $16\,384^2$ collocation points in a steady-state regime dominated by the direct cascade of scalar variance. We show that mixed structure functions, coupling velocity increments with scalar differences, develop well-defined scaling ranges, highlighting the role of anomalous fluxes of all the scalar moments. However, the clean multiscaling properties of SQG transport are blurred when considering velocity and scalar fields separately. In particular, the usual (unmixed) structure functions do no follow any power-law scaling in any range of scales, neither for the velocity nor for the scalar increments. This specific form of the intermittency phenomenon reflects the specific kinematic properties of SQG turbulence, involving the interplay between long-range interactions, structures and geometry. Revealing the multiscaling in single-field statistics requires us to resort to generalised notions of scale invariance, such as extended self-similarity and a specific form of refined self-similarity. Our findings emphasise the fundamental entanglement of scalar and velocity fields in SQG turbulence: they evolve hand in hand and any attempt to isolate them destroys scaling in its usual sense. This perspective sheds new lights on the discrepancies in spectra and structure functions that have been repeatedly observed in SQG numerics for the past 20 years.

We study the mechanics of evaporation and precipitate formation in pure and bacteria-laden sessile whole blood droplets in the context of disease diagnostics. Using experimental and theoretical analysis, we show that the evaporation process has three stages based on evaporation rate. In the first stage, edge evaporation results in a gelated contact line along the periphery through a sol–gel phase transition. The intermediate stage consists of a gelated front propagating radially inwards due to capillary flow and droplet height regression in pinned mode, forming a wet-gel phase. We unearthed that the gelation of the entire droplet occurs in the second stage, and the wet-gel formed contains trace amounts of water. In the final slowest stage, the wet gel transforms into a dry gel, leading to desiccation-induced stress forming diverse crack patterns in the precipitate. Slow evaporation in the final stage is quantitatively measured using evaporation of trace water and associated transient delamination of the precipitate. Using the axisymmetric lubrication approximation, we compute the transient droplet height profile and the erythrocytes concentration for the first two stages of evaporation. We show that the precipitate thickness profile computed from the theoretical analysis conforms to the optical profilometry measurements. We show that the drop evaporation rate and final dried residue pattern do not change appreciably within the parameter variation of the bacterial concentration typically found in bacterial infection of living organisms. However, at exceedingly high bacterial concentrations, the cracks formed in the coronal region deviate from the typical radial cracks found in lower concentrations.

This study compares turbulent channel flows over elastic walls with those over rough walls, to explore the role of the dynamic change of shape of the wall in turbulence. The comparison is made meaningful by generating rough walls from instantaneous configurations of elastic cases. The aim of this comparison is to individually understand the role of fluid–structure interaction effects and the role of wall shape/undulations in determining the overall physics of flow near elastic walls. With an increase in the compliance of the wall, qualitatively similar trends for many of the effects produced by a rough wall are also seen in the elastic wall. However, specific features can be observed for the elastic-wall cases only, arising from the mutual interaction between the solid and fluid, leading to a further increase in drag. To understand them, we look at the turbulent structures, which exhibit clear differences across the various configurations: roughness induces only a slight reduction of streamwise coherency, resulting in a situation qualitatively similar to what is found in classical turbulent channel flows, whereas elasticity causes the emergence of a novel dominant spanwise coherency. Additionally, we explored the effect of vertical disturbances on elastic-wall dynamics by comparing with permeable walls having similar (average) wall-normal velocity fluctuations at the interface. The permeable walls were found to have minimal similarities to elastic walls. Overall, we can state that the wall motion caused by the complex fluid–structure interaction contributes significantly to the flow and must be considered when modelling it. In particular, we highlight the emergence of strong wall-normal fluctuations near the wall, which result in strong ejection events, an attribute not observed for rigid walls.

The simulation of turbulent flow requires many degrees of freedom to resolve all the relevant time and length scales. However, due to the dissipative nature of the Navier–Stokes equations, the long-term dynamics is expected to lie on a finite-dimensional invariant manifold with fewer degrees of freedom. In this study, we build low-dimensional data-driven models of pressure-driven flow through a circular pipe. We impose the ‘shift-and-reflect’ symmetry to study the system in a minimal computational cell (e.g. the smallest domain size that sustains turbulence) at a Reynolds number of 2500. We build these models by using autoencoders to parametrise the manifold coordinates and neural ordinary differential equation to describe their time evolution. Direct numerical simulations (DNSs) typically require of the order of $\mathcal{O}(10^5)$ degrees of freedom, while our data-driven framework enables the construction of models with fewer than 20 degrees of freedom. Remarkably, these reduced-order models effectively capture crucial features of the flow, including the streak breakdown. In short-time tracking, these models accurately track the true trajectory for one Lyapunov time, as well as the leading Lyapunov exponent, while at long-times, they successfully capture key aspects of the dynamics such as Reynolds stresses and energy balance. The model can quantitatively capture key characteristics of the flow, including the streak breakdown and regeneration cycle. Additionally, we report new exact coherent states found in the DNS with the aid of these low-dimensional models. This approach leads to the discovery of seventeen previously unknown solutions within the turbulent pipe flow system, notably featuring relative periodic orbits characterised by the longest reported periods for such flow conditions.

Particle suspensions at the interface of turbulent liquids are governed by the balance of capillary attraction, strain-induced drag and lubrication. Here, we extend previous findings, obtained for small particles whose capillary interactions are dominated by quadrupolar-mode deformation of the interface, to larger spherical and disc-shaped particles experiencing monopole-dominant capillarity. By combining pair-approach experiments, two-dimensional turbulent flow realizations and particle imaging, we demonstrate that particles experiencing monopole-dominant attraction exhibit enhanced clustering compared with their quadrupole-dominant counterparts. We introduce an interaction scale defined by balancing viscous drag and capillary attraction, which is compared with the particle size and interparticle distance. This allows us to map the clustering behaviour onto a parameter space solely defined by those characteristic length scales. This yields a unified framework able to predict the tendency to cluster (and the concentration threshold for those clusters to percolate) in a vast array of fluid–particle systems.

Turbulence accounts for most of the energy losses associated with the pumping of fluids in pipes. Pulsatile drivings can reduce the drag and energy consumption required to supply a desired mass flux, when compared with steady driving. However, not all pulsation waveforms yield reductions. Here, we compute drag- and energy-optimal driving waveforms using direct numerical simulations and a gradient-free black-box optimisation framework. Specifically, we show that Bayesian optimisation is vastly superior to ordinary gradient-based methods in terms of computational efficiency and robustness, due to its ability to deal with noisy objective functions, as they naturally arise from the finite-time averaging of turbulent flows. We identify optimal waveforms for three Reynolds numbers and two Womersley numbers. At a Reynolds number of $8600$ and a Womersley number of 10, optimal waveforms reduce total energy consumption by 22 % and drag by 37 %. These reductions are rooted in the suppression of turbulence prior to the acceleration phase, the resulting delay in turbulence onset, and the radial localisation of turbulent kinetic energy and production towards the pipe centre. Our results pinpoint that the predominant, steady operation mode of pumping fluids through pipes is far from optimal.

The present experiments investigated the combustion dynamics of single and coaxial laminar diffusion flames within a closed cylindrical acoustic waveguide, focusing on their response to acoustic forcing at a pressure antinode. Nine alternative fuel injectors were used to examine the effect of injector jet diameter and configuration, tube wall thickness, annular-to-inner area and velocity ratio, and jet Reynolds number (below 100) on flame behaviour under different applied frequencies and pressure perturbation amplitudes. Fundamental flame–acoustic coupling phenomena were identified, all of which involved symmetric flame perturbations. These included sustained oscillatory combustion (SOC), multi-frequency periodic liftoff and reattachment (PLOR), permanent flame lift-off (PFLO) with low-level oscillations, and flame blowoff (BO). The phase lag between acoustic forcing and flame response was quantified, providing valuable insights into the coupling dynamics and transition behaviours. Findings revealed how various geometrical and flow characteristics could affect flame stability and resistance to blowoff, even under similar acoustic forcing conditions. Analysis of high-speed spatiotemporal visible imaging using proper orthogonal decomposition (POD) uncovered additional distinct phase portraits and spectral signatures associated with instability transitions, which, coupled with specific dynamical characteristics, enabled new insights into the relevance of injector geometrical characteristics and flow conditions in addressing acoustically coupled combustion instabilities.

The growth of small perturbations in isotropic turbulence is studied using massive ensembles of direct numerical simulations. These ensembles capture the evolution of the ensemble-averaged flow field and the ensemble variance in the fully nonlinear regime of perturbation growth. Evolution equations for these two fields are constructed by applying the ensemble average operator to the Navier–Stokes equations and used to study uncertainty growth in scale and physical space. It is shown that uncertainty growth is described by a flux of energy from the ensemble-averaged flow to the ensemble variance. This flux is formally equivalent to the subgrid scale (SGS) energy fluxes of the turbulence cascade, and can be interpreted as an inverse uncertainty cascade from small to large scales. In the absence of information sources (measurements), the uncertainty cascade is unsteady and leads to the progressive filtering of the small scales in the ensemble-averaged flow, a process that represents the loss of predictability due to chaos. Similar to the kinetic energy cascade, the uncertainty cascade displays an inertial range with a constant average uncertainty flux, which is bounded from below by the average kinetic energy dissipation. Locally in space, uncertainty fluxes differ from the SGS energy fluxes at the same scale, but both have similar statistics and are significantly correlated with each other in space. This suggests that uncertainty propagation is partly connected to the energy cascade and that they share similar mechanisms. These findings open avenues to model uncertainty propagation in turbulence following an approach similar to the SGS models in large-eddy simulations. This is relevant not only to efficiently assess the reliability and accuracy of turbulence forecasts, but also to design uncertainty-robust reconstruction techniques for data assimilation or SGS modelling.

Research on water wave metamaterials based on local resonance has advanced rapidly. However, their application to floating structures for controlling surface gravity waves remains underexplored. In this work, we introduce the floating metaplate, a periodic array of resonators on a floating plate that leverages locally resonant bandgaps to effectively manipulate surface gravity waves. We employ the eigenfunction matching method combined with Bloch’s theorem to solve the wave–structure interaction problem and obtain the band structure of the floating metaplate. An effective model based on averaging is developed, which agrees well with the results of numerical simulation, elucidating the mechanism of bandgap formation. Both frequency- and time-domain simulations demonstrate the floating metaplate’s strong wave attenuation capabilities. Furthermore, by incorporating a gradient in the resonant frequencies of the resonators, we achieve the rainbow trapping effect, where waves of different frequencies are reflected at distinct locations. This enables the design of a broadband wave reflector with a tuneable operation frequency range. Our findings may lead to promising applications in coastal protection, wave energy harvesting and the design of resilient offshore renewable energy systems.

Mass dispersion in oscillatory flows is closely tied to various environmental and biological processes, differing markedly from dispersion in steady flows due to the periodic expansion and contraction of particle patches. In this study, we investigate the Taylor–Aris dispersion of active particles in laminar oscillatory flows between parallel plates. Two complementary approaches are employed: a two-time-variable expansion of the Smoluchowski equation is used to facilitate Aris’ method of moments for the pre-asymptotic dispersion, while the generalised Taylor dispersion theory is extended to capture phase-dependent periodic drift and dispersivity in the long-time asymptotic limit. Applying both frameworks, we find that spherical non-gyrotactic swimmers can exhibit greater or lesser diffusivity than passive solutes in purely oscillatory flows, depending on the oscillation frequency. This behaviour arises primarily from the disruption of cross-streamline migration governed by Jeffery orbits. When a steady component is superimposed, oscillation induces a non-monotonic dual effect on diffusivity. We further examine two well-studied shear-related accumulation mechanisms, arising from gyrotaxis and elongation. Although these accumulation effects are less pronounced than in steady flows due to flow unsteadiness, gyrotactic swimmers respond more strongly to the unsteady shear profile, significantly modifying their drift and dispersivity. This work offers new insights into the dispersion of active particles in oscillatory flows, and also provides a foundation for studying periodic active dispersion beyond the oscillatory flow, such as periodic variations in shape and swimming speed.

Flutter in lightweight airfoils under unsteady flows presents a critical challenge in aeroelastic stability and control. This study uncovers phase-dependent effects that drive the onset and suppression of flutter in a freely pitching airfoil at low Reynolds number. By introducing targeted impulsive stiffness perturbations, we identify critical phases that trigger instability. Using phase-sensitivity functions, energy-transfer metrics and dynamic mode decomposition, we show that flutter arises from phase lock-on between structural and fluid modes. Leveraging this insight, we design an energy-optimal, phase-based control strategy that applies transient heaving motions to disrupt synchronisation and arrest unstable growth. This minimal, time-localised control suppresses subharmonic amplification and restores stable periodic motion.

In biology, cells undergo deformations under the action of flow caused by the fluid surrounding them. These flows lead to shape changes and instabilities that have been explored in detail for single component vesicles. However, cell membranes are often multicomponent in nature, made up of multiple phospholipids and cholesterol mixtures that give rise to interesting thermodynamics and fluid mechanics. Our work analyses shear flow around a multicomponent vesicle using a small-deformation theory based on vector and scalar spherical harmonics. We set up the problem by laying out the governing momentum equations and the traction balance arising from the phase separation and bending. These equations are solved along with a Cahn–Hilliard equation that governs the coarsening dynamics of the phospholipid–cholesterol mixture. We provide a detailed analysis of the vesicle dynamics (e.g. tumbling, breathing, tank-treading and swinging/phase-treading) in two regimes – when flow is faster than coarsening dynamics (Péclet number ${\textit{Pe}} \gg 1$) and when the two time scales are comparable ($\textit{Pe} \sim O(1)$) – and provide a discussion on when these behaviours occur. The analysis aims to provide an experimentalist with important insights pertaining to the phase separation dynamics and their effect on the deformation dynamics of a vesicle.

In this investigation, the effect of Ekman pumping on a quasi-geostrophic (QG) system is explored via the vertical buoyancy flux. The vertical buoyancy flux is the quantity in QG flows that is responsible for the adiabatic transfer between kinetic energy (KE) and available potential energy (APE), as well as the slow-time evolution of the mean buoyancy. Ekman pumping (or suction) is a phenomenon that arises through conservation of mass at no-slip boundaries of rotating fluid systems. Three-dimensional QG numerical simulations are run with and without Ekman pumping at the bottom boundary, as well as with and without a realistic stratification profile. Through theory and numerical experiment, it is shown that Ekman pumping drives a conversion of energy from APE to KE at small scales, and from KE to APE at large scales, even in the absence of a mean isopycnal slope. It is also shown that Ekman pumping affects the mean buoyancy by slightly weakening the stratification near the bottom boundary.

Large numbers of relative periodic orbits (RPOs) have been found recently in doubly periodic, two-dimensional Kolmogorov flow at moderate Reynolds numbers ${\textit{Re}} \in \{40, 100\}$. While these solutions lead to robust statistical reconstructions at the ${\textit{Re}}$ values where they were obtained, it is unclear how their dynamical importance changes with ${\textit{Re}}$. Arclength continuation on this library of solutions reveals that large numbers of RPOs quickly become dynamically irrelevant, reaching dissipation values either much larger or smaller than the values typical of the turbulent attractor at high ${\textit{Re}}$. The scaling of the high-dissipation RPOs is shown to be consistent with a direct connection to solutions of the unforced Euler equation, and is observed for a wide variety of states beyond the ‘unimodal’ solutions considered in previous work (Kim & Okamoto, Nonlinearity vol. 28, 2015, p. 3219). However, the weakly dissipative states have properties indicating a connection to exact solutions of a forced Euler equation. The dynamical irrelevance of many solutions leads to poor statistical reconstruction at higher ${\textit{Re}}$, raising serious questions for the future use of RPOs for estimating probability densities. Motivated by the Euler connection of some of our RPOs, we also show that many of these states can be well described by exact relative periodic solutions in a system of point vortices. The point vortex RPOs are converged via gradient-based optimisation of a scalar loss function which (i) matches the dynamics of the point vortices to the turbulent vortex cores and (ii) insists the point vortex evolution is itself time-periodic.