This paper addresses the stability of plane Couette flow in the presence of strong density and viscosity stratifications. It demonstrates the existence of a generalised inflection point that satisfies the generalised Fjørtoft criterion of instability when a minimum of kinematic viscosity is present in the base flow. The characteristic scales associated with this minimum are identified as the primary controlling parameters of the associated instability, regardless of the type of stratification. To support this finding, analytical stability models are derived in the long-wave approximation using piecewise linear base flows. Numerical stability calculations are carried out to validate these models and to provide further information on the production of disturbance vorticity. All instabilities are interpreted as arising from the interaction between two vorticity waves. Depending on the type of stratification, these two waves are produced by different physical mechanisms. When both strong density and viscosity stratifications are present, we show that they result from the concurrent action of shear and inertial baroclinic effects. The stability models developed for simple fluid models ultimately shed light on a recently observed unstable mode in supercritical fluids (Ren et al., J. Fluid Mech., vol. 871, 2019, pp. 831–864), providing a quantitative prediction of the stability diagram and identifying the dominant mechanisms at play. Furthermore, our study suggests that the minimum of kinematic viscosity reached at the Widom line in these fluids is the leading cause of their instability. The existence of similar instabilities in different fluids and flows (e.g. miscible fluids) is finally discussed.

A study of the physics of separating and reattaching flows around bodies with sharp edges is reported. Data from direct numerical simulations of the flow around a rectangular cylinder with aspect ratio 5 at different Reynolds numbers are used. The flow is decomposed into multiple interacting flow phenomena such as the laminar boundary layer in the front face, the separated shear layer, the flow impingement at reattachment, the reverse boundary layer within the recirculating bubble and the near- and far-wake flow. A detailed analysis of the physics of these phenomena is provided, including the slow modulation induced by large-scale instabilities related with vortex shedding. The entrainment phenomena acting along the separated shear layer and their unbalance between its inner and outer sides are recognised as fundamental mechanisms determining the tendency of the flow to reattach and the overall fluxes of momentum and heat. The behaviour of entrainment is found to be strictly related with the shear-layer velocity difference that in turn is determined by the behaviour of the reverse boundary layer and by its strength in counteract adverse pressure gradients. The physical understanding of the compound role played by these and all the other mechanisms composing the flow, poses the basis for the formulation of theoretical frameworks able to unify all these interacting phenomena. Finally, the present work provides access to high-fidelity flow statistics of relevance for benchmark activities on bluff bodies with sharp edges.

We investigate the linear instability of two-layer stratified shear flows in a sloping two-dimensional channel, subject to non-zero longitudinal gravitational forces. We reveal three previously unknown instabilities, distinct from the well-known Kelvin–Helmholtz instability and Holmboe wave instability, in that they have longer wavelengths (of the order of 10 to $10^3$ shear-layer depths) and often slower growth rates. Importantly, they can grow in background flows with gradient Richardson number $\gg 1$, which offers a new mechanism to sustain turbulence and mixing in strongly stratified flows. These instabilities are shown to be generic and relatively insensitive to Reynolds number, Prandtl number, base flow profile and boundary conditions. The nonlinear evolution of these instabilities is investigated through a forced direct numerical simulation, in which the background momentum and density are sustained. The growth of long unstable waves in background flows initially stable to short wave causes a decrease in the local gradient Richardson number. This leads to local nonlinear processes that result in small-scale overturns resembling Kelvin–Helmholtz billows. Our results establish a new energy exchange pathway, where the mean kinetic energy of a strongly stratified flow is extracted by primary unstable long waves and secondary short waves, and subsequently dissipated into internal energy.

Theoretical analysis and numerical results have shown that frequency lock-in in vortex-induced vibration (VIV) is caused by the instability of the structural mode rather than a resonant response to external excitations. However, there is a lack of experimental evidence supporting relevant theoretical research findings. This study investigates VIV suppression with a passive modal controller (PMC) for a circular cylinder at Reynolds numbers $Re = 60$ and $Re = 40$, using experiments to distinguish the effects of stable and unstable wake modes. Comparative analysis before and after the implementation of the PMC reveals significant reduction in the vibration amplitude and the disappearance of the lock-in phenomenon at $Re = 60$. The vibration frequency closely follows the vortex shedding frequency after control, while dynamic mode decomposition of the flow field indicates that the wake mode is dominant. For $Re = 40$, the vibration is eliminated and the flow becomes steady. Additionally, the root loci of the coupled system are investigated before and after the PMC implementation via linear stability analysis. The results indicate that the PMC can alter the dynamic characteristics of the original system, causing the structural mode and PMC mode to couple when approaching the PMC frequency. Then, the interaction typically improves the stability of the structural mode. Finally, a parametric study is conducted in the experiment, as well as a linear stability analysis. The study provides experimental evidence that stability control of the structural mode is the key to suppressing VIV and eliminating the lock-in phenomenon.

The dynamics of fluid-conveying pipelines with different shapes has received extensive research attention. Significant wall shear stress and flow separation occur when the fluid flows through pipelines with various curvatures. These phenomena trigger pipeline vibration, the generation of mechanical and hydrodynamic noise, damage, and even the rupture of the pipeline. However, previous studies have not considered the mechanism of internal pipeline flow to eliminate flow separation and the generation of secondary flow inside bent pipelines by redirecting and manipulating the flow. To steer the fluid flow, a ‘hydrodynamic transformation strategy’ based on the metamaterial technology is proposed for the first time in this work; through this strategy, the fluid in pipelines can be made to flow along trajectories that are always parallel to the central axis of the bent pipelines. Interestingly, this innovative method can effectively eliminate the elbow-induced secondary flow and prevent the generation of a pressure gradient toward the pipeline wall. Using the soft lithography technology or the three-dimensional printing technology, the hydrodynamic metamaterial microstructure required to manipulate the fluid flow path in actual engineering applications can be achieved. Our work paves the way for developing new approaches for controlling the flow characteristics and reducing the turbulence intensity of the fluid flowing in pipelines with elbows and corners.

An experimental investigation identifying the effects of surface roughness on the drag coefficient ($C_{D}$) of freely rolling spheres is reported. Although lubrication theory predicts an infinite drag force for an ideally smooth sphere in contact with a smooth wall, finite drag coefficients are obtained in experiments. It is proposed that surface roughness provides a finite effective gap ($G$) between the sphere and panel, resulting in a finite drag force while also allowing physical contact between the sphere and plane. The measured surface roughnesses of both the sphere and panel are combined to give a total relative roughness ($\xi$). The measured $C_{D}$ increases with decreasing $\xi$, in agreement with analytical predictions. Furthermore, the measured $C_{D}$ is also in good agreement with the combined analytical and numerical predictions for a smooth sphere and wall, with a gap approximately equal to the root-mean-square roughness ($R_q$). The accuracy of these predictions decreases for low mean Reynolds numbers ($\overline {Re}$), due to the existence of multiple scales of surface roughness that are not effectively captured by $R_{q}$. Experimental flow visualisations have been used to identify critical flow transitions that have been previously predicted numerically. Path tracking of spheres rolling on two panels with different surface roughnesses indicates that surface roughness does not significantly affect the sphere path or oscillations. Analysis of sphere Strouhal number ($St$) highlights that wake shedding and sphere oscillations are coupled at low $\overline {Re}$ but with increasing $\overline {Re}$, the influence of wake shedding on the sphere path diminishes.

Two-dimensional oblique detonation wave (ODW) propagations in partially prevaporized n-heptane sprays are numerically simulated with a skeletal chemical mechanism. The influences of the droplet diameter and total equivalence on oblique detonation are considered. The initiation length is found to increase first and then decrease with increasing initial droplet diameter, and the effect of droplet size is maximized when the initial droplet diameter is approximately $10\ \mathrm {\mu } {\rm m}$. As the initial droplet diameter varies, unsteady and steady ODWs are observed. In the cases of unsteady ODWs, temperature gradients and non-uniform distributions of the reactant mixture due to droplet evaporation lead to formation of unsteady detonation propagation, therefore leading to fluctuations in the initiation length. The fluctuations in initiation length decrease as the pre-evaporation gas equivalence ratio increases for the unsteady cases. The results further suggest that the relationship between the evaporation layer thickness along the streamline and the corresponding theoretical initiation length can be used to identify an unsteady or steady ODW in cases with large droplets that evaporate behind an oblique shock wave or ODW under the effects of different initial droplet diameters.

This paper presents a novel Hamiltonian formulation of the isotropic Navier–Stokes problem based on a minimum-action principle derived from the principle of least squares. This formulation uses the velocities $u_{i}(x_{j},t)$ and pressure $p(x_{j},t)$ as the field quantities to be varied, along with canonically conjugate momenta deduced from the analysis. From these, a conserved Hamiltonian functional $H^{*}$ satisfying Hamilton's canonical equations is constructed, and the associated Hamilton–Jacobi equation is formulated for both compressible and incompressible flows. This Hamilton–Jacobi equation reduces the problem of finding four separate field quantities ($u_{i}$,$p$) to that of finding a single scalar functional in those fields – Hamilton's principal functional ${S}^{*}[u_{i},p,t]$. Moreover, the transformation theory of Hamilton and Jacobi now provides a prescribed recipe for solving the Navier–Stokes problem: find ${S}^{*}$. If an analytical expression for ${S}^{*}$ can be obtained, it will lead via canonical transformation to a new set of fields which are simply equal to their initial values, giving analytical expressions for the original velocity and pressure fields. Failing that, if one can only show that a complete solution to this Hamilton–Jacobi equation does or does not exist, that will also resolve the question of existence of solutions. The method employed here is not specific to the Navier–Stokes problem or even to classical mechanics, and can be applied to any traditionally non-Hamiltonian problem.

We introduce a numerical strategy to study the evolution of two-dimensional water waves in the presence of a plunging jet. The free-surface Navier–Stokes solution is obtained with a finite, but small, viscosity. We observe the formation of a surface boundary layer where the vorticity is localised. We highlight convergence to the inviscid solution. The effects of dissipation on the development of a singularity at the tip of the wave is also investigated by characterising the vorticity boundary layer appearing near the interface.

Shear significantly influences turbulence in the energy-containing range of shear-dominated flows, and the longitudinal structure functions do not have a universal form as they do in homogeneous isotropic turbulence. Despite this, the relative scaling of structure functions exhibits universal sub-Gaussian behaviour in shear-dominated flows, in particular for turbulent boundary layers, channels and Taylor–Couette flows. Our investigation of a turbulent vertical buoyancy layer at $Pr = 0.71$ using direct numerical simulation shows this universality even in moderate-Reynolds-number buoyancy-driven but shear-dominated boundary layers. It is demonstrated that the universality is related to the energy density of the eddies, which attains a hierarchical equilibrium in the energy-containing range of shear-dominated turbulence. We conjecture that the universal sub-Gaussian behaviour of the energy density of the energy-containing range, which was considered to be non-trivial in prior studies, is related to the universal anomalous scaling exponents of the inertial subrange turbulence. Based on this conjecture, we propose a hypothesis that relates large-scale eddies and the intermittent dissipation field in shear-dominated turbulence, highlighting a relationship between large and small scales. A phenomenological model is also developed to predict the scaling, which is verified using data from a turbulent boundary layer, half-channel and vertical buoyancy layer at friction Reynolds numbers spanning four orders of magnitude. Excellent agreement is observed.

A new model of the steady boundary layer flow around a rotating sphere is developed that includes the widely observed collision and subsequent eruption of boundary layers at the equator. This is derived following the Segalini & Garrett (J. Fluid Mech., vol. 818, 2017, pp. 288–318) asymptotic approach for large Reynolds numbers but replacing the Smith & Duck (Q. J. Mech. Appl. Maths, vol. 30, issue 2, 1977, pp. 143–156) correction with a higher-order version of the Stewartson (Grenzschichtforschung/Boundary Layer Research, 1958, pp. 59–71. Springer) model of the equatorial flow. The Stewartson model is then numerically solved, for the first time, via a geometric multigrid method that solves the steady planar Navier–Stokes equations in streamfunction-vorticity form on large rectangular domains in a quick and efficient manner. The results are then compared with a direct numerical simulation of the full unsteady problem using the Semtex software package where it is found that there is broad qualitative agreement, namely the separation and reattachment of the boundary layer at the equator. However, the presence of unobserved behaviour such as a large area of reverse flow seen at lower Reynolds numbers than those observed in other studies, and that the absolute error increases with Reynolds number suggest the model needs improvement to better capture the physical dynamics.

The importance of shear-induced lift, i.e. Saffman lift, is investigated in two-way coupling point-particle direct numerical simulation of turbulent channel flow laden with heavy particles. A new criterion is established for the critical particle diameter normalized by wall viscous length, and the effect of Saffman lift on particle motion can be neglected when the particle diameter is smaller than the critical value. The new criterion agrees well with the numerical results. We further study the influence of Saffman lift on particle motion and turbulence modulation. The results suggest that the Saffman lift will increase near-wall particle velocity and reduce particle accumulation therein. By integrating the particle motion equation, it is demonstrated that the mean particle streamwise velocity is dependent on the particle wall-normal velocity. This theoretically explains why the presence of a wall-normal force alters the streamwise velocity of particles. A theoretical profile of mean particle streamwise velocity at different Stokes numbers is then obtained by introducing the diffusive gradient hypothesis and Prandtl mixing length model. In addition, the present results reveal that the Saffman lift tends to augment the level of turbulence attenuation. The spectral analysis shows that the Saffman lift reduces the turbulent energy of small scales in the buffer region and large structures in the logarithmic region. When the Saffman lift is absent, near-wall particle streaks can be observed in the buffer region, which tend to separate fluid streaks. By contrast, the Saffman lift tends to prevent the formation of near-wall particle streaks.

Compressible direct numerical simulations are employed to elucidate the low-wavenumber behaviour of wall-pressure fluctuations in turbulent channel flow and the effect of flow Mach number in the nearly incompressible regime. Simulations are conducted at bulk Mach numbers 0.4, 0.2 and 0.1, and friction Reynolds number 180. In addition to the convective ridge that is virtually Mach-number-independent, acoustic ridges, whose magnitudes are orders of magnitude lower, are identified in the two-dimensional wavenumber–frequency spectrum. At lower frequencies, the acoustic ridges represent propagating longitudinal and oblique waves that match the theoretical predictions of two-dimensional duct modes with a uniform mean flow. They decay with decreasing Mach number but remain distinctly identifiable even at Mach 0.1. At high frequencies, in contrast, no propagating waves are found, and the spectral level in the supersonic wavenumber range is broadly elevated and increases with decreasing Mach number.

The nexus between turbulence, particle interaction and interfacial tension is virtually unexplored, despite being highly relevant to a wealth of industrial and environmental settings. Here we investigate it by conducting experiments on non-Brownian spherical particles at the interface of turbulent liquid layers. The latter are electromagnetically stirred in a quasi-two-dimensional apparatus, while the particles are individually tracked. By systematically varying interfacial conditions, turbulence intensity, particle size and concentration from dilute to dense, we map the system behaviour over a wide parameter space. We reveal how the dynamics is governed by the balance of drag, capillarity and lubrication. Based on their scaling, we propose a phase diagram comprising three distinct regimes, characterized by widely different levels of clustering and fluctuating energy of the particles. This is quantitatively confirmed by the experimental results.