The motion of a bubble of negligible viscosity, such as air, forced down a tube filled with a viscous fluid which wets the walls of the tube has become a classic of the fluid dynamical literature. The differential motion of the bubble and the fluid are determined by the thin film which surrounds the bubble, whose shape and thickness are set by the interplay between gradients in surface tension and viscous shear stresses. Bretherton (J. Fluid Mech., vol. 10, issue 2, 1961, p. 166) provided a first, clear mathematical analysis in the lubrication limit coupled with carefully constructed experimental confirmation of the thin films deposited by a bubble moving in the confining geometry of the capillary tube. Its lasting impact has been not only in the migration of bubbles, but in a host of related fluid dynamical, industrial, biological and environmental processes for which thin lubricating films on the sometimes convoluted geometries of complex microstructures, such as porous media, determine the large-scale behaviour.

The entrainment hypothesis states that the mean inflow velocity across the boundary of a turbulent flow is proportional to a characteristic velocity of the flow. Proposed by G. I. Taylor approximately 80 years ago, it is still a common model of turbulence closure widely used in environmental engineering and geophysical fluid mechanics. Although it is a very simple concept and mathematical model, it has proven to be able to predict the entrainment in a variety of geophysical flows, e.g. convective clouds and plumes from erupting volcanoes in the atmosphere; dense water overflows and turbidity currents in the ocean; magma injection in a magma chamber in the interior of the Earth, to name just a few. In a seminal paper, Turner (J. Fluid Mech., vol. 173, 1986, pp. 431–471) presents a variety of laboratory and geophysical flows to illustrate the success of the entrainment hypothesis and discusses why such a simple hypothesis works so well even when the original assumptions are no longer valid.

The paper by Castaing et al. (J. Fluid Mech., vol. 204, 1989, pp. 1–30) on turbulent Rayleigh–Bénard convection has been one of the most impactful papers on the subject – not by giving the right and complete answers but by developing versatile concepts and by asking the right questions, namely: (i) What is the overall flow organization? (ii) What is the dependence of the Nusselt number ${\textit {Nu}}$ (the dimensionless heat transport) on the Rayleigh number ${\textit {Ra}}$ (the thermal driving strength)? (iii) What is the ultimate state of turbulence for extremely large ${\textit {Ra}}$? Thanks to Castaing et al. having asked the right questions, the field has made tremendous progress over the last 35 years.

The featured article ‘Break-up of a falling drop containing dispersed particles’ (Nitsche and Batchelor, J. Fluid Mech., 1997, vol. 340, pp. 161–175) is G. K. Batchelor's last published paper with his former postdoctoral associate J. M. Nitsche. The objective of the study was to investigate the randomness of the velocities of interacting rigid particles falling under gravity through a viscous fluid at a small Reynolds number and its consequence for the breakup of a falling cloud of particles. The study focused on a quintessential problem of the collective dynamics of interacting particles and has been an inspiration for subsequent work.

Boundary layers are present in many natural and industrial fluid flows. The concept of boundary layers can be traced back to Leonardo da Vinci's paintings of pipe flow, where he was aware of a higher velocity away from the walls. During the 19th century, the physics of boundary conditions had been extensively debated, and the well-known Maxwell–Navier slip length was proposed in 1823. In most cases, the no-slip boundary condition is valid at a fluid–solid interface. However, with the advancement of measurement techniques, slip lengths ranging from nanometre to micrometre scales were experimentally measured, raising questions regarding the applicability of the no-slip condition. In 2003, Lauga & Stone (J. Fluid Mech., vol. 489, 2003, pp. 55–77) proposed a simple model to elucidate the effect of surface heterogeneities on the slip length, elegantly bridging the microscopic structure of the wall-boundary conditions to the macroscopic effective slip length.

A long-standing issue in pipe flow physics is whether the friction of the fluid follows a logarithmic or an algebraic decay. In 2005, McKeon et al. (J. Fluid Mech., vol. 538, 2005, pp. 429–443) published a detailed analysis of new measurements in the Princeton facility, and apparently settled the debate by showing that ‘the log is the law’. Almost 20 years later, no better data are presently available to reinforce their statement. Still, the story may not be totally over, and this is bad news for mathematicians who were hoping to get a long awaited final answer to one of their most elusive questions.

Since its publication in 2010, the paper by Schmid (J. Fluid Mech., vol. 656, 2010, pp. 5–28) has wielded considerable influence, an impact we examine here. That seminal work introduced dynamic mode decomposition, a method for performing flow-field spectral analysis of snapshot sequences of data. As a data-driven approach aimed at uncovering spatial and temporal patterns or modes within datasets, its applicability has extended far beyond fluid mechanics, reaching into a wide array of fields.

The nanoscale is the new frontier of fluid dynamics and its phenomenology can echo at the macroscale as in the canonical example of drop impact on a planar substrate. Unprecedented advances in measurement technology have recently equipped fluid dynamicists with the ability to probe nanoscale effects. The paper by Li et al. (J. Fluid Mech., vol. 785, 2015, R2) uses ultrafast imaging at the hundreds of nanoseconds scale to resolve the first contact between the drop and the substrate and thereby reveal the effect of prescribed nano-roughness on contact line motion.

Losses due to wake interactions between wind turbines can significantly reduce the power output of wind farms. The possibility of active flow control by wake deflection downstream of yawed horizontal-axis wind turbines has motivated research on the fluid mechanics involved. We summarize the findings of a wind tunnel study (Bastankhah & Porté-Agel, J. Fluid Mech., vol. 806, 2016, pp. 506–541) of the flow associated with a yawed model wind turbine, and the insights and modelling developments that have followed this important study.

In recent years, microfluidic systems have underpinned a wealth of biotechnology applications and proposed solutions for complex problems, including the sorting and enrichment of deformable particle suspensions. Motivated by such applications of microfluidic systems, Lu et al. (J. Fluid Mech., vol. 923, 2021, A11) present a three-dimensional computational study of a train of deformable capsules flowing through a branched microchannel. Insights into the intricacies of the underlying complex fluid–structure interactions between the suspended capsules and the surrounding fluid can inform experimental scenarios whereby strong capsule interactions are avoided, facilitating precise operating control of microfluidic devices for sorting and enrichment.

Microscopic irregularity (roughness) of bounding surfaces affects macroscopic dynamics of fluid flows. Its effect on bulk flow is usually quantified empirically by means of a roughness coefficient. A new approach, which treats rough surfaces (e.g. parallel plates) as random fields whose statistical properties can be inferred from measurements, is presented. The mapping of a random flow domain onto its deterministic counterpart, and the subsequent stochastic averaging of the transformed Stokes equations, yield expressions for the effective viscosity and roughness coefficient in terms of the statistical characteristics of the irregular geometry of the boundaries. The analytical nature of the solutions allows one to handle surface roughness characterized by short correlation lengths, a challenging feature for numerical stochastic simulations.

We present herein the derivation of a lubrication-mediated (LM) quasi-potential model for droplet rebounds off deep liquid baths, assuming the presence of a persistent dynamic air layer which acts as a lubricating pressure transfer. We then present numerical simulations of the LM model for axisymmetric rebounds of solid spheres and compare quantitatively to current results in the literature, including experimental data in the low-speed impact regime. In this regime the LM model has the advantage of being far more computationally tractable than direct numerical simulation (DNS) and is also able to provide detailed behaviour within the micro-metric thin lubrication region. The LM system has an interesting mathematical structure, with the lubrication layer providing a free-boundary elliptic problem mediating the drop and bath free-boundary evolutionary equations.

Turbulent flow induced by elastorotational instability in viscoelastic Taylor–Couette flow (TCF) with Keplerian rotation is analogous to a turbulent accretion disk destabilized by magnetorotational instability. We examine this novel viscoelastic Keplerian turbulence via direct numerical simulations (DNS) for the shear Reynolds number ($Re$) ranging from $10^2$ to $10^4$. The observed characteristic flow structure consists of penetrating streamwise vortices with axial length scales much smaller than the gap width, distinct from the classic centrifugally induced Taylor vortices, which have axial lengths of the gap width. These intriguing vortices persist for the wide $Re$ range considered and give rise to intriguing scaling behaviour in key flow quantities. Specifically, the characteristic axial length of the penetrating vortices is shown to scale as $Re^{-0.22}$; the angular momentum transport scales as $Re^{0.42}$; the kinetic and elastic boundary-layer thicknesses based on angular velocity and hoop stress near the inner cylinder wall scale as $Re^{-0.48}$ and $Re^{-0.49}$, respectively. This implies that the viscoelastic Keplerian turbulence belongs to the classical turbulent regime of TCF with the Prandtl–Blasius-type boundary layer. Furthermore, we present an analytical relation between the viscous and elastic dissipation rates of kinetic energy and the angular momentum transport and in turn demonstrate its validity using our DNS data. This study has paved the way for future research to explore astrophysics-related Keplerian turbulence and angular momentum transport via the scaling relations of the analogous TCF of dilute polymeric solutions.

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.