We study the dispersion of bubble swarms rising in initially quiescent water using three-dimensional Lagrangian tracking of deformable bubbles and tracer particles in an octagonal bubble column. Two different bubble sizes (3.5 mm and 4.4 mm) and moderate gas volume fractions ($0.52\,\%{-}1.20\,\%$) are considered. First, we compare the dispersion inside bubble swarms with that for single-bubble cases, and find that the horizontal mean squared displacement (MSD) in the swarm cases exhibits oscillations around the asymptotic scaling predicted for a diffusive regime. This occurs due to wake-induced bubble motion; however, the oscillatory behaviour is heavily damped compared to the single-bubble cases due to the presence of bubble-induced turbulence (BIT) and bubble–bubble interactions in the swarm. The vertical MSD in bubble swarms is nearly an order of magnitude faster than in the single-bubble cases, due to the much higher vertical fluctuating bubble velocities in the swarms. We also investigate tracer dispersion in BIT, and find that concerning the time to transition away from the ballistic regime, larger bubbles with a higher gas void fraction transition earlier than tracers, consistent with Mathai et al. (2018, Phys. Rev. Lett., vol. 121, 054501). However, for bubble swarms with smaller bubbles and a lower gas void fraction, they transition at the same time. This differing behaviour is due to the turbulence being more well-mixed for the larger bubble case, whereas for the smaller bubble case, the tracer dispersion is highly dependent on the wake fluctuations generated by the oscillating motion of nearby bubbles.

Studying rotating convection under geo- and astrophysically relevant conditions has proven to be extremely difficult. For the rotating Rayleigh–Bénard system, van Kan et al. (J. Fluid Mech., vol. 1010, 2025,A42)have now been able to massively extend the parameter space accessible by direct numerical simulations. Their progress relies on a rescaling of the governing Boussinesq equations, which vastly improves numerical conditioning (Julien et al., arXiv:2410.02702). This opens the door for investigating previously inaccessible dynamical regimes and bridges the gap to the asymptotic branch of rapidly rotating convection.

Previous studies claimed that the non-monotonic effects of wettability came mainly from the heterogeneity of geometries or flow conditions on multiphase displacements in porous media. For macroscopic homogeneous porous media, without permeability contrast or obvious preferential flow pathways, most pore-scale evidence showed a monotonic trend of the wettability effect. However, this work reports transitions from monotonic to non-monotonic wettability effects when the dimension of the model system rises from two-dimensional (2-D) to three-dimensional (3-D), validated by both the network modelling and the microfluidic experiments. The mechanisms linking the pore-scale events to macroscopic displacement patterns have been analysed through direct simulations. For 2-D porous media, the monotonic effect of wettability comes from the consistent transition pattern for the full range of capillary numbers $Ca$, where the capillary fingering mode transitions to the compact displacement mode as the contact angle $\theta$ decreases. Yet, it is indicated that the 3-D porous geometries, even though homogeneous without permeability contrast or obvious preferential flow pathways, introduce a different $Ca$–$\theta$ phase diagram with new pore-scale events, such as the coupling of capillary fingering with snap-off during strong drainage, and frequent snap-off events during strong imbibition. These events depend strongly on geometric confinements and capillary numbers, leading to the non-monotonicity of wettability effects. Our findings provide new insights into the multiphase displacement dependent on wettability in various natural porous media and offer design principles for engineering artificial porous media to achieve desired immiscible displacement behaviours.

Ultra-thin liquid sheets generated by impinging two liquid jets are crucial high-repetition-rate targets for laser ion acceleration and ultra-fast physics, and serve widely as barrier-free samples for structural biochemistry. The impact of liquid viscosity on sheet thickness should be comprehended fully to exploit its potential. Here, we demonstrate experimentally that viscosity significantly influences thickness distribution, while surface tension primarily governs shape. We propose a thickness model based on momentum exchange and mass transport within the radial flow, which agrees well with the experiments. These results provide deeper insights into the behaviour of liquid sheets and enable accurate thickness control for various applications, including atomization nozzles and laser-driven particle sources.

In this study, we obtain the continuum equations of Arctic sea ice motion starting from the dynamics of a single floe and show that the rheology that emerges from floe–floe interactions is viscous – as conjectured by Reed and Campbell (J. Geophys. Res., vol. 67 (1), 1962, pp. 281–297). The motion of the floe is principally driven by the wind and ocean currents and by inelastic collisions with the neighbouring floes. A mean-field representation of these collisions is developed, neglecting any changes in the floe thickness due to thermal growth and mechanical deformation. This mean-field representation depends on the state of the ice cover, and is expressed in terms of ice concentration and mean thickness. The resulting Langevin equation for the floe velocity, or the corresponding kinetic equation (Kramers–Chandrasekhar equation (KCE)) for its probability density, provides a complete description of the floe’s motion. We then use the floe-scale dynamics to obtain a continuum description of sea ice motion through a Chapman–Enskog analysis of the KCE. The local equilibrium solution to the kinetic equation is found to be the Laplace distribution, in qualitative agreement with observations. Our approach also allows us to establish the dependence of pressure and shear viscosity of the ice cover on ice concentration and mean thickness. Lastly, we show that our results resolve a conflict associated with the choice of the value of shear viscosity in previous idealised numerical studies of Arctic sea ice motion.

Stochastic resonance (SR) is universal phenomenon, where noise amplifies a weak periodic signal in bistable nonlinear systems, with wide applications in biology, climate science, engineering etc., although in fluid dynamics it remains underexplored. Recently, we unexpectedly found SR above non-modal elastic instability onset in an inertialess viscoelastic channel flow, where it emerges on the top of a chaotic streamwise velocity power spectrum $E_u$ due to its interaction with white-noise spanwise velocity power spectrum $E_w$ and weak elastic waves. These three conditions necessary for SR emergence differ from those required for the classical SR emergence mentioned above. Here, we consider SR in an inertialess viscoelastic channel flow with a smoothed inlet causing order of magnitude lower noise intensity than in our former studies. Our observations reveal that SR appears at the same conditions mentioned above, where SR is found just upon the instability onset in a lower subrange of a transition regime, in contrast, here, SR persists across all flow regimes – transition, elastic turbulence and drag reduction. Furthermore, we provide experimental evidence that SR, presented by a sharp peak in $E_u$, manifests as either a standing or propagating wave in the $x$-direction, with a rather uniform amplitude of streamwise velocity fluctuations and zero propagation velocity in the $z$-direction. These findings reveal a new mechanism underpinning the transition to a chaotic channel flow of viscoelastic fluids and establish SR as a robust framework for understanding complex flow dynamics. This work opens new avenues for exploring SR in other nonlinear systems and practical applications such as mixing enhancement and flow control in industrial and biological contexts.

The accurate quantification of wall-shear stress dynamics is of substantial importance for various applications in fundamental and applied research, spanning areas from human health to aircraft design and optimization. Despite significant progress in experimental measurement techniques and postprocessing algorithms, temporally resolved wall-shear stress fields with adequate spatial resolution and within a suitable spatial domain remain an elusive goal. Furthermore, there is a systematic lack of universal models that can accurately replicate the instantaneous wall-shear stress dynamics in numerical simulations of multiscale systems where direct numerical simulations (DNSs) are prohibitively expensive. To address these gaps, we introduce a deep learning architecture that ingests wall-parallel streamwise velocity fields at $y^+ \approx 3.9 \sqrt {Re_\tau }$ of turbulent wall-bounded flows and outputs the corresponding two-dimensional streamwise wall-shear stress fields with identical spatial resolution and domain size. From a physical perspective, our framework acts as a surrogate model encapsulating the various mechanisms through which highly energetic outer-layer flow structures influence the governing wall-shear stress dynamics. The network is trained in a supervised fashion on a unified dataset comprising DNSs of statistically one-dimensional turbulent channel and spatially developing turbulent boundary layer flows at friction Reynolds numbers ranging from $390$ to $1500$. We demonstrate a zero-shot applicability to experimental velocity fields obtained from particle image velocimetry measurements and verify the physical accuracy of the wall-shear stress estimates with synchronized wall-shear stress measurements using the micro-pillar shear-stress sensor for Reynolds numbers up to $2000$. In summary, the presented framework lays the groundwork for extracting inaccessible experimental wall-shear stress information from readily available velocity measurements and thus, facilitates advancements in a variety of experimental applications.

Experiments are carried out in a smooth-wall turbulent boundary layer (TBL) ($\textit{Re}_\tau \geq 3500$) subjected to different pressure gradient (PG) histories. Oil-film interferometry is used to measure the skin friction evolution over the entire history while wide-field particle image velocimetry captures the mean flow field. This data are used to demonstrate the influence of PG history on skin friction as well as other integral quantities such as displacement ($\delta ^*$) and momentum thickness ($\theta$). Based on observations from the data, a new set of ordinary differential equations are proposed to model the streamwise evolution of a TBL subjected to different PG histories. The model is calibrated using a limited number of experimental cases and its utility is demonstrated on other cases. Moreover, the model is applied to data from large-eddy simulations of flows in adverse PG conditions (Bobke et al. 2017, J. Fluid Mech. 820, 667–692). The model is subsequently used to identify the impact of PG history length on the boundary layer. This can also be interpreted as determining the spatial frequency response of the boundary layer to PG disturbances. Results suggest that short spatial variations in PGs primarily affect a small portion of the TBL evolution, whereas longer-lasting ones have a more extensive impact.

This study utilises large-eddy simulation with the actuator line model to examine the effects of the tip speed ratio (TSR) on the wake-meandering characteristics of a wind turbine in uniform and turbulent inflows. It is shown that as the TSR grows, the onset position of the wake meandering moves closer to the rotor, and the magnitude of wake oscillation is stronger. This aligns with previous work showing that a higher TSR can accelerate the instability and breakdown of tip vortices. Without a nacelle, the Strouhal number of the wake meandering is found to be independent of the TSR under both the uniform and turbulent inflows. However, with a relatively large nacelle, the Strouhal number first increases and then decreases with TSR. Therefore, the current discovery elucidates the crucial role of the nacelle and clarifies the origin of the TSR dependence of the Strouhal number in wake meandering. In addition, the characteristic frequency of the wake meandering under the turbulent inflow is much smaller than that under the uniform inflow, because of the significant influence of the freestream turbulence. Furthermore, the proper orthogonal decomposition (POD) and spectral POD (SPOD) methods are employed to study the spatiotemporal characteristics of the meandering wake and its TSR dependence. It is found that the tip and root vortices are the prominent wake structures under the uniform inflow, whereas more complex multiscale structures from the interaction between the freestream turbulence and tip/root vortices exist under the turbulent inflow. Moreover, an amplitude modulation phenomenon of the POD time coefficients at the optimal TSR is observed in the uniform inflow case. Finally, a reduced-order model is constructed for predicting the wake dynamics by combining the SPOD and the ‘sparse identification of nonlinear dynamics’ algorithm with high accuracy and interpretability.

We present a theoretical framework and validation for manipulating instability growth in shock-accelerated dual-layer material systems, which feature a light–heavy interface followed by two sequential heavy–light interfaces. An analytical model is first developed to predict perturbation evolution at the two heavy–light interfaces, explicitly incorporating the effects of reverberating waves within the dual-layer structure. The model identifies five distinct control regimes for instability modulation. Shock-tube experiments and numerical simulations are designed to validate these regimes, successfully realising all five predicted states. Notably, the selective growth stagnation of a perturbation at either the upstream or downstream heavy–light interface is realised numerically by tuning the initial separation distances of the three interfaces. This work elucidates the critical role of the wave dynamics in governing interface evolution of a shocked dual layer, offering insights for mitigating hydrodynamic instabilities in practical scenarios such as inertial confinement fusion.

We consider two-dimensional (2-D) free surface gravity waves in prismatic channels, including bathymetric variations uniquely in the transverse direction. Starting from the Saint-Venant equations (shallow-water equations) we derive a one-dimensional transverse averaged model describing dispersive effects related solely to variations of the channel topography. These effects have been demonstrated in Chassagne et al. 2019 J. Fluid Mech. 870, 595–616 to be predominant in the propagation of bores with Froude numbers below a critical value of approximately 1.15. The model proposed is fully nonlinear, Galilean invariant, and admits a variational formulation under natural assumptions about the channel geometry. It is endowed with an exact energy conservation law, and admits exact travelling-wave solutions. Our model generalises and improves the linear equations proposed by Chassagne et al. 2019 J. Fluid Mech. 870, 595–616, as well as in Quezada de Luna and Ketcheson, 2021 J. Fluid Mech. 917, A45. The system is recast in two useful forms appropriate for its numerical approximations, whose properties are discussed. Numerical results allow the verification of the implementation of these formulations against analytical solutions, and validation of our model against fully 2-D nonlinear shallow-water simulations, as well as the famous experiments by Treske 1994 J. Hyd. Res. 32, 355–370.

Nonlinear dynamical systems often allow for multiple statistically stationary states for the same values of the control parameters. Herein, we introduce a framework that selectively eliminates specific nonlinear triad interactions, thereby suppressing emergence of a particular state, and enabling the emergence of another. The methodology is applied to yield the multiple convection-roll states in two-dimensional planar Rayleigh–Bénard convection (e.g. Wang et al., 2020, Phys. Rev. Lett., vol. 125, 074501) in the turbulent regime. The intrusive framework presented here is based on the observation that the characteristic wavenumber associated with the mean horizontal size of the convection rolls mediates triadic scale interactions resulting in both kinetic energy and temperature variance cascades that are dominant energy transfer processes in a statistically stationary state. Suppression of these cascades mediated by a candidate wavenumber hinders the formation of the convection rolls at that wavenumber. Consequently, convection rolls are formed at another candidate wavenumber which is allowed to mediate energy to establish the cascade processes. In case no stable convection-roll states are possible, this technique does not tend to yield any convection rolls, making it a suitable method for discovering multiple states. Whereas in previous investigations the signature of different states in the initial condition in simulations yielded the multiple states, the present method alleviates such dependence of the emergence of multiple states on initial conditions. It is also demonstrated that accurate predictions of statistical quantities, such as Nusselt number and volume-averaged Reynolds numbers, can also be obtained using this method. The convection-roll states yielded using this technique may be used as initial conditions for direct simulations quickly converging to the target roll state without taking long convergence routes involving state transitions. Additionally, because only one state can possibly emerge in each simulation, this technique can empirically designate each of the multiple states with respect to their stability.

Ventilated cavities in the wake of a two-dimensional bluff body are studied experimentally via time-resolved X-ray densitometry. With a systematic variation of flow velocity and gas injection rate, expressed as Froude number ($\textit{Fr}$) and ventilation coefficient ($C_{qs}$), four cavities with different closure types are identified. A regime map governed by $\textit{Fr}$ and $C_{qs}$ is constructed to estimate flow conditions associated with each cavity closure type. Each closure exhibits a different gas ejection mechanism, which in turn dictates the cavity geometry and the pressure in the cavity. Three-dimensional cavity closure is seen to exist for the supercavities at low $\textit{Fr}$. However, closure is nominally two-dimensional for supercavities at higher $\textit{Fr}$. At low $C_{qs}$, cavity closure is seen to be wake-dominated, while supercavities are seen to have interfacial perturbation near the closure at higher $C_{qs}$, irrespective of $\textit{Fr}$. With the measured gas fraction, a gas balance analysis is performed to quantify the gas ejection rate at the transitional cavity closure during its formation. For a range of $\textit{Fr}$, the transitional cavity closure is seen to be characterised by re-entrant flow, whose intensity depends on the flow inertia, dictating the gas ejection rates. Two different ventilation strategies were employed to systematically investigate the formation and maintenance gas fluxes. The interaction of wake and gas injection is suspected to dominate the cavity formation process and not the maintenance, resulting in ventilation hysteresis. Consequently, the ventilation gas flux required to maintain the supercavity is significantly less than the gas flux required to form the supercavity.

The jet from a model-scale, internally mixed nozzle produced a loud howling when operated at jet Mach numbers between 0.80 and 1.00. Discrete tones dominated the noise radiated to the far field and powerful oscillations were present in the jet. To explain these observations, this paper leverages a blend of experimental acoustic and flow measurements and modal analyses thereof via the spectral proper orthogonal decomposition, computational fluid dynamics simulations and local, linear stability analyses of vortex-sheet models for the flow inside the nozzle. This blend of experiments, computations and theory makes clear the cause of the howling, what sets its characteristic frequency and how it may be suppressed. The flow around a small-radius, convex bend just upstream of the final-nozzle exit led to a pocket of locally supersonic flow that was terminated by a shock. The shock was strong enough to separate the boundary layer, but neither the attached nor separated states were stable. A periodic, shock-induced separation of the boundary layer resulted, and this shock-wave/boundary-layer interaction coupled with a natural acoustic mode of the nozzle’s interior in a feedback phenomenon of sorts. Acoustic tones and large flow oscillations were produced at the associated natural frequency of the nozzle’s interior.

Cavitation bubble pulsation and liquid jet loads are the main causes of hydraulic machinery erosion. Methods to weaken the load influences have always been hot topics of related research. In this work, a method of attaching a viscous layer to a rigid wall is investigated in order to reduce cavitation pulsations and liquid jet loads, using both numerical simulations and experiments. A multiphase flow model incorporating viscous effects has been developed using the Eulerian finite element method (EFEM), and experimental methods of a laser-induced bubble near the viscous layer attached on a rigid wall have been carefully designed. The effects of the initial bubble–wall distance, the thickness of the viscous layer, and the viscosity on bubble pulsation, migration and wall pressure load are investigated. The results show that the bubble migration distance, the normalised thickness of the oil layer and the wall load generally decrease with the initial bubble–wall distance or the oil-layer parameters. Quantitative analysis reveals that when the initial bubble–wall distance remains unchanged, there exists a demarcation line for the comparison of the bubble period and the reference period (the bubble period without viscous layer under the same initial bubble–wall distance), and a logarithmic relationship is observed that $\delta \propto \log_{10} \mu ^*$, where $\delta =h/R_{max}$ is the thickness of the viscous layer h normalised by the maximum bubble radius $R_{max}$, $\mu ^* = \mu /({R_{max }}\sqrt {{\rho }{{\mathop {P}\nolimits } _{{atm}}}})$ is the dynamic viscosity $\mu$ normalised by water density $ \rho $ and atmospheric pressure $P_{atm}$. The results of this paper can provide technical support for related studies of hydraulic cavitation erosion.

This research investigates the spanwise oscillation patterns of turbulent non-premixed flames in a tandem configuration, using both experimental methods and large eddy simulations under cross-airflow conditions. Based on the heat release rate (17.43–34.86 kW) and the burner size (0.15 $\times$ 0.15 m), the flame behaves like both a buoyancy-controlled fire (such as a pool fire) and, due to cross-wind effects, a forced flow-controlled fire. The underlying fire dynamics was modelled by varying the spacing between the square diffusion burners, cross-wind velocity and heat release rate. Two flapping modes, the oscillating and bifurcating modes, were observed in the wake of the downstream diffusion flame. This behaviour depends on the wake of the upstream diffusion flame. As the backflow of the upstream flame moved downstream, the maximum flame width of the downstream flame became broader. The flapping amplitude decreased with a stronger cross-wind. Furthermore, the computational fluid dynamics simulation was performed by FireFOAM based on OpenFOAM v2006 2020 to investigate the flapping mechanism. The simulation captured both modes well. Disagreement of the flapping period on the left and right sides results in the oscillating mode, while an agreement of the flapping period results in the bifurcating mode. Finally, the scaling law expressed the dimensionless maximum flame width with the proposed set of basic dimensional parameters, following observations and interpretation by simulations. The results help prevent the potential hazards of this type of basic fire scenario and are fundamentally significant for studying wind-induced multiple fires.

The rupture of a liquid film, where a thin liquid layer between two other fluids breaks and forms holes, commonly occurs in both natural phenomena and industrial applications. The post-rupture dynamics, from initial hole formation to the complete collapse of the film, are crucial because they govern droplet formation, which plays a significant role in many applications such as disease transmission, aerosol formation, spray drying nanodrugs, oil spill remediation, inkjet printing and spray coating. While single-hole rupture has been extensively studied, the dynamics of multiple-hole ruptures, especially the interactions between neighbouring holes, are less well understood. Here, this study reveals that when two holes ‘meet’ on a curved film, the film evolves into a spinning twisted ribbon before breaking into droplets, distinctly different from what occurs on flat films. We explain the formation and evolution of the spinning twisted ribbon, including its geometry, orbits, corrugations and ligaments, and compare the experimental observations with models. We compare and contrast this phenomena with its counterpart on planar films. While our experiments are based on the multiple-hole ruptures in corona splash, the underlying principles are likely applicable to other systems. This study sheds light on understanding and controlling droplet formation in multiple-hole rupture, improving public health, climate science and various industrial applications.

Supersonic free jets are extensively employed across a range of applications, especially in high-tech industries such as semiconductor processing and aerospace propulsion. Due to the difficulties involved in flow measurement, previous research on supersonic free jets has primarily focused on investigating near-field shockwave structures, with quantitative experimental analysis of the far-field zone being relatively scarce. However, physical understanding of the far-field flow, particularly post-shockwave energy dissipation, holds significant importance for the application and utilisation of these jets in vacuum environments. Therefore, this study aims to provide a robust experimental foundation for a rarefied supersonic free jet through the analysis of the flow field in both the near- and far-field zones. Nanometre-sized tracer particles and molecules were utilised to measure the rarefied supersonic jet flow field using particle image velocimetry and acetone molecular tagging velocimetry, respectively. The experiments revealed that in rarefied conditions, the supersonic jet exhibits a one-barrel shockwave structure in the near field, and after passing the Mach disk, a long annular viscous layer develops downstream. Experimental data on the jet velocity profile and width demonstrated a transition to a laminar flow regime in the far-field zone. This transition aligns with the theoretically inferred flow regimes based on the complex Reynolds number. The velocity profile and potential core length of the laminar flow regime could be modelled using a bi-modal distribution, which represents the summation of symmetric Gaussian distributions.

We extend the perceived velocity gradient defined by a group of particles that was previously used to investigate the Lagrangian statistics of fluid turbulence to the study of inertial particle dynamics. Using data from direct numerical simulations, we observe the correlation between the strong compression in the particle phase and the instantaneous local fluid compression. Furthermore, the Lagrangian nature of the particle velocity gradient defined in this way allows an investigation of its evolution along particle trajectories, including the process after the caustic event, or the blow-up of the particle velocity gradient. Observations reveal that, for particles with Stokes number in the range $St \lesssim 1$, inertial particles experience the maximum compression by local fluid before the caustic event. Interestingly, data analyses show that, while the post-caustic process is mainly the relaxation of the particle motion and the particle relaxation time is the relevant time scale for the dynamics, the pre-caustic dynamics is controlled by the fluid–particle interaction and the proper time scale is determined by both the Kolmogorov time and the particle relaxation time.

Contactless manipulation of small objects is essential for biomedical and chemical applications, such as cell analysis, assisted fertilisation and precision chemistry. Established methods, including optical, acoustic and magnetic tweezers, are now complemented by flow control techniques that use flow-induced motion to enable precise and versatile manipulation. However, trapping multiple particles in fluid remains a challenge. This study introduces a novel control algorithm capable of steering multiple particles in flow. The system uses rotating disks to generate flow fields that transport particles to precise locations. Disk rotations are governed by a feedback control policy based on the optimising a discrete loss framework, which combines fluid dynamics equations with path objectives into a single loss function. Our experiments, conducted in both simulations and with the physical device, demonstrate the capability of the approach to transport two beads simultaneously to predefined locations, advancing robust contactless particle manipulation for biomedical applications.