A novel theoretical model for bubble dynamics is established that simultaneously accounts for the liquid compressibility, phase transition, oscillation, migration, ambient flow field, etc. The bubble dynamics equations are presented in a unified and concise mathematical form, with clear physical meanings and extensibility. The bubble oscillation equation can be simplified to the Keller–Miksis equation by neglecting the effects of phase transition and bubble migration. The present theoretical model effectively captures the experimental results for bubbles generated in free fields, near free surfaces, adjacent to rigid walls, and in the vicinity of other bubbles. Based on the present theory, we explore the effect of the bubble content by changing the vapour proportion inside the cavitation bubble for an initial high-pressure bubble. It is found that the energy loss of the bubble shows a consistent increase with increasing Mach number and initial vapour proportion. However, the radiated pressure peak by the bubble at the collapse stage increases with decreasing Mach number and increasing vapour proportion. The energy analyses of the bubble reveal that the presence of vapour inside the bubble not only directly contributes to the energy loss of the bubble through phase transition but also intensifies the bubble collapse, which leads to greater radiation of energy into the surrounding flow field due to the fluid compressibility.

Stratified wakes past an isolated conical seamount are simulated at a Froude number of $Fr = 0.15$ and Rossby numbers of $Ro = 0.15$, 0.75 and $\infty$. The wakes exhibit a Kármán vortex street, unlike their unstratified, non-rotating counterpart. Vortex structures are studied in terms of large-scale global modes, as well as spatially localised vortex evolution, with a focus on rotation effects. The global modes are extracted by spectral proper orthogonal decomposition (SPOD). For all three studied $Ro$ ranging from mesoscale, submesoscale and non-rotating cases, the frequency of the SPOD modes at different heights remains coupled as a global constant. However, the shape of the SPOD modes changes from slanted ‘tongues’ at zero rotation ($Ro=\infty$) to tall hill-height columns at strong rotation ($Ro=0.15$). A novel method for vortex centre tracking shows that, in all three cases, the vortices at different heights advect uniformly at approximately $0.9U_{\infty }$ beyond the near wake, consistent with the lack of variability of the global modes. Under system rotation, cyclonic vortices and anticyclonic vortices (AVs) present considerable asymmetry, especially at $Ro = 0.75$. The vorticity distribution as well as the stability of AVs are tracked downstream using statistics conditioned to the identified vortex centres. At $Ro=0.75$, intense AVs with relative vorticity up to $\omega _z/f_{c}=-2.4$ (where $\omega_z$ is the vertical vorticity and $f_c$ is the Coriolis frequency) are seen with small regions of instability and they maintain large $\omega _z/f_{c}$ magnitude in the far wake. Recent stability analysis that accounts for stratification and viscosity is found to improve on earlier criteria and show that these intense AVs are stable.

A stream of free-falling particles from a rectangular hopper, hereafter called a ‘curtain’, was characterised systematically using a well-resolved, non-intrusive optical shadowgraphic method, to reveal both an additional axial region and an additional dilute region distributed laterally on either side of the curtain, relative to those identified previously. The effects of particle size and hopper outlet thickness on the evolution of the particle curtain were separately isolated, whilst measuring particle mass flow rate. The curtains were characterised into four distinct axial regions, namely a near-field expansion region near to the hopper exit, a neck zone where the curtain contracts, a region of intermediate-field expansion and a far field with particles reaching terminal velocity. The initial expansion half-angle, $2.3^{\circ } \le \alpha \le 4.8^{\circ }$, was found to be insensitive to particle size, but to increase with hopper outlet thickness. The ‘trough’ in the neck zone was deduced to be caused by a pressure gradient driven by particle acceleration. The curtain expansion rate at the intermediate field was found to increase with a decrease in particle size and hopper outlet thickness. The outermost dilute-particle region was deduced to be caused by collisions, induced by gradients in the velocity profile near to the hopper exit. New dimensionless analysis reveals that the dynamics of curtains can be characterised broadly into two regimes, one in which the aerodynamics is dominant and the other where it is weak. Curtain transmittance was found to scale with the Froude number, highlighting the importance of particle momentum.

Unlike in solids, heat transfer in fluids can be greatly enhanced due to the presence of convection. Under gravity, an unevenly distributed temperature field results in differences in buoyancy, driving fluid motion that is seen in Rayleigh–Bénard convection (RBC). In RBC, the overall heat flux is found to have a power-law dependence on the imposed temperature difference, with enhanced heat transfer much beyond thermal conduction. In a bounded domain of fluid such as a cube, how RBC responds to thermal perturbations from the vertical sidewall is not clear. Will sidewall heating or cooling modify flow circulation and heat transfer? We address these questions experimentally by adding heat to one side of the RBC. Through careful flow, temperature and heat flux measurements, the effects of adding side heating to RBC are examined and analysed, where a further enhancement of flow circulation and heat transfer is observed. Our results also point to a direct and simple control of the classical RBC system, allowing further manipulation and control of thermal convection through sidewall conditions.

The analysis of Scott (J. Fluid Mech., vol. 741, 2014, pp. 316–349) is implemented numerically. Decaying turbulence is confined to a channel between two infinite, parallel, rotating walls. The Rossby and Ekman numbers are supposed small, the former condition making nonlinearity small, while the latter allows the turbulence to persist for the many rotational periods needed for the small nonlinearity to be effective. The flow is expressed as a combination of inertial waveguide modes, indexed by a two-dimensional wave vector $\boldsymbol{k}$ and an integer n. The $n = 0$ modes form a two-dimensional component of the flow, whereas the remainder is the wave component, on which attention is focused in this article. Assuming statistical axisymmetry and homogeneity in directions parallel to the walls, the second-order moments of the mode amplitudes yield a spectral matrix ${A_{nm}}(k,t)$ (where $k = |\boldsymbol{k} |$), of which the diagonal elements describe the distribution of energy over different modes. Wave-turbulence analysis provides an equation governing the time evolution of ${A_{nn}}$, $n \ne 0$, the wave spectra, which forms the basis for the present work. The initial distribution of energy is Gaussian and depends on a parameter $\varXi$, the initial spectral width. The problem has two other parameters, ${\beta _w}$ and ${\beta _v}$, which correspond to two distinct viscous dissipative mechanisms: wall damping due to boundary layers and volumetric damping by viscous effects throughout the flow. Results obtained by numerical solution include the time evolution of the total wave energy, E, and the detailed description of its distribution over k and n provided by ${A_{nn}}(k)$.

The reflection of a shock pulse at a liquid–gas interface occurs in many applications, from lithotripsy to underwater explosions and additive manufacturing. In linear theory, reflection and transmission at an interface depend only on the impedance difference, but this does not hold for a nonlinear pulse. This work develops an analytical framework for computing the reflection and transmission coefficients for an impulsive shock wave at a liquid–gas interface. The problem is treated analytically by considering idealised pulses and solving a series of consecutive Riemann problems. These correspond to the initial interaction with the interface and important subsequent wave interactions that enable a complete description of the process to be obtained. Comparisons with numerical and existing analytical approaches are made for the case of a water–air interface. In the acoustic limit, the method produces results identical to those of linear acoustic theory. As the pulse strength increases, the proposed method agrees well with numerical simulation results, whereas existing analytical methods that consider only the interface fail. We detail how a reflecting pulse can put water into tension without any incident negative pressure. It is further shown that the magnitude of the reflection coefficient decreases with increasing incident shock pressure, and the reflected pulse widens. Reflections of pulses with positive and negative pressures temporarily create negative pressure regions with greater magnitude than the incident pulse. Finally, we consider non-idealised waves. Comparisons with simulations show that the reflection characteristics can be explained qualitatively using the analytical method, and the reflection coefficients are predicted accurately.

Turbulent convection in the interiors of the Sun and the Earth occurs at high Rayleigh numbers $Ra$, low Prandtl numbers $Pr$, and different levels of rotation rates. To understand the combined effects better, we study rotating turbulent convection for $Pr = 0.021$ (for which some laboratory data corresponding to liquid metals are available), and varying Rossby numbers $Ro$, using direct numerical simulations in a slender cylinder of aspect ratio 0.1; this confinement allows us to attain high enough Rayleigh numbers. We are motivated by the earlier finding in the absence of rotation that heat transport at high enough $Ra$ is similar between confined and extended domains. We make comparisons with higher aspect ratio data where possible. We study the effects of rotation on the global transport of heat and momentum as well as flow structures (a) for increasing rotation at a few fixed values of $Ra$, and (b) for increasing $Ra$ (up to $10^{10}$) at the fixed, low Ekman number $1.45 \times 10^{-6}$. We compare the results with those from unity $Pr$ simulations for the same range of $Ra$ and $Ro$, and with the non-rotating case over the same range of $Ra$ and low $Pr$. We find that the effects of rotation diminish with increasing $Ra$. These results and comparison studies suggest that for high enough $Ra$, rotation alters convective flows in a similar manner for small and large aspect ratios, so useful insights on the effects of high thermal forcing on convection can be obtained by considering slender domains.

Turbulent flows over rough surfaces can be encountered in a wide range of engineering applications. Despite the progress made after several decades of studies, the prediction of drag and roughness function from the surface geometrical parameters remains an open question. Several methods have shown encouraging results. However, they lack generality and present some scatter in the data. In this paper we propose a new parameter, the effective distribution ($ED$), which lays foundation on the effective slope with some changes to take into account the sheltering effect of large roughness elements and the drag induced by pinnacles higher than the average roughness elements. To develop this new correlation between geometrical features of the wall and the drag, we performed a set of simulations of the turbulent flow over a rough surface made of triangular elements varying their height and spatial distribution. The $ED$ correlates quite well both with the drag and the roughness function for a wide range of cases having different mean roughness height, skewness and kurtosis. To further validate the $ED$, and assessing how it can be generalized to real rough wall, an irregular wall made from the superposition of random sinusoidal function was considered. Results were consistent with the correlation here presented.

The development of thermal boundary layers and plume near a section-triangular roof under different isothermal heating conditions has been the focus of numerous numerical studies. However, flow transition in this type of flow has never been observed experimentally. Here, phase-shifting interferometry and thermistor measurements are employed to experimentally observe and quantify the flow transitions in a buoyancy-driven flow over an isothermal section-triangular roof. Visualisation of temperature contours is conducted across a wide range of Rayleigh numbers from laminar at 103 to chaotic state at 4 × 106. Power spectral density of the temperature measurements reveals the type of bifurcations developing as the Rayleigh number is increased. This flow transition is characterised as a complex bifurcation route with the presence of two fundamental frequencies, a low and a high frequency. We found that the thermal stratification in the environment plays a significant role in the flow transition. The spatial development of flow is also quantitatively and qualitatively described. In addition to clarifying flow transition in experiments, the work demonstrates the implementation of phase-shifting interferometry and punctual temperature measurements for characterisation of near-field flow over a heated surface.

The flow field of a bluff body, a circular disk, that moves horizontally in a stratified environment is studied using large-eddy simulations. Five levels of stratification (body Froude numbers of ${{Fr}} = 0.5, 1, 1.5, 2$ and $5$) are simulated at Reynolds number of ${{Re}} = 5000$ and Prandtl number of $Pr =1$. A higher ${{Re}} = 50\,000$ database at ${{Fr}} = 2, 10$ and $Pr =1$ is also examined for comparison. The wavelengths and amplitudes of steady lee waves are compared with a linear-theory analysis. Excellent agreement is found over the entire range of ${{Fr}}$ if an ‘equivalent body’ that includes the separation region is employed for the linear theory. For asymptotically large distances, the velocity amplitude varies theoretically as ${{Fr}}^{-1}$ but a correction owing to the dependence of the separation zone on ${{Fr}}$ is needed. The wake waves propagate in a narrow band of angles with the vertical, and have a wavelength that increases with increasing ${{Fr}}$. The envelope of wake waves, demarcated using buoyancy variance, exhibits self-similar behaviour. The higher ${{Re}}$ results are consistent with the buoyancy effects exhibited at the lower ${{Re}}$. The wake wave energy is larger at ${{Re}} = 50\,000$. Nevertheless, independent of ${{Fr}}$ and ${{Re}}$, the ratio of the wake wave potential energy to the wake turbulent energy increases to approximately 0.6–0.7 in the non-equilibrium stage showing their energetic importance besides suggesting universality in this statistic. There is a crossover of energetic dominance of lee waves at ${{Fr}} <2$ to wake-wave dominance at ${{Fr}} \approx 5$.

This study investigates the interaction between a freely rising, deformable bubble and a freely settling particle of the same size due to gravity. Initially, an in-line configuration is considered while varying the Bond, Galilei and Archimedes numbers. The study shows that as the bubble and particle approach each other, a liquid film forms between them that undergoes drainage. The formation of the liquid film leads to dissipation of kinetic energy, and for sufficiently large bubble velocities, particle flotation takes place. Increasing the Bond number causes the bubble to deform more severely, which may allow the particle to pass through the bubble as it ruptures. This work also considers an offset configuration, which shows that the bubble slides away from the particle, affecting its settling trajectory.

We investigate the heterogeneous cavitation phenomenon in water when a spherical surface is abruptly separated from a nearby flat substrate, at a distance of approximately 10 nm. By tracking the surface separation using Newton ring positions, and capturing the bubble evolution with a high-speed camera on a microscope, we compare our experimental findings with hydrodynamic predictions at low Reynolds numbers. Upon upward movement of the spherical surface, the resulting bubble develops branched fingers through the Saffman–Taylor instability. Simultaneously, negative liquid pressures in the range $\sim$10 atm are observed. These large tension values occasionally lead to secondary nucleation events. The bubble sizes satisfy a predicted Family–Vicsek scaling law where the bubble area is proportional to the inverse bubble lifetime. The fact that creeping flow cavitation bubbles are more short-lived the larger they are separates them from bubbles that are governed by inertial dynamics.

The penetration of a spherical vortex into turbulence is studied theoretically and experimentally. The characteristics of the vortex are first analysed from an integral perspective that reconciles the far-field dipolar flow with the near-field source flow. The influence of entrainment on the vortex drag force is elucidated, extending the Maxworthy (J. Fluid Mech., vol. 81, 1977, pp. 465–495) model to account for turbulent entrainment into the vortex movement and vortex penetration into an evolving turbulent field. The physics are explored numerically using a spherical vortex (initial radius $R_0$, speed $U_{v0}$), characterised by a Reynolds number $Re_0(=2R_0U_{v0}/\nu$, where $\nu$ is the kinematic viscosity) of 2000, moving into decaying homogeneous turbulence (root-mean-square $u_0$, integral scale $L$) with turbulent intensity $I_t=u_0/U_{v0}$. When the turbulence is absent ($I_t=0$), a wake volume flux leads to a reduction of vortex impulse that causes the vortex to slow down. In the presence of turbulence ($I_t> 0$), the loss of vortical material is enhanced and the vortex speed decreases until it is comparable to the local turbulent intensity and quickly fragments, penetrating a distance that scales as $I_t^{-1}$. In the experimental study, a vortex ($Re_0\sim 1490\unicode{x2013}5660$) propagating into a statistically steady, spatially varying turbulent field ($I_{ve}=0.02$ to 0.98). The penetration distance is observed to scale with the inverse of the turbulent intensity. Incorporating the spatially and temporally varying turbulent fields into the integral model gives a good agreement with the predicted trend of the vortex penetration distance with turbulent intensity and insight into its dependence on the structure of the turbulence.

In the present study, we propose a Reynolds analogy model for compressible wall turbulence. This model is demonstrated to be able to alleviate the defects of the generalized Reynolds analogy model (GRA) (Zhang et al., J. Fluid Mech., vol. 739, 2014, pp. 392–420), and maintains its success in describing the mean velocity–temperature relation. Furthermore, the present model is superior to the GRA in depicting the relationship between their fluctuating fields and also bridges the gap between the phenomenological model and the mathematical representation of the Reynolds analogy. The key points of the present model are validated by analysing the data of compressible wall-bounded turbulence with different Mach numbers, Reynolds numbers and wall thermal conditions.

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.

Microorganism motility often takes place within complex, viscoelastic fluid environments, e.g. sperm in cervicovaginal mucus and bacteria in biofilms. In such complex fluids, strains and stresses generated by the microorganism are stored and relax across a spectrum of length and time scales and the complex fluid can be driven out of its linear response regime. Phenomena not possible in viscous media thereby arise from feedback between the swimmer and the complex fluid, making swimming efficiency co-dependent on the propulsion mechanism and fluid properties. Here, we parameterize a flagellar motor and filament properties together with elastic relaxation and nonlinear shear-thinning properties of the fluid in a computational immersed boundary model. We then explore swimming efficiency, defined as a particular flow rate divided by the torque required to spin the motor, over this parameter space. Our findings indicate that motor efficiency (measured by the volumetric flow rate) can be boosted or degraded by relatively moderate or strong shear thinning of the viscoelastic environment.

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.

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.

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.