The study presents observations on the interaction of double-blade propeller tip vortices with a smooth-wall turbulent boundary layer (TBL). The wall-bounded helicoidal vortices from the propeller modify the velocity profiles and turbulence statistics. The effects of two different tip clearances, $\epsilon = 0.1\delta _0$ and $0.5\delta _0$, at a matched thrust, are explored with particle image velocimetry to understand the dynamics of tip-vortex formation within the logarithmic and wake regions of the boundary layer. The measurements are performed with $\lambda =U_{tip}/U_{\infty }$ in the range 5.3–5.9, and a blade passing frequency ($\,f_{prop}$) of the same order of the boundary-layer time scale ($\,f_{TBL}$). Observations indicate a reduction in the extent of the log region and an enhancement of the wake parameter $\varPi$, mirroring the behaviour seen in TBLs under adverse pressure gradient conditions. Notably, the slipstream most contracted region exhibits a significant reduction in the skin friction coefficient $C_f$ and an amplification of the velocity fluctuation statistics across the entire boundary layer. At a clearance of $\epsilon = 0.1\delta _0$, there is evidence of the formation of paired coherent wall-bounded structures. The presence of the wall decreases the amplitude of both periodic and stochastic fluctuations obtained with a phase-locked triple decomposition. An exception is observed behind the propeller for the stochastic fluctuations of the wall-normal component of the flow, which become amplified as the blades move away from the wall. This leads to the creation of a more intense phase-locked two-point spatial coherence than that observed in fluctuations aligned with the streamwise direction. Furthermore, results reveal that reduced tip clearances lead to higher viscous dissipation and more active energy exchange between the mean flow and organized motions.

The motion of a sphere freely rising or falling in a 5d (d is the diameter of the sphere) square tube was numerically studied for the sphere-to-fluid density ratio ranging from 0.1 to 2.3 (0.1 ≤ ρs/ρ ≤ 2.3, ρs is the density of spheres and ρ the fluid density) and Galileo number from 140 to 230 (140 ≤ Ga ≤ 230). We report that Hopf bifurcation occurs at Gacrit ≈ 157, where both the heavy and light spheres lose stability. The helical motion is widely seen for all spheres at Ga > 160 resulting from a double-threaded vortex interacting with the tube walls, which becomes irregular at Ga ≥ 190 where heavy spheres act differently from their counterparts; that is, heavy spheres change their helical directions alternately while light spheres exhibit helical trajectories with jaggedness in connection with the shedding of the double-threaded vortices. This is because of the difference in inertia between the heavy and light spheres. We also checked the oscillation periods for the helical motion of the spheres. They show opposite variations with ρs/ρ for the two types of spheres. Light spheres (ρs/ρ ≤ 0.7) reach a zigzagging regime at Ga ≥ 200 where a vortex loop (hairpin-like vortical structure) is formed which may develop into a vortex ring downstream at small ρs/ρ. This might be the first time a transition from the helical motion to the zigzagging motion for heavy spheres (ρs/ρ ≥ 1.8) has been reported. Finally, we examined the dependence of both the terminal Reynolds number and the drag coefficient of the spheres on the Galileo number.

This study introduces vector autoregression (VAR) as a linear procedure that can be used for synthesizing turbulence time series over an entire plane, allowing them to be imposed as an efficient turbulent inflow condition in simulations requiring stationary and cross-correlated turbulence time series. VAR is a statistical tool for modelling and prediction of multivariate time series through capturing linear correlations between multiple time series. A Fourier-based proper orthogonal decomposition (POD) is performed on the two-dimensional (2-D) velocity slices from a precursor simulation of a turbulent boundary layer at a momentum thickness-based Reynolds number, $Re_{\theta }=790$. A subset of the most energetic structures in space are then extracted, followed by applying a VAR model to their complex time coefficients. It is observed that VAR models constructed using time coefficients of 5 and 30 most energetic POD modes per wavenumber (corresponding to $66\,\%$ and $97\,\%$ of turbulent kinetic energy, respectively) are able to make accurate predictions of the evolution of the velocity field at $Re_{\theta }=790$ for infinite time. Moreover, the 2-D velocity fields from the POD–VAR when used as a turbulent inflow condition, gave a short development distance when compared with other common inflow methods. Since the VAR model can produce an infinite number of velocity planes in time, this enables reaching statistical stationarity without having to run an extremely long precursor simulation or applying ad hoc methods such as periodic time series.

We explore the application of the reference map technique, originally developed for Eulerian simulation of solid mechanics, in Lagrangian kinematics of turbulent flows. Unlike traditional methods based on explicit particle tracking, the reference map facilitates the calculation of flow maps and gradients without the need for particles. This is achieved through an Eulerian update of the reference map, which records the take-off positions of fluid particles. This approach is found to be mathematically equivalent to the work of Leung (J. Comput. Phys., vol. 230, issue 9, 2011, pp. 3500–3524), who computed the flow map of simple two-dimensional flows using an Eulerian approach. We discuss important modifications necessary for its first application to complex three-dimensional turbulent flows, including the conservative, low-dissipation update of the flow map and the treatment of periodic boundary conditions. We first demonstrate the accuracy of finite-time Lyapunov exponent (FTLE) calculations based on the reference map against the standard particle-based approach in a two-dimensional Taylor–Green vortex. Then we apply it to turbulent channel flow at $Re_\tau =180$, where Lagrangian coherent structures identified as ridges of the backward-time FTLE are found to bound vortical regions of flow, consistent with Eulerian coherent structures from the $Q$-criterion. The reference map also proves suitable for material surface tracking despite not explicitly tracking particles. This capability can provide valuable insights into the Lagrangian landscape of turbulent momentum transport, complementing Eulerian velocity field analysis. The evolution of initially wall-normal material surfaces in the viscous sublayer, buffer layer and log layer sheds light on the Reynolds stress-generating events from a Lagrangian perspective.

Compound flows consist of two or more parallel compressible streams in a duct and their theoretical treatment has gained attention for the analysis and modelling of ejectors. Recent works have shown that these flows can experience choking upstream of the geometric throat. While it is well known that friction can push the sonic section downstream of the throat, no mechanism has been identified yet to explain its displacement in the opposite direction. This study extends the existing compound flow theory and proposes a one-dimensional (1-D) model, including friction between the streams and the duct walls. The model captures the upstream and downstream displacements of the sonic section. Through an analytical investigation of the singularity at the sonic section, it is demonstrated that friction between the streams is the primary driver of upstream displacement. The 1-D formulation is validated against axisymmetric Reynolds averaged Navier–Stokes simulations of a compound nozzle for various inlet pressures and geometries. The effect of friction is investigated using an inviscid simulation for the isentropic case and viscous simulations with both slip and no-slip conditions at the wall. The proposed extension accurately captures the displacement of the sonic section, offering a new tool for in-depth analysis and modelling of internal compound flows.

We investigate the effect of high wind speeds on the breakup mechanisms that govern the formation of a spray from nozzles that form liquid sheets, which subsequently break up. The fragmentation mechanism of liquid sheets from spray nozzles has recently been described in detail under quiescent conditions. With high wind speeds, measurements of the droplet size distribution reveal two rather than one characteristic drop sizes, suggesting the existence of two distinct breakup mechanisms. High-speed images of the spray are used to identify these two mechanisms. We show that the smaller droplets result from the breakup of ‘bags’ formed in the spray sheet by the wind, while the larger droplets result from the breakup of the remaining perforated sheet. Based on the two mechanisms, a probability density function is constructed and fitted to the measured droplet size distributions. We show that the spray sheet destabilises due to the Rayleigh–Taylor instability induced by the airflow, and that the experimentally observable breakup length and size of the holes blown in the sheet are predicted by the fastest growing wavenumber. From this, a theoretical prediction for the droplet size from bag breakup and remaining sheet breakup is derived.

The fingers known as bubbles (spikes) resulting from the penetration of light (heavy) fluids into heavy (light) fluids are significant large-scale features of Richtmyer–Meshkov instability (RMI). Through shock-tube experiments, we study finger collisions in light fluid layers under reshock conditions. Four unperturbed fluid layers with varying thicknesses are created to analyse the motion of waves and interfaces during finger collisions. The wave dynamics, sensitive to initial layer thicknesses, are characterized by a one-dimensional theory. Eight perturbed fluid layers, with four thicknesses and two interface phase combinations, are generated to explore the finger collision mechanism. It is shown that after reshock, the initial in-phase and anti-phase cases undergo spike–bubble rear-end collisions (SBCs) and spike–spike head-on collisions (SSCs), respectively. Compared with SBCs, SSCs significantly suppress spike growth, leading to the attenuation of perturbation growth, especially for larger thicknesses. As the initial thickness decreases, an SSC impedes the downstream interface from reversing its phase, resulting in abnormal RMI, thereby reducing the SSC's effectiveness in attenuating growth. The effects of rarefaction waves enhance both interfaces’ amplitudes and the whole layer's thickness, diminishing the intensity of finger collisions, while the second reshock exerts an opposing influence. Linear and nonlinear models, incorporating the influence of reshocks and rarefaction waves, are developed to predict the interface perturbation growth before and after finger collisions.

We study numerically the flow around a spherical droplet set fixed in a linear shear flow with moderate shear rates ($Sr\leq 0.5$, $Sr$ being the ratio between the velocity difference across the drop and the relative velocity) over a wide range of external Reynolds numbers ($0.1<{{Re}}\leq 250$, ${{Re}}$ based on the slip velocity and the viscosity of the external fluid) and drop-to-fluid viscosity ratios ($0.01\leq \mu ^\ast \leq 100$). The flow structure, the vorticity field and their intrinsic connection with the lift force are analysed. Specifically, the results on lift force are compared with the low-${{Re}}$ solution derived for droplets of arbitrary $\mu ^\ast$, as well as prior data at finite ${{Re}}$ available in both the clean-bubble limit ($\mu ^\ast \to 0$) and the solid-sphere limit ($\mu ^\ast \to \infty$). Notably, at ${{Re}}=O(100)$, the lift force exhibits a non-monotonic transition from $\mu ^\ast \to 0$ to $\mu ^\ast \to \infty$, peaking at $\mu ^\ast \approx 1$. This behaviour is related to an internal three-dimensional flow bifurcation also occurring under uniform-flow conditions, which makes the flow to evolve from axisymmetric to biplanar symmetric. This flow bifurcation occurs at low-but-finite $\mu ^\ast$ when the internal Reynolds number (${{Re}}^i$, based on the viscosity of the internal fluid) exceeds approximately 300. In the presence of shear, the corresponding imperfect bifurcation enhances the extensional rate of the flow in the wake. Consequently, the streamwise vortices generated behind the droplet can be more intense compared with those behind a clean bubble. Given the close relation between the lift and these vortices, a droplet with ${{Re}}=O(100)$ and $\mu ^\ast \approx 1$ typically experiences a greater lift force than that in the inviscid limit.

We theoretically and experimentally study gravity currents of a Newtonian fluid advancing in a two-dimensional, infinite and saturated porous domain over a horizontal impermeable bed. The driving force is due to the density difference between the denser flowing fluid and the lighter, immobile ambient fluid. The current is taken to be in the Darcy–Forchheimer regime, where a term quadratic in the seepage velocity accounts for inertial contributions to the resistance. The volume of fluid of the current varies as a function of time as $\sim T^{\gamma }$, where the exponent parameterizes the case of constant volume subject to dam break ($\gamma =0$), of constant ($\gamma =1$), waning ($\gamma <1$) and waxing inflow rate ($\gamma >1$). The nonlinear governing equations, developed within the lubrication theory, admit self-similar solutions for some combinations of the parameters involved and for two limiting conditions of low and high local Forchheimer number, a dimensionless quantity involving the local slope of the current profile. Another parameter $N$ expresses the relative importance of the nonlinear term in Darcy–Forchheimer's law; values of $N$ in practical applications may vary in a large interval around unity, e.g. $N\in [10^{-5},10^{2}]$; in our experiments, $N\in [2.8,64]$. Sixteen experiments with three different grain sizes of the porous medium and different inflow rates corroborate the theory: the experimental nose speed and current profiles are in good agreement with the theory. Moreover, the asymptotic behaviour of the self-similar solutions is in excellent agreement with the numerical results of the direct integration of the full problem, confirming the validity of a relatively simple one-dimensional model.

We present new constrained and free-swimming experiments and simulations in the inertial regime, with Reynolds number $\mbox{Re} = O(10^4)$, of a pair of two-dimensional and three-dimensional pitching hydrofoils interacting in a minimal school. The hydrofoils have an out-of-phase synchronisation, and they are varied through in-line, staggered and side-by-side formations within the two-dimensional interaction plane. It is discovered that there is a two-dimensionally stable equilibrium point for a side-by-side formation. This formation is super-stable, meaning that hydrodynamic forces will passively maintain this formation even under external perturbations, and the school as a whole has no net forces acting on it that cause it to drift to one side or the other. Previously discovered one-dimensionally stable equilibria driven by wake vortex interactions are shown to be, in fact, two-dimensionally unstable, at least for an out-of-phase synchronisation. Additionally, it is discovered that a trailing-edge vortex mechanism provides the restorative force to stabilise a side-by-side formation. The stable equilibrium is further verified by experiments and simulations for freely swimming foils where dynamic recoil motions are present. When constrained, swimmers in compact side-by-side formations experience collective efficiency and thrust increases up to 40 % and 100 %, respectively, whereas slightly staggered formations output an even higher efficiency improvement of 84 %, with an 87 % increase in thrust. Freely swimming foils in a stable side-by-side formation show efficiency and speed enhancements of up to 9 % and 15 %, respectively. These newfound schooling performance and stability characteristics suggest that fluid-mediated equilibria may play a role in the control strategies of schooling fish and fish-inspired robots.

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.

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.

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

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.

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.

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.