In Rayleigh–Bénard convection and Taylor–Couette flow cellular patterns emerge at the onset of instability and persist as large-scale coherent structures in the turbulent regime. Their long-term dynamics has been thoroughly characterised and modelled for the case of turbulent convection, whereas turbulent Taylor rolls have received much less attention. Here we present direct numerical simulations of axisymmetric Taylor–Couette flow in the corotating regime and reveal a transition to spatio–temporal chaos as the system size increases. Beyond this transition, Taylor rolls suddenly undergo erratic drifts evolving on a very slow time scale. We estimate an effective diffusion coefficient for the drift and compare the dynamics with analogous motions in Rayleigh–Bénard convection and Poiseuille flow, suggesting that this spontaneous diffusive displacement of large coherent structures is common among different types of wall-bounded turbulent flows.

Active fluids encompass a wide range of non-equilibrium fluids, in which the self-propulsion or rotation of their units can give rise to large-scale spontaneous flows. Despite the diversity of active fluids, they are commonly viscoelastic. Therefore, we develop a hydrodynamic model of isotropic active liquids by accounting for their viscoelasticity. Specifically, we incorporate an active stress term into a general viscoelastic liquid model to study the spontaneous flow states and their transitions in two-dimensional channel, annulus and disk geometries. We have discovered rich spontaneous flow states in a channel as a function of activity and Weissenberg number, including unidirectional flow, travelling-wave and vortex-roll states. The Weissenberg number acts against activity by suppressing the spontaneous flow. In an annulus confinement, we find that a net flow can be generated only if the aspect ratio of the annulus is not too large nor too small, akin to some three-dimensional active-flow phenomena. In a disk geometry, we observe a periodic chirality switching of a single vortex flow, resembling the bacteria-based active fluid experiments. The two phenomena reproduced in our model differ in Weissenberg number and frictional coefficient. As such, our active viscoelastic model offers a unified framework to elucidate diverse active liquids, uncover their connections and highlight the universality of dynamic active-flow patterns.

We investigate the sliding dynamics of a millimetre-sized particle trapped in a horizontal soap film. Once released, the particle moves toward the centre of the film in damped oscillations. We study experimentally and model the forces acting on the particle, and evidence the key role of the mass of the film on the shape of the film and particle dynamics. Not only is the gravitational distortion of the film measurable, it completely determines the force responsible for the motion of the particle – the catenoid-like deformation induced by the particle has negligible effect on the dynamics. Surprisingly, this is expected for all film sizes as long as the particle radius remains much smaller than the film width. We also measure the friction force, and show that ambient air and the film contribute almost equally to the friction. The theoretical model that we propose predicts exactly the friction coefficient as long as inertial effects can be neglected in air (for the smallest and slowest particles). The fit between theory and experiments sets an upper boundary $\eta _s \leqslant 10^{-8}$ Pa s m for the surface viscosity, in excellent agreement with recent interfacial microrheology measurements.

This article explores how a submerged elastic plate, clamped at one edge, interacts with water waves. Submerged elastic plates have been considered as potentially effective design elements in the development of wave energy harvesters but their behaviour in a wave field remains largely unexplored, especially experimentally. Positioned at a fixed depth in a wave tank, the flexible plate demonstrates significant wave reflection capabilities, a characteristic absent in rigid plates of identical dimensions. The experiments thus reveal that plate motion is crucial for wave reflection. Sufficiently steep waves are shown to induce a change in the mean position of the plate, with the trailing edge reaching the free surface in some cases. This configuration change is found to be particularly efficient to break water waves. These findings contribute to understanding the potential of elastic plates for wave energy harvesting and wave attenuation scenarios.

A collection of secondary instability calculations in streaky boundary layers is presented. The data are retrieved from well-resolved numerical simulations of boundary layers forced by free-stream turbulence (FST), considering different geometries and FST conditions. The stability calculations are performed before streak breakdown, taking place at various $Rey_x$ the Reynolds number based on the streamwise coordinate. Despite the rich streak population of various sizes, it is found that breaking streaks have similar aspect ratios, independently of the streamwise position where they appear. This suggests that wider streaks will break down further downstream than thinner ones, making the appearance of secondary instabilities somewhat independent of the streak’s wavelength. Moreover, the large difference in the integral length scale among the simulations suggests that this aspect ratio is also independent of the FST scales. An explanation for this behaviour is provided by showing that these breaking streaks are in the range of perturbations that can experience maximum transient growth according to optimal disturbance theory. This could explain why, at a given streamwise position, there is a narrow spanwise wavelength range where streak breakdown is more likely to occur.

Skin-friction drag reduction (DR) in a turbulent boundary layer (TBL) using plasma-generated streamwise vortices (PGSVs) is governed by plasma-induced spanwise wall-jet velocity $W$, the distance $L$ between the positive electrodes of two adjacent plasma actuators (PAs) and the friction Reynolds number $Re_\tau$. It is found experimentally that DR increases logarithmically with the growing maximum spanwise mean velocity $\overline {W}_{max}^+$ but decreases with rising $L^+$ and $Re_\tau$, where superscript ‘+’ denotes normalization by the inner scales. It is further found from theoretical and empirical scaling analyses that the dimensionless drag variation $\Delta F = g_1 (\overline {W}_{max}^+, L^+, {Re_\tau })$ may be reduced to $\Delta F = g_2 (\xi )$, where $g_1$ and $g_2$ are different functions and the scaling factor $\xi = [k_{2} \log _{10} (k_{1} \overline {W}_{max }^{+} ) ] / (L^{+} Re_{\tau } )$ ($k_{2}$ and $k_{1}$ are constants) is physically the circulation of the PGSVs. Discussion is conducted based on $\Delta F = g_2 (\xi )$, which provides important insight into the physics of TBL control based on PAs.

Elastoinertial turbulence (EIT) is a chaotic state that emerges in the flows of dilute polymer solutions. Direct numerical simulation (DNS) of EIT is highly computationally expensive due to the need to resolve the multiscale nature of the system. While DNS of two-dimensional (2-D) EIT typically requires $O(10^6)$ degrees of freedom, we demonstrate here that a data-driven modelling framework allows for the construction of an accurate model with 50 degrees of freedom. We achieve a low-dimensional representation of the full state by first applying a viscoelastic variant of proper orthogonal decomposition to DNS results, and then using an autoencoder. The dynamics of this low-dimensional representation is learned using the neural ordinary differential equation (NODE) method, which approximates the vector field for the reduced dynamics as a neural network. The resulting low-dimensional data-driven model effectively captures short-time dynamics over the span of one correlation time, as well as long-time dynamics, particularly the self-similar, nested travelling wave structure of 2-D EIT in the parameter range considered.

Explaining fast magnetic reconnection in electrically conducting plasmas has been a theoretical challenge in plasma physics since its first description by Eugene N. Parker. In recent years, the observed reconnection rate has been shown by numerical simulations to be explained by the plasmoid instability that appears in highly conductive plasmas. In this work, by studying numerically the Orszag–Tang vortex, we show that the plasmoid instability is very sensitive to the numerical resolution used. It is shown that well-resolved runs display no plasmoid instability even at Lundquist numbers as large as $5\times 10^5$ achieved at resolutions of $32\,768^2$ grid points. On the contrary, in simulations that are under-resolved below a threshold, the plasmoid instability manifests itself with the formation of larger plasmoids the larger the under-resolving is. The present results thus emphasize the importance of performing convergence tests in numerical simulations and suggest that further investigations on the nonlinear evolution of the plasmoid instability are required.

Electro-osmotic flow (EOF) in nanochannels exhibits a puzzling non-monotonic dependence on salt concentration, which contrasts with observations in microchannels and remains not fully understood. In this work, we address this phenomenon through a theoretical investigation of EOF in $\mathrm{pH}$-regulated channels. New analytical approximations for electrostatic potential, EOF profile and electro-osmotic mobility beyond the Debye–Hückel limit are derived through asymptotic analysis. Our findings reveal that the surface electrostatic potential is independent of the channel size only when the half-channel size exceeds the Gouy–Chapman length. In contrast, surface ionization and net charge distribution play more crucial roles in EOF at the nanoscale, as they govern both the magnitude and the spatial distribution of the Coulomb driving force. As salt concentration increases, EOF velocity initially rises due to enhanced surface ionization, followed by a decline attributed to increased wall shear stress. This work provides key insights for EOF applications in nanofluidics and biomedical devices, and deepens the understanding of electrokinetic phenomena influenced by $\mathrm{pH}$-regulation effects.

This study introduces a novel approach to investigate the Reynolds analogy in complex flow scenarios. It is shown that the total mechanical energy $\mathit {B}$, viz. the sum of kinetic energy and pressure work, and the field $\Gamma =\theta ^2/2$ (where $\theta$ is the transported passive scalar) are governed by two equations that are similar in form, when time-averaged for statistically stationary flows. For fully developed channel flows the integral energy balance links the mean bulk velocity and scalar with the volume averages of the respective dissipation rates, allowing the assessment of the Reynolds analogy in terms of the dissipation fields. This approach is tested on direct numerical simulation data of rough-wall turbulent channel flow at two different roughness Reynolds numbers, namely $k^+=15$ and $k^+=90$. For a unit Prandtl number, the same qualitative behaviour is observed for the mean wall-normal distributions of the budget-equation terms of $B$ and $\Gamma$, the latter being larger than the corresponding terms in the mechanical-energy budget. The Reynolds decomposition of the flow into temporal mean and stochastic parts reveals that roughness primarily affects the mean-flow dissipation. For the $k^+=90$ case, the analysis shows that attached-flow and high-shear regions dominate the integral mean scalar and momentum transfer and exhibit the greatest differences between the mean mechanical and scalar dissipation rates. In contrast, well-mixed regions, sheltered by large roughness elements, contribute similarly and minimally to the integral scalar and momentum transfer.

Inspired by laboratory experiments showing internal waves generated by a plume impinging upon a stratified fluid layer (Ansong & Sutherland. 2010 J. Fluid Mech. 648, 405–434), we perform large eddy simulations in three dimensions to examine the structure and source of internal waves emanating from the top of a plume that rises vertically into stratification whose strength ranges over two orders of magnitude between different simulations. Provided the plume is sufficiently energetic to penetrate into the stratified layer, internal waves are generated with frequencies in a relatively narrow band moderately smaller than the buoyancy frequency. Through adaptations of ray theory including viscosity and use of dynamic mode decomposition, we show that the waves originate from within the turbulent flow rather than at the turbulent/non-turbulent interface between the fountain top and the surrounding stratified fluid.

We report the first shock-tube experiments on Richtmyer–Meshkov instability at a single-mode light–heavy interface accelerated by a strong shock wave with Mach number higher than 3.0. Under the proximity effect of the transmitted shock and its induced secondary compression effect, the interface profile is markedly different from that in weakly compressible flows. For the first time, the validity of the compressible linear theory and the failure of the impulsive model in predicting the linear amplitude evolution in highly compressible flows are verified through experiments. Existing nonlinear and modal models fail to accurately describe the perturbation evolution, as they do not account for the shock proximity and secondary compression effects on interface evolution. The shock proximity effect manifests mainly in the early stages when the transmitted shock remains close to the interface, while the effect of secondary compression manifests primarily at the period when interactions of transverse shocks occur at the bubble tips. Based on these findings, we propose an empirical model capable of predicting the bubble evolution in highly compressible flows.

The dynamics of evolving fluid films in the viscous Stokes limit is relevant to various applications, such as the modelling of lipid bilayers in cells. While the governing equations were formulated by Scriven (1960), solving for the flow of a deformable viscous surface with arbitrary shape and topology has remained a challenge. In this study, we present a straightforward discrete model based on variational principles to address this long-standing problem. We replace the classical equations, which are expressed with tensor calculus in local coordinates, with a simple coordinate-free, differential-geometric formulation. The formulation provides a fundamental understanding of the underlying mechanics and translates directly to discretization. We construct a discrete analogue of the system using Onsager's variational principle, which, in a smooth context, governs the flow of a viscous medium. In the discrete setting, instead of term-wise discretizing the coordinate-based Stokes equations, we construct a discrete Rayleighian for the system and derive the discrete Stokes equations via the variational principle. This approach results in a stable, structure-preserving variational integrator that solves the system on general manifolds.

Super-resolution of turbulence is a term used to describe the prediction of high-resolution snapshots of a flow from coarse-grained observations. This is typically accomplished with a deep neural network and training usually requires a dataset of high-resolution images. An approach is presented here in which robust super-resolution can be performed without access to high-resolution reference data, as might be expected in an experiment. The training procedure is similar to data assimilation, wherein the model learns to predict an initial condition that leads to accurate coarse-grained predictions at later times, while only being shown coarse-grained observations. Implementation of the approach requires the use of a fully differentiable flow solver in the training loop to allow for time-marching of predictions. A range of models are trained on data generated from forced, two-dimensional turbulence. The networks have reconstruction errors which are similar to those obtained with ‘standard’ super-resolution approaches using high-resolution data. Furthermore, the methods are comparable to the performance of standard data assimilation for state estimation on individual trajectories, outperforming these variational approaches at initial time and remaining robust when unrolled in time where performance of the standard data-assimilation algorithm improves.

We develop a simple model which describes the repeated injection and extraction of hydrogen in a permeable water-saturated rock which has the form of an anticline. We demonstrate that the flow is controlled by the dimensionless ratio of the square of the buoyancy speed to the product of the two-dimensional volume injection rate and the injection–extraction frequency, and we explore the cases in which this ratio is large and small. Over the first few cycles, the volume of hydrogen in the system gradually builds up since during the extraction phase, some of the water eventually reaches the extraction well, and in our model the system ceases to extract fluid for the remainder of this extraction phase. After many cycles, there is sufficient hydrogen in the system that a quasi-equilibrium state develops in which the mass of fluid injected matches the mass extracted over the course of a cycle. We show that in this equilibrium, the ratio between the mass of gas remaining in the aquifer at the end of the extraction phase, known as the cushion gas, to the mass of gas injected, known as the working gas, decreases if either the flow rate or frequency of the cycles decrease or the buoyancy speed increases, leading to more efficient storage.

We study experimentally the onset of Faraday waves near the end walls of a rectangular vessel containing two stably stratified fluid layers, subject to horizontal oscillations. These subharmonic waves (SWs) are excited, because the horizontal inertial forcing drives a harmonic propagating wave which displaces the interface in the vertical direction at the end walls. We find that the onset of SWs is regulated by a balance between capillary and viscous forces, where the rate of damping is set by the Stokes layer thickness at the wall rather than the wavelength of the SWs. We model the onset of SWs with a weakly damped Mathieu equation and find that the dimensional critical acceleration scales as $\nu _m^{1/2} \omega ^{3/2}$, where $\nu _m$ is the mean viscosity and $\omega$ is the frequency of forcing, in excellent agreement with the experiment over a wide range of parameters.

We present a theory that quantifies the interplay between intrapore and interpore flow variabilities and their impact on hydrodynamic dispersion. The theory reveals that porous media with varying levels of structural disorder exhibit notable differences in interpore flow variability, characterised by the flux-weighted probability density function (PDF), $\hat {\psi }_\tau (\tau ) \sim \tau ^{-\theta -2}$, for advection times $\tau$ through conduits. These differences result in varying relative strengths of interpore and intrapore flow variabilities, leading to distinct scaling behaviours of the hydrodynamic dispersion coefficient $D_L$, normalised by the molecular diffusion coefficient $D_m$, with respect to the Péclet number $Pe$. Specifically, when $\hat {\psi }_\tau (\tau )$ exhibits a broad distribution of $\tau$ with $\theta$ in the range of $(0, 1)$, the dispersion undergoes a transition from power-law scaling, $D_L/D_m \sim Pe^{2-\theta }$, to linear scaling, $D_L/D_m \sim Pe$, and eventually to logarithmic scaling, $D_L/D_m \sim Pe\ln (Pe)$, as $Pe$ increases. Conversely, when $\tau$ is narrowly distributed or when $\theta$ exceeds 1, dispersion consistently follows a logarithmic scaling, $D_L/D_m \sim Pe\ln (Pe)$. The power-law and linear scaling occur when interpore variability predominates over intrapore variability, while logarithmic scaling arises under the opposite condition. These theoretical predictions are supported by experimental data and network simulations across a broad spectrum of porous media.

Empirical evidence is provided that within the inertial sublayer (i.e. logarithmic region) of adiabatic turbulent flows over smooth walls, the skewness of the vertical-velocity component $S_w$ displays universal behaviour, being a positive constant and constrained within the range $S_w \approx 0.1\unicode{x2013}0.16$, regardless of flow configuration and Reynolds number. A theoretical model is then proposed to explain this behaviour, including the observed range of variations of $S_w$. The proposed model clarifies why $S_w$ cannot be predicted from down-gradient closure approximations routinely employed in large-scale meteorological and climate models. The proposed model also offers an alternative and implementable approach for such large-scale models.

Capsules, which are potentially active fluid droplets enclosed in a thin elastic membrane, experience large deformations when placed in suspension. The induced fluid–structure interaction stresses can potentially lead to rupture of the capsule membrane. While numerous experimental studies have focused on the rheological behaviour of capsules until rupture, there remains a gap in understanding the evolution of their mechanical properties and the underlying mechanisms of damage and breakup under flow. We here investigate the damage and rupture of bioartificial microcapsules made of ovalbumin reticulated with terephthaloyl chloride and placed in simple shear flow. We characterize damage by identifying how the surface shear modulus of the capsule membrane changes over time. Rupture is then characterized by comparing the number and size distribution of capsules before and after exposure to shear, while varying the shear rates and time during which capsules are sheared. Our findings reveal how the percentage of ruptured capsules increases with their size, exposure time to shear, and the ratio of viscous to elastic forces at rupture.

The velocity gradient tensor can be decomposed into normal straining, pure shearing and rigid rotation tensors, each with distinct symmetry and normality properties. We partition the strength of turbulent velocity gradients based on the relative contributions of these constituents in several canonical flows. These flows include forced isotropic turbulence, turbulent channels and turbulent boundary layers. For forced isotropic turbulence, the partitioning is in excellent agreement with previous results. For wall-bounded turbulence, the partitioning collapses onto the isotropic partitioning far from the wall, where the mean shearing is relatively weak. By contrast, the near-wall partitioning is dominated by shearing. Between these two regimes, the partitioning collapses well at sufficiently high friction Reynolds numbers and its variations in the buffer layer and the log-law region can be reasonably modelled as a function of the mean shearing strength. Altogether, our results highlight the expressivity and broad applicability of the velocity gradient partitioning as advantages for turbulence modelling.