Predicting the magnitude of the slip velocity of non-tracer particles with respect to the surrounding fluid is crucial to address both fundamental and practical questions involving dispersed turbulent flows. Here we derive an analytical model to predict the slip velocity of spherical particles in homogeneous isotropic turbulence. We modulate the particle equation of motion according to the inertial filtering framework, and obtain closed-form expressions for the mean slip velocity magnitude as a function of the governing parameters. These are compared against laboratory measurements and direct numerical simulations, demonstrating close agreement for both light and heavy particles, both smaller and larger than the Kolmogorov scales. The predictive value of the model and its implications are discussed, as well as the range of validity of the underlying assumptions.

Direct numerical simulation of the three-dimensional (3-D) wake transition of a heated square cylinder subjected to horizontal cross-flow is performed in the presence of buoyancy. In order to capture the effects of large-scale heating, a non-Oberbeck–Boussinesq model is utilized, which includes the governing equations for compressible gas flow. All computations are performed at low free stream Mach number $M=0.1$ using air (free stream Prandtl number, $Pr=0.71$) as the working fluid. The 3-D instability modes A and B, which correspond to free stream Reynolds numbers of 180 and 250, are observed with longer and shorter spanwise wavelengths, respectively, and the onset of three-dimensionality is triggered at a Reynolds number of 173. In the presence of buoyancy, baroclinic vorticity production in the near-wake plays an important role for streamwise vorticity generation. The chaotic wake of the Mode-A instability bifurcates into periodic and quasiperiodic wakes at various heating levels, expressed by the overheat ratio, $\varepsilon =(T_w-T_\infty )/T_\infty$, where $T_w$ and $T_\infty$ are the temperature of the cylinder surface and the ambient air, respectively. At low heating ($\varepsilon =0.2$), the 3-D Mode-A instability is suppressed leading to a two-dimensional wake flow. Further increase in heating, again brings back the three-dimensionality in the wake through Mode-E instability. The variation of thermophysical properties and the effective Reynolds number with increase in heating level around the cylinder is examined. It is shown that the effect of thermophysical properties competes with the baroclinic streamwise vorticity generation at higher levels of heating ($\varepsilon \geqslant 0.4$) to control the 3-D modes and wake dynamics.

Buoyancy-driven bubbly flows play pivotal roles in various scenarios, such as the oxygenation and mixing in the upper ocean and the reaction kinetics in chemical and bio-reactors. This work focuses on the convective flow induced by the localised release of large air bubbles ($D_b=3.7$ mm, ${Re}_b=950$) in a water tank, exploring the resulting flow and the transition from laminar to disturbed states as a function of the Rayleigh number ranging from $3\times 10^3$ to $2\times 10^5$. At low ${Ra}$ the flow is smooth and laminar with weak temporal oscillations, while a highly disturbed state appears above a critical value ${Ra}_c$. A theoretical analysis is presented that links the mean flow circulation to the Rayleigh number. Through an experimental investigation, utilising three-dimensional particle tracking velocimetry and flow visualisation, we confirm the theory presented, and characterise the laminar to disturbed transition in the system. These findings not only enhance our fundamental understanding of buoyancy-driven convective flows but also hold significant implications for practical applications, particularly in the optimisation of bio-reactor design and other industrial processes reliant on controlled convective dynamics.

We use the Korteweg–de Vries (KdV) equation, supplemented with several forcing/friction terms, to describe the evolution of wind-driven water wave packets in shallow water. The forcing/friction terms describe wind-wave growth due to critical level instability in the air, wave decay due to laminar friction in the water at the air–water interface, wave stress in the air near the interface induced by a turbulent wind and wave decay due to a turbulent bottom boundary layer. The outcome is a modified KdV–Burgers equation that can be a stable or unstable model depending on the forcing/friction parameters. To analyse the evolution of water wave packets, we adapt the Whitham modulation theory for a slowly varying periodic wave train with an emphasis on the solitary wave train limit. The main outcome is the predicted growth and decay rates due to the forcing/friction terms. Numerical simulations using a Fourier spectral method are performed to validate the theory for various cases of initial wave amplitudes and growth and/or decay parameter ranges. The results from the modulation theory agree well with these simulations. In most cases we examined, many solitary waves are generated, suggesting the formation of a soliton gas.

The effectiveness of utilizing heating patterns as a drag-reduction tool in sloping channels is analysed. The usefulness of heating is judged by determining the pressure gradient required to maintain the same flow rate as in the isothermal case. The key to reducing pressure loss is the formation of separation bubbles, although these bubbles are washed away at relatively large Reynolds numbers. The bubbles reduce the direct contact between the stream and the side walls, thereby reducing the friction experienced by the flow. Moreover, the fluid inside the bubbles tends to rotate, a motion provoked by longitudinal temperature gradients. This rotation also seems to reduce the resistance. On the other hand, the existence of the bubbles tends to obstruct the stream, increasing the flow resistance. In general, channels oriented close to horizontal experience a relatively small pressure loss, but this loss grows markedly as the channel inclines towards the vertical. When modest heating is applied, the pressure loss is approximately proportional to the square of the associated Rayleigh number. It is also shown that if the heating wavelength is too short or too long, the heating loses its effectiveness. In certain circumstances, it turns out that the theoretical pressure-gradient reduction achieved by judicious heating is so large that it exceeds the pressure gradient required to drive the flow in the isothermal problem. The conclusion is that in these instances, a pressure gradient of the opposite sign must be applied to prevent flow acceleration.

The present experimental study shows that a nozzle with optimal flexibility can enhance the impulse and entrainment of a pulsed jet. Near the nozzle exit, vortex rings emanating from the flexible nozzle move faster because of the timely release of the elastic energy (stored during the expansion) to the jet, which is maximized at the structural stiffness that needs to be optimally tuned to the jet acceleration. The total circulation, hydrodynamic impulse and entrained fluid volume are enhanced substantially. Interestingly, we find that the same condition for optimal flexibility to maximize the hydrodynamic impulse and circulation of the primary vortex ring of the continuous jet (Choi & Park, J. Fluid Mech., vol. 949, 2022, A39) holds universally for the pulsed jet, indicating that abrupt jet termination is irrelevant to the impulse augmentation mechanism. Compared to the rigid counterpart, increments of the impulse (${\sim }400\,\%$) and entrainment (${\sim }220\,\%$) of a pulsed jet in the present study are considerably larger than those ($200\,\%$ and $50\,\%$, respectively) in a continuous jet from previous studies, which is attributed to the significant suppression of negative pressure at the nozzle exit by the collapsing motion of the flexible nozzle in the phase with the jet-driven upstream propagation of the surface wave on the nozzle. This universal mechanism provides a guideline for a novel jet propulsor using a flexible nozzle, for example, for small-scale underwater robots.

High-frequency observation data, including all three components of instantaneous fluctuating velocity, temperature, as well as particulate matter 10 ($PM_{10}$), collected from the unstable atmospheric surface layer at $z/L = -0.11$ and $-$0.12, $L$ being the Obukhov length, during sand and dust storms (SDS), were used to explore the scaling of vertical coherence and the logarithmic energy profile for wall-attached eddies. The present results demonstrate good agreement with the self-similar range of the wall-attached features for velocity and temperature components, as well as for $PM_{10}$ at lower heights ($z<15$ m) during SDS. Following the idea depicted by Davenport (Q. J. R. Meteorol., vol. 372, 1961, pp. 194–211), an empirically derived transfer kernel comprises implicit filtering via a scale-dependent gain and phase, parametrically defined as $|H_L^2(f)|=\exp (c_1-c_2\delta /\lambda _x)$, where $c_1$ and $c_2$ are parameters, $\delta$ is the boundary layer thickness and $\lambda _x$ is the streamwise wavelength. Linear coherence spectrum analysis is applied as a filter to separate the coherent and incoherent portions. After this separation procedure, the turbulence intensity decay for wall-attached eddies is described in a log–linear manner, which also identifies how the scaling parameter differs between the measured components. These findings present abundant features of wall-attached eddies during SDS which further are used to improve/enrich existing near-wall models.

Two-dimensional turbulence transfers its energy towards the lowest mode in the domain, but domain geometry exerts a powerful control. On the sphere, with its three axes of rotational symmetry, angular momentum conservation prevents energy from entering the three lowest modes – those corresponding to the spherical harmonics $Y_1^0$ and $Y_1^{\pm 1}$ – because the amplitudes of these three modes are proportional to the three conserved components of the angular momentum vector. Non-spherical ellipsoids partly or completely break the rotational symmetry corresponding to angular momentum conservation. The flow on spheroids, which have only one axis of rotational symmetry, conserves only a single component of angular momentum. If the axis of symmetry is taken to be the $z$-axis, then only the $z$-component of angular momentum is conserved. Energy can flow into the other two lowest modes. The general triaxial ellipsoid breaks all rotational symmetries, thus angular momentum is not conserved, and energy can flow into any mode. We describe numerical experiments that confirm these predictions.

In this work we investigate the spatio-temporal nature of various coherent modes present in a wind turbine wake using a combination of new particle image velocimetry experiments and data from Biswas & Buxton (J. Fluid Mech., vol. 979, 2024, A34). A multiscale triple decomposition of the acquired velocity field is sought to extract the coherent modes and, thereafter, the energy exchanges to and from them are studied using the multiscale triple decomposed coherent kinetic energy budgets developed by Baj & Buxton (Phys. Rev. Fluids, vol. 2, 2017, 114607). Different frequencies forming the tip vortex system (such as the blade passing frequency, turbine's rotational frequency and their harmonics) are found to be energised by different sources such as production from the mean flow or nonlinear triadic interaction or both, similar to the primary, secondary or the mixed modes discussed in Biswas et al. (J. Fluid Mech., vol. 941, 2022, A36). The tip vortex system forms a complex network of nonlinear triadic energy transfers, the nature and the magnitudes of which depend on the tip speed ratio ($\lambda$). Contrastingly, the modes associated with the sheddings from the nacelle or tower and wake meandering are found to be primarily energised by the mean flow. We show that the tip vortex system exchanges energy with the mean flow primarily through the turbine's rotational frequency. In fact, the system transfers energy back to the mean flow through the turbine's rotational frequency at some distance downstream marking the onset location of wake recovery ($x_{wr}$). Here $x_{wr}$ is shown to reduce with $\lambda$ due to stronger interaction and earlier merging of the tip vortices at a higher $\lambda$.

Drop collision with a solid particle is a ubiquitous phenomenon in a wide range of applications, including rain, spray coating, cooling or cleaning, particle encapsulation, inkjet printing, and additive manufacturing. Understanding the dynamics of drop collision is essential for optimizing these processes. In this study, we present a comprehensive experimental and analytical investigation of non-axisymmetric as well as axisymmetric drop impact on a solid particle. We use a high-speed video system to visualize the drop profile during the impact, and measure the drop height and spreading diameter for different liquid viscosities, ratios of the target to drop diameters, offsets, and various other impact parameters. We then develop a theoretical model for drop spreading on a solid spherical particle that relies on the formulation of a remote asymptotic solution for the inviscid flows, generated by non-axisymmetric drop impact. Next, the viscous effects in a thin viscous boundary layer are considered, which allows the formulation of an expression for the residual lamella thickness and maximum spreading. The theoretically predicted evolution of the lamella thickness, the residual film thickness, and the maximum spreading angle agree well with the experimental data presented in this work and the literature. Finally, we present a novel approach for in situ measurement of liquid viscosity, drop impact viscometry, at high shear rates via a single drop impact experiment, with potential application in industries where non-Newtonian drops play a major role, such as pesticide spraying, paint droplet spreading, blood drop impact and fuel injectors.

Three-dimensional non-rotating odd viscous liquids give rise to Taylor columns and support axisymmetric inertial-like waves (J. Fluid Mech., vol. 973, 2023, A30). When an odd viscous liquid is subjected to rigid-body rotation however, there arise in addition a plethora of other phenomena that need to be clarified. In this paper, we show that three-dimensional incompressible or two-dimensional compressible odd viscous liquids, rotating rigidly with angular velocity $\varOmega$, give rise to both oscillatory and evanescent inertial-like waves or a combination thereof (which we call of mixed type) that can be non-axisymmetric. By evanescent, we mean that along the radial direction, typically when moving away from a solid boundary, the velocity field decreases exponentially. These waves precess in a prograde or retrograde manner with respect to the rotating frame. The oscillatory and evanescent waves resemble respectively the body and wall-modes observed in (non-odd) rotating Rayleigh–Bénard convection (J. Fluid Mech., vol. 248, 1993, pp. 583–604). We show that the three types of waves (wall, body or mixed) can be classified with respect to pairs of planar wavenumbers $\kappa$ which are complex, real or a combination, respectively. Experimentally, by observing the precession rate of the patterns, it would be possible to determine the largely unknown values of the odd viscosity coefficients. This formulation recovers as special cases recent studies of equatorial or topological waves in two-dimensional odd viscous liquids which provided examples of the bulk–interface correspondence at frequencies $\omega <2\varOmega$. We finally point out that the two- and three-dimensional problems are formally equivalent. Their difference then lies in the way data propagate along characteristic rays in three dimensions, which we demonstrate by classifying the resulting Poincaré–Cartan equations.

Dynamic stall at low Reynolds numbers, $Re \sim O(10^4)$, exhibits complex flow physics with co-existing laminar, transitional and turbulent flow regions. Current state-of-the-art stall onset criteria use parameters that rely on flow properties integrated around the leading edge. These include the leading edge suction parameter or $LESP$ (Ramesh et al., J. Fluid Mech., vol. 751, 2014, pp. 500–538) and boundary enstrophy flux or $BEF$ (Sudharsan et al., J. Fluid Mech., vol. 935, 2022, A10), which have been found to be effective for predicting stall onset at moderate to high $Re$. However, low-$Re$ flows feature strong vortex-shedding events occurring across the entire airfoil surface, including regions away from the leading edge, altering the flow field and influencing the onset of stall. In the present work, the ability of these stall criteria to effectively capture and localize these vortex shedding events in space and time is investigated. High-resolution large-eddy simulations for an SD7003 airfoil undergoing a constant-rate, pitch-up motion at two $Re$ (10 000 and 60 000) and two pitch rates reveal a rich variety of unsteady flow phenomena, including instabilities, transition, vortex formation, merging and shedding, which are described in detail. While stall onset is reflected in both $LESP$ and $BEF$, local vortex-shedding events are identified only by the $BEF$. Therefore, $BEF$ can be used to identify both dynamic stall onset and local vortex-shedding events in space and time.

The dynamics of turbulent kinetic energy (TKE), turbulence dissipation rate (TDR) and turbulence production rate (TPR) are explored in fully developed turbulent channel flow using direct numerical simulations up to $\textit {Re}_\tau \approx 2000$ with minimal computational box for large-scale structures. Time correlation analysis based on volume-averaged TKE and TDR shows a well-defined average time lag, as in periodic/homogeneous turbulence, which, unlike periodic/homogeneous turbulence, appears to be Reynolds-number-dependent. On the basis of a spatio-temporal correlation analysis, we show that plane-averaged TKE fluctuations in the near-equilibrium region are transported towards both the core and near-wall regions, and are positively correlated with plane-averaged TDR fluctuations there with combined wall-distance and time lags. In the path towards the core region, the wall-distance lag is very close to the time lag multiplied by the friction velocity. The path towards the near-wall region has a wide spread of time lags, which increases with Reynolds number. The spatio-temporal correlation paths both towards the core and towards the wall are reproduced when the reference plane TKE is conditionally averaged on either ejections or sweeps, and are in fact stronger in correlation values in the case of ejections, which are better organised than sweeps. While volume-averaged TPR evidently precedes volume-averaged TKE, a more complex picture of non-local space–time correlations between reference plane TKE and TPR is revealed. A mechanistic model is proposed to elucidate these correlations between TKE and TPR through the interaction between the mean shear and the Reynolds shear stress.

Cavity evolution in granular media is crucial in explosion-driven gas–particle flows, particularly in many engineering applications. In this study, a theoretical model was first proposed to describe the cavity evolution in granular media by extending the classical Rayleigh–Plesset model. A closed equation set comprising the radius, pressure and gas leak-off velocity equations was built by considering the gas expansion and non-Darcy gas-penetration effects. Both centrally symmetric and non-centrally symmetric cases of gas injection into granular media were investigated. Especially for modelling the non-symmetric scenario, the radius and gas leak-off velocity equations were proposed in each radial direction with angle $\theta$, and then the pressure equation was built up based on the integral gas leak-off along the cavity outline. Through non-dimensionalizing the theoretical equations, four key dimensionless numbers $\varPi_1,\ \varPi_4$ were obtained to characterize the response time of cavity expansion and the intensity of non-Darcy effects for both cases. This allowed us to determine a scaling law of effective cavity radius $R_{eff}^*=\sqrt {2\varPi _2/(7{\rm \pi} )}t^{*1/2}$ and the critical time $t_{cr}^*=\sqrt {12.5/\varPi _1}$ for two-dimensional cavity evolution. Additionally, the necessity of incorporating non-Darcy effects was ascertained under conditions of $\varPi _4>400$. The findings demonstrate that the proposed theoretical equations effectively predict the cavity evolution results under various operational conditions ($0.7<\varPi _1<7\times 10^2, 3<\varPi _4<1.1\times 10^3$), as validated by refined Euler–Lagrange numerical simulations.