We experimentally investigate the rotational dynamics of neutrally buoyant flat bodies of revolution (spheroids, disks and rings with different cross-sectional shapes) in shear flows. In the Stokes regime, the axis of revolution of these rigid particles moves in one of a family of closed periodic Jeffery orbits. Inertia is able to lift the orbit degeneracy and induces drift among several rotations towards limiting stable orbits. Furthermore, permanent alignment can be achieved for disks and rings with triangular cross-sectional shapes, provided the inertia is sufficiently high. The bifurcations between the different dynamics are compared with those predicted by small-inertia asymptotic theories and numerical simulations.

Fingering instabilities readily occur if a less viscous fluid displaces a more viscous fluid in a narrow gap due to the action of destabilising viscous forces. If the fluids are miscible, the instability can be suppressed in the limit of large advection as complicated flow structures are formed across the gap. Using a fluid to displace a monolayer of non-colloidal particles suspended in the same fluid, Luo et al. (2025 J. Fluid Mech. vol. 1011, A48) suppress the formation of the cross-gap structures and identify a new fingering mechanism which instead relies on long-range dipolar disturbance flows generated by the particle confinement.

Dynamics of a spherical particle and the suspending low-Reynolds-number fluid confined between two concentric spherical walls were studied numerically. We calculated the particle’s hydrodynamic mobilities at various locations in the confined space. It was observed that the mobility is largest near the middle of confined space along the radial direction, and decays as the particle becomes closer to no-slip walls. At a certain confinement level, the maximal mobility occurs near the inner wall. We also calculated the drift velocity of the particle perpendicular to the external force. The magnitude of the drift velocity normalised by the velocity along the external force was found to depend on particle location and the confinement level; it is observed that the maximal drift velocity occurs near the wall. Fluid vortices in the confined space induced by particle motion were observed and analysed. In addition, we studied particle trajectories in the flow when the walls rotate at constant angular velocities. The externally applied force, rotation-induced flow and centrifugal/centripetal force, and particle–wall interaction lead to various modes of particle motion. This work lays the foundation to understand and manipulate particulate transport in microfluidic applications such as intracellular transport and encapsulation technologies.

Large-scale spanwise motions in shock wave–turbulent boundary-layer interactions over a $ 25^{\circ }$ compression ramp at Mach 2.95 are investigated using large-eddy simulations. Spectral proper orthogonal decomposition (SPOD) identifies coherent structures characterised by low-frequency features and a large-scale spanwise wavelength of $ O(15\delta _{0})$, where $ \delta _{0}$ is the incoming boundary-layer thickness. The dominant frequency is at least one order of magnitude lower than that of the shock motions. These large-scale spanwise structures are excited near the shock foot and are sustained along the separation shock. Global stability analysis (GSA) is then employed to investigate the potential mechanisms driving these structures. The GSA identifies a stationary three-dimensional (3-D) mode at a wavelength of $ 15\delta _{0}$ with a similar perturbation field, particularly near the separation shock. Good agreement is achieved between the leading SPOD mode and the 3-D GSA mode both qualitatively and quantitatively, which indicates that global instability is primarily responsible for the large-scale spanwise structures surrounding the shock. The reconstructed turbulent separation bubble (TSB) using the 3-D global mode manifests as spanwise undulations, which directly induce the spanwise rippling of the separation shock. Furthermore, the coupled TSB motions in the streamwise and spanwise directions are examined. The TSB oscillates in the streamwise direction while simultaneously exhibiting spanwise undulations. The filtered wall-pressure signals indicate the dominant role of the streamwise motions.

This paper explores the role of barodiffusion in the dynamics of gas bubble growth in highly viscous gas-saturated magma subjected to instant decompression. A mathematical model describing the growth of a single isolated bubble is formulated in terms of the modified Rayleigh–Plesset equation coupled with the mass transfer and material balance equations. The model simultaneously takes into account both dynamic and diffusion mechanisms, including the effect of barodiffusion caused by emergence of a large pressure gradient in the liquid, which, in turn, is associated with formation of a diffusion boundary layer around the bubble. An analytical solution of the problem is found, the construction of which is based on the existence of a quasi-stationary state of the bubble growth process. It is shown that barodiffusion manifests itself at the initial and transient stages and under certain conditions can play a paramount role.

Natural fliers and marine swimmers twist and turn their lifting or control surfaces to manipulate the unsteady forces experienced in air and water. The passive deformation of such surfaces has been investigated by several researchers, but the aspect of controlled deformation has received comparatively less attention. In this paper, we experimentally measure the forces and the flow fields of a flat-plate wing (aspect ratio (AR) = 3), translating at a constant Reynolds number (Re) of 10 000, with a dynamically twisting span. We show that the unsteady forces can be dependably estimated by a three-dimensional discrete vortex model. In this model, we account for the leading-edge separation with the help of the leading-edge-suction parameter. Experiments are conducted for two angles of attack (AoAs), $5^\circ$ and $15^\circ$. In addition, two rates of twisting are implemented where part of the leading edge, closer to the tip region, is twisted away from the incoming flow, increasing the effective AoA. The results show that twisting away from the flow augments the lift forces in all cases, although the rate of increase of lift is higher for the highest twist rate. The act of twisting causes an increase in effective AoA beyond the static stall angle in the AoA $=15^\circ$ case. This is highlighted by a distinct dip in the force data following the initial rise after twisting is activated. The increase in effective AoA from the reference case (without twisting) causes separation of the flow below the mid-span. This, in turn, creates higher levels of vorticity in those regions and results in a leading-edge vortex with increased cross-section and strength when compared with the reference case without twisting. Finally, we apply force partitioning and reveal that dynamic twisting leads to a localised increase in vorticity-induced forces along the twisted part of the span, which is approximately twice that of the untwisted case.

Turbulent flows are strongly chaotic and unpredictable, with a Lyapunov exponent that increases with the Reynolds number. Here, we study the chaoticity of the surface quasi-geostrophic system, a two-dimensional model for geophysical flows that displays a direct cascade similar to that of three-dimensional turbulence. Using high-resolution direct numerical simulations, we investigate the dependence of the Lyapunov exponent on the Reynolds number and find an anomalous scaling exponent larger than that predicted by dimensional arguments. We also study the finite-time fluctuation of the Lyapunov exponent by computing the Cramér function associated with its probability distribution. We find that the Cramér function attains a self-similar form at large $\textit{Re}$.

This study presents a numerical investigation of wall-mounted tandem flexible plates with unequal lengths in a laminar boundary layer flow, examining both two-dimensional (2-D) and three-dimensional (3-D) configurations. Key parameters influencing the system include the plate’s bending stiffness ($K$), Reynolds number (${Re}$) and length ratio ($L^*$). Five motion modes are identified: dual collapse (DC), flapping collapse (FC), dual flapping (DF), static flapping (SF) and dual static (DS). A phase diagram in the ($K,L^*$) space is constructed to illustrate their regimes. We focus on DF and SF modes, which significantly amplify oscillations in the downstream plate – critical for energy harvesting. These amplification mechanisms are classified into externally driven and self-induced modes, with the self-induced mechanism, which maximises the downstream plate’s amplitude, being the main focus of our study. A rigid–flexible (RF) configuration is introduced by setting the upstream plate as rigid, showing enhanced performance at high ${Re}$, with oscillation amplitudes up to 100 % larger than the isolated flexible (IF) plate configuration. A relation is developed to explain these results, relating oscillation amplitude to trailing-edge velocity, oscillation frequency and chord length. Force analysis reveals that the RF configuration outperforms both IF and flexible–flexible (FF) configurations. Unlike frequency lock-in, the RF configuration exhibits frequency unlocking, following a $-2/3$ scaling law between the Strouhal number ($St$) and ${Re}$. Results from the 3-D RF configuration confirm that the 2-D model remains applicable, with the self-induced amplification mechanism persisting in 3-D scenarios. These findings enhance understanding of fluid–structure interactions, and offer valuable insights for designing efficient energy harvesting systems.

This paper considers the propagation, arrest and recession of a planar hydraulic fracture in a porous elastic medium whose footprint is constrained to a growing or shrinking rectangular region with a constant height. Hydraulic fractures with large aspect ratio rectangular footprints are frequently referred to as PKN fractures in recognition of the original researchers (Perkins & Kern 1961 J. Petrol. Tech. 13, 937–949) and (Nordgren 1972 J. Petrol Technol. 1972, 306–314) who first analyzed models of such fracture geometries. We investigate the one-dimensional non-local PKN approximation to a fully planar rectangular hydraulic fracture model in a three-dimensional elastic medium. By analysing the tip behaviour of the non-local PKN model, a transformation procedure is established to render the asymptotic equations for the dynamics of the steady semi-infinite PKN and plane strain models formally identical, which implies that all the existing multiscale plane strain asymptotes can be converted directly to the PKN case by making use of this transformation. Using this transformation, it is shown that the appropriate PKN asymptotes for the average aperture $\bar {w}$ with distance $\hat {x}$ to the fracture front are $\bar {w}\sim \hat {x}^{1/2},\ \hat {x}^{5/8}\ {\textrm{and}}\, \ \hat {x}^{2/3}$ in the toughness, leak-off and viscous modes of propagation, respectively; as well as the linear elastic fracture mechanics tip asymptote $\bar {w}\sim \hat {x}^{1/2}$ for arrest, which transitions to the linear asymptote tip $\bar {w}\sim \hat {x}$ for a fracture driven to recede due to fluid leak-off. Both the arrest and recession tip asymptotes share the intermediate leak-off asymptote $\bar {w}\sim \hat {x}^{3/4}$. A scaling analysis yields the arrest time, length and aperture as functions of a dimensionless injection-cessation time $\omega$. An asymptotic analysis of the non-local PKN model is used to establish the fundamental decoupling between dynamics and kinematics, which leads to the emergence of a similarity solution – termed the sunset solution – close to the time of collapse of the fracture. The multiscale PKN numerical solutions agree well with those for a fully planar multiscale rectangular hydraulic fracture model in a three-dimensional elastic medium. The scaling laws and the emergence of the sunset solution are confirmed by the PKN numerical model. The sunset solution also emerges in the fully planar numerical model and persists beyond the collapse time of the PKN model, by which time its footprints have separated from the upper and lower constraining sedimentary layer boundaries and have assumed self-similar elliptic shapes that shrink as they approach collapse.

In the rapidly rotating limit, we derive a balanced set of reduced equations governing the strongly nonlinear development of the convective wall-mode instability in the interior of a general container. The model illustrates that wall-mode convection is a multiscale phenomenon where the dynamics of the bulk interior diagnostically determine the small-scale dynamics within Stewartson boundary layers at the sidewalls. The sidewall boundary layers feedback on the interior via a nonlinear lateral heat-flux boundary condition, providing a closed system. Outside the asymptotically thin boundary layer, the convective modes connect to a dynamical interior that maintains scales set by the domain geometry. In many ways, the final system of equations resembles boundary-forced planetary geostrophic baroclinic dynamics coupled with barotropic quasi-geostrophic vorticity. The reduced system contains the results from previous linear instability theory but captured in an elementary fashion, providing a new avenue for investigating wall-mode convection in the strongly nonlinear regime. We also derive the dominant Ekman-flux correction to the onset Rayleigh number for large Taylor number, ${\textit {Ra}} \approx 31.8 \,{\textit{Ta}}^{1/2} - 4.43 \,{\textit{Ta}}^{5/12} + {\mathcal{O}}({\textit{Ta}}^{1/3})$ for no-slip boundaries. However, we find that the linear onset in a finite cylinder differs noticeably compared with a Cartesian channel. We demonstrate some of the reduced model’s nonlinear dynamics with numerical simulations in a cylindrical container.

This study employs a direct numerical simulation method to investigate the wake pattern evolutions of flows past an insulated spheroid and provides expressions of force and torque coefficients influenced by a streamwise magnetic field in an incompressible, conducting, viscous fluid. A total of 1150 cases are examined covering a parameter range of Reynolds number $50 \leqslant \textit{Re} \leqslant 250$, aspect ratio $1.5 \leqslant \beta \leqslant 6$, inclination angle $0^\circ \leqslant \theta \leqslant 90^\circ$, and interaction parameter $0 \leqslant N \leqslant 10$, where $\beta$ and $N$, respectively, reflect the anisotropy of the spheroid and the strength of magnetic field. Nine wake patterns are classified based on wake structure features and summarised in three maps of regimes according to the inclination angle. The transition mechanisms among these wake patterns are also investigated under the influence of a streamwise magnetic field. Furthermore, expressions for drag, lift and torque coefficients are derived with the help of three fundamental physical criteria. Results indicate that the force and torque expressions give a good prediction within the present parameter space $\{\textit{Re}, \beta , \theta , N\}$.

In conventional hypersonic wind tunnels, tunnel noise is dominated by acoustic radiation from turbulent nozzle-wall boundary layers, which can directly influence the boundary-layer transition (BLT) over the model in the test section. To offer new insights into BLT in conventional ground facilities, direct numerical simulations (DNS) were performed to simulate the receptivity and transition processes of a Mach 8 boundary layer over a nearly sharp $7^\circ$ half-angle cone, with transition triggered by tunnel-like broadband free-stream acoustic disturbances radiated from the nozzle wall of the Sandia hypersonic wind tunnel at Mach 8 (Sandia HWT-8). The DNS captured all the stages of the transition to turbulence caused by tunnel noise, including the passage of broadband free-stream noise through the shock wave, the receptivity process leading to the generation of Mack’s second-mode waves, their nonlinear growth to saturation, the laminar breakdown to turbulence and the post-transitional, fully turbulent flow. The transition location predicted by DNS compared well with that of Pate’s theory and was also consistent with the locations of peak pressure fluctuations as measured in the Sandia HWT-8 facility. The computed skin friction and Stanton number distributions in the initial breakdown region showed an overshoot compared with the turbulent predictions by the van Driest II theory. The wall-pressure spectra in both the transitional and turbulent regions of the cone compared well with those measured in the Sandia HWT-8. The second-mode breakdown amplitude $A_{max}$ predicted by the DNS was also consistent with sharp-cone measurements from multiple conventional wind tunnels.