Droplet impact on oscillating substrates is important for both natural and industrial processes. Recognizing the importance of the dynamics that arises from the interplay between droplet transport and substrate motion, in this work, we present an experimental investigation of the spreading of a droplet impacting a sinusoidally oscillating hydrophobic substrate. We focus particularly on the maximum spread of droplets as a function of various parameters of substrate oscillation. We first quantify the maximum spreading diameter attained by the droplets as a function of frequency, amplitude of vibration, and phase at the impact for various impact velocities. We highlight that there can be two stages of spreading. Stage I, which is observed at all impact conditions, is controlled by the droplet inertia and affected by the substrate oscillation. For certain conditions, a Stage II spreading is also observed, which occurs during the retraction process of Stage I due to additional energies imparted by the substrate oscillation. Subsequently, we derive scaling analyses to predict the maximum spreading diameters and the time for this maximum spread for both Stage I and Stage II. Furthermore, we identify the necessary condition for Stage II spreading to be greater than Stage I spreading. The results will enable optimization of the parameters in applications where substrate oscillation is used to control the droplet spread, and thus heat and mass transfer between the droplet and the substrate.

Particulate flows at moderate particle Reynolds numbers are important in critical engineering and geological applications. This experimental study explores neutrally buoyant suspensions in an outer-rotating coaxial rheometer for solid fractions, $\phi$, from 0.1 to 0.5, and particle Reynolds number, $Re$, from 0.5 to 800, covering laminar, transitional and turbulent regimes; $Re$ is defined in terms of the square of the particle diameter and the shear rate. For $0.1 < \phi < 0.4$ and $0.5 < Re <10$, the direct torque measurements normalised by the laminar flow torque, $M/M_{lam}$, are independent of $Re$, but depend on $\phi$. For the same range of $\phi$ and for $10< Re<100$, the normalised torques depend on both $\phi$ and $Re$, and show an increasing dependence on $Re$. As $Re$ increases, the flow transitions to turbulence. Small particles delay the turbulent transition for $\phi \leqslant 0.3$, while large particles augment the transition. A modified Reynolds number, $Re^\prime$, that depends linearly on the particle diameter and the maximum velocity, $U_{o}$, is introduced for both laminar and turbulent flows and shows a better correlation of the results as compared with $Re$. For $\phi = 50\,\%$, the normalised torque minus the torque at zero rotational speed is nearly independent of $Re^\prime$. Rheological models based on $Re^\prime$ and the Krieger–Dougherty relative viscosity are proposed in the laminar regime for $10< Re^\prime <500$; in the turbulent regime, a correlation is proposed in terms of $Re^\prime$ and $\phi$ for $1000< Re^\prime < 6000$.

We consider two-dimensional flows above topography, revisiting the selective decay (or minimum enstrophy) hypothesis of Bretherton and Haidvogel. We derive a ‘condensed branch’ of solutions to the variational problem where a domain-scale condensate coexists with a flow at the (smaller) scale of the topography. The condensate arises through a supercritical bifurcation as the conserved energy of the initial condition exceeds a threshold value, a prediction that we quantitatively validate using direct numerical simulations. We then consider the forced–dissipative case, showing how weak forcing and dissipation select a single dissipative state out of the continuum of solutions to the energy-conserving system predicted by selective decay. As the forcing strength increases, the condensate arises through a supercritical bifurcation for topographic-scale forcing and through a subcritical bifurcation for domain-scale forcing, both predictions being quantitatively validated by direct numerical simulations. This method provides a way of determining the equilibrated state of forced–dissipative flows based on variational approaches to the associated energy-conserving system, such as the statistical mechanics of two-dimensional flows or selective decay.

In this work, formal asymptotic solutions of a problem for linear water waves in a bounded basin are constructed. The solutions have the form of asymptotic quasimodes and are used for the description of standing water waves localised near the shoreline. Such short-wavelength quasimodes exist only for a discrete set of frequencies, which are determined by means of a quantisation-type condition. Some numerical results are also addressed.

We use $123$ three-dimensional direct numerical simulations to study fingering convection in non-rotating spherical shells. We investigate the scaling behaviour of the flow length scale, the non-dimensional heat and compositional fluxes $Nu$ and $Sh$ and the mean convective velocity over the fingering convection instability domain defined by $1 \leq R_\rho < Le$, $R_\rho$ being the ratio of density perturbations of thermal and compositional origins and $Le$ the Lewis number. We show that the chemical boundary layers are marginally unstable and adhere to the laminar Prandtl–Blasius model, hence explaining the asymmetry between the inner and outer spherical shell boundary layers. We develop scaling laws for two asymptotic regimes close to the two edges of the instability domain, namely $R_\rho \lesssim Le$ and $R_\rho \gtrsim 1$. For the former, we develop novel power laws of a small parameter $\epsilon$ measuring the distance to onset, which differ from theoretical laws published to date in Cartesian geometry. For the latter, we find that the Sherwood number $Sh$ gradually approaches a scaling $Sh\sim Ra_\xi ^{1/3}$ when $Ra_\xi \gg 1$; and that the Péclet number accordingly follows $Pe \sim Ra_\xi ^{2/3} |Ra_T|^{-1/4}$, $Ra_T$ and $Ra_{\xi}$ being the thermal and chemical Rayleigh numbers. When the Reynolds number exceeds a few tens, we report on a secondary instability which takes the form of large-scale toroidal jets which span the entire spherical domain. Jets distort the fingers, resulting in Reynolds stress correlations, which in turn feed the jet growth until saturation. This nonlinear phenomenon can yield relaxation oscillation cycles.

Analysing the dynamics of velocity gradients is useful for understanding various nonlinear turbulence processes. This work focuses on how the vibrational non-equilibrium of the constituent molecules in a gaseous medium affects the dynamics of the velocity gradient and the pressure Hessian tensors. We first derive the exact evolution equation of the pressure-Hessian tensor in the presence of the vibrational non-equilibrium process. Subsequently, we perform several direct numerical simulations of compressible isotropic turbulence, including the vibrational relaxation process therein. Using flow fields extracted from these simulations, we conduct several parametric studies over a range of the Damköhler number (ratio of the relevant fluid time scale to that of the mean vibrational relaxation process) and the initial ratio of the vibrational temperature to the mean local temperature. We find that a variation in the initial Damköhler number does influence the evolution of the pressure-Hessian and the velocity gradient tensors. As the vibrational relaxation process becomes more rapid (an increase in the value of the initial Damköhler number), it causes a decrease in the strength of the pressure-Hessian tensor and simultaneous suppression of dilatational fluctuations in the flow field. On the other hand, a variation in the initial value of the ratio of the vibrational temperature to the local temperature does not seem to affect the pressure-Hessian or the velocity gradient tensor. These findings are expected to aid in the development of closure models for the pressure-Hessian tensor in compressible flows under vibrational non-equilibrium conditions.

An ultrasonic phased array system is introduced to study the three-dimensional (3-D) movement of a single bubble in a GaInSn alloy under a transverse magnetic field (MF), which is verified by bubble experiments in water. The 3-D motion trajectories of individual bubbles in the GaInSn are obtained under a horizontal MF. As the MF becomes stronger, the bubble successively oscillates in random directions (R mode), a direction perpendicular to the MF (V mode), a direction parallel to the MF (P mode) and finally it rises straight (S mode). The significant anisotropy of the oscillation directions at a moderate MF intensity may be due to the anisotropy of the vortex structure around the bubble. Furthermore, the oscillation amplitude gradually declines with increasing MF intensity until the bubble trajectory finally becomes a straight line. Our measurements allow us to specify the characteristic regions for the observed bubble modes in the $N-Eo-Re$ parameter space (N is the magnetic interaction parameter, Eo is the Eötvös number and Re is the Reynolds number). In addition, more detailed characteristics of bubble terminal velocity are revealed, showing that the bubble velocities are closely related to the motion modes. The increase in bubble velocity at a moderate MF intensity is caused by the weakening oscillation. At a high strength, the MF monotonically suppresses the rise velocity of the bubble with a fixed scaling law.

We explored the settling dynamics of vertically aligned particles in a quiescent, stratified two-layer fluid using particle tracking velocimetry. Glass spheres of $d=4\,{\rm mm}$ diameter were released at frequencies of 4, 6 and 8 Hz near the free surface, traversing through an upper ethanol layer ($H_1$), where H is height or layer thickess, varying from $10d$ to $40d$ and a lower oil layer. Results reveal pronounced lateral particle motion in the ethanol layer, attributed to a higher Galileo number ($Ga = 976$, ratio of buoyancy–gravity to viscous effects), compared with the less active behaviour in the oil layer ($Ga = 16$). The ensemble vertical velocity of particles exhibited a minimum just past the density interface, becoming more pronounced with increasing $H_1$, and suggesting that enhanced entrainment from ethanol to oil resulted in an additional buoyancy force. This produced distinct patterns of particle acceleration near the density interface, which were marked by significant deceleration, indicating substantial resistance to particle motion. An increased drag coefficient occurred for $H_1/d = 40$ compared with a single particle settling in oil; drag reduced as the particle-release frequency ($\,f_p$) increased, likely due to enhanced particle interactions at closer proximity. Particle pair dispersions, lateral ($R^2_L$) and vertical ($R^2_z$), were modulated by $H_1$, initial separation $r_0$ and $f_p$. The $R^2_L$ dispersion displayed ballistic scaling initially, Taylor scaling for $r_0 < H_1$ and Richardson scaling for $r_0 > H_1$. In contrast, $R^2_z$ followed a $R^2_z \sim t^{5.5}$ scaling under $r_0 < H_1$. Both $R^2_L$ and $R^2_z$ plateaued at a distance from the interface, depending on $H_1$ and $f_p$.

Cowl-induced incident shock wave/boundary layer interactions (ISWBLIs) under the influence of shoulder expansion represent one of the dominant phenomena in supersonic inlets. To provide a more comprehensive understanding of how an expansion corner affects the ISWBLI, a detailed experimental and analytical study is performed in a Mach 2.73 flow in this work. Pressure measurement, schlieren photography and surface oil-flow visualisation are used to record flow features, including the pressure distribution, separation extent and surface-flow topological structures. Our results reveal three types of ISWBLIs influenced by the expansion corner. When the shock intensity is weak, the separation is small scale with the expansion waves emanating from the expansion corner. This is the first type of expansion-corner-affected ISWBLI (EC-ISWBLI). When the incident shock wave is strong, large-scale separation occurs, accompanied by the disappearance of expansion waves, forming the second type of EC-SWBLI. The expansion corner induces a ‘lock-in’ effect in which the separation onset is consistently locked near the expansion corner regardless of the incident shock intensity and impingement position. The third type of EC-ISWBLI occurs when the shock is sufficiently strong and the impingement point is close to the expansion corner. In this interaction, the ‘lock-in’ effect ceases to manifest. Moreover, a shock polar-incorporating inviscid model is employed to elucidate the shock patterns. Two criteria are established by combining free interaction theory with this model. The first criterion provides valuable insights into the evolution of separations with a minimal overall pressure rise and the second criterion determines the threshold for the occurrence of the ‘lock-in’ effect.

A simple continuum model is proposed for discrete granular avalanches, as observed in slowly rotating drums. Subject to granular transfer across the basal interface, the model evolves the transient surface inclination and mid-slope velocity profile of the avalanches, from failure to arrest. This is done using new boundary conditions that allow entrainment or detrainment contingent on the basal shear rate. For entrainment to occur, the basal shear stress must overcome the erosion resistance of the jammed deposit, while for detrainment to take place the basal shear rate must vanish. The resulting avalanche dynamics is controlled by a single dimensionless number, the ratio of excess slope at failure to excess resistance to erosion, that can be determined from inclination history data. This number provides a measure of slope brittleness, from which the peak flow rate and arrest inclination can be predicted. Analytical techniques are used to derive model solutions, and clarify the resulting behaviour. Finally, the model is tested by comparing solutions with laboratory experiments and discrete particle simulations.

In the present study, we performed direct numerical simulations for a hypersonic turbulent boundary layer over the windward side of a lifting body, the HyTRV model, at Mach number $6$ and attack angle 2$^{\circ }$ to investigate the global and local turbulent features, and evaluate its difference from canonical turbulent boundary layers. By scrutinizing the instantaneous and averaged flow fields, we found that the transverse curvature on the windward side of the HyTRV model induces the transverse opposing pressure gradients that push the flow on both sides towards the windward symmetry plane, yielding significant effects of the azimuthal inhomogeneity and large-scale cross-stream circulations, moderate and azimuthal independent influences of adverse pressure gradient, and negligible impact of the mean flow three-dimensionality. Further inspecting the local turbulent statistics, we identified that the mean and fluctuating velocity become increasingly similar to the highly decelerated turbulent boundary layers over flat plates in that the mean velocity deficit is enhanced, and the outer layer Reynolds stresses are amplified as it approaches the windward symmetry plane, and prove to be self-similar under the scaling of Wei & Knopp (J. Fluid Mech., vol. 958, 2023, A9) for adverse-pressure-gradient turbulent boundary layers. Conditionally averaged Reynolds stresses based on strong sweeping and ejection events demonstrated that the Kelvin–Helmholtz instability of the strong embedded shear layer induced by the large-scale cross-stream circulations is responsible for the turbulence amplification in the outer layer. The strong Reynolds analogy that relates the mean velocity and temperature was refined to incorporate the non-canonical effects, showing considerable improvements in the accuracy of such a formula. On the other hand, the temperature fluctuations are still transported passively, as indicated by their resemblance to the velocity. The conclusions obtained in the present study provide potentially profitable information for turbulent modelling modification for the accurate predictions of skin friction and wall heat transfer.

In this paper, direct numerical simulations in hypersonic turbulent boundary layers over a $24^{\circ }$ compression ramp at Mach 6.0 are performed. The wall skin friction and its spanwise non-homogeneity in the interaction region are analysed via the spectral analysis and drag decomposition method. On the compression ramp, the premultiplied spanwise energy spectrum of wall shear stress $\tau _{w}$ reveals two energetic spanwise length scales. One occurs in the region of $x/\delta _{ref}=0\unicode{x2013}3$ ($x=0$ lies in the compression corner; $\delta _{ref}$ is the boundary layer thickness upstream of the interaction region) and is consistent with that of the large-scale streamwise vortices, indicating that the fluctuation intensity of $\tau _{w}$ is associated with the Görtler-type structures. The other one is observed downstream of $x/\delta _{ref}=3.0$ and corresponds to the regenerated elongated streaky structures. The fluctuation intensity of $\tau _{w}$ peaks at $x/\delta _{ref}=3.0$, where both the above energetic length scales are observed. The drag decomposition method proposed by Li et al. (J. Fluid Mech., vol. 875, 2019, pp. 101–123) is extended to include the effects of spanwise non-homogeneity so that it can be used in the interaction region where the mean flow field and the mean skin friction $C_f$ exhibit an obvious spanwise heterogeneity. The results reveal that, in the upstream turbulent boundary layer, the drag contribution arising from the spanwise heterogeneity can be neglected, while this value on the compression ramp is up to 20.7 % of $C_f$, resulting from the Görtler-type vortices. With the aid of the drag decomposition method, it is found that the main flow features that contribute positively to the amplification of $C_f$ and its rapid increase on the compression ramp includes: the density increase across the shock, the high mean shear stress and turbulence amplification around the detached shear layer and the Favre-averaged downward velocity towards the ramp wall. Compared with the spanwise-averaged value, $C_f$ and its components at the spanwise station where the downwash and upwash of the Görtler-type vortices occur reveal a spanwise variation exceeding 10 %.

When a viscous fluid partially fills a Hele-Shaw channel, and is pushed by a pressure difference, the fluid interface is unstable due to the Saffman–Taylor instability. We consider the evolution of a fluid region of finite extent, bounded between two interfaces, in the limit that the interfaces are close, that is, when the fluid region is a thin liquid filament separating two gases of different pressure. In this limit, we derive a second-order ‘thin-filament’ model that describes the normal velocity of the filament centreline, and evolution of the filament thickness, as functions of the thickness, centreline curvature and their derivatives. We show that the second-order terms in this model, that include the effect of transverse flow along the filament, are necessary to regularise the instability. Numerical simulation of the thin-filament model is shown to be in accordance with level-set computations of the complete two-interface model. Solutions ultimately evolve to form a bubble of rapidly increasing radius and decreasing thickness.

We investigate the effects of bacterial activity on the mixing and transport properties of a passive scalar in time-periodic flows in experiments and in a simple model. We focus on the interactions between swimming Escherichia coli and the Lagrangian coherent structures (LCSs) of the flow, which are computed from experimentally measured velocity fields. Experiments show that such interactions are non-trivial and can lead to transport barriers through which the scalar flux is significantly reduced. Using the Poincaré map, we show that these transport barriers coincide with the outermost members of elliptic LCSs known as Lagrangian vortex boundaries. Numerical simulations further show that elliptic LCSs can repel elongated swimmers and lead to swimmer depletion within Lagrangian coherent vortices. A simple mechanism shows that such depletion is due to the preferential alignment of elongated swimmers with the tangents of elliptic LCSs. Our results provide insights into understanding the transport of micro-organisms in complex flows with dynamical topological features from a Lagrangian viewpoint.

The turbulent boundary layer (TBL) is a widely existing flow phenomenon in nature and engineering applications. Its strong mixing effect can achieve more sufficient material mixing, heat transport, etc. The understanding of the entrainment process and mechanism of irrotational fluids entering the turbulent region can be promoted by studying the geometric and dynamic characteristics of turbulent${/}$non-turbulent interfaces (TNTI). In compressible flow, it is unclear whether the properties of TNTI will change and whether the entrainment will show different features due to the influence of compressibility. Based on the direct numerical simulation results of supersonic compressible plate TBLs with Mach number of 2.9, the geometric and dynamic characteristics of TNTI are investigated in this paper. The interface is identified by the enstrophy method, and the height, thickness, fractal dimension, enstrophy transportation and entrainment characteristics of the interface are investigated. It is found that for the enstrophy transportation in a TBL, the contribution of compressibility-related terms accounts for approximately 13.4 % of the total enstrophy transportation, which tends to transfer the enstrophy of turbulence near the interface to both directions vertical to the interface. This promotes the expansion of the turbulent region towards the non-turbulent region, and the mean height, thickness and entrainment velocity are increased by approximately 3.7 %, 7.0 % and 8.5 %, respectively, while the fractal dimension is basically unaffected. Different from the incompressible flow, the contribution of the compressibility-related terms to the entrainment velocity is independent of the local curvature, and the intense entrainment process is more likely to occur on a highly curved concave surface.

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.

We investigate two-dimensional (2-D) axisymmetric flow in toroidal geometry, with a focus on a transition between 2-D three-component flow and 2-D two-component flow. This latter flow state allows a self-organization of the system to a quiescent dynamics, characterized by long-living coherent structures. When these large-scale structures orient in the azimuthal direction, the radial transport is reduced. Such a transition, if it can be triggered in toroidally confined fusion plasmas, is beneficial for the generation of zonal flows and should consequently result in a flow field beneficial for confinement.

The transition route and bifurcations of the buoyant flows developing on a heated horizontal circular surface are elaborated using direct numerical simulations and direct stability analysis. A series of bifurcations, as a function of Rayleigh numbers ($Ra$) ranging from $10^6$ to $6.0\times 10^7$, are found on the route to chaos of the flows at $Pr=7$. When $Ra<1.0\times 10^3$, the buoyant flows above the heated horizontal surface are dominated by conduction, because of which the distinct thermal boundary layer and plume are not present. At $Ra=1.1\times 10^6$, a Hopf bifurcation occurs, resulting in the flow transition from a steady state to a periodic puffing state. As $Ra$ increases further, the flow enters a periodic rotating state at $Ra=1.9\times 10^6$, which is a unique state that was rarely discussed in the literature. These critical transitions, leaving from a steady state and subsequently entering a series of periodic states (puffing, rotating, flapping and period-doubling) and finally leading to chaos, are diagnosed using two-dimensional Fourier transforms. Moreover, direct stability analysis is conducted by introducing random numerical perturbations into the boundary condition of the surface heating. We find that when the state of a flow is in the vicinity of critical values (e.g. $Ra=2.0\times 10^6$), the flow is conditionally unstable to perturbations, and it can bifurcate from the rotating state to the flapping state in advance. However, for relatively stable flow states, such as at $Ra=1.5\times 10^6$, the flow remains in its periodic puffing state even though it is being perturbed.

We present experiments on oscillating hydrofoils undergoing combined heaving and pitching motions, paying particular attention to connections between propulsive efficiency and coherent wake features extracted using modal analysis. Time-averaged forces and particle image velocimetry measurements of the flow field downstream of the foil are presented for a Reynolds number of $Re=11\times 10^3$ and Strouhal numbers in the range $St=0.16\unicode{x2013}0.35$. These conditions produce 2S and 2P wake patterns, as well as a near-momentumless wake structure. A triple decomposition using the optimized dynamic mode decomposition method is employed to identify dominant modal components (or coherent structures) in the wake. These structures can be connected to wake instabilities predicted using spatial stability analyses. Examining the modal components of the wake provides insightful explanations into the transition from drag to thrust production, and conditions that lead to peak propulsive efficiency. In particular, we find modes that correspond to the primary vortex development in the wakes. Other modal components capture elements of bluff body shedding at Strouhal numbers below the optimum for peak propulsive efficiency and characteristics of separation for Strouhal numbers higher than the optimum.

The bubbly shock-driven partial cavitation in an axisymmetric venturi is studied with time-resolved two-dimensional X-ray densitometry. The bubbly shock waves are characterised using the vapour fraction and pressure changes across it, propagation velocity, and Mach number. The sharp changes in vapour fraction measured with X-ray densitometry, combined with high-frequency dynamic pressure measurements, reveal that the interaction of the pressure wave with the vapour cavity dictates the shedding dynamics. At the lowest cavitation number ($\sigma \sim 0.47$), the condensation shock front is the predominant shedding mechanism. However, as $\sigma$ increases ($\sigma \sim 0.78$), we observe an upstream travelling pressure discontinuity that changes into a condensation shock as it approaches the venturi throat. This coincides with the increasing strength of the bubbly shock wave as it propagates upstream, manifested by the increasing velocity of the shock front and the pressure rise across it. Consequently, the Mach number of the shock front increases and surpasses the critical value 1, favouring condensation shocks. Further, at higher $\sigma$ (${\sim }0.84\unicode{x2013}0.9$), both the re-entrant jet and pressure wave can cause cavity detachment. However, at such $\sigma$, the pressure wave likely remains subsonic. Hence cavity condensation is not favoured readily. This leads to the re-entrant jet causing the cavity detachment at higher $\sigma$. The shock front is accelerated as it propagates upstream through the variable cross-section of the venturi. This enhances its strength, favouring cavity condensation and eventual shedding. These observations explain the existence of shock fronts in an axisymmetric venturi for a large range of $\sigma$.