We recently lost a valuable colleague and friend, Paul Steen, and the text below is a brief overview of his life.

It has been reported that a fully localized turbulent band in channel flow becomes sustained when the Reynolds number is above a threshold. Here we show evidence that turbulent bands are of a transient nature instead. When the band length is controlled to be fixed, the lifetime of turbulent bands appears to be stochastic and exponentially distributed, a sign of a memoryless transient nature. Besides increasing with the Reynolds number, the mean lifetime also strongly increases with the band length. Given that the band length always changes over time in real channel flow, this size dependence may translate into a time dependence, which needs to be taken into account when clarifying the relationship between channel flow transition and the directed percolation universality class.

The transition to dripping in the gravity-driven flow of a liquid film under an inclined plate is investigated at zero Reynolds number. Computations are carried out on a periodic domain assuming either a fixed fluid volume or a fixed flow rate for a hierarchy of models: two lubrication models with either linearised curvature or full curvature (the LCM and FCM, respectively), and the full equations of Stokes flow. Of particular interest is the breakdown of travelling-wave solutions as the plate inclination angle is increased. For any fixed volume, the LCM reaches the horizontal state where it attains a cosine-shaped profile. For sufficiently small volume, the FCM and Stokes solutions attain a weak Young–Laplace equilibrium profile, the approach to which is described by an asymptotic analysis generalising that of Kalliadasis & Chang (J. Fluid Mech., vol. 261, 1994, pp. 135–168) for the LCM. For large volumes, the bifurcation curves for the FCM and Stokes model have a turning point so that the fully inverted state is never reached. For fixed flow rate, the LCM blows up at a critical angle that is well predicted by asymptotic analysis. The bifurcation curve for the FCM either has a turning point or else reaches a point at which the surface profile has an infinite slope singularity, indicating the onset of multi-valuedness. The latter is confirmed by the Stokes model, which can be continued to obtain overturning surface profiles. Overall, the thin-film models either provide an accurate prediction for dripping onset or else supply an upper bound on the critical inclination angle.

We consider a theoretical model for the settling of rod-shaped particles of a dilute, initially homogeneous, suspension in rapid rotation. The particle Reynolds number and the particle Taylor number of the detailed flow around the particles are assumed small, representing a relevant limit for an industrial centrifugal separation process. By applying a statistical approach using the Fokker–Planck equation, and neglecting particle–particle interactions, we obtain an explicit, analytical solution for the time dependent, spatially uniform particle orientation distribution function. Not only does the volume fraction in the bulk of the suspension decrease with time due to the divergent centrifugal field, as similarly described in the literature for suspensions of spherical particles, the orientation of the rod particles also changes with time from an initially uniform distribution to one where the particles tend to align with a plane perpendicular to the axis of rotation. The corresponding particle trajectories, as also influenced by first-order effects from the Coriolis acceleration and gyroscopic effects, are obtained numerically for different initial particle orientation angles.

Convection from a buoyancy source distributed over a vertical wall has diverse applications, from the natural ventilation of buildings to the melting of marine-terminating glaciers which impacts on future sea level. A key challenge involves determining how the rate and mechanisms of turbulent heat transfer should be extrapolated across a range of scales. Ke et al. (J. Fluid Mech., vol. 964, 2023, A24) explore transitions in the turbulent flow dynamics using direct numerical simulation of a convective boundary layer at a heated vertical wall. A classical regime of heat transfer, consistent with previous laboratory experiments, gives way with increasing accumulation of buoyancy to an ultimate regime with enhanced heat transfer. The key to this transition lies in a near-wall sublayer, with a switch from laminar buoyancy-driven dynamics to a sublayer dominated by turbulence and shear instability from the mean flow.

An important problem in passive scalar transport is to parametrize the effect of a fluctuating component of the flow, in order to overcome a limited resolution. A local effective diffusivity is one such parametrization, and over the years there have been many different suggestions for ‘closures’ that relate the advective flux to gradients of the mean concentration. Souza et al. (J. Fluid Mech., 2023, in press) introduce a stochastic framework where the local effective diffusivity is replaced by an exact effective diffusivity operator. By computing this operator for various examples, they quantify deviations from the local approximation, which can suggest areas of improvement and novel closure models.

The refraction of surface gravity waves by currents leads to spatial modulations in the wave field and, in particular, in the significant wave height. We examine this phenomenon in the case of waves scattered by a localised current feature, assuming (i) the smallness of the ratio between current velocity and wave group speed, and (ii) a swell-like, highly directional wave spectrum. We apply matched asymptotics to the equation governing the conservation of wave action in the four-dimensional position–wavenumber space. The resulting explicit formulas show that the modulations in wave action and significant wave height past the localised current are controlled by the vorticity of the current integrated along the primary direction of the swell. We assess the asymptotic predictions against numerical simulations using WAVEWATCH III for a Gaussian vortex. We also consider vortex dipoles to demonstrate the possibility of ‘vortex cloaking’ whereby certain currents have (asymptotically) no impact on the significant wave height. We discuss the role of the ratio of the two small parameters characterising assumptions (i) and (ii) above, and show that caustics are significant only for unrealistically large values of this ratio, corresponding to unrealistically narrow directional spectra.

Despite the widespread occurrence of pendant drops in nature, there is still a lack of combined studies on their dynamic and static stability. This study focuses on the dynamic and static stability of elongated drops with either a free or pinned contact line on a plane. We first examine static stability for both axisymmetric and non-axisymmetric perturbations subject to volume or pressure constraints. The stability limits for volume and pressure disturbances (axisymmetric) correspond to the maximum volume and pressure of the drops, respectively. Drops with free contact lines are marginally stable to non-axisymmetric perturbations because of their horizontal translational invariance, whereas pinned drops are stable. The linear dynamic stability is then investigated numerically through a boundary element model, restricted to volume disturbances. Results show that when the stability limit is reached, the first zonal mode has a zero frequency, suggesting that the thresholds for static and dynamic stability are essentially equivalent. Furthermore, natural frequencies experience sharp changes as the stability limit is approached. Another zero frequency mode associated with the horizontal motion of the centre of mass is also revealed by the numerical results, reflecting the horizontal translational invariance of drops with free contact lines. Finally, the frequency spectrum modified by gravity is explored, resulting in the identification of five gravity-induced frequency shift patterns. The frequency shifts break the spectral degeneracy for hemispherical drops with free contact lines, leading to various spectral orderings according to polar and azimuthal wavenumbers.

The secondary motion caused by turbulence anisotropy is one of the crucial factors for determining the size of corner-flow separation in a side-wall interference flow field. Therefore, through a wall-resolved large-eddy simulation (LES) of a side-wall interference flow field, this study investigates the effects of the secondary motion on the corner-flow separation and explores the turbulence modelling that can reproduce the secondary flow motion. The momentum transport analysis using the LES results shows that the secondary vortex has twofold effects on delaying the corner-flow separation: the convective transport of the streamwise momentum towards the corner, and the enhanced production of turbulence by increasing the shear. Also, the vorticity transport analysis reconfirms that the secondary motion is caused primarily by turbulence anisotropy in the outer layer of the turbulent boundary layer. Furthermore, a quadratic constitutive relation (QCR) is proposed based on the analysis of the relationship between the Reynolds stress and velocity gradient. The proposed QCR consists of two quadratic terms and three constant parameters. The a priori analysis using the LES data shows that the proposed QCR represents the anisotropy of the Reynolds stress overall better than the existing QCR. Reynolds-averaged Navier–Stokes simulation using the proposed QCR with the Spalart–Allmaras turbulence model shows improvements in the prediction of the corner-flow separation compared to the results obtained by the existing QCR with the same turbulence model.

Inclining a fluid layer subjected to a temperature gradient introduces a profusion of fascinating patterns and regimes. Previous experimental and computational studies form the starting point for an extensive numerical bifurcation study by Reetz & Schneider (J. Fluid Mech., vol. 898, 2020, A22) and Reetz et al. (J. Fluid Mech., vol. 898, 2020, A23). Intricate trajectories passing through multiple steady and periodic states organize the dynamics. The consequences for chaotic patterns in large geometries is discussed.

This study presents a buoyancy-driven stability analysis in a three-dimensional inclined porous medium with a capillary transition zone that is formed between a non-wetting and an underlying wetting phase. In this two-phase, two-component, partially miscible system, a solute from a non-wetting phase diffuses into a porous layer saturated with a wetting-phase fluid, creating a dense diffusive boundary layer beneath an established capillary transition zone. Transient concentration and gravity-driven velocity fields are derived for the wetting phase while the saturation field remains fixed. Linear stability analysis with the quasi-steady-state approximation is employed to determine the onset of solutal convective instability for buoyancy-dominant, in-transition and capillary-dominant systems. The analysis of the problem leads to a differential eigenvalue problem composed of a system of three complex-valued equations that are numerically solved to determine the critical times, critical wavenumbers and neutral stability curves as a function of inclination angle for different Bond numbers. The layer inclination is shown to play an essential role in the stability of the problem, where the gravity-driven flow removes solute concentrations in the diffusive boundary layer. The results indicate that the horizontal porous layer exhibits the fastest onset of instability, and longitudinal rolls are always more unstable than oblique and transverse rolls. The inclination angle has a more substantial impact on stabilizing the diffusive boundary layer in the buoyancy-dominant than in the capillary-dominant systems. Furthermore, for both buoyancy-dominant and capillary-dominant systems, the critical times and wavenumbers vary exponentially with inclination angle ≤ 60° and follow the Stirling model.

We experimentally investigate the transport of liquid by a toroidal bubble that rises vertically and penetrates a horizontal interface between two immiscible liquids. Experiments are conducted with various strengths of vortex circulation in the bubble, and with different liquid densities and viscosities. In contrast to a spherical bubble, a rising toroidal bubble carries a great amount of the lower liquid by virtue of the self-induction of circulating flow. The lower liquid is entrained by the toroidal bubble and forms an ellipsoidal body enclosing the bubble after it penetrates the interface. The downward net force acting on the ellipsoidal body results in the radial contraction of the bubble, reducing the volume of the entrained lower liquid. As the entrained volume decreases, the nearby upper liquid eventually pierces the ellipsoidal body, making the bubble unstable. At this instant, the net force acting on the ellipsoidal body approaches zero, and the volume ratio of the entrained lower liquid and bubble converges to a specific value. For smaller vortex circulation and larger density difference between the liquids, the volume of the entrained lower liquid within the ellipsoidal body becomes smaller and the travel distance of the bubble from the initial interface until it becomes unstable decreases. The effective Froude number, which accounts for both the inertial effect of vortex circulation and the gravitational effect of liquid density difference, is found to characterise the temporal changes in the ring radius, propagation speed and entrained volume.

We study the gravity-induced collisions of charged spheres of dielectric materials dispersed in a gaseous medium. When the gap thickness between the surfaces of two spheres is shorter than the mean free path of the surrounding fluid medium, continuum assumptions for the hydrodynamics interactions are no longer valid, and the non-continuum lubrication interactions result in surface-to-surface contact in finite time. Two like-charged dielectric spheres attract each other at close separations for a wide range of size and charge ratio values. We use trajectory analysis to calculate the collision rate and, thus, explore the role of electrostatic interactions in the collision dynamics of a pair of like-charged dielectric spheres. We present the modifications of pair trajectories due to electrostatic forces and show how collision efficiencies vary with the non-dimensional parameter capturing the relative strength of the electrostatic force to gravity as well as the charge ratio and size ratio.

In this study, the effects of antagonistic muscle actuation on the propulsion of a bilaminar-structure fish fin ray were investigated using a two-dimensional computational flow–structure interaction (FSI) model. The structure and material properties of the model were based on the realistic biological data of the sunfish fin. The effect of muscle actuation was modelled using root displacement offset between the two hemitrichs. Parametric FSI simulations were conducted by assuming a sinusoidal function of the offset over a cycle and varying the amplitude and phase difference between the actuations and pitching/plunging motions. The results show that the phase of muscle actuation is a critical factor affecting its effects. Three performance regions can be identified with different phase ranges, including a thrust-favour region, an efficiency-favour region and a thrust-efficiency-unfavour region. In each region, the relationships among the root actuations, fin-ray kinematics, vortex dynamics and resulting performance are studied and discussed. Furthermore, a strong positive correlation between the trailing–leading amplitude ratio and thrust coefficient as well as a negative relationship between the efficiency and angle of attack at the centre of mass of the fin ray are observed.

The hydrodynamic influence of surface texture on static surfaces ranges from large drag penalties (roughness) to potential performance benefits (shark-like skin). Although it is of wide-ranging research interest, the impact of roughness on flapping systems has received limited attention. In this work, we explore the effect of roughness on the unsteady performance of a harmonically pitching foil through experiments using foils with different surface roughness, at a fixed Strouhal number and within the Reynolds number ($Re$) range of $17\,000\unicode{x2013}33\,000$. The foils’ surface roughness is altered by changing the distribution of spherical-cap-shaped elements over the propulsor area. We find that the addition of surface roughness does not improve the performance compared with a smooth surface over the $Re$ range considered. The analysis of the flow fields shows near-identical wakes regardless of the foil's surface roughness. The performance reduction mainly occurs due to an increase in profile drag. However, we find that the drag penalty due to roughness is reduced from $76\,\%$ for a static foil to $16\,\%$ for a flapping foil at the same mean angle of attack, with the strongest decrease measured at the highest $Re$. Our findings highlight that the effect of roughness on dynamic systems is very different than that on static systems; thereby, it cannot be estimated by only using information obtained from static cases. This also indicates that the performance of unsteady, flapping systems is more robust to the changes in surface roughness.

The average wave drag in unsteady motion is studied experimentally with force measurements. Towing hulls of size $L$ at sinusoidal speed, the mean drag is measured for different amplitudes and frequencies of the fluctuating velocity, as well as different Froude numbers $\mathcal {F}_0$ associated with the mean velocity $V_0$ ($\mathcal {F}_0 = V_0/\sqrt {gL}$). The wave drag is reported to be either increased or decreased by velocity fluctuations depending on $\mathcal {F}_0$. For small fluctuation amplitudes, this drag change is proportional to the square of the amplitude. The effect is maximized for a resonance frequency identified as the Wehausen frequency, which scales as $\sqrt {g/L}$ times the inverse of the Froude number. All these results are rationalized by developing an extension to Havelock's theory.

In this paper, a high-level Green–Naghdi (HLGN) model for large-amplitude internal waves in a two-layer fluid system, where the upper-fluid layer is of finite depth and the lower-fluid layer is of infinite depth, is developed under the rigid-lid free-surface approximation. The equations of the present HLGN model follow Euler's equations under the sole assumption that the horizontal and vertical velocity distributions along the vertical column are presented by known shape functions for each layer. The linear dispersion relations of the HLGN model for different levels are presented and compared with those obtained by other strongly nonlinear models for deep water, including the fully nonlinear models that include the dispersion effects $O(\mu )$ (where $\mu$ is the ratio of the upper-fluid layer depth to a typical wavelength) derived by Choi & Camassa (Phys. Rev. Lett., vol. 77, 1996, pp. 1759–1762) and $O(\mu ^2)$ derived by Debsarma et al. (J. Fluid Mech., vol. 654, 2010, pp. 281–303). It is shown that the HLGN model has a wider application range than other models. Solutions of travelling large-amplitude internal solitary waves in the absence and presence of background shear-current are then investigated by using the HLGN model. For the no-current cases, results obtained by the HLGN model show better agreement with Euler's solution on wave profile, velocity profile at the maximum interface displacement and wave speed compared with those obtained by other models. For the background shear-current cases, results obtained by the HLGN model also show good agreement with those obtained by solving the Dubreil-Jacotin–Long equation.

The tangential strain rate in premixed flames impacts significantly the flame surface area generation and thus the combustion process. Studies on incompressible isotropic turbulence have revealed that the mean tangential strain rate at material and iso-scalar surfaces is positive and exhibits a universal value when normalized by the Kolmogorov time. This is associated with the preferential alignment of the surface normal with the most compressive principal strain rate. The present study investigates such effects in premixed hydrogen and iso-octane flame kernels using direct numerical simulations. It is shown that the normalized mean tangential strain rate of the investigated flames has a very similar value compared with the incompressible flows. However, in the reaction zone, the flame surface normal aligns preferentially with the most extensive principal strain rate. Furthermore, this alignment depends on the reaction progress variable and the Lewis number, while the tangential strain rate remains independent of these parameters. Such counter-intuitive behaviour is systematically investigated by decomposing the effects of dilatation and residual solenoidal turbulence. It is found that the solenoidal turbulence influences significantly the tangential strain rate. A general effect of turbulence on the tangential strain rate is identified, which is consistent with incompressible flows and independent of the Lewis number and the reaction progress variable. This is a remarkable finding indicating that models of the tangential strain rate developed based on incompressible flows apply also to premixed flames with different Lewis numbers, and, for the modelling, only the solenoidal turbulence should be considered.