In this paper, we study the linear stability of a two-dimensional shear flow of an air layer overriding a water layer of finite depth. The air layer is considered to be of an infinite extent with an exponential velocity profile. Three different background conditions are considered in the finite-depth water layer: a quiescent background, a linear velocity profile and a quadratic velocity profile. It is known that the cases of the quiescent water layer and the linear velocity profile allow for analytical treatment. We further provide an analytical solution for the case of the quadratic velocity field: we specifically consider a flow-reversal profile, although the result could be generalized to other quadratic profiles as well. The role of water layer depth on the growth rate of the Miles and rippling instabilities is studied in each of the three cases. Using asymptotic analysis, with the air–water density ratio being a small parameter, we obtain an analytical expression for the growth rate of the Miles mode and discuss the condition for the existence of a long-wave cutoff for these profiles. We provide analytical expressions for the stability boundary in the parameter space of inverse squared Froude number and wavenumber. In scenarios where a long-wave cutoff does not exist, we have carried out a long-wave asymptotic study to obtain the growth rate behaviour in that regime.

We consider steady, fully-developed flows of deformable, inelastic grains driven by gravity between identical bumpy walls. Using constitutive relations from extended kinetic theory (EKT) for the erodible bed near the centreline and the collisional flow between the surfaces of the bed and the walls, we calculate the fields of mean velocity, fluctuation velocity and solid volume fraction across the chute. We consider both situations in which the solid volume fraction at and near the centreline is high enough to form a bed and when it is not. We compare results predicted by EKT with recent discrete element simulations results, and obtain very good agreement.

The linear stability of the two-dimensional Taylor–Green vortices, which is a spatially periodic array of vortices, in rotating stratified fluids is investigated by local and modal stability analysis. Five types of instability appear in general: the pure hyperbolic instability, the strato-hyperbolic instability, the rotational-hyperbolic instability, the centrifugal instability and the elliptic instability. The condition for each instability and the estimate of the growth rate, which are useful in interpreting numerical results, are obtained in the framework of local stability analysis. Realizability of an instability is introduced to predict whether an unstable mode corresponding to an unstable region found in the local stability analysis exists at finite Reynolds numbers. In the absence of stratification, the pure hyperbolic instability is dominant for weak rotation; it is stabilized for strong rotation. For strong anti-cyclonic rotation, the elliptic instability or the centrifugal instability becomes dominant depending on the parameter values; further stronger rotation stabilizes both instabilities. For strong cyclonic rotation, the rotational-hyperbolic instability or the elliptic instability becomes dominant, although the growth rate is smaller than the anti-cyclonic cases. Strong stratification changes the stability properties. The strato-hyperbolic instability occurs for weak rotation. The rotational-hyperbolic instability and the elliptic instability are weakened under cyclonic rotation, while the latter survives and extends the unstable range under anti-cyclonic rotation. The pure hyperbolic instability and the centrifugal instability are less affected by stratification. The mode structures of each instability are in good agreement with the corresponding solution to local stability equations, confirming the physical mechanism of the instability.

In turbulence near a free surface, strong vortices attach to the surface, creating surface imprints visible as nearly circular ‘dimples’. By studying these imprints in direct numerical simulation (DNS) data we make two observations. First, the imprints of surface-attached vortices can be very effectively distinguished from other turbulent surface features using two physical features: they are nearly circular in shape, and persist for a long time compared with other pertinent time scales. Secondly, the instantaneous number of surface dimples from surface-attached vortices in an area, $N(t)$, is intimately related to its mean-square surface divergence, $\beta ^2(t)$. We develop a simple and physically transparent computer vision procedure which, using the properties of low eccentricity and longevity, detects and tracks vortices from their surface features only, with sensitivity and accuracy of $90\,\%$ or better. We compare $N(t)$ and $\beta ^2(t)$, finding a normalised cross-correlation of $0.90$, with changes in $N$ lagging around $0.8T_\infty$ behind those in $\beta ^2$ ($T_\infty$ is an integral time scale), confirming the common observation that vortices are spawned by strong upwelling events where $\beta ^2$ is large. These findings suggest that the rate of mass flux across the surface, being closely related to surface divergence, can be estimated remotely in some natural flows using visible free-surface dimples as proxy.

Drainage and spreading processes in thin liquid films have received considerable attention in the past decades. Yet, our understanding of three-dimensional cases remains sparse, with only a few studies focusing on flat and axisymmetric substrates. Here, we exploit differential geometry to understand the drainage and spreading of thin films on curved substrates, under the assumption of negligible surface tension and hydrostatic gravity effects. We develop a solution for the drainage on a local maximum of a generic substrate. We then investigate the role of geometry in defining the spatial thickness distribution via an asymptotic expansion in the vicinity of the maximum. Spheroids with a much larger (respectively smaller) height than the equatorial radius are characterized by an increasing (respectively decreasing) coating thickness when moving away from the pole. These thickness variations result from a competition between the variations of the substrate's slope and mean curvature. The coating of a torus presents larger thicknesses and a faster spreading on the inner region than on the outer region, owing to the different curvatures in these two regions. In the case of an ellipsoid with three different axes, spatial modulations in the drainage solution are observed as a consequence of a faster drainage along the short principal axis, faithfully reproduced by a three-dimensional asymptotic solution. Leveraging the conservation of mass, an analytical solution for the average spreading front is obtained. The solutions are in agreement with numerical simulations and experimental measurements obtained from the coating of a curing polymer on diverse substrates.

Floating rather stably above the convective mantle, continents exert a thermal blanket effect on the mantle by locally accumulating heat and altering the flow structure, which in turn affect continent motion. This thermal-mechanic feedback has been captured by a simplified model―a thermally insulating plate floating freely over a bottom-heated convecting fluid. Previous studies have revealed that, as plate size increases, plate motion transitions from long-term stagnancy over cold downwelling to short-term stagnancy over hot upwelling, termed stagnant modes (SMs) I and II, respectively, and eventually to a unidirectionally moving mode (UMM) with no stagnancy. During SM I/II, the adjacent hot/cold plumes are observed to move toward/away from the plate. The present study quantifies these plume motions through analysis on a partially insulated loop model established for SM I and SM II, respectively. The model predictions are well verified by results of numerical simulation. Moreover, the dominant frequency of plate speed is shown to increase linearly with plate size, with an abrupt jump as the coupling mode transitions from SM I to SM II. This suggests that the enhanced thermal blanket effect with plate size enables an increasingly fast feedback between the slowdown of a plate and the strengthening of hot upwelling which boosts plate motion. When plate motion is boosted before it stops, the UMM emerges. The predicted critical plate size for the UMM is well validated by simulation results. Geological implications of the results for continent motion, continent–mantle coupling and particularly the migration of plume and subduction sites are discussed.

This study presents the findings of a wind tunnel experiment investigating the behaviour of micrometric inertial particles with Stokes numbers around unity in the turbulent wake of a stationary porous disk. Various concentrations $\varPhi _{v}\in ([6-19] \times 10^{-6})$ of poly-disperse water droplets (average diameter 40–50 $\mathrm {\mu }$m) are compared with sub-inertial tracer particles. Hot-wire anemometry, phase Doppler interferometry and particle image velocimetry were implemented in the near- and far-wake regions to study the complex dynamics of such particles. Quadrant analysis is used to explore the shear effects of the particle wake interaction. Turbulence statistics and particle size distributions reveal distinct differences in the structure of the wake when inertial particles are present in the flow. Additionally, there are different structures in the near and far wake regions and structures change with particle volume fraction.

This study discusses turbulent suspension flows of non-Brownian, non-colloidal, neutrally buoyant and rigid spherical particles in a Newtonian fluid over porous media with particles too large to penetrate and move through the porous layer. We consider suspension flows with the solid volume fraction ${{\varPhi _b}}$ ranging from 0 to 0.2, and different wall permeabilities, while porosity is constant at 0.6. Direct numerical simulations with an immersed boundary method are employed to resolve the particles and flow phase, with the volume-averaged Navier–Stokes equations modelling the flow within the porous layer. The results show that in the presence of particles in the free-flow region, the mean velocity and the concentration profiles are altered with increasing porous layer permeability because of the variations in the slip velocity and wall-normal fluctuations at the suspension-porous interface. Furthermore, we show that variations in the stress condition at the interface significantly affect the particle near-wall dynamics and migration toward the channel core, thereby inducing large modulations of the overall flow drag. At the highest volume fraction investigated here, ${{\varPhi _b}}= 0.2$, the velocity fluctuations and the Reynolds shear stress are found to decrease, and the overall drag increases due to the increase in the particle-induced stresses.

Following our previous work on periodic ray paths (He et al., J. Fluid Mech., vol. 939, 2022, A3), we study asymptotically and numerically the structure of internal shear layers for very small Ekman numbers in a three-dimensional spherical shell and in a two-dimensional cylindrical annulus when the rays converge towards an attractor. We first show that the asymptotic solution obtained by propagating the self-similar solution generated at the critical latitude on the librating inner core describes the main features of the numerical solution. The internal shear layer structure and the scaling for its width and velocity amplitude in $E^{1/3}$ and $E^{1/12}$, respectively, are recovered. The amplitude of the asymptotic solution is shown to decrease to $E^{1/6}$ when it reaches the attractor, as is also observed numerically. However, some discrepancies are observed close to the particular attractors along which the phase of the wave beam remains constant. Another asymptotic solution close to those attractors is then constructed using the model of Ogilvie (J. Fluid Mech., vol. 543, 2005, pp. 19–44). The solution obtained for the velocity has an $O(E^{1/6})$ amplitude, but a self-similar structure different from the critical-latitude solution. It also depends on the Ekman pumping at the contact points of the attractor with the boundaries. We demonstrate that it reproduces correctly the numerical solution. Surprisingly, the numerical solution close to an attractor with phase shift (that is, an attractor touching the axis in three or two dimensions with a symmetric forcing) is also found to be $O(E^{1/6})$, but its amplitude is much weaker. However, its asymptotic structure remains a mystery.

In this paper, we conduct a systematic study of the instability of a boundary layer over a rotating cone that is inserting into a supersonic stream with zero angle of attack. The base flow is obtained by solving the compressible boundary-layer equations using a marching scheme, whose accuracy is confirmed by comparing with the full Navier–Stokes solution. Setting the oncoming Mach number and the semi-apex angle to be 3 and 7$^\circ$, respectively, the instability characteristics for different rotating rates ($\bar \varOmega$, defined as the ratio of the rotating speed of the cone to the axial velocity) and Reynolds numbers ($R$) are revealed. For a rather weak rotation, $\bar \varOmega \ll 1$, only the modified Mack mode (MMM) exists, which is an extension of the supersonic Mack mode in a quasi-two-dimensional boundary layer to a rotation configuration. Further increase of $\bar \varOmega$ leads to the appearance of a cross-flow mode (CFM), coexisting with the MMM but in the quasi-zero frequency band. The unstable zones of the MMM and CFM merge together, and so they are referred to as the type-I instability. When $\bar \varOmega$ is increased to an $O(1)$ level, an additional unstable zone emerges, which is referred to as the type-II instability to be distinguished from the aforementioned type-I instability. The type-II instability appears as a centrifugal mode (CM) when $R$ is less than a certain value, but appears as a new CFM for higher Reynolds numbers. The unstable zone of the type-II CM enlarges as $\bar \varOmega$ increases. The vortex structures of these types of instability modes are compared, and their large-$R$ behaviours are also discussed.

The impact of a drop train, a series of identical liquid drops separated by a constant distance, on a liquid pool initially generates a long slender cavity. However, the cavity soon collapses and turns into a shallow funnel. Here we theoretically model the dynamic profile of the elongated cavity and the steady shape of the funnel, which are then shown to agree well with experiment. When the liquid inertia plays a dominant role, the cavity assumes a parabolic profile that depends only on the drop diameter and the centre-to-centre spacing of adjacent drops. We consider the capillary forces as well as the drop impact force to obtain the shape of the funnel that persists once the elongated cavity collapses. Our study allows for predicting the interfacial morphology by the impact of a drop train and designing impact conditions useful for semiconductor cleaning processes.

A yield stress is added to Taylor's (1952, Proc. Royal Soc. A, 211, 225-239) model of a microscopic organism with a wavy cylindrical tail swimming through a viscous fluid. Viscoplastic slender-body theory is employed for the task, generalising existing results for Bingham fluid to the Herschel–Bulkley constitutive model. Numerical solutions are provided over a range of the two key parameters of the problem: the wave amplitude relative to the wavelength, and a Bingham number which describes the strength of the yield stress. Numerical solutions are supplemented with discussions of various limits of the problem in which analytical progress is possible. If the wave amplitude is sufficiently small, the yield stress of the material inevitably dominates the flow; the resulting ‘plastic locomotion’ results in swimming speeds that depend strongly on the swimming gait, and can, in some cases, even be negative. Conversely, when the yield stress is large, swimming becomes possible at the wave speed, with the swimmer sliding or burrowing along its centreline with a relatively high efficiency.

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths, vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn., vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids, vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. (Phys. Fluids, vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.