A liquid drop impacting a rigid substrate undergoes deformation and spreading due to normal reaction forces, which are counteracted by surface tension. On a non-wetting substrate, the drop subsequently retracts and takes off. Our recent work (Zhang et al., Phys. Rev. Lett., vol. 129, 2022, 104501) revealed two peaks in the temporal evolution of the normal force $F(t)$ – one at impact and another at jump-off. The second peak coincides with a Worthington jet formation, which vanishes at high viscosities due to increased viscous dissipation affecting flow focusing. In this article, using experiments, direct numerical simulations and scaling arguments, we characterize both the peak amplitude $F_1$ at impact and the one at takeoff ($F_2$) and elucidate their dependency on the control parameters: the Weber number $We$ (dimensionless impact kinetic energy) and the Ohnesorge number $Oh$ (dimensionless viscosity). The first peak amplitude $F_1$ and the time $t_1$ to reach it depend on inertial time scales for low viscosity liquids, remaining nearly constant for viscosities up to 100 times that of water. For high viscosity liquids, we balance the rate of change in kinetic energy with viscous dissipation to obtain new scaling laws: $F_1/F_\rho \sim \sqrt {Oh}$ and $t_1/\tau _\rho \sim 1/\sqrt {Oh}$, where $F_\rho$ and $\tau _\rho$ are the inertial force and time scales, respectively, which are consistent with our data. The time $t_2$ at which the amplitude $F_2$ appears is set by the inertiocapillary time scale $\tau _\gamma$, independent of both the viscosity and the impact velocity of the drop. However, these properties dictate the magnitude of this amplitude.

We provide scaling relations for the Nusselt number $Nu$ and the friction coefficient $C_{S}$ in sheared Rayleigh–Bénard convection, i.e. in Rayleigh–Bénard flow with Couette- or Poiseuille-type shear forcing, by extending the Grossmann & Lohse (J. Fluid Mech., vol. 407, 2000, pp. 27–56, Phys. Rev. Lett., vol. 86, 2001, pp. 3316–3319, Phys. Rev. E, vol. 66, 2002, 016305, Phys. Fluids, vol. 16, 2004, pp. 4462–4472) theory to sheared thermal convection. The control parameters for these systems are the Rayleigh number $Ra$, the Prandtl number $Pr$ and the Reynolds number $Re_S$ that characterises the strength of the imposed shear. By direct numerical simulations and theoretical considerations, we show that, in turbulent Rayleigh–Bénard convection, the friction coefficients associated with the applied shear and the shear generated by the large-scale convection rolls are both well described by Prandtl's (Ergeb. Aerodyn. Vers. Gött., vol. 4, 1932, pp. 18–29) logarithmic friction law, suggesting some kind of universality between purely shear-driven flows and thermal convection. These scaling relations hold well for $10^6 \leq Ra \leq 10^8$, $0.5 \leq Pr \leq 5.0$, and $0 \leq Re_S \leq 10^4$.

We employ direct numerical simulations to investigate the heat transfer and flow structures in turbulent Rayleigh–Bénard convection in both cylindrical cells and laterally periodic domains, spanning an unprecedentedly wide range of aspect ratios $0.075 \leqslant \varGamma \leqslant 32$. We focus on Prandtl number ${Pr}=1$ and Rayleigh numbers ${{Ra}}=2\times 10^7$ and ${{Ra}}=10^8$. In both cases, with increasing aspect ratio, the heat transfer first increases, then reaches a maximum (which is more pronounced for the cylindrical case due to confinement effects), and then slightly goes down again before it finally saturates at the large aspect ratio limit, which is achieved already at $\varGamma \approx 4$. Already for $\varGamma \gtrsim 0.75$, the heat transfers in both cylindrical and laterally periodic domains become identical. The large-$\varGamma$ limit for the volume-integrated Reynolds number and the boundary layer thicknesses are also reached at $\varGamma \approx 4$. However, while the integral flow properties converge at $\varGamma \approx 4$, the confinement of a cylindrical domain impacts the temperature and velocity variance distributions up to $\varGamma \approx 16$, as thermal superstructures cannot form close to the sidewall.

The paper by Castaing et al. (J. Fluid Mech., vol. 204, 1989, pp. 1–30) on turbulent Rayleigh–Bénard convection has been one of the most impactful papers on the subject – not by giving the right and complete answers but by developing versatile concepts and by asking the right questions, namely: (i) What is the overall flow organization? (ii) What is the dependence of the Nusselt number ${\textit {Nu}}$ (the dimensionless heat transport) on the Rayleigh number ${\textit {Ra}}$ (the thermal driving strength)? (iii) What is the ultimate state of turbulence for extremely large ${\textit {Ra}}$? Thanks to Castaing et al. having asked the right questions, the field has made tremendous progress over the last 35 years.

Turbulent shear flows driven by a combination of a pressure gradient and buoyancy forcing are investigated using direct numerical simulations. Specifically, we consider the set-up of a differentially heated vertical channel subject to a Poiseuille-like horizontal pressure gradient. We explore the response of the system to its three control parameters: the Grashof number $Gr$, the Prandtl number $Pr$, and the Reynolds number $Re$ of the pressure-driven flow. From these input parameters, the relative strength of buoyancy driving to the pressure gradient can be quantified by the Richardson number $Ri=Gr/Re^2$. We compare the response of the mixed vertical convection configuration to that of mixed Rayleigh–Bénard convection, and find a nearly identical behaviour, including an increase in wall friction at higher $Gr$, and a drop in the heat flux relative to natural convection for $Ri=O(1)$. This closely matched response is despite vastly different flow structures in the systems. No large-scale organisation is visible in visualisations of mixed vertical convection – an observation that is confirmed quantitatively by spectral analysis. This analysis, combined with a statistical description of the wall heat flux, highlights how moderate shear suppresses the growth of small-scale plumes and reduces the likelihood of extreme events in the local wall heat flux. Vice versa, starting from a pure shear flow, the addition of thermal driving enhances the drag due to the emission of thermal plumes.

Melting and solidification in periodically time-modulated thermal convection are relevant for numerous natural and engineering systems, for example, glacial melting under periodic sun radiation and latent thermal energy storage under periodically pulsating heating. It is highly relevant for the estimation of melt rate and melt efficiency management. However, even the dynamics of a solid–liquid interface shape subjected to a simple sinusoidal heating has not yet been investigated in detail. In this paper, we offer a better understanding of the modulation frequency dependence of the melting and solidification front. We numerically investigate periodic melting and solidification in turbulent convective flow with the solid above and the melted liquid below, and sinusoidal heating at the bottom plate with the mean temperature equal to the melting temperature. We investigate how the periodic heating can prevent the full solidification, and the resulting flow structures and the quasi-equilibrium interface height. We further study the dependence on the heating modulation frequency. As the frequency decreases, we found two distinct regimes, which are ‘partially solid’ and ‘fully solid’. In the fully solid regime, the liquid freezes completely, and the effect of the modulation is limited. In the partially solid regime, the solid partially melts, and a steady or unsteady solid–liquid interface forms depending on the frequency. The interface height can be derived based on the energy balance through the interface. In the partially solid regime, the interface height oscillates periodically, following the frequency of modulation. Here, we propose a perturbation approach that can predict the dependency of the oscillation amplitude on the modulation frequency.

Frozen water might appear opaque since gas bubbles can get trapped in the ice during the freezing process. They nucleate and then grow near the advancing solidification front, due to the formation of a gas supersaturation region in its vicinity. A delicate interplay between the rate of mass transfer and the rate of freezing dictates the final shapes and sizes of the entrapped gas bubbles. In this work, we experimentally and numerically investigate the initial growth of such gas bubbles that nucleate and grow near the advancing ice front. We show that the initial growth of these bubbles is governed by diffusion and is enhanced due to a combination of the presence of the background gas concentration gradient and the motion of the approaching front. Additionally, we recast the problem into that of mass transfer to a moving spherical object in a homogeneous concentration field, finding good agreement between our experimental data and the existing scaling relations for that latter problem. Lastly, we address how fluid flow around the bubble might further affect this growth and qualitatively explore this through numerical simulations.

We consider the process of convective dissolution in a homogeneous and isotropic porous medium. The flow is unstable due to the presence of a solute that induces a density difference responsible for driving the flow. The mixing dynamics is thus driven by a Rayleigh–Taylor instability at the pore scale. We investigate the flow at the scale of the pores using Hele-Shaw type experiment with bead packs, two-dimensional direct numerical simulations and physical models. Experiments and simulations have been specifically designed to mimic the same flow conditions, namely matching porosities, high Schmidt numbers and linear dependency of fluid density with solute concentration. In addition, the solid obstacles of the medium are impermeable to fluid and solute. We characterise the evolution of the flow via the mixing length, which quantifies the extension of the mixing region and grows linearly in time. The flow structure, analysed via the centreline mean wavelength, is observed to grow in agreement with theoretical predictions. Finally, we analyse the dissolution dynamics of the system, quantified through the mean scalar dissipation, and three mixing regimes are observed. Initially, the evolution is controlled by diffusion, which produces solute mixing across the initial horizontal interface. Then, when the interfacial diffusive layer is sufficiently thick, it becomes unstable, forming finger-like structures and driving the system into a convection-dominated phase. Finally, when the fingers have grown sufficiently to touch the horizontal boundaries of the domain, the mixing reduces dramatically due to the absence of fresh unmixed fluid. With the aid of simple physical models, we explain the physics of the results obtained numerically and experimentally. The solute evolution presents a self-similar behaviour, and it is controlled by different length scales in each stage of the mixing process, namely the length scale of diffusion, the pore size and the domain height.

Evaporation of multi-component liquid mixtures in confined geometries, such as capillaries, is crucial in applications such as microfluidics, two-phase cooling devices and inkjet printing. Predicting the behaviour of such systems becomes challenging because evaporation triggers complex spatio-temporal changes in the composition of the mixture. These changes in composition, in turn, affect evaporation. In the present work, we study the evaporation of aqueous glycerol solutions contained as a liquid column in a capillary tube. Experiments and direct numerical simulations show three evaporation regimes characterised by different temporal evolutions of the normalised mass transfer rate (or Sherwood number $Sh$), namely $Sh (\tilde{t} ) = 1$, $Sh \sim 1/\sqrt {\tilde{t} }$ and $Sh \sim \exp (-\tilde{t} )$, where $\tilde {t}$ is a normalised time. We present a simplistic analytical model that shows that the evaporation dynamics can be expressed by the classical relation $Sh = \exp ( \tilde{t} )\,\mathrm {erfc} ( \sqrt {\tilde{t} })$. For small and medium $\tilde{t}$, this expression results in the first and second of the three observed scaling regimes, respectively. This analytical model is formulated in the limit of pure diffusion and when the penetration depth $\delta (t)$ of the diffusion front is much smaller than the length $L(t)$ of the liquid column. When $\delta \approx L$, finite-length effects lead to $Sh \sim \exp (-\tilde{t} )$, i.e. the third regime. Finally, we extend our analytical model to incorporate the effect of advection and determine the conditions under which this effect is important. Our results provide fundamental insights into the physics of selective evaporation from a multi-component liquid column.

Direct numerical simulations are utilised to investigate mass-transfer processes at gas-evolving electrodes that experience successive formation and detachment of bubbles. The gas–liquid interface is modelled employing an immersed boundary method. We simulate the growth phase of the bubbles followed by their departure from the electrode surface in order to study the mixing induced by these processes. We find that the growth of the bubbles switches from a diffusion-limited mode at low to moderate fractional bubble coverages of the electrode to a reaction-limited growth dynamics at high coverages. Furthermore, our results indicate that the net transport within the system is governed by the effective buoyancy driving induced by the rising bubbles and that mechanisms commonly subsumed under the term ‘microconvection’ do not significantly affect the mass transport. Consequently, the resulting gas transport for different bubble sizes, current densities and electrode coverages can be collapsed onto one single curve and only depends on an effective Grashof number. The same holds for the mixing of the electrolyte when additionally taking the effect of surface blockage by attached bubbles into account. For the gas transport to the bubble, we find that the relevant Sherwood numbers also collapse onto a single curve when accounting for the driving force of bubble growth, incorporated in an effective Jakob number. Finally, linking the hydrogen transfer rates at the electrode and the bubble interface, an approximate correlation for the gas-evolution efficiency has been established. Taken together, these findings enable us to deduce parametrisations for all response parameters of the systems.

Rotational effects are commonly neglected when considering the dynamics of freely rising or settling isotropic particles. Here, we demonstrate that particle rotations play an important role for rising as well as for settling cylinders in situations when mass eccentricity, and thereby a new pendulum time scale, is introduced to the system. We employ two-dimensional simulations to study the motion of a single cylinder in a quiescent unbounded incompressible Newtonian fluid. This allows us to vary the Galileo number, density ratio, relative moment of inertia (MOI) and centre-of-mass (COM) offset systematically and beyond what is feasible experimentally. For certain buoyant density ratios, the particle dynamics exhibits a resonance mode, during which the coupling via the Magnus lift force causes a positive feedback between translational and rotational motions. This mode results in vastly different trajectories with significantly larger rotational and translational amplitudes and an increase of the drag coefficient easily exceeding a factor two. We propose a simple model that captures how the occurrence of the COM offset induced resonance regime varies, depending on the other input parameters, specifically the density ratio, the Galileo number and the relative MOI. Remarkably, depending on the input parameters, resonance can be observed for COM offsets as small as a few per cent of the particle diameter, showing that the particle dynamics can be highly sensitive to this parameter.

For dissolving active oil droplets in an ambient liquid, it is generally assumed that the Marangoni effect results in repulsive interactions, while the buoyancy effects caused by the density difference between the droplets, diffusing product and the ambient fluid are usually neglected. However, it has been observed in recent experiments that active droplets can form clusters due to buoyancy-driven convection (Krüger et al., Eur. Phys. J. E, vol. 39, 2016, pp. 1–9). In this study we numerically analyse the buoyancy effect, in addition to the propulsion caused by Marangoni flow (with its strength characterized by the Péclet number $Pe$). The buoyancy effects have their origin in (i) the density difference between the droplet and the ambient liquid, which is characterized by the Galileo number $Ga$; and (ii) the density difference between the diffusing product (i.e. filled micelles) and the ambient liquid, which can be quantified by a solutal Rayleigh number $Ra$. We analyse how the attracting and repulsing behaviour of neighbouring droplets depends on the control parameters $Pe$, $Ga$ and $Ra$. We find that while the Marangoni effect leads to the well-known repulsion between the interacting droplets, the buoyancy effect of the reaction product leads to buoyancy-driven attraction. At sufficiently large $Ra$, even collisions between the droplets can take place. Our study on the effect of $Ga$ further shows that with increasing $Ga$, the collision becomes delayed. Moreover, we derive that the attracting velocity of the droplets, which is characterized by a Reynolds number $Re_d$, is proportional to $Ra^{1/4}/( \ell /R)$, where $\ell /R$ is the distance between the neighbouring droplets normalized by the droplet radius. Finally, we numerically obtain the repulsive velocity of the droplets, characterized by a Reynolds number $Re_{rep}$, which is proportional to $PeRa^{-0.38}$. The balance of attractive and repulsive effect leads to $Pe\sim Ra^{0.63}$, which agrees well with the transition curve between the regimes with and without collision.

Iceberg melting is a critical factor for climate change. However, the shape of an iceberg is an often neglected aspect of its melting process. Our study investigates the influence of different ice shapes and ambient flow velocities on melt rates by conducting direct numerical simulations of a simplified system of bluff body flow. Our study focuses on the ellipsoidal shape, with the aspect ratio as the control parameter. We found the shape plays a crucial role in the melting process, resulting in significant variations in the melt rate between different shapes. Without flow, the optimal shape for a minimal melt rate is the disk (two-dimensional) or sphere (three-dimensional), due to the minimal surface area. However, as the ambient flow velocity increases, the optimal shape changes with the aspect ratio. We find that ice with an elliptical shape (when the long axis is aligned with the flow direction) can melt up to 10 % slower than a circular shape when exposed to flowing water. Following the approach considered by Huang et al. (J. Fluid Mech., vol. 765, 2015, R3) for dissolving bodies, we provide a quantitative theoretical explanation for this optimal shape, based on the combined contributions from both surface-area effects and convective-heat-transfer effects. Our findings provide insight into the interplay between phase transitions and ambient flows, contributing to our understanding of the iceberg melting process and highlighting the need to consider the aspect-ratio effect in estimates of iceberg melt rates.

‘Storm oil’ – nearly water-insoluble oil poured into the ocean and acting as a surfactant – has been used since ancient times to smooth the waves on the ocean. It was first scientifically described by Benjamin Franklin (Phil. Trans. R. Soc. Lond., vol. 64, 1774, pp. 445–460). In a recent paper, by combining highly controlled experiments in a wave tank and direct numerical simulations, Erinin et al. (J. Fluid Mech., vol. 972, 2023, R5) have now beautifully revealed the strong effect of soluble surfactants on the dynamics of plunging breakers. Remarkably, it is not the change in surface tension which mainly matters, but the surface tension gradient which emerges through compression and dilation of the plunging breaker surface.

Two-dimensional horizontally periodic Rayleigh–Bénard convection between stress-free boundaries displays two distinct types of states, depending on the initial conditions. Roll states are composed of pairs of counter-rotating convection rolls. Windy states are dominated by strong horizontal wind (also called zonal flow) that is vertically sheared, precludes convection rolls and suppresses heat transport. Windy states occur only when the Rayleigh number $Ra$ is sufficiently above the onset of convection. At intermediate $Ra$ values, windy states can be induced by suitable initial conditions, but they undergo a transition to roll states after finite lifetimes. At larger $Ra$ values, where windy states have been observed for the full duration of simulations, it is unknown whether they represent chaotic attractors or only metastable states that would eventually undergo a transition to roll states. We study this question using direct numerical simulations of a fluid with a Prandtl number of 10 in a layer whose horizontal period is eight times its height. At each of seven $Ra$ values between $9\times 10^6$ and $2.25\times 10^7$ we have carried out 200 or more simulations, all from initial conditions leading to windy convection with finite lifetimes. The lifetime statistics at each $Ra$ indicate a memoryless process with survival probability decreasing exponentially in time. The mean lifetimes grow with $Ra$ approximately as $Ra^4$. This analysis provides no $Ra$ value at which windy convection becomes stable; it might remain metastable at larger $Ra$ with extremely long lifetimes.

The efficiency of water electrolysis is significantly impacted by the generation of micro- and nanobubbles on the electrodes. Here molecular dynamics simulations are used to investigate the dynamics of single electrolytic nanobubbles on nanoelectrodes. The simulations reveal that, depending on the value of current, nucleated nanobubbles either grow to an equilibrium state or grow unlimitedly and then detach. To account for these findings, the stability theory for surface nanobubbles is generalized by incorporating the electrolytic gas influx at the nanobubble's contact line and adopting a real gas law, leading to accurate predictions for the numerically observed transient growth and stationary states of the nanobubbles. With this theory, the minimum current for bubble detachment can also be derived analytically. In the detachment regime, the radius of the nanobubble first increases with time (t) as $R\propto t^{1/2}$ and then as $R\propto t^{1/3}$, up to bubble detachment.

The presence of salt in seawater strongly affects the melt rate and the shape evolution of ice, both of utmost relevance in ice–ocean interactions and thus for the climate. To get a better quantitative understanding of the physical mechanics at play in ice melting in salty water, we numerically investigate the lateral melting of an ice block in stably stratified saline water. The developing ice shape from our numerical results shows good agreement with the experiments and theory from Huppert & Turner (J. Fluid Mech., vol. 100, 1980, pp. 367–384). Furthermore, we find that the melt rate of ice depends non-monotonically on the mean ambient salinity: it first decreases for increasing salt concentration until a local minimum is attained, and then increases again. This non-monotonic behaviour of the ice melt rate is due to the competition among salinity-driven buoyancy, temperature-driven buoyancy and salinity-induced stratification. We develop a theoretical model based on the force balance which gives a prediction of the salt concentration for which the melt rate is minimal, and is consistent with our data. Our findings give insight into the interplay between phase transitions and double-diffusive convective flows.

Our understanding of the material organization of complex fluid flows has benefited recently from mathematical developments in the theory of objective coherent structures. These methods have provided a wealth of approaches that identify transport barriers in three-dimensional (3-D) turbulent flows. Specifically, theoretical advances have been incorporated into numerical algorithms that extract the most influential advective, diffusive and active barriers to transport from data sets in a frame-indifferent fashion. To date, however, there has been very limited investigation into these objectively defined transport barriers in 3-D unsteady flows with complicated spatiotemporal dynamics. Similarly, no systematic comparison of advective, diffusive and active barriers has been carried out in a 3-D flow with both thermally driven and mechanically modified structures. In our study, we utilize simulations of turbulent rotating Rayleigh–Bénard convection to uncover the interplay between advective transport barriers (Lagrangian coherent structures), material barriers to diffusive heat transport, and objective Eulerian barriers to momentum transport. For a range of (inverse) Rossby numbers, we identify each type of barrier and find intriguing relationships between momentum and heat transport that can be related to changes in the relative influence of mechanical and thermal forces. Further connections between bulk behaviours and structure-specific behaviours are also developed.

When an immiscible oil drop is immersed in a stably stratified ethanol–water mixture, the Marangoni flow on the surface of the drop can experience an oscillatory instability, so that the drop undergoes a transition from levitating to bouncing. The onset of the instability and its mechanisms have been studied previously (Li et al., Phys. Rev. Lett., vol. 126, 2021, 124502; Li et al., J. Fluid Mech., vol. 932, 2022, A11), yet the bouncing motion of the drop itself, which is a completely different problem, has not yet been investigated. Here we study how the bouncing characteristics (jumping height, rising and sinking time) depend on the control parameters (drop radius, stratification strength, drop viscosity). We first record experimentally the bouncing trajectories of drops of different viscosities in different stratifications. Then a simplified dynamical analysis is performed to get the scaling relations of the jumping height and the rising and sinking times. The rising and sinking time scales are found to depend on the drag coefficient $C_D^S$ of the drop in the stratified liquid, which is determined empirically for the current parameter space (Zhang et al., J. Fluid Mech., vol. 875, 2019, 622–656). For low-viscosity (5 cSt) oil drops, the results on the drag coefficient match those from the literature (Yick et al., J. Fluid Mech., vol. 632, 2009, pp. 49–68; Candelier et al., J. Fluid Mech., vol. 749, 2014, pp. 184–200). For high-viscosity (100 cSt) oil drops, the parameter space had not been explored and the drag coefficients are not readily available. Numerical simulations are therefore performed to provide external verification for the drag coefficients, which well match with the experimental results.

Ambient air cushions the impact of drops on solid substrates, an effect usually revealed by the entrainment of a bubble, trapped as the air squeezed under the drop drains and liquid–solid contact occurs. The presence of air becomes evident for impacts on very smooth surfaces, where the gas film can be sustained, allowing drops to bounce without wetting the substrate. In such a non-wetting situation, Mandre & Brenner (J. Fluid Mech., vol. 690, 2012, p. 148) numerically and theoretically evidenced that two physical mechanisms can act to prevent contact: surface tension and nonlinear advection. However, the advection dominated regime has remained hidden in experiments as liquid–solid contact prevents rebounds being realised at sufficiently large impact velocities. By performing impacts on superheated surfaces, in the so-called dynamical Leidenfrost regime (Tran et al., Phys. Rev. Lett., vol. 108, issue 3, 2012, p. 036101), we enable drop rebound at higher impact velocities, allowing us to reveal this regime. Using high-speed total internal reflection, we measure the minimal gas film thickness under impacting drops, and provide evidence for the transition from the surface tension to the nonlinear inertia dominated regime. We rationalise our measurements through scaling relationships derived by coupling the liquid and gas dynamics, in the presence of evaporation.