In homogeneous and isotropic turbulence, the relative contributions of different physical mechanisms to the energy cascade can be quantified by an exact decomposition of the energy flux (Johnson, Phys. Rev. Lett., vol. 124, 2020, 104501; J. Fluid Mech., vol. 922, 2021, A3). We extend the formalism to the transfer of kinetic helicity across scales, important in the presence of large-scale mirror‐breaking mechanisms, to identify physical processes resulting in helicity transfer and quantify their contributions to the mean flux in the inertial range. All subfluxes transfer helicity from large to small scales. Approximately 50 % of the mean flux is due to the scale-local vortex flattening and vortex twisting. We derive a new exact relation between these effects, similar to the Betchov relation for the energy flux, revealing that the mean contribution of the former is three times larger than that of the latter. Multi-scale effects account for the remaining 50 % of the mean flux, with approximate equipartition between multi-scale vortex flattening, twisting and entangling.

In this note we provide an analytical proof of the zero helicity condition for systems governed by the Gross–Pitaevskii equation (GPE). The proof is based on the hydrodynamic interpretation of the GPE, and the direct use of Noether's theorem by applying Kleinert's multi-valued gauge theory. As a by-product we also demonstrate the conservation and quantization of the circulation for the GPE.

The transition to turbulence in conduits is among the longest-standing problems in fluid mechanics. Challenges in producing or saving energy hinge on understanding promotion or suppression of turbulence. While a global picture based on an intrinsically 3-D subcritical mechanism is emerging for 3-D turbulence, subcritical turbulence is yet to even be observed when flows approach two dimensions, e.g. under intense rotation or magnetic fields. Here, stability analysis and direct numerical simulations demonstrate a subcritical quasi-two-dimensional (quasi-2-D) transition from laminar flow to turbulence, via a radically different 2-D mechanism to the 3-D case, driven by nonlinear Tollmien–Schlichting waves. This alternative scenario calls for a new line of thought on the transition to turbulence and should inspire new strategies to control transition in rotating devices and nuclear fusion reactor blankets.

We have performed particle-resolved direct numerical simulations of many heavy non-spherical particles settling under gravity in the dilute regime. The particles are oblate spheroids of aspect ratio $1.5$ and density ratio $1.5$. Two Galileo numbers are considered, namely $111$ and $152$, for which a single oblate spheroid follows a steady vertical and a steady oblique path, respectively. In both cases, a strongly inhomogeneous spatial distribution of the disperse phase in the form of columnar clusters is observed, with a significantly enhanced average settling velocity as a consequence. Thus, in contrast to previous results for spheres, the qualitative difference in the single-particle regime does not result in a qualitatively different behaviour of the many-particle cases. In addition, we have carried out an analysis of pairwise interactions of particles in the well-known drafting–kissing–tumbling set-up, for oblate spheroids of aspect ratio $1.5$ and for spheres. We have varied systematically the relative initial position between the particle pair and we have considered free-to-rotate particles and rotationally locked ones. We have found that the region of attraction for both particle shapes, with and without rotation, is very similar. However, significant differences occur during the drafting and tumbling phases. In particular, free-to-rotate spheres present longer drafting phases and separate quickly after the collision. Spheroids remain close to each other for longer times after the collision, and free-to-rotate ones experience two or more collision events. Therefore, we have observed a shape-induced increase in the interaction time which might explain the increased tendency to cluster of the many-particle cases.

Ambient air cushions the impact of drops on solid substrates, an effect usually revealed by the entrainment of a bubble, trapped as the air squeezed under the drop drains and liquid–solid contact occurs. The presence of air becomes evident for impacts on very smooth surfaces, where the gas film can be sustained, allowing drops to bounce without wetting the substrate. In such a non-wetting situation, Mandre & Brenner (J. Fluid Mech., vol. 690, 2012, p. 148) numerically and theoretically evidenced that two physical mechanisms can act to prevent contact: surface tension and nonlinear advection. However, the advection dominated regime has remained hidden in experiments as liquid–solid contact prevents rebounds being realised at sufficiently large impact velocities. By performing impacts on superheated surfaces, in the so-called dynamical Leidenfrost regime (Tran et al., Phys. Rev. Lett., vol. 108, issue 3, 2012, p. 036101), we enable drop rebound at higher impact velocities, allowing us to reveal this regime. Using high-speed total internal reflection, we measure the minimal gas film thickness under impacting drops, and provide evidence for the transition from the surface tension to the nonlinear inertia dominated regime. We rationalise our measurements through scaling relationships derived by coupling the liquid and gas dynamics, in the presence of evaporation.

This work proposes a machine-learning or artificial intelligence (AI) control of a low-drag Ahmed body with a rear slant angle φ = 35° with a view to finding strategies for efficient drag reduction (DR). The Reynolds number Re investigated is 1.7 × 105 based on the square root of the body cross-sectional area. The control system comprises of five independently operated arrays of steady microjets blowing along the edges of the rear window and vertical base, twenty-six pressure taps on the rear end of the body and a controller based on an ant colony algorithm for unsupervised learning of a near-optimal control law. The cost function is designed such that both DR and control power input are considered. The learning process of the AI control discovers forcing that produces a DR up to 18 %, corresponding to a drag coefficient reduction of 0.06, greatly exceeding any previously reported DR for this body. Furthermore, the discovered forcings may provide alternative solutions, i.e. a tremendously increased control efficiency given a small sacrifice in DR. Extensive flow measurements performed with and without control indicate significant alterations in the flow structure around the body, such as flow separation over the rear window, recirculation bubbles and C-pillar vortices, which are linked to the pressure rise on the window and base. The physical mechanism for DR is unveiled, along with a conceptual model for the altered flow structure under the optimum control or biggest DR. This mechanism is further compared with that under the highest control efficiency.

The behaviour of granular column collapses is associated with the dynamics of geohazards, such as debris flows, landslides and pyroclastic flows, yet their underlying physics is still not well understood. In this paper, we explore granular column collapses using the sphero-polyhedral discrete element method, where the system contains two types of particles with different frictional properties. We impose three different mixing ratios and multiple different particle frictional coefficients, which lead to different run-out distances and deposition heights. Based on our previous work and a simple mixture theory, we propose a new effective initial aspect ratio for the bi-frictional granular mixture, which helps unify the description of the relative run-out distances. We analyse the kinematics of bi-frictional granular column collapses and find that deviations from classical power-law scaling in both the dimensionless terminal time and the dimensionless time when the system reaches the maximum kinetic energy may result from differences in the initial solid fraction and initial structures. To clarify the influence of initial states, we further decrease the initial solid fraction of granular column collapses, and propose a trial function to quantitatively describe its influence. Due to the utilization of a simple mixture theory of contact occurrence probability, this study can be associated with the friction-dependent rheology of granular systems and friction-induced granular segregations, and further generalized to applications with multiple species of particles in various natural and engineering mixtures.

Through direct numerical simulations we investigate the effects of spatiotemporal intermittency as a result of stable stratification in surface heated stratified open channel flow. By adapting the density inversion criterion method of Portwood et al. [J. Fluid Mech., vol. 807, 2016, R2] for our flow, we demonstrate that the flow may be robustly separated into regions of active turbulence for which $Re_B \gtrsim {O}(10)$ and surrounding quiescent fluid, where $Re_B$ is the buoyancy Reynolds number. The intermittency in the flow spontaneously manifests as a deformed horizontal interface between the upper quiescent and lower turbulent flow, characterised by vigorous mixing from ‘overturning’ shear instabilities. The resulting vertical intermittency profile is accurately predicted by a local Monin–Obukhov length normalised by viscous wall units $\varLambda ^+$ such that the flow displays intermittency within the parameter range of $2.5 \lesssim \varLambda ^+ \lesssim 260$. By considering conditional averages of the ‘turbulent’ and ‘quiescent’ flow separately, we find the ‘turbulent’ flow within this region to be described by constant critical gradient Richardson and turbulent Froude numbers of $Ri_{g,c} \approx 0.2$ and $Fr_c \approx 0.3$. We find that the turbulent flow continues to display a $\varGamma \sim Fr^{-1}$ relationship when $Fr < Fr_c$, whereas the quiescent flow shows no correlation between $\varGamma$ and $Fr$, where $\varGamma$ is the flux coefficient. Hence, we demonstrate directly that for our flow, the emergence of an asymptotic ‘saturated’ $\varGamma$ regime in the limit of a low ‘global’ $Fr$ occurs due to intermittency and increasing contributions to measurements of $\varGamma$ from the quiescent flow.

Turbulence modulation by finite-size particles in homogeneous isotropic turbulence (HIT) has been investigated numerically and experimentally in many studies, but its controlling parameters are not fully clear. In this work, four non-dimensional parameters governing the turbulent modulation by non-settling particles, i.e. $Re_\lambda$ of the background HIT, the particle-to-fluid density ratio $\rho _p/\rho _f$, the relative particle size $d_p/\eta$ and the particle volume fraction $\phi _v$, are identified through dimensional analysis. Then, a parameterization study is conducted based on results from fully resolved direct numerical simulations to investigate the influence of the above non-dimensional parameters on the modulation of turbulent kinetic energy (TKE) and viscous dissipation rate. Empirical models that quantitatively predict the modulation of TKE and dissipation rate are then developed by fitting in the simulation results. These models are used to examine the turbulence modulation results reported in the literature. The model predictions and the data points of TKE modulation show reasonable agreement, but the model predicting the modulation of dissipation rate needs further deliberation as the credibility of the available data points is currently difficult to assess. The generality and the physics behind these empirical models also require further investigation.

Kinetic models of polyatomic gas typically account for the internal degrees of freedom at the level of the two-particle distribution function. However, close to the hydrodynamic limit, the internal (rotational) degrees of freedom tend to be well represented just by rotational kinetic energy density. We account for the rotational energy by augmenting the ellipsoidal statistical Bhatnagar–Gross–Krook (ES–BGK) model, an extension of the BGK model, at the level of the single-particle distribution function with an advection–diffusion–relaxation equation for the rotational energy. This reduced model respects the $H$ theorem and recovers the compressible hydrodynamics for polyatomic gases as its macroscopic limit. As required for a polyatomic gas model, this extension of the ES–BGK model not only has the correct specific heat ratio but also allows for three independent tunable transport coefficients: thermal conductivity, shear viscosity and bulk viscosity. We illustrate the effectiveness of the model via a lattice Boltzmann method implementation.

It is known that, apart from the growth of instability modes, the non-modal (algebraic) mechanism of linear growth plays an important role in near-wall flows. In unbounded flows, including submerged jets, the theoretical analysis of the non-modal growth mechanism started only in the last decade; this mechanism has not yet been identified in experiments. In the present work, experiments were conducted on the excitation of a non-modal ‘lift-up’ growth mechanism. Special wavy structures (deflectors) were introduced into a laminar round submerged jet, which excited a roller-like transverse motion. Based on experimental results, we definitely identify the non-modal ‘lift-up’ growth mechanism of introduced disturbances. The development of perturbations in the experiment qualitatively corresponds to the theoretically calculated optimal perturbations. The features of the transition to turbulence caused by non-modal growth are considered.

Magneto-gravity-elliptic instability is addressed here considering an unbounded strained vortex (with constant vorticity $2\varOmega$ and with ellipticity parameter $\varepsilon$) of a perfectly conducting fluid subjected to a uniform axial magnetic field (with Alfvén velocity scaled from the basic magnetic field $B)$ and an axial stratification (with a constant Brunt–Väisälä frequency $N$). Such a simple model allows us to formulate the stability problem as a system of equations for disturbances in terms of Lagrangian Fourier (or Kelvin) modes which is universal for wavelengths of the perturbation sufficiently small with respect to the scale of variation of the basic velocity gradients. It can model localised patches of elliptic streamlines which often appear in some astrophysical flows (stars, planets and accretion discs) that are tidally deformed through gravitational interaction with other bodies. In the limit case where the flow streamlines are exactly circular ($\varepsilon =0),$ there are fast and slow magneto-inertia-gravity waves with frequencies $\omega _{1,2}$ and $\omega _{3,4},$ respectively. Under the effect of finite ellipticity, the resonant cases of these waves, $\omega _i-\omega _j=n\varOmega$ $(i\ne j)$ ($n$ being an integer), can become destabilising. The maximal growth rate of the subharmonic instability (related to the resonance of order $n=2)$ is determined by extending the asymptotic method by Lebovitz & Zweibel (Astrophys. J., vol. 609, 2004, pp. 301–312). The domains of the $(k_0B/\varOmega, N/\varOmega )$ plane for which this instability operates are identified ($1/k_0$ being a characteristic length scale). We demonstrate that the $N\rightarrow 0$ limit is, in fact, singular (discontinuous). The axial stable stratification enhances the subharmonic instability related to the resonance between two slow modes because, at large magnetic field strengths, its maximal growth rate is twice that found in the case without stratification.

Although superhydrophobic surfaces (SHSs) show promise for drag reduction applications, their performance can be compromised by traces of surfactant, which generate Marangoni stresses that increase drag. This question is addressed for soluble surfactant in a three-dimensional laminar channel flow, with periodic SHSs made of long finite-length longitudinal grooves located on both walls. We assume that bulk diffusion is sufficiently strong for cross-channel concentration gradients to be small. Exploiting long-wave theory and accounting for the difference between the rapid transverse and slower longitudinal Marangoni flows, we derive a one-dimensional model for surfactant transport from the full three-dimensional transport equations. Our one-dimensional model allows us to predict the drag reduction and surfactant distribution across the parameter space. The system exhibits multiple regimes, involving competition between Marangoni effects, bulk and interfacial diffusion, bulk and interfacial advection, shear dispersion and surfactant exchange between the bulk and the interface. We map out asymptotic regions in the high-dimensional parameter space, and derive explicit closed-form approximations of the drag reduction, without any fitting or empirical parameters. The physics underpinning the drag reduction effect and the negative effect of surfactant is discussed through analysis of the velocity field and surfactant concentrations, which show both uniform and non-uniform stress distributions. Our theoretical predictions of the drag reduction compare well with results from the literature solving numerically the full three-dimensional transport problem. Our atlas of maps provides a comprehensive analytical guide for designing surfactant-contaminated channels with SHSs, to maximise the drag reduction in applications.

We propose a physical model to predict indirect noise generated by the acceleration of compositional inhomogeneities in nozzles with viscous dissipation (non-isentropic nozzles). First, we derive the quasi-one-dimensional equations from the conservation laws of multicomponent flows. Second, we validate the proposed model with the experimental data available in the literature for binary mixtures of four gases. Third, we calculate the nozzle transfer functions for different Helmholtz numbers and friction factors, in both subsonic and supersonic flows with/without shock waves. We show that friction and dissipation have a significant effect on the generation of indirect noise, for which the physical mechanism is identified and explained. Fourth, we find a semi-analytical solution with path integrals, which provides an asymptotic expansion with respect to the Helmholtz number. Fifth, we introduce the compositional-noise scaling factor, which is applied to quickly estimate compositional noise from the knowledge of only one single-component gas transfer function. The approximation error is less than $1\,\%$. The proposed low-order model provides accurate estimates of the transfer functions and physical insight into indirect noise for multicomponent gases. This opens up new possibilities to accurately predict, and understand, sound generation in gas turbines.

Water stuck in the ear is a common problem during showering, swimming or other water activities. Having water trapped in the ear canal for a long time can lead to ear infections and possibly result in hearing loss. A common strategy for emptying water from the ear canal is to shake the head, where high acceleration helps remove the water. In this present study, we rationalize the underlying mechanism of water ejection/removal from the ear canal by performing experiments and developing a stability theory. From the experiments, we measure the critical acceleration to remove the trapped water inside different sizes of canals. Our theoretical model, modified from the Rayleigh–Taylor instability, can explain the critical acceleration observed in experiments, which strongly depends on the radius of the ear canal. The resulting critical acceleration tends to increase, especially in smaller ear canals, which indicates that shaking heads for water removal can be more laborious and potentially threatening to children due to their small size of the ear canal compared with adults.

The presence of large-scale coherent structures can significantly impact the dynamics of a turbulent flow field and the behaviour of a flame stabilized in that flow. The goal of this study is to analyse how increasing free-stream turbulence can change the response of the flow to longitudinal acoustic excitation of varying amplitudes. We study the flow in the wake of a cylindrical bluff body at both non-reacting and reacting conditions, as the presence of a flame can significantly alter the global stability of the flow. The frequency of longitudinal acoustic excitation is set to match the natural frequency of anti-symmetric vortex shedding for this configuration and we vary the free-stream turbulence using perforated plates upstream of the bluff body. The results show that varying the level of free-stream turbulence can influence not only the amplitude of the coherent flow response, but also the symmetry of vortex shedding in the presence of longitudinal acoustic excitation. Increasing the turbulence intensity can fundamentally change the structure of the time-averaged flow and can directly impact the coherent flow response in two ways. First, increasing turbulence intensity can enhance the amplitude of the natural anti-symmetric vortex shedding mode in the wake. Second, increasing turbulence intensity weakens the symmetric response of the flow to the longitudinal acoustic excitation. In the non-reacting and reacting conditions, both symmetric and anti-symmetric modes are present and are characterized using a spectral proper orthogonal decomposition. We see evidence of interaction between the symmetric and anti-symmetric modes, which leads to an interference pattern in the coherent vorticity response in the shear layers. We conclude by presenting a conceptual model for the influence that turbulence has on these flows.

It is known that the water entry of a body with a recessed, cupped nose can suppress the splash and air cavity typically observed for solid body entry (Mathai, Govardhan & Arakeri, Appl. Phys. Lett., vol. 106, 2015, 064101). However, the interplay between the captive gas in the cup, the cavity and the splash is quite subtle and has not been thoroughly explored. Here we study the cavity and splash dynamics associated with the vertical water entry of cups and find a variety of regimes over a range of Weber numbers ($We_D$) and dimensionless cup depths. Our parameter space spans a transition between slow-developing cavities with long closure times (low $We_D$) to fast-sealing cavities (high $We_D$). An important dynamic event is the evacuation of trapped gas from within the cup, which drives the ensuing cavity and splash behaviour. Through modelling, we predict the conditions for which the evacuating gas inflates a cavity that opens to the atmosphere versus inflating a submerged cavity that suppresses air entrainment from above the surface. We also compare our cup water entry findings to the impact phenomena observed for flat disks, which entrap gas on the front surface similar to cups. In doing so, we reveal the sensitivity of disk splash and cavity behaviour to impact angle, and show that disks share a common regime with cups, in which a thin splash quickly seals on the body. We deduce the mechanisms by which increasing cup depth delays the cavity seal time in this regime. These findings reveal that cups may in fact promote or suppress cavity growth, depending on the cup depth and impact conditions.

Dust storms are typical dispersed two-phase atmospheric turbulence involving electrified charged dust particles. Previous observations have demonstrated that clean-air atmospheric turbulence is strongly intermittent. However, the intermittency of the wind velocity, concentration of dust particles with a diameter smaller than $10\ \mathrm {\mu }{\rm m}$ (PM10) and electric fields, known as multifield intermittency, has not been reported or characterized yet. Here, we quantify the small-scale multifield intermittency of dust storms using datasets obtained from the Qingtu Lake Observation Array and a wavelet-based data analysis technique. The results indicate that the probability density functions of the multifield increments are scale dependent, and the scaling exponents of the multifield structure functions exhibit anomalous scaling, suggesting that the multiple fields in dust storms are also highly intermittent. Specifically, the wind velocity during dust storms appears to be more intermittent as compared with clean-air conditions. Among the multiple fields, the small-scale intermittency is strongest for PM10 dust concentration, moderate for electric fields and weakest for wind velocity. Furthermore, the anomalous scaling of multiple fields is well described by the hierarchical structure theory of turbulence. It is theoretically predicted that the wind velocity displays a one-dimensional filamentary structure, while the PM10 dust concentration and electric fields display two-dimensional sheet-like structures. Finally, after removing the coherent components of the observed time series by the proposed wavelet conditioning statistics, Kolmogorov linear scaling is recovered for the multiple fields, suggesting that small-scale multifield intermittency is caused by the presence of small-scale coherent structures.

Flow around a tethered model of a swimming batoid fish is studied by using the wall-modelled large-eddy simulation in conjunction with the immersed boundary method. A Reynolds number ($Re$) up to 148 000 is chosen, and it is comparable to that of a medium-sized aquatic animal in cruising swimming state. At such a high $Re$, we provide, to the best of our knowledge, the first evidence of hairpin vortical (HV) structures near the body surface using three-dimensional high-fidelity flow field data. It is observed that such small-scale vortical structures are mainly formed through two mechanisms: the leading-edge vortex (LEV)–secondary filament–HV and LEV–HV transformations in different regions. The HVs create strong fluctuations in the pressure distribution and frequency spectrum. Simulations are also conducted at $Re=1480$ and 14 800 to reveal the effect of Reynolds number. Variations of the flow separation behaviour and local pressure with $Re$ are presented. Our results indicate that low-$Re$ simulations are meaningful when the focus is on the force variation tendency, whereas high-$Re$ simulations are needed when concerning flow fluctuations and turbulence mechanisms.

Eccentric compound drops, which are ubiquitous in many naturally inspired and engineering systems, can migrate under the sole presence of a uniform electric field, unlike the case of isolated single drops. Here, we report the migration of eccentric compound drops under a uniform electric field, imposed parallel to the line of centres of the constituting drops, by developing an approximate analytical model that applies to low Reynolds number limits under negligible droplet deformation, following axisymmetric considerations. In contrast to the sole influence of the electrohydrodynamic forces that has thus far been established to be emphatic for the eccentric configuration, here we report the additional effects induced by the dielectrophoretic forces to result in decisive manipulation in the drop migration. We show that the relative velocity between the inner and outer drops, which is a function of the eccentricity itself, dictates the dynamical evolution of the eccentricity variation under the competing electrohydrodynamic and dielectrophoretic interactions. This brings out four distinct regimes of the migration characteristics of the two drops based on their relative electro-physical properties. Our results reveal that an increase in eccentricity and the size ratio of the inner and outer droplets may induce monotonic or non-monotonic variation in the drop velocities, depending on the operating regime. We show how the interplay of various properties holds the control of selectively increasing or suppressing the eccentricity with time. These findings open up various avenues of electrically manipulative motion of encapsulated fluidic phases in various applications encompassing engineering and biology.