We study the behaviour of the streamwise velocity variance in turbulent wall-bounded flows using a direct numerical simulation (DNS) database of pipe flow up to friction Reynolds number ${{Re}}_{\tau } \approx 12000$. The analysis of the spanwise spectra in the viscous near-wall region strongly hints to the presence of an overlap layer between the inner- and the outer-scaled spectral ranges, featuring a $k_{\theta }^{-1+\alpha }$ decay (with $k_{\theta }$ the wavenumber in the azimuthal direction, and $\alpha \approx 0.18$), hence shallower than suggested by the classical formulation of the attached-eddy model. The key implication is that the contribution to the streamwise velocity variance $(\langle{u}^2\rangle)$ from the largest scales of motion (superstructures) slowly declines as ${{Re}}_{\tau }^{-\alpha }$, and the integrated inner-scaled variance follows a defect power law of the type $\langle u^2 \rangle ^+ = A - B \, {{Re}}_{\tau }^{-\alpha }$, with constants $A$ and $B$ depending on $y^+$. The DNS data very well support this behaviour, which implies that strict wall scaling is restored in the infinite-Reynolds-number limit. The extrapolated limit distribution of the streamwise velocity variance features a buffer-layer peak value of $\langle u^2 \rangle ^+ \approx 12.1$, and an additional outer peak with larger magnitude. The analysis of the velocity spectra also suggests a similar behaviour of the dissipation rate of the streamwise velocity variance at the wall, which is expected to attain a limiting value of approximately $0.28$, hence slightly exceeding the value $0.25$ which was assumed in previous analyses (Chen & Sreenivasan, J. Fluid Mech., vol. 908, 2021, R3). We have found evidence suggesting that the reduced near-wall influence of wall-attached eddies is likely linked to the formation of underlying turbulent Stokes layers.

In 2017, Brosseau & Vlahovska (Phys. Rev. Lett, vol. 119, no. 3, 2017, p. 034501) found that, in a strong electric field, a weakly conductive, low-viscosity droplet immersed in a highly conductive, high-viscosity medium formed a lens shape, and liquid rings continuously detached from its equatorial plane and subsequently broke up into satellite droplets. This fascinating multiphase electrohydrodynamic (EHD) phenomenon is known as droplet equatorial streaming. In this paper, based on the unified lattice Boltzmann method framework proposed by Luo et al. (Phil. Trans. R. Soc. A Math. Phys. Engng Sci, vol. 379, no. 2208, 2021, p. 20200397), a novel lattice Boltzmann (LB) model is constructed for multiphase EHD by coupling the Allen–Cahn type of multiphase LB model and two new LB equations to solve the Poisson equation of the electric field and the conservation equation of the surface charge. Using the proposed LB model, we successfully reproduced, for the first time, the complete process of droplet equatorial streaming, including the continuous ejection and breakup of liquid rings on the equatorial plane. In addition, it is found that, under conditions of high electric field strength or significant electrical conductivity contrast, droplets exhibit fingering equatorial streaming that was unknown before. A power-law relationship is discovered for droplet total charge evolution and a theoretical model is then proposed to describe the droplet radius and height over time. The breakup of liquid rings is found to be dominated by capillary instability, while the breakup of liquid fingers is governed by the end-pinching mechanism. Finally, a phase diagram is constructed for fingering equatorial streaming and ring equatorial streaming, and a criterion equation is established for the phase boundary.

In particle-laden turbulent wall flows, lift forces can influence the near-wall turbulence. This has been observed recently in particle-resolved simulations, which, however, are too expensive to be used in upscaled models. Instead, point-particle simulations have been the method of choice to simulate the dynamics of these flows during the last decades. While this approach is simpler, cheaper and physically sound for small inertial particles in turbulence, some issues remain. In the present work, we address challenges associated with lift force modelling in turbulent wall flows and the impact of lift forces in the near-wall flow. We performed direct numerical simulations of small inertial point particles in turbulent channel flow for fixed Stokes number and mass loading while varying the particle size. Our results show that the particle dynamics in the buffer region, causing the apparent particle-to-fluid slip velocity to vanish, raises major challenges for modelling lift forces accurately. While our results confirm that lift forces have little influence on particle dynamics for sufficiently small particle sizes, for inner-scaled diameters of order one and beyond, lift forces become quite important near the wall. The different particle dynamics under lift forces results in the modulation of streamwise momentum transport in the near-wall region. We analyse this lift-induced turbulence modulation for different lift force models, and the results indicate that realistic models are critical for particle-modelled simulations to correctly predict turbulence modulation by particles in the near-wall region.

The orientational dynamics of a spherical magnetic particle in linear shear flow subjected to an oscillating magnetic field in the flow plane is analysed in the viscous limit. The shear is in the $X$–$Y$ plane, the magnetic field is in the $X$ direction and the vorticity is perpendicular to the flow in the $Z$ direction. The relevant dimensionless groups are $\omega ^\ast$, the ratio of the frequency of the magnetic field and the strain rate, and $\varSigma$, the ratio of the magnetic and hydrodynamic torques. As $\varSigma$ is decreased, there is a transition from in-plane rotation, where the rotation is in the flow ($X$–$Y$) plane, to out-of-plane rotation, where the orientation vector is not necessarily in the $X$–$Y$ plane and the dynamics depends on the initial orientation. The particle rotation is phase-locked for in-plane rotation with discrete odd rotation number (number of rotations in one period of magnetic field oscillation), while the orbits are quasi-periodic with non-integer rotation number for out-of-plane rotation. For $\varSigma \gg 1$, regions of odd rotation number $n_o$ are bound by the lines $8 (n_o-1) \varSigma \omega ^\ast = 1$ and $8 (n_o+1) \varSigma \omega ^\ast = 1$, and there are discontinuous changes in the rotation number and mean and root-mean-square torque at these lines. For $\varSigma \ll 1$, the domains of in-plane rotation of finite width in the $\omega ^\ast$–$\varSigma$ plane extend into downward cusps at $\omega ^\ast = {1}/{2 n_o}$. The orbits are quasi-periodic between these domains, where the rotation is out of plane.

Coherent small-amplitude unsteadiness of the shock wave and the separation region over a canonical double cone flow, termed in literature as oscillation-type unsteadiness, is experimentally studied at Mach 6. The double cone model is defined by three non-dimensional geometric parameters: fore- and aft-cone angles ($\theta _1$ and $\theta _2$), and ratio of the conical slant lengths ($\varLambda$). Previous studies of oscillations have been qualitative in nature, and mostly restricted to a special case of the cone model with fixed $\theta _1 = 0^\circ$ and $\theta _2 = 90^\circ$ (referred to as the spike-cylinder model), where $\varLambda$ becomes the sole governing parameter. In the present effort we investigate the self-sustained flow oscillations in the $\theta _1$-$\varLambda$ parameter space for fixed $\theta _2 = 90^\circ$ using high-speed schlieren visualisation. The experiments reveal two distinct subtypes of oscillations, characterised by the motion (or lack thereof) of the separation point on the fore-cone surface. The global time scale associated with flow oscillation is extracted using spectral proper orthogonal decomposition. The non-dimensional frequency (Strouhal number) of oscillation is seen to exhibit distinct scaling for the two oscillation subtypes. The relationship observed between the local flow properties, instability of the shear layer, and geometric constraints on the flow suggests that an aeroacoustic feedback mechanism sustains the oscillations. Based on this understanding, a simple model with no empiricism is developed for the Strouhal number. The model predictions are found to match well with experimental measurements. The model provides helpful physical insight into the nature of the self-sustained flow oscillations over a double cone at high speeds.

Emptying–filling boxes have been studied in a wide range of configurations for decades, but the flow created in the box by two plumes rising from sources of arbitrary strength and elevation was previously unsolved. Guided by experiments and simplified analytical modelling, we reveal a rich array of two- and three-layer stratifications across seven possible flow regimes. The governing equations for these regimes show how the prevailing regime and stratification properties vary with three key parameters: the relative strength of the plumes, the height difference between their sources and a parameter characterising the resistance of the box to emptying. We observe and explain new behaviours not described in previous studies that are crucial to understanding emptying–filling boxes with multiple plumes. In particular, we demonstrate that the oft-assumed premise that $n$ plumes leads to a stratification with $n+1$ layers is not necessarily true, even in the absence of mixing. Two emptying–filling box models are developed: an analytical model addressing all combinations of the governing parameters and an extended model for three-layer stratifications that incorporates two mixing processes observed in the experiments. The predictions of these two models are in generally excellent agreement with measurements from the experimental campaign covering 69 combinations of the governing parameters. This study improves our understanding of emptying–filling boxes and could facilitate improvements to natural ventilation building design, as demonstrated by an example scenario in which occupants feel cooler upon the addition of a second source of heat.

The non-Oberbeck–Boussinesq effects on the stability of a vertical natural convection boundary layer are investigated using the linearised disturbance equations for air flows up to a temperature difference of $\Delta T=100\,{\rm K}$. Based on the linear stability results, the neutral curve is shown to be sensitive to the choice of reference temperature. When evaluated using the film temperature $T_f$, a lower film Grashof number is required to trigger the linear instability for larger $\Delta T$. The relative contributions of shear and buoyant production to the perturbation kinetic energy budget reveals that the marginally unstable modes are amplified based on different mechanisms: for lower wavenumbers at relatively small Grashof number, the instability is driven by buoyancy; whereas for higher wavenumbers and larger Grashof number, the flow becomes unstable due to a shear instability. The use of reference temperature is found to scale the shear- and buoyant-driven instabilities differently so that no single reference temperature definition would collapse the neutral curves. The linear stability result further demonstrates that at a given Grashof number a higher temperature difference would give a larger amplification rate of the perturbation, which then leads to an earlier onset of the nonlinearities when evaluated at $T_f$. Finally, by comparing the amplification rates obtained from direct numerical simulation and the linear stability results, the extent of the linear regime is determined for $\Delta T = 100\,{\rm K}$.

In this study we propose a novel data-driven reduced-order model for complex dynamics, including nonlinear, multi-attractor, multi-frequency and multiscale behaviours. The starting point is a fully automatable cluster-based network model (CNM) (Li et al., J. Fluid Mech., vol. 906, 2021, A21) that kinematically coarse grains the state with clusters and dynamically predicts the transitions in a network model. In the proposed dynamics-augmented CNM (dCNM) the prediction error is reduced with trajectory-based clustering using the same number of centroids. The dCNM is first exemplified for the Lorenz system and then demonstrated for the three-dimensional sphere wake featuring periodic, quasi-periodic and chaotic flow regimes. For both plants, the dCNM significantly outperforms the CNM in resolving the multi-frequency and multiscale dynamics. This increased prediction accuracy is obtained by stratification of the state space aligned with the direction of the trajectories. Thus, the dCNM has numerous potential applications to a large spectrum of shear flows, even for complex dynamics.

The evolution of the water-entry cavity affects the impact load and the motion of the body. This paper adopts the Eulerian finite element method for multiphase flow for simulations of the high-speed water-entry process. The accuracy and convergence of the numerical method are verified by comparing it with the experimental data and the results of the transient cavity dynamics theory. Based on the results, the representative characteristics of the cavity are discussed from the perspective of the cavity cross-section. It is found that the asymmetry of the cavity expansion and contraction durations is related to the motion of the free surface and the closure of the cavity. The uplift of the free surface suppresses cavity expansion, while the jet generated from free surface closure accelerates cavity contraction. The duration of the contraction of the cavity near the free surface is shorter than the expansion duration due to the change in the velocity distribution caused by the free surface motion. The necking phenomenon during deep closure leads to an increase in the internal pressure of the cavity, prolonging cavity contraction near the deep closure area. This work provides new insights into the cavity dynamics in high-speed water entry.

While the problem governing Stokes flow about a single particle that is subject to an external force is ill posed in two dimensions (the ‘Stokes paradox’), the related problem of two mutually repellent particles is well posed. Motivated by self-assembly phenomena in thin viscous membranes, we consider this problem in the limit of remote particles. Such limits are typically handled in the literature using reflection techniques, which provide successive approximations to the mutual hydrodynamic interactions. Since their starting point is a single particle in an unbounded fluid domain, these techniques are futile in the present two-dimensional problem. We show how this apparent contradiction is resolved via use of singular perturbations. We obtain a two-term approximation for the velocity acquired by circular disks, considering both rigid and free particle surfaces. We also illustrate our perturbation scheme for elliptic disks, deriving a renormalised single-particle velocity. The utility of our asymptotic scheme is illustrated in the general problem of hydrodynamic interaction between a cluster of remote disks.

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$.

We explore predictions of two models of one-dimensional capillary rise in rigid and partially saturated porous media. One is an existing one from the literature and the second is a free-boundary model based on Richards’ equation with two moving boundaries of the evolving partially saturated region. Both models involve the specification of saturation-dependent functions for local capillary pressure and permeability and connect to classical models for saturated porous media. Existing capillary-rise experiments show two notable regimes: (i) an early-time regime typically well-described by classical capillary-rise theory in a fully saturated porous media, and (ii) a long-time regime that has anomalous dynamics in which the capillary-rise height may scale with a non-classical power law in time or have more complicated dynamics. We demonstrate that the predictions of both models compare well with experimental capillary-rise data over early- and long-time regimes gathered from three independent studies in the literature. The model predictions also shed light on recent scaling laws that relate the capillary pressure and permeability of the partially saturated media to the capillary-rise height. We use these models to probe computationally observed permeability relationships to capillary-rise height. We demonstrate that a recently proposed permeability scaling for the anomalous capillary-rise regime is indeed realized and is particularly apparent in the lower portion of the partially saturated media. For our free-boundary model we also compute capillary pressure measures and show that these reveal the linear relation between the capillary pressure and capillary-rise height expected for a capillarity–gravity balance in the upper portion of the partially saturated porous media.

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.

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.

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.

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.

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.

Although the roll-streak (R-S) is fundamentally involved in the dynamics of wall turbulence, the physical mechanism responsible for its formation and maintenance remains controversial. In this work we investigate the dynamics maintaining the R-S in turbulent Poiseuille flow at $R=1650$. Spanwise collocation is used to remove spanwise displacement of the streaks and associated flow components, which isolates the streamwise-mean flow R-S component and the second-order statistics of the streamwise-varying fluctuations that are collocated with the R-S. This partition of the dynamics into streamwise-mean and fluctuation components facilitates exploiting insights gained from the analytic characterization of turbulence in the second-order statistical state dynamics (SSD), referred to as S3T, and its closely associated restricted nonlinear dynamics (RNL) approximation. Symmetry of the statistics about the streak centreline permits separation of the fluctuations into sinuous and varicose components. The Reynolds stress forcing induced by the sinuous and varicose fluctuations acting on the R-S is shown to reinforce low- and high-speed streaks, respectively. This targeted reinforcement of streaks by the Reynolds stresses occurs continuously as the fluctuation field is strained by the streamwise-mean streak and not intermittently as would be associated with streak-breakdown events. The Reynolds stresses maintaining the streamwise-mean roll arise primarily from the dominant proper orthogonal decomposition (POD) modes of the fluctuations, which can be identified with the time average structure of optimal perturbations growing on the streak. These results are consistent with a universal process of R-S growth and maintenance in turbulent shear flow arising from roll forcing generated by straining turbulent fluctuations, which was identified using the S3T SSD.

We investigate by direct numerical simulations the fluid–solid interaction of non-dilute suspensions of spherical particles moving in triperiodic turbulence, at the relatively large Reynolds number of $Re_\lambda \approx 400$. The solid-to-fluid density ratio is varied between $1.3$ and $100$, the particle diameter $D$ is in the range $16 \le D/\eta \le 123$ ($\eta$ is the Kolmogorov scale) and the volume fraction of the suspension is $0.079$. Turbulence is sustained using the Arnold–Beltrami–Childress cellular-flow forcing. The influence of the solid phase on the largest and energetic scales of the flow changes with the size and density of the particles. Light and large particles modulate all scales in an isotropic way, while heavier and smaller particles modulate the largest scales of the flow towards an anisotropic state. Smaller scales are isotropic and homogeneous for all cases. The mechanism driving the energy transfer across scales changes with the size and the density of the particles. For large and light particles the energy transfer is only marginally influenced by the fluid–solid interaction. For small and heavy particles, instead, the classical energy cascade is subdominant at all scales, and the energy transfer is essentially driven by the fluid–solid coupling. The influence of the solid phase on the flow intermittency is also discussed. Besides, the collective motion of the particles and their preferential location in relation to properties of the carrier flow are analysed. The solid phase exhibits moderate clustering; for large particles the level of clustering decreases with their density, while for small particles it is maximum for intermediate values.

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.