In this study, direct numerical simulation of the particle dispersion and turbulence modulation in a sonic transverse jet injected into a supersonic cross-flow with a Mach number of 2 was carried out with the Eulerian–Lagrangian point-particle method. One single-phase case and two particle-laden cases with different particle diameters were simulated. The jet and particle trajectories, the dispersion characteristics of particles, and the modulation effect of particles on the flow were investigated systematically. It was found that large particles primarily accumulate around shear layer structures situated on the windward side of the jet trajectory. In contrast, small particles exhibit radial transport, accessing both upstream and downstream recirculation zones. Moreover, small particles disperse extensively within the boundary layer and large-scale shear layers, evidently influenced by the streamwise vortices. The particles increase the mean wall-normal velocity near the wall in the wake region of the transverse jet, while reducing the mean streamwise and wall-normal velocities in outer regions. Particles significantly alter the flow velocity adjacent to shock fronts. In particular, the turbulent fluctuations near the windward barrel shock and bow shock are reduced, while those around the leeward barrel shock are increased. An upward displacement of the bow shock in the wall-normal direction is also observed due to particles. In the regions away from the shocks, small particles tend to amplify the Reynolds stress, while large particles attenuate the turbulent kinetic energy.

This work aims to perform a parametric study on a round supersonic jet with a design Mach number Md = 1.8, which is manipulated using a single steady radial minijet with a view to enhancing its mixing. Four control parameters are examined, i.e. the mass flow rate ratio Cm and diameter ratio d/D of the minijet to main jet, and exit pressure ratio Pe/Pa and fully expanded jet Mach number Mj, where Pe and Pa are the nozzle exit and atmospheric pressures, respectively. Extensive pressure and schlieren flow visualization measurements are conducted on the natural and manipulated jets. The supersonic jet core length Lc/D exhibits a strong dependence on the four control parameters. Careful scaling analysis of experimental data reveals that Lc/D = f1(Cm, d/D, Pe/Pa, Mj) may be reduced to Lc/D = f2(ξ), where f1 and f2 are different functions. The scaling factor $\xi = J({d_i}/{D_j})/(\gamma M_j^2{P_e}/{P_a})$ is physically the penetration depth of the minijet into the main jet, where $J({d_i}/{D_j})$ is the square root of the momentum ratio of the minijet to main jet (di and Dj are the fully expanded diameters of d and D, respectively), γ is the specific heat ratio and $\gamma M_j^2{P_e}/{P_a}$ is the non-dimensional exit pressure ratio. Important physical insight may be gained from this scaling law into the optimal choice of control parameters such as d/D and Pe/Pa for practical applications. It has been found for the first time that the minijet may induce a street of quasi-periodical coherent structures once Cm exceeds a certain level for a given ${P_e}/{P_a}$. Its predominant dimensionless frequency Ste (≡ feDj/Uj) scales with a factor $\zeta = J({d_i}/{D_j})\; \sqrt {\gamma M_j^2{P_e}/{P_a}} $, which is physically the ratio of the minijet momentum thrust to the ambient pressure thrust. The formation mechanism of the street and its role in enhancing jet mixing are also discussed.

We consider the initial ‘slumping phase’ of a lock-release gravity current (GC) on a down slope with focus on particle-driven (turbidity) flows, in the inertia–buoyancy (large Reynolds number) and Boussinesq regime. We use a two-layer shallow-water (SW) model for the depth-averaged variables, and compare the predictions with previously published experimental data. In particular, we analyse the empirical conclusion of Gadal et al. (J. Fluid Mech., vol. 974, 2023, A4) that the slumping displays a constant speed for a significant range of slopes and particle-sedimentation speeds. We emphasize the physical definition of the slumping phase (stage): the adjustment process during which (a) the fluid in the lock is set into motion by the dam break, then (b) forms a tail from the backwall to the nose. We focus on the question of if and when the propagation speed $u_N$ of the nose (front) of the GC is constant during this process (there is consensus that a significant deceleration of $u_N$ appears in the post-slumping stage.) The SW theory predicts correctly the adjustment of the flow field during the slumping stage, but indicates that a constant $u_N$ appears only for the classical case ($\gamma =E=c_D=\beta =0$) where $\gamma, E, c_D, \beta$ are the slope, entrainment and drag coefficients, and the scaled particle settling speed for a particle-driven GC. However, since $\gamma, E, c_D, \beta$ are typically small, the change of $u_N$ during the slumping phase is also small in many cases of interest. The interaction between the various driving and hindering mechanisms is elucidated. We show that, in a system with a horizontal (open) top (typical laboratory experiments), the height of the ambient increases along the slope, and this compensates for buoyancy loss due to particle sedimentation. We point out the need for further experimental and simulation studies for a better understanding of the slumping phase and transition to the next phases, and further assessment/improvement of the SW predictions.

This paper investigates the amplification and propagation of swirl fluctuations in turbulent swirling flows using resolvent analysis. Swirl fluctuations have been repeatedly observed in acoustically excited swirl flows and play a significant role in triggering thermoacoustic instabilities in swirl-stabilized flames. While recent research on simplified rotating laminar base flows suggests that the linear inertial-wave mechanism is a key driver of swirl fluctuations, it remains unclear whether this applies to the fully turbulent regime and whether a linear method is sufficient for modelling. To address this issue, a turbulent swirling pipe flow is considered using large-eddy simulations and phase-locked particle image velocimetry, which are combined with mean-field resolvent analysis. A sound agreement between the empirical and physics-based modes is found in terms of shape and propagation velocity. The latter is particularly important for thermoacoustic time-lag models. The comparison with a generic rotating pipe flow shows that the observed swirl fluctuations are indeed driven by a linear inertial wave mechanism. The resolvent framework is, then, exploited to further investigate the coupling and amplification mechanisms in detail. It is discovered that the combined effects of inertia and strong shear lead to very high amplification rates of the swirl fluctuations, explaining the high potential of these structures to trigger combustion instabilities. The study further demonstrates the capability of the resolvent to reveal the driving mechanisms of flow response structures in highly complex turbulent flows, and it opens the path for efficient physics-based optimization to prevent combustion instabilities.

This paper presents a numerical study on the flow around two tandem circular cylinders beneath a free surface at a Reynolds number of $180$. The free-surface effects on the wake dynamics and hydrodynamic forces are investigated through a parametric study, covering a parameter space of gap ratios from $0.20$ to $2.00$, spacing ratios from $1.50$ to $4.00$ and Froude numbers from $0.2$ to $0.8$. A jet-like flow accompanied by a shear layer of positive vorticity separating from the free surface is formed in the wake at small gap ratios, which significantly alters the wake pattern through its dynamic behaviours. At shallow submergence depths, the three-dimensional wake transitions from mode B to mode A as the distance between the cylinders increases. As submergence depth increases, the wavy deformation of the primary vortex cores disappears in the wake, and the flow transitions to a two-dimensional state. Higher Froude numbers can extend the effect of the free surface to deeper submergence depths. The critical spacing ratio tends to be larger at higher Froude numbers. Furthermore, the free-surface deformation is examined. The free-surface profile typically comprises a hydraulic jump immediately ahead of the upstream cylinder, trapped waves in the vicinity of the two tandem cylinders and well-defined travelling waves on the downstream side. The frequencies of the waves cluster around the vortex shedding frequency, indicating a close association between the generation of waves and the vortex shedding process.

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative $\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$ and residual $\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$ components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in $\boldsymbol {\varPhi }_{I}(\boldsymbol {x},t)$, whereas $\boldsymbol {\varPhi }_{R}(\boldsymbol {x},t)$ offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

We analyse the motion of a flagellated bacterium in a two-fluid medium using slender body theory. The two-fluid model is useful for describing a body moving through a complex fluid with a microstructure whose length scale is comparable to the characteristic scale of the body. This is true for bacterial motion in biological fluids (entangled polymer solutions), where the entanglement results in a porous microstructure with typical pore diameters comparable to or larger than the flagellar bundle diameter, but smaller than the diameter of the bacterial head. Thus, the polymer and solvent satisfy different boundary conditions on the flagellar bundle and move with different velocities close to it. This gives rise to a screening length $L_B$ within which the fluids exchange momentum and the relative velocity between the two fluids decays. In this work, both the solvent and polymer of the two-fluid medium are modelled as Newtonian fluids with different viscosities $\mu _s$ and $\mu _p$ (viscosity ratio $\lambda = \mu _p/\mu _s$), thereby capturing the effects solely introduced by the microstructure of the complex fluid. From our calculations, we observe an increased drag anisotropy for a rigid, slender flagellar bundle moving through this two-fluid medium, resulting in an enhanced swimming velocity of the organism. The results are sensitive to the interaction between the bundle and the polymer, and we discuss two physical scenarios corresponding to two types of interaction. Our model provides an explanation for the experimentally observed enhancement of swimming velocity of bacteria in entangled polymer solutions and motivates further experimental investigations.

Metastructures composed of a closely spaced plate array have been widely used in bespoke manipulation of waves in contexts of acoustics, electromagnetics, elasticity and water waves. This paper focuses on wave scattering by discrete plate array metastructures of arbitrary cross-sections, including isolated vertical metacylinders, periodic arrays and horizontal surface-piercing metacylinders. A suitable transform-based method has been applied to each problem to reduce the influence of barriers in a two-dimensional problem to a set of points in a one-dimensional wave equation wherein the solution is constructed using a corresponding Green's function. A key difference from the existing work is the use of an exact description of the plate array rather than an effective medium approximation, enabling the exploration of wave frequencies above resonance where homogenisation models fail but where the most intriguing physical findings are unravelled. The new findings are particularly notable for graded plate array metastructures that produce a dense spectrum of resonant frequencies, leading to broadband ‘rainbow reflection’ effects. This study provides new ideas for the design of structures for the bespoke control of waves with the potential for innovative solutions to coastal protection schemes or wave energy converters.

The combined effects of heater position, mean flow parameters and flame models on thermoacoustic instability in a one-dimensional Rijke tube are studied systematically by classic linear stability analysis (LSA) and lattice Boltzmann method (LBM) simulation. In the former, the stability range of the linear flame model under low Mach number assumption is solved analytically, while in the more general case, it is obtained by numerically solving the dispersion relation. Both the linear and nonlinear flame model cases are studied using the LBM with a spectral multiple-relaxation-time collision model and a newly developed heat source term. With the linear flame model, the LBM is in good agreement with LSA in predicting the transition point and growth rates, while with the nonlinear flame model, LBM simulations are consistent with solutions of limit cycle theory in the fully developed state. These results demonstrate the applicability of the LBM in solving complex thermoacoustic problems.

A gas bubble sitting at a liquid–gas interface can burst following the rupture of the thin liquid film separating it from the ambient, owing to the large surface energy of the resultant cavity. This bursting bubble forms capillary waves, a Worthington jet and subsequent droplets for a Newtonian liquid medium. However, rheological properties of the liquid medium like elastoviscoplasticity can greatly affect these dynamics. Using direct numerical simulations, this study exemplifies how the complex interplay between elasticity (in terms of elastic stress relaxation) and yield stress influences the transient interfacial phenomenon of bursting bubbles. We investigate how bursting dynamics depends on capillary, elastic and yield stresses by exploring the parameter space of the Deborah number ${{\textit {De}}}$ (dimensionless relaxation time of elastic stresses) and the plastocapillary number $\mathcal {J}$ (dimensionless yield-stress of the medium), delineating four distinct characteristic behaviours. Overall, we observe a non-monotonic effect of elastic stress relaxation on the jet development while plasticity of the elastoviscoplastic (EVP) medium is shown to affect primarily the jet evolution only at faster relaxation times (low ${{\textit {De}}}$). The role of elastic stresses on jet development is elucidated with the support of energy budgets identifying different modes of energy transfer within the EVP medium. The effects of elasticity on the initial progression of capillary waves and droplet formation are also studied. In passing, we study the effects of solvent–polymer viscosity ratio on bursting dynamics and show that polymer viscosity can increase the jet thickness apart from reducing the maximum height of the jet.

This study is dedicated to achieving efficient active noise control in a supersonic underexpanded planar jet, utilizing control parameters informed by resolvent analysis. The baseline supersonic underexpanded jet exhibits complex wave structures and substantial high-amplitude noise radiations. To perform the active control, unsteady blowing and suction are applied along the nozzle inner wall close to the exit. Employing both standard and acoustic resolvent analyses, a suitable frequency and spanwise wavenumber range for the blowing and suction is identified. Within this range, the control forcing can be significantly amplified in the near field, effectively altering the original sound-producing energetic structure while minimizing far-field amplification to prevent excessive noise. A series of large-eddy simulations are further conducted to validate the control efficiency, demonstrating an over 10 dB reduction in upstream-propagated screech noise. It is identified that the present unsteady control proves more effective than steady control at the same momentum coefficient. The controlled jet flow indicates that the shock structures become more stable, and the stronger the streamwise amplification of the forcing, the more likely it is to modify the mean flow characteristics, which is beneficial for reducing far-field noise radiation. Spectral proper orthogonal decomposition analysis of the controlled flow confirms that the control redistributes energy to higher forcing frequencies and suppresses large-scale antisymmetric and symmetric modes related to screech and its harmonics. The findings of this study highlight the potential of resolvent-guided control techniques in reducing noise in supersonic underexpanded jets and provide a detailed understanding of the inherent mechanisms for effective noise reduction through active control strategies.

A model of imbibition dynamics in a channel of flattened triangular cross-section is presented, taking into account the liquid film flow in the corners of the channel. The quasi-analytical solutions are derived on the basis of a lubrication approximation. The analysis encompasses two imbibition scenarios corresponding to a constant flow rate or constant pressure imposed in the wetting fluid at the inlet of the channel. In the former case, the process starts with a liquid film flow regime in the corners that is followed by a bulk and corner film flow regime characterised by a triple point advancing (far) ahead of the bulk meniscus after its entrance in the channel. In the latter case, the occurrence of the bulk and corner film flow regime is conditioned by an imposed pressure yielding a capillary pressure at the inlet smaller than a threshold capillary pressure. Above this threshold, the liquid film regime remains. For both imbibition scenarios under concern, important features are highlighted, including (i) the time scalings of the dynamics of both the triple point and apex of the bulk meniscus (when it exists), (ii) the contrast in the positions of these two points showing that the classical Washburn approach, which neglects the effect of the corner films, overpredicts the dynamics of the bulk meniscus. The important consequence is an early wetting fluid breakthrough at the channel outlet much before the bulk meniscus arrival. Comparisons with experimental data available in the literature are provided, validating the approach proposed in this work.

Understanding interfacial instability in a coflow system has relevance in the effective manipulation of small objects in microfluidic applications. We experimentally elucidate interfacial instability in stratified coflow systems of Newtonian and viscoelastic fluid streams in microfluidic confinements. By performing a linear stability analysis, we derive equations that describe the complex wave speed and the dispersion relationship between wavenumber and angular frequency, thus categorizing the behaviour of the systems into two main regimes: stable (with a flat interface) and unstable (with either a wavy interface or droplet formation). We characterize the regimes in terms of the capillary numbers of the phases in a comprehensive regime plot. We decipher the dependence of interfacial instability on fluidic parameters by decoupling the physics into viscous and elastic components. Remarkably, our findings reveal that elastic stratification can both stabilize and destabilize the flow, depending on the fluid and flow parameters. We also examine droplet formation, which is important for microfluidic applications. Our findings suggest that adjusting the viscous and elastic properties of the fluids can control the transition between wavy and droplet-forming unstable regimes. Our investigation uncovers the physics behind the instability involved in interfacial flows of Newtonian and viscoelastic fluids in general, and the unexplored behaviour of interfacial waves in stratified liquid systems. The present study can lead to a better understanding of the manipulation of small objects and production of droplets in microfluidic coflow systems.

The self-generated magnetic field in three-dimensional (3-D) single-mode ablative Rayleigh–Taylor instability (ARTI) relevant to the acceleration phase of a direct-drive inertial confinement fusion (ICF) implosion is investigated. It is found that stronger magnetic fields up to a few thousand teslas can be generated by 3-D ARTI rather than by its two-dimensional (2-D) counterpart. The Nernst effects significantly alter the magnetic field convection and amplify the magnetic fields. The magnetic field of thousands of teslas yields the Hall parameter of the order of unity, leading to profound magnetized heat flux modification. While the magnetic field significantly accelerates the bubble growth in the short-wavelength 2-D modes through modifying the heat fluxes, the magnetic field mostly accelerates the spike growth but has little influence on the bubble growth in 3-D ARTI. The accelerated growth of spikes in 3-D ARTI is expected to enhance material mixing and degrade ICF implosion performance. This work is focused on a regime relevant to direct-drive ICF parameters at the National Ignition Facility, and it also covers a range of key parameters that are relevant to other ICF designs and hydrodynamic/astrophysical scenarios.

In this study we investigate the sedimentation of prolate spheroids in a quiescent fluid by means of the particle-resolved direct numerical simulation. With the increase of the particle volume fraction $\phi$ from $0.1\,\%$ to $10\,\%$, we observe a non-monotonic variation of the mean settling velocity of particles, $\langle V_s \rangle$. By virtue of the Voronoi analysis, we find that the degree of particle clustering is highest when $\langle V_s \rangle$ reaches the local maximum at $\phi =1\,\%$. Under the swarm effect, clustered particles are found to preferentially sample downward fluid flows in the wake regions, leading to the enhancement of the settling speed. As for lower or higher volume fractions, the tendency of particle clustering and the preferential sampling of downward flows are attenuated. The hindrance effect becomes predominant when the volume fraction exceeds 5 % and reduces $\langle V_s \rangle$ to less than the isolated settling velocity. Particle orientation plays a minor role in the mean settling velocity, although individual prolate particles still tend to settle faster in suspensions when they deviate more from the broad-side-on alignment. Moreover, we also demonstrate that particles are prone to form column-like microstructures in dilute suspensions under the effect of wake-induced hydrodynamic attractions. The radial distribution function is higher at a lower volume fraction. As a result, the collision rate scaled by the particle number density decreases with the increasing volume fraction. By contrast, as another contribution to the particle collision rate, the relative radial velocity for nearby particles shows a minor degree of variation due to the lubrication effect.

We study the generation, nonlinear development and secondary instability of unsteady Görtler vortices and streaks in compressible boundary layers exposed to free-stream vortical disturbances and evolving over concave, flat and convex walls. The formation and evolution of the disturbances are governed by the compressible nonlinear boundary-region equations, supplemented by initial and boundary conditions that characterise the impact of the free-stream disturbances on the boundary layer. Computations are performed for parameters typical of flows over high-pressure turbine blades, where the Görtler number, a measure of the curvature effects, and the disturbance Reynolds number, a measure of the nonlinear effects, are order-one quantities. At moderate intensities of the free-stream disturbances, increasing the Görtler number renders the boundary layer more unstable, while increasing the Mach number or the frequency stabilises the flow. As the free-stream disturbances become more intense, vortices over concave surfaces no longer develop into the characteristic mushroom-shaped structures, while the flow over convex surfaces is destabilised. An occurrence map identifies Görtler vortices or streaks for different levels of free-stream disturbances and Görtler numbers. Our calculations capture well the experimental measurements of the enhanced skin friction and wall-heat transfer over turbine-blade pressure surfaces. The time-averaged wall-heat transfer modulations, termed hot fingers, are elongated in the streamwise direction and their spanwise wavelength is half of the characteristic wavelength of the free-stream disturbances. Nonlinearly saturated disturbances are unstable to secondary high-frequency modes, whose growth rate increases with the Görtler number. A new varicose even mode is reported, which may promote transition to turbulence at the stem of nonlinear streaks.

The coherent vortical structures in turbulent flow through a strong 16 : 1 three-dimensional contraction are studied using time-resolved volumetric measurements. Visualization using the vorticity magnitude criterion shows the emergence of long, stretched cylindrical vortices aligned with the mean flow. This alignment is quantified by probability density functions (p.d.f.s) of the direction cosines. We propose two measures to quantify the alignment, the peak height in the probability and a coefficient from the moment of the p.d.f., both of which reaffirm the strong streamwise alignment. The root mean square streamwise vorticity grows within the contraction to become 4.5 times larger than the transverse component, at the downstream location where the contraction ratio is $C=11$. The characteristic vortices become as long as the measurement volume, or more than 4 times the integral scale at the entrance to the contraction. We also characterize the vorticity enhancement along individual vortices, measuring 65 % strengthening over the distance where $C$ goes from 4 to 11. The prevalence of these coherent structures is estimated from 700 000 measured volumes, showing that near the outlet, it is more likely to have one or two of these structures present than none.