This study explores precession-driven flows in a non-axisymmetric ellipsoid spinning around its medium axis. Using an experimental approach, we focus on two aspects of the flow: the base flow of uniform vorticity and the development of fluid instabilities. In contrast to a preceding paper (J. Fluid. Mech., vol. 932, 2022, A24), where the ellipsoid rotated around its shortest axis, we do not observe bi-stability or hysteresis of the base flow, but a continuous transition from small to large differential rotation and tilt of the fluid rotation axis. We then use the model developed by Noir & Cébron (J. Fluid. Mech., vol. 737, 2013, pp. 412–439) to numerically determine regions in the parameter space of axial and equatorial deformations for which bi-stability may exist. Concerning fluid instabilities, we use three independent observations to track the onset of both boundary layer and parametric instabilities. Our results clearly show the presence of a parametric instability, yet the exact nature of the underlying mechanism (conical shear layer instability, shear instability and elliptical instability) is not unambiguously identified. A coexisting boundary layer instability, although unlikely, cannot be ruled out based on our experimental data. To make further progress on this topic, a new generation of experiments at significantly lower Ekman numbers (ratio of rotation and viscous time scales) is clearly needed.

In confined systems, the entrapment of a gas volume with an equivalent spherical diameter greater than the dimension of the channel can form extended bubbles that obstruct fluid circuits and compromise performance. Notably, in sealed vertical tubes, buoyant long bubbles cannot rise if the inner tube radius is below a critical value near the capillary length. This critical threshold for steady ascent is determined by geometric constraints related to the matching of the upper cap shape with the lubricating film surrounding the elongated part of the bubble. Developing strategies to overcome this threshold and release stuck bubbles is essential for applications involving narrow liquid channels. Effective strategies involve modifying the matching conditions with an external force field to facilitate bubble ascent. However, it is unclear how changes in acceleration conditions affect the motion onset of buoyancy-driven long bubbles. This study investigates the mobility of elongated bubbles in sealed tubes with an inner radius near the critical value inhibiting bubble motion in a vertical setting. Two strategies are explored to tune bubble motion, leveraging variations in axial and transversal accelerations: tube rotation around its axis and tube inclination relative to gravity. By revising the geometrical constraints of the simple vertical setting, the study predicts new thresholds based on rotational speed and tilt angle, respectively, providing forecasts for the bubble rising velocity under modified apparent gravity. Experimental measurements of motion threshold and rising velocity compare well with theoretical developments, thus suggesting practical approaches to control and tune bubble motion in confined environments.

This study explores heat and turbulent modulation in three-dimensional multiphase Rayleigh–Bénard convection using direct numerical simulations. Two immiscible fluids with identical reference density undergo systematic variations in dispersed-phase volume fractions, $0.0 \leq \varPhi \leq 0.5$, and ratios of dynamic viscosity, $\lambda _{\mu }$, and thermal diffusivity, $\lambda _{\alpha }$, within the range $[0.1\unicode{x2013}10]$. The Rayleigh, Prandtl, Weber and Froude numbers are held constant at $10^8$, $4$, $6000$ and $1$, respectively. Initially, when both fluids share the same properties, a 10 % Nusselt number increase is observed at the highest volume fractions. In this case, despite a reduction in turbulent kinetic energy, droplets enhance energy transfer to smaller scales, smaller than those of single-phase flow, promoting local mixing. By varying viscosity ratios, while maintaining a constant Rayleigh number based on the average mixture properties, the global heat transfer rises by approximately 25 % at $\varPhi =0.2$ and $\lambda _{\mu }=10$. This is attributed to increased small-scale mixing and turbulence in the less viscous carrier phase. In addition, a dispersed phase with higher thermal diffusivity results in a 50 % reduction in the Nusselt number compared with the single-phase counterpart, owing to faster heat conduction and reduced droplet presence near walls. The study also addresses droplet-size distributions, confirming two distinct ranges dominated by coalescence and breakup with different scaling laws.

We describe the rising trajectory of bubbles in isotropic turbulence and quantify the slowdown of the mean rise velocity of bubbles with sizes within the inertial subrange. We perform direct numerical simulations of bubbles, for a wide range of turbulence intensity, bubble inertia and deformability, with systematic comparison with the corresponding quiescent case, with Reynolds number at the Taylor microscale from 38 to 77. Turbulent fluctuations randomise the rising trajectory and cause a reduction of the mean rise velocity $\tilde {w}_b$ compared with the rise velocity in quiescent flow $w_b$. The decrease in mean rise velocity of bubbles $\tilde {w}_b/w_b$ is shown to be primarily a function of the ratio of the turbulence intensity and the buoyancy forces, described by the Froude number $Fr=u'/\sqrt {gd}$, where $u'$ is the root-mean-square velocity fluctuations, $g$ is gravity and $d$ is the bubble diameter. The bubble inertia, characterised by the ratio of inertial to viscous forces (Galileo number), and the bubble deformability, characterised by the ratio of buoyancy forces to surface tension (Bond number), modulate the rise trajectory and velocity in quiescent fluid. The slowdown of these bubbles in the inertial subrange is not due to preferential sampling, as is the case with sub-Kolmogorov bubbles. Instead, it is caused by the nonlinear drag–velocity relationship, where velocity fluctuations lead to an increased average drag. For $Fr > 0.5$, we confirm the scaling $\tilde {w}_b / w_b \propto 1 / Fr$, as proposed previously by Ruth et al. (J. Fluid Mech., vol. 924, 2021, p. A2), over a wide range of bubble inertia and deformability.

We studied the transport and deposition behaviour of point particles in Rayleigh–Bénard convection cells subjected to Couette-type wall shear. Direct numerical simulations (DNSs) are performed for Rayleigh number ($Ra$) in the range $10^{7} \leq Ra \leq ~10^9$ with a fixed Prandtl number $Pr = 0.71$, while the wall-shear Reynolds number ($Re_w$) is in the range $0 \leq Re_w \leq ~12\,000$. With the increase of $Re_w$, the large-scale rolls expanded horizontally, evolving into zonal flow in two-dimensional simulations or streamwise-oriented rolls in three-dimensional simulations. We observed that, for particles with a small Stokes number ($St$), they either circulated within the large-scale rolls when buoyancy dominated or drifted near the walls when shear dominated. For medium $St$ particles, pronounced spatial inhomogeneity and preferential concentration were observed regardless of the prevailing flow state. For large $St$ particles, the turbulent flow structure had a minor influence on the particles’ motion; although clustering still occurred, wall shear had a negligible influence compared with that for medium $St$ particles. We then presented the settling curves to quantify the particle deposition ratio on the walls. Our DNS results aligned well with previous theoretical predictions, which state that small $St$ particles settle with an exponential deposition ratio and large $St$ particles settle with a linear deposition ratio. For medium $St$ particles, where complex particle–turbulence interaction emerges, we developed a new model describing the settling process with an initial linear stage followed by a nonlinear stage. Unknown parameters in our model can be determined either by fitting the settling curves or using empirical relations. Compared with DNS results, our model also accurately predicts the average residence time across a wide range of $St$ for various $Re_w$.

Diffusion-driven flow is a boundary layer flow arising from the interplay of gravity and diffusion in density-stratified fluids when a gravitational field is non-parallel to an impermeable solid boundary. This study investigates diffusion-driven flow within a nonlinearly density-stratified fluid confined between two tilted parallel walls. We introduce an asymptotic expansion inspired by the centre manifold theory, where quantities are expanded in terms of derivatives of the cross-sectional averaged stratified scalar (such as salinity or temperature). This technique provides accurate approximations for velocity, density and pressure fields. Furthermore, we derive an evolution equation describing the cross-sectional averaged stratified scalar. This equation takes the form of the traditional diffusion equation but replaces the constant diffusion coefficient with a positive-definite function dependent on the solution's derivative. Numerical simulations validate the accuracy of our approximations. Our investigation of the effective equation reveals that the density profile depends on a non-dimensional parameter denoted as $\gamma$ representing the flow strength. In the large $\gamma$ limit, the system is approximated by a diffusion process with an augmented diffusion coefficient of $1+\cot ^{2}\theta$, where $\theta$ signifies the inclination angle of the channel domain. This parameter regime is where diffusion-driven flow exhibits its strongest mixing ability. Conversely, in the small $\gamma$ regime, the density field behaves like pure diffusion with distorted isopycnals. Lastly, we show that the classical thin film equation aligns with the results obtained using the proposed expansion in the small $\gamma$ regime but fails to accurately describe the dynamics of the density field for large $\gamma$.

Aqueous foams coarsen with time due to gas diffusion through the liquid between the bubbles. The mean bubble size grows, and small bubbles vanish. However, coarsening is little understood for foams with an intermediate liquid content, particularly in the presence of surfactant-induced attractive forces between the bubbles, measured by the interface contact angle where thin films meet the bulk liquid. Rigorous bubble growth laws have yet to be developed, and the evolution of bulk foam properties is unclear. We present a quasistatic numerical model for coarsening in two-dimensional wet foams, focusing on growth laws and related bubble properties. The deformation of bubble interfaces is modelled using a finite-element approach, and the gas flow through both films and Plateau borders is approximated. We give results for disordered two-dimensional wet foams with $256$ to $1024$ bubbles, at liquid fractions from $2\,\%$ to $25\,\%$, beyond the zero-contact-angle unjamming transition, and with contact angles up to $10^\circ$. Simple analytical models for the bubble pressures, film lengths and coarsening growth rates are developed to aid interpretation. If the contact angle is non-zero, we find that a prediction of the coarsening rate approaches a non-zero value as the liquid fraction is increased. We also find that an individual bubble's effective number of neighbours determines whether it grows or shrinks to a good approximation.

We examine the mechanisms responsible for the onset of the three-dimensional mode B instability in the wake behind a circular cylinder. We show that it is possible to explicitly account for the stabilising effect of spanwise viscous diffusion and then demonstrate that the remaining mechanisms involved in this short-wavelength instability are preserved in the limit of zero wavelength. Using the resulting simplified equations, we show that perturbations in different fluid particles interact only through the in-plane viscous diffusion which turns out to have a destabilising effect. We also show that in the presence of viscous diffusion, the closed trajectories which had been conjectured to play a crucial role in the onset of the mode B instability are not actually a prerequisite for the growth of mode B type perturbations. We combine these observations to identify the three essential ingredients for the development of the mode B instability: (i) the amplification of perturbations in the braid regions due to the stretching mechanism; and the spreading of perturbations through (ii) viscous diffusion, and (iii) cross-flow advection which transports fluid between the two braid regions on either side of the cylinder. Finally, we develop a simple criterion that allows the prediction of the regions where three-dimensional short-wavelength perturbations are amplified by the stretching mechanism. The approach used in our study is general and has the potential to give insights into the onset of three-dimensionality via short-wavelength instabilities in other flows.

Wave-assisted propulsion (WAP) systems directly convert wave energy into thrust using elastically mounted hydrofoils. The wave conditions as well as the design of the hydrofoil drive the fluid–structure interaction of the hydrofoil and, consequently, its performance. We employ simulations using a sharp-interface immersed boundary method to examine the effect of three key parameters on the flow physics, the fluid–structure interaction as well as thrust performance of these systems – the stiffness of the torsional spring, the location of the pitch axis and the Strouhal number. We demonstrate the utility of ‘maps’ of energy exchange between the flow and the hydrofoil system, as a way to understand and predict these characteristics. The force-partitioning method (FPM) is used to decompose the pressure forces into interpretable components and to quantify the mechanisms associated with thrust generation. Based on the results from FPM, a phenomenological model for the thrust generated by the WAP foil is presented. The parameters associated with this model are estimated based on data from over 450 distinct simulations. The predictions of the model are compared with the simulations and the use of this model for guiding WAP design is discussed.