The present study aims to demonstrate the application of the oblique temperature gradient (OTG) to manipulate the buoyancy instabilities and resulting mixing. To accomplish this, we consider a model consisting of a liquid layer supported by a perfectly conducting bottom wall from below and the upper surface is exposed to ambient inert gas in a shallow container. The stability analysis employs the pseudospectral method and perturbation energy budget approach. The imposed OTG consists of a negative vertical temperature gradient (VTG), which leads to unstable density stratification, and a negative horizontal temperature gradient (HTG) component responsible for the Hadley circulation, which can be utilised to control instabilities. Without the HTG, a Rayleigh–Bénard mode of instability exists due to the imposed VTG termed as a VTG mode. The imposed HTG induces a linear VTG term and a quintic polynomial VTG term in the bottom wall-normal coordinate $y$. The induced linear VTG term reinforces the imposed VTG, while the induced quintic polynomial VTG term opposes the imposed VTG. For the vanishing Biot number ($Bi$), the induced quintic VTG term stabilises the VTG mode for the Prandtl number, $Pr=7$. However, for $Pr=0.01$, the Reynolds stresses associated with the Hadley circulation destabilise the VTG mode. For non-zero $Bi$, new streamwise and spanwise instability modes arise for $Pr=7$ due to the induced linear VTG term. Physical arguments show that the unstable density stratification near the gas–liquid interface caused by the synergistic effect of the imposed VTG and the induced linear VTG term lead to the existence of the predicted new modes. The present study thus shows that the strength of the HTG can be utilised to control the instability modes and critical parameters, thereby demonstrating the utility of the OTG in manipulating the buoyancy instabilities.

We study the energy stability of pressure-driven laminar magnetohydrodynamic flow in a rectangular duct with a transverse homogeneous magnetic field and electrically insulating walls. For sufficiently strong fields, the laminar velocity distribution has a uniform core and convex Hartmann and Shercliff boundary layers on the walls perpendicular and parallel to the magnetic field. The problem is discretized by a double expansion in Chebyshev polynomials in the cross-stream coordinates. The linear eigenvalue problem for the critical Reynolds number depends on the streamwise wavenumber, Hartmann number and the aspect ratio. We consider the limits of small and large aspect ratios in order to compare with stability models based on one-dimensional base flows. For large aspect ratios, we find good numerical agreement with results based on the quasi-two-dimensional approximation. The lift-up mechanism dominates in the limit of a zero streamwise wavenumber and provides a linear dependence between the critical Reynolds and Hartmann numbers in the duct. As the aspect ratio is reduced away from unity, the duct results converge to Orr's original energy stability result for spanwise uniform perturbations imposed on the plane Poiseuille base flow. We also examine different possible symmetries of eigenmodes as well as the purely hydrodynamic case in the duct geometry.

A turbulent two-vortex system (T-2VS) is obtained by inserting analytical model wake vortices into very weak homogeneous isotropic turbulence (HIT) and by evolving them in time using large-eddy simulation until a turbulent state at statistical equilibrium is reached. The T-2VS is characterised as follows: circulation distribution of the vortices; energy of the mean and fluctuating fields; energy dissipation rate. It is also verified that essentially the same T-2VS is obtained when varying the initial model or initial HIT perturbation. A wall-resolved simulation of the T-2VS further interacting with a smooth ground is then performed at $Re_\varGamma = 2 \times 10^5$; this is $10 \times$ higher than in previous works, which allows us to better capture the high Reynolds number behaviour. The high release height of the T-2VS also ensures a physically correct approach to the ground. The results are compared with the literature and also to what is obtained for the case of non-turbulent vortices interacting with the same ground at the same Reynolds number. The flow topologies are discussed, and significant differences are highlighted regarding the separation of the boundary layer generated at the ground, and the way this secondary vorticity interacts with the primary vortices and makes them decay. The vortex trajectories are also measured, together with their circulation distribution and global circulation evolution, and the differences are discussed.

Bubble growth, departure and sliding in low-pressure flow boiling has received considerable attention in the past. However, most applications of boiling heat transfer rely on high-pressure flow boiling, for which very little is known, as experimental data are scarce and very difficult to obtain. In this work, we conduct an experiment using high-resolution optical techniques. By combining backlit shadowgraphy and phase-detection imaging, we track bubble shape and physical footprint with high spatial ($6\,\mathrm {\mu }{\rm m}$) and temporal ($33\,\mathrm {\mu }{\rm s}$) resolutions, as well as bubble size and position as bubbles nucleate and slide on top of the heated surface. We show that at pressures above 1 MPa bubbles retain a spherical shape throughout the growth and sliding process. We analytically derive non-dimensional numbers to correlate bubble velocity and liquid velocity throughout the turbulent boundary layer and predict the sliding of bubbles on the surface, solely from physical properties and the bubble growth rate. We also show that these non-dimensional solutions can be leveraged to formulate elementary criteria that predict the effect of pressure and flow rate on bubble departure diameter and growth time.

The long- and finite-wavelength instabilities of weakly viscoelastic film on an oscillating inclined plane are investigated. By using the Chebyshev series solution with the Floquet theory, the combined effects of viscoelasticity and forcing amplitude on instability are described when the inclined plane oscillates in streamwise and wall-normal directions. For long-wavelength instability, the solution to the eigenvalue problem is obtained analytically by the asymptotic expansion method. Results show that the Floquet exponent is independent of the wall-normal oscillation amplitude. The effects of inclined angle, gravity and surface tension on the stability of viscoelastic liquid film are also discussed. For finite-wavelength instability, numerical results corresponding to the wall-normal oscillation disclose that with the increase of viscoelasticity, the unstable gravitational boundary moves to a higher wavenumber, and the critical amplitudes of subharmonic and harmonic instabilities are reduced. The neutral curve of gravity instability for streamwise oscillatory flow is divided into two parts, and a stable bandwidth is formed for a large value of the viscoelastic parameter. Besides, a new oscillatory mode is identified at small angles of inclination, which will be enhanced if the effect of viscoelasticity is incorporated.

Rain drops form in clouds by collision of submillimetric droplets falling under gravity: larger drops fall faster than smaller ones and collect them on their path. The puzzling stability of fogs and non-precipitating warm clouds with respect to this avalanche mechanism has been a longstanding problem. How can droplets of diameter around 10 $\mathrm {\mu }$m have a low collision probability, inhibiting the cascade towards larger and larger drops? Here we review the dynamical mechanisms that have been proposed in the literature and quantitatively investigate the frequency of drop collisions induced by Brownian diffusion, electrostatics and gravity, using an open-source Monte Carlo code taking all of them into account. Inertia dominates over aerodynamic forces for large drops, when the Stokes number is larger than $1$. Thermal diffusion dominates over aerodynamic forces for small drops, when the Péclet number is smaller than $1$. We show that there exists a range of size (typically 3–30 $\mathrm {\mu }$m for water drops in air) where neither inertia nor Brownian diffusion are significant, leading to a gap in the collision rate. The effect is particularly important, due to the lubrication film forming between the drops immediately before collision, and secondarily to the long-range aerodynamic interaction. Two different mechanisms regularise the divergence of the lubrication force at vanishing separation: the transition to a non-continuum regime in the lubrication film, when the separation is comparable to the mean free path of air, and the induction of a flow inside the drops due to shear at their surfaces. In the gap between inertia-dominated and diffusion-dominated regimes, dipole–dipole electrostatic interactions becomes the major effect controlling the efficiency of drop collisions.

Electrohydrodynamic (EHD)-induced droplet emission is an efficient method for the production of micron- and submicron-sized droplets in technological applications. Existing studies propose several scaling laws to determine the size of the emitted droplet. However, they have usually focused on the tip streaming phenomena of a droplet when subjected to a uniform electric field. In most applications, a non-uniform distribution of the electric field is created owing to the nozzle-to-plate configuration. Here, we employ an arbitrary Lagrangian–Eulerian method to demonstrate the mechanism of the first droplet emission from an electrified liquid meniscus with a fixed volume hanging at the nozzle tip. The critical condition when tip streaming occurs is determined using our numerical results. A phase diagram in terms of the electric field and initial liquid volume is presented to obtain the commonly used jetting mode. The effects of the liquid volume, electric field strength and electrical conductivity of the liquid on the processes of jet formation and breakup are further investigated. We find a particularly non-monotonic dependence of the size of the emitted droplet on the electrical conductivity. These findings could be useful for generating microdroplets and improving injection frequency in EHD printing technology.

Turbulence in gravity-driven film flow is usually discussed in terms of three-dimensional solitary-wave pulses as they are frequently observed in flows along smooth walls. Here, we show that free-surface turbulence arises in films along rippled substrates, as they are commonly employed in process engineering applications, already at rather low Reynolds numbers from the irregular break-up of the solitary-wave fronts. Short waves in the capillary regime replace the broken solitary waves beyond a certain Reynolds number. The crossover coincides with the occurrence of steady three-dimensional surface patterns upstream, which suppress travelling waves. The waves show power spectral densities with power-law exponents typical for weak capillary-wave turbulence. With increasing Reynolds number, the steepness of the power law declines to lower levels.

The present study investigates the modal stability of the steady incompressible flow inside a toroidal pipe for values of the curvature $\delta$ (ratio between pipe and torus radii) approaching zero, i.e. the limit of a straight pipe. The global neutral stability curve for $10^{-7} \leq \delta \leq ~10^{-2}$ is traced using a continuation algorithm. Two different families of unstable eigenmodes are identified. For curvatures below $1.5 \times 10^{-6}$, the critical Reynolds number ${{Re}}_{cr}$ is proportional to $\delta ^{-1/2}$. Hence, the critical Dean number is constant, ${{De}}_{cr} = 2\,{{Re}}_{cr}\,\sqrt {\delta } \approx 113$. This behaviour confirms that the Hagen–Poiseuille flow is stable to infinitesimal perturbations for any Reynolds number and suggests that a continuous transition from the curved to the straight pipe takes place as far as it regards the stability properties. For low values of the curvature, an approximate self-similar solution for the steady base flow can be obtained at a fixed Dean number. Exploiting the proposed semi-analytic scaling in the stability analysis provides satisfactory results.

Surfactant-like impurities are omnipresent in multiphase emulsions and may substantially affect the motion of small droplets by altering their interfacial properties. Usually these surfactants are soluble in the bulk and undergo adsorption–desorption onto the interface which modifies their surface concentration and hence their overall influence on droplet motion. Yet, the impact of the bulk solubility and transport of surfactants on droplet dynamics, especially in the presence of bounding walls, remains poorly understood. As such, in this article, we assess the impact of bulk soluble surfactants on the settling of a spherical drop towards a plane wall. We consider coupled bulk and interfacial transport of surfactants, mediated by adsorption–desorption processes and construct a semi-analytical framework for arbitrary values of ‘bulk interaction parameter’, which dictates the strength of adsorption–desorption kinetics compared with bulk diffusion. Our results indicate that while mass exchange between the bulk and the interface can remobilize the drop, a finite bulk diffusion rate restricts this process and therefore slows down the drop. This also results in bulk concentration depletion near the south pole and accumulation near the north pole, the extent of which becomes strongly asymmetric with an enhanced intensity of depletion, as the drop approaches the wall. Presence of the wall and bulk solubility are found to aid each other towards remobilizing the drop by aptly modifying the interfacial concentration. Our results may provide fundamental insights into the kinetics of surfactant-laden drops, with potential applications in food and pharmaceutical industries, separation processes, etc.

We study self-similar viscous fingering for the case of divergent flow within a wedge-shaped Hele-Shaw cell. Previous authors have conjectured the existence of a countably infinite number of selected solutions, each distinguished by a different value of the relative finger angle. Interestingly, the associated solution branches have been posited to merge and disappear in pairs as the surface tension decreases. For the first time, we demonstrate how the selection mechanism can be derived based on exponential asymptotics. Asymptotic predictions of the finger-to-wedge angle are additionally given for different sized wedges and surface-tension values. The merging of solution branches is explained; this feature is qualitatively different to the case of classic Saffman–Taylor viscous fingering in a parallel channel configuration. Moreover, because the asymptotic framework does not highly depend on specifics of the wedge geometry, the proposed theory for branch merging in our self-similar problem likely relates much more widely to tip-splitting instabilities in time-dependent flows in circular and other geometries, where the viscous fingers destabilise and divide in two.

Wave turbulence provides a conceptual framework for weakly nonlinear interactions in dispersive media. Dating from five decades ago, applications of wave turbulence theory to oceanic internal waves assigned a leading-order role to interactions characterized by a near equivalence between the group velocity of high-frequency internal waves with the phase velocity of near-inertial waves. This scale-separated interaction leads to a Fokker–Planck (generalized diffusion) equation. More recently, starting four decades ago, this scale-separated paradigm has been investigated using ray tracing methods. These ray methods characterize spectral transport of energy by counting the amplitude and net velocity of wave packets in phase space past a high-wavenumber gate prior to ‘breaking’. This explicitly advective characterization is based on an intuitive assignment and lacks theoretical underpinning. When one takes an estimate of the net spectral drift from the wave turbulence derivation and makes the corresponding assessment, one obtains a prediction of spectral transport that is an order of magnitude larger than either observations or reported ray tracing estimates. Motivated by this contradiction, we report two parallel derivations for transport equations describing the refraction of high-frequency internal waves in a sea of random inertial waves. The first uses standard wave turbulence techniques and the second is an ensemble-averaged packet transport equation characterized by the dispersion of wave packets about a mean drift in the spectral domain. The ensemble-averaged transport equation for ray tracing differs in that it contains the intuitively motivated advective term. We conclude that the aforementioned contradiction between theory, numerics and observations needs to be taken at face value and present a pathway for resolving this contradiction.

A theory is presented for wave-driven propulsion of floating bodies driven into oscillation at the fluid interface. By coupling the equations of motion of the body to a quasipotential flow model of the fluid, we derive expressions for the drift speed and propulsive thrust of the body which in turn are shown to be consistent with global momentum conservation. We explore the efficacy of our model in describing the motion of SurferBot (Rhee et al., Bioinspir. Biomim., vol. 17, issue 5, 2022), demonstrating close agreement with the experimentally determined drift speed and oscillatory dynamics. The efficiency of wave-driven propulsion is then computed as a function of driving oscillation frequency and the forcing location, revealing optimal values for both of these parameters which await confirmation in experiments. A comparison with other modes of locomotion and applications of our model with competitive water sports is discussed in conclusion.