Magneto-gravity-precessional instability, which results from the excitation of resonant magneto-inertia-gravity (MIG) waves by a background shear generic to precessional flows, is addressed here. Two simple background precession flows, that of Kerswell (1993 Geophys. Astrophys. Fluid Dyn. vol 72, no. 1–4, pp. 107–144), and that of Mahalov (1993 Phys. Fluids A: Fluid Dyn. vol. 5, no. 4, pp. 891–900), are considered. We analytically perform an asymptotic analysis to order ${ O}(\varepsilon ),$ where $\varepsilon$ denotes the Poincaré number, i.e. the precession parameter, and determine the maximum growth rate of the destabilizing subharmonic resonances of MIG waves: that between two fast modes, that between two slow modes and that between a fast mode and a slow mode (mixed modes). The domains of the $(K_0 B_0/\varOmega _0, N/\varOmega _0)\hbox{-}$plane for which this instability operates are identified, where $1/K_0$ denotes a characteristic length scale, $B_0$ is the unperturbed Alfvén velocity, $\varOmega _0$ is the rotation rate and $N$ denotes the Brunt–Väisälä frequency. We demonstrate that the $N\rightarrow 0$ limit is, in fact, singular (discontinuous). At large $K_0B_0/\varOmega _0,$ stable stratification acts to suppress the destabilizing resonance between two fast modes as well as that between two slow modes, whereas it revives the destabilizing resonance between a fast mode and a slow mode provided $N\lt \varOmega _0,$ because, without stratification, the maximal growth rate of this instability approaches zero as $K_0B_0/\varOmega _0\rightarrow +\infty .$ This would be relevant for the generation of the mean electromotive force, and hence, the $\alpha \hbox{-}$effect in helical magnetized precessional flows under weak stable stratification. Diffusive effects on the instability is considered in the simple case where the magnetic and thermal Prandtl numbers are both equal to one.

Interfaces subjected to strong time-periodic horizontal accelerations exhibit striking patterns known as frozen waves. In this study, we experimentally and numerically investigate the formation of such structures in immiscible fluids under high-frequency forcing. In the inertial regime – characterised by large Reynolds and Weber numbers, where viscous and surface tension effects become negligible – we demonstrate that the amplitude of frozen waves scales proportionally with the square of the forcing velocity. These results are consistent with vibro-equilibria theory and extend the theoretical framework proposed by Gréa & Briard (2019 Phys. Rev. Fluids 4, 064608) to immiscible fluids with large density contrasts. Furthermore, we examine the influence of both Reynolds and Weber numbers, not only in the onset of secondary Faraday instabilities – which drive the transition of frozen wave patterns toward a homogenised turbulent state – but also in selecting the dominant wavelength in the final saturated regime.

This study considers the global instability of unidirectional flows through single, and double, bifurcation models using linear stability and direct numerical simulation (DNS). The motivation is respiratory flows, so we consider flow in both directions, through two geometries. We identify conditions (quantified by the Reynolds number, ${Re}=U^*D/\nu$, where $U^*$ is the peak centreline velocity, $D$ is the primary pipe diameter and $\nu$ is the kinematic viscosity) where temporal fluctuations occur using DNS. We calculate the linear stability of the steady flows, identifying the critical Reynolds number and leading unstable modes. For flows from single to double pipe, the critical Reynolds number is dependent on the number of bifurcations in the domain, but the mode structures are similar, with growth observed in regions dominated by longitudinal vortices formed by the centrifugal imbalance of flows passing through curved bifurcations. Flows in the opposite direction, from double to single pipe, also depend on the number of bifurcations in the domain. The flow through the double-bifurcation case undergoes two spatial symmetry-breaking bifurcations, altering the mode structure and critical Reynolds number. In all cases, the critical Reynolds number closely matches with temporal fluctuations observed from DNS, suggesting transition is the result of a linear instability, similar to other curved geometries like toroidal and helical pipes. We compare the frequencies of the modes with the frequencies observed from DNS, finding a close match during both initial and saturated flows. These results are important for understanding respiratory flows where turbulent mixing and streaming contribute to gas transport.

Vertically bounded, horizontally propagating internal waves may become unstable through triad resonant instability, in which two sibling waves in background noise draw energy from a parent internal tide. If the background stratification is uniform, then the condition for pure resonance between the parent and sibling wave frequencies and horizontal and vertical wavenumbers can be found semi-analytically from the roots of a polynomial expression. In non-uniform stratification, determining the frequencies and horizontal wavenumbers for which resonance occurs is less straightforward. We develop a theory for near-resonant excitation of a pair of sibling waves from a low-mode internal wave in which the proximity to pure resonance is characterised by the discrepancy between the forced sibling wave frequencies and the natural frequency of these modes. Knowing this discrepancy can be used methodically to determine pure resonance conditions. This inviscid theory is compared with numerical simulations of effectively inviscid waves. For comparison with laboratory experiments, the theory is adapted to include viscous effects both in the bulk of the fluid and at the side walls of the tank. We find that our theoretical predictions for frequencies and wavenumbers of the fastest growing sibling waves are generally consistent between theory, simulations and experiments, though theory overpredicts the growth rate observed in experiments. In all cases, the growth rate of sibling waves decreases with decreasing parent wave frequency, becoming negligibly small in experiments if the parent wave has frequency less than $\approx 0.7$ of the buoyancy frequency at the surface.

Understanding how bubbles on a substrate respond to ultrasound is crucial for applications from industrial cleaning to biomedical treatments. Under ultrasonic excitation, bubbles can undergo shape deformations due to Faraday instability, periodically producing high-speed jets that may cause damage. While recent studies have begun to elucidate this behaviour for free bubbles, the dynamics of wall-attached bubbles is still largely unexplored. In particular, the selection and evolution of non-spherical modes in these bounded systems have not previously been resolved in three dimensions, and the resulting jetting dynamics has yet to be compared with that observed in free bubbles. In this study, we investigate individual micrometric air bubbles in contact with a rigid substrate and subjected to ultrasound. We introduce a novel dual-view imaging technique that combines top-view bright-field microscopy with side-view phase-contrast X-ray imaging, enabling visualisation of bubble shape evolution from two orthogonal perspectives. This set-up reveals the progression of bubble shape through four distinct dynamic regimes: purely spherical oscillations, onset of harmonic axisymmetric meniscus waves, emergence of half-harmonic axisymmetric Faraday waves and the superposition of half-harmonic sectoral Faraday waves. This stepwise evolution contrasts with the behaviour of free bubbles, which exhibit their ultimate Faraday wave pattern immediately upon instability onset. For the substrate chosen, the resulting shape-mode spectrum appears to be degenerate and exhibits a continuous range of shape mode degrees, in line with our theoretical predictions derived from kinematic arguments. While free bubbles also display a degenerate spectrum, their shape mode degrees remain discrete, constrained by the bubble spherical periodicity. Experimentally measured ultrasound pressure thresholds for the onset of Faraday instability agree well with classical interface stability theory, modified to incorporate the effects of a rigid boundary. Complementary three-dimensional boundary element simulations of bubble shape evolution align closely with experimental observations, validating this method’s predictive capability. Finally, we determine the acceleration threshold at which shape mode lobes initiate cyclic jetting. Unlike free bubbles, jetting in wall-attached bubbles consistently emerges from the side not restricted by the substrate.

This study explores the Faraday instability as a mechanism to enhance heat transfer in two-phase systems by exciting interfacial waves through resonance. The approach is particularly applicable to reduced-gravity environments where buoyancy-driven convection is ineffective. A reduced-order model, based on a weighted residual integral boundary layer method, is used to predict interfacial dynamics and heat flux under vertical oscillations with a stabilising thermal gradient. The model employs long-wave and one-way coupling approximations to simplify the governing equations. Linear stability theory informs the oscillation parameters for subsequent nonlinear simulations, which are then qualitatively compared against experiments conducted under Earth’s gravity. Experimental results show up to a 4.5-fold enhancement in heat transfer over pure conduction. Key findings include: (i) reduced gravity lowers interfacial stability, promoting mixing and heat transfer; and (ii) oscillation-induced instability significantly improves heat transport under Earth’s gravity. Theoretical predictions qualitatively validate experimental trends in wavelength-dependent enhancement of heat transfer. Quantitative discrepancies between model and experiment are rationalised by model assumptions, such as neglecting higher-order inertial terms, idealised boundary conditions, and simplified interface dynamics. These limitations lead to underprediction of interface deflection and heat flux. Nevertheless, the study underscores the value of Faraday instability as a means to boost heat transfer in reduced gravity, with implications for thermal management in space applications.

A generalised multiparameter model for linear modal stability and sensitivity analysis is developed. The stability and sensitivity equations are derived from a generalised vector-form governing equation comprised of multiple dimensionless parameters that represent different physical forces affecting the system’s stability. By introducing adjoint variables and constructing the Lagrangian identity, a differential relationship between the eigenvalue of the perturbation mode and dimensionless parameters is determined and defined as the global sensitivity gradient. It provides the constraint that must be satisfied for changes in different dimensionless parameters along the isoeigenvalue curve, which aids in the fast computation of the neutral curve. Moreover, the global sensitivity gradient can directly and intuitively evaluate the competitive relationship among the influences of various parameters on system instability. Based on the global sensitivity gradient, an optimal stability control strategy for transitioning from an unstable state to a stable state is discussed. Additionally, the relative sensitivity function is also introduced to investigate the influence of relative parameter variations on instability. To demonstrate the effectiveness of this method, three applications are presented: two-dimensional flow around a circular cylinder with a single dimensionless parameter Re; three-dimensional axisymmetric magnetohydrodynamic (MHD) flow around a sphere with two parameters Re and $N$; and two-dimensional MHD mixed convection with three parameters Re, ${\textit{Gr}}$ and $\textit{Ha}$.

Parametric oscillations of an interface separating two fluid phases create nonlinear surface waves, called Faraday waves, which organise into simple patterns, such as squares and hexagons, as well as complex structures, such as double hexagonal and superlattice patterns. In this work, we study the influence of surfactant-induced Marangoni stresses on the formation and transition of Faraday-wave patterns. We use a control parameter, $B$, that assesses the relative importance of Marangoni stresses as compared with the surface-wave dynamics. Our results show that the threshold acceleration required to destabilise a surfactant-covered interface through vibration increases with increasing $B$. For a surfactant-free interface, a square-wave pattern is observed. As $B$ is incremented, we report transitions from squares to asymmetric squares, weakly wavy stripes and ultimately to ridges and hills. These hills are a consequence of the bidirectional Marangoni stresses at the neck of the ridges. The mechanisms underlying the pattern transitions and the formation of exotic ridges and hills are discussed.

Understanding wave kinematics is crucial for analysing the thermodynamic effects of sloshing, which can lead to pressure drops in non-isothermal cryogenic fuel tanks. In the research reported here, Faraday waves in a horizontal circular tank (partially filled with water) under vertical excitation are investigated. The tank geometry is referred to as a horizontal circular tank throughout, with its circular face oriented perpendicular to the horizontal plane. Firstly, this paper addresses the eigenvalue problem through linear potential flow theory, in order to provide theoretical evidence of Faraday waves in horizontal circular tanks, the impact the density ratio has on the eigenvalues is then considered. Secondly, an experimental investigation testing multiple liquid fill levels is conducted. A soft-spring nonlinear response is demonstrated throughout the parameter space. The results showed larger sloshing amplitudes for low fill levels and smaller sloshing amplitudes for high fill levels. Asymmetry between anti-nodes at the container sidewalls and through the tank centreline are evident for low fill levels. Moreover, the sloshing wave amplitude at which breaking waves occur is smaller for high fill level conditions. Finally, period tripling was observed for all fill levels tested, confirming nonlinear mode interactions before the onset to wave breaking.

We study experimentally the onset of Faraday waves near the end walls of a rectangular vessel containing two stably stratified fluid layers, subject to horizontal oscillations. These subharmonic waves (SWs) are excited, because the horizontal inertial forcing drives a harmonic propagating wave which displaces the interface in the vertical direction at the end walls. We find that the onset of SWs is regulated by a balance between capillary and viscous forces, where the rate of damping is set by the Stokes layer thickness at the wall rather than the wavelength of the SWs. We model the onset of SWs with a weakly damped Mathieu equation and find that the dimensional critical acceleration scales as $\nu _m^{1/2} \omega ^{3/2}$, where $\nu _m$ is the mean viscosity and $\omega$ is the frequency of forcing, in excellent agreement with the experiment over a wide range of parameters.

We report the dynamics of a droplet levitated in a multi-emitter, single-axis acoustic levitator. The deformation and atomisation behaviour of the droplet in the acoustic field displays myriad complex phenomena, in a series of events. These include the primary breakup of the droplet, wherein it exhibits stable levitation, deformation, sheet formation and equatorial atomisation, followed by its secondary breakup, which could be of various types such as umbrella, bag, bubble or multi-stage breakup. A large number of tiny atomised droplets, formed as a result of the primary and secondary breakup, can remain levitated in the acoustic levitator and exhibit aggregation and coalescence. The visualisation of the interfacial instabilities on the surface of the liquid sheet using both side- and top-view imaging is presented. An approximate size distribution of the atomised droplets is also provided. The stable levitation of the droplet is due to a balance of acoustic and gravitational forces while the resulting ellipsoidal shape of the droplet is a consequence of the balance of the deforming acoustic force and the restoring capillary force. Stronger acoustic forces can no longer be balanced by capillary forces, resulting in a highly flattened droplet, with a thin liquid sheet at the edge (equatorial region). The thinning of the sheet is caused by the differential acceleration induced by the increasing pressure difference between the poles and the equator as the sheet deforms. When the sheet thickness reduces to of the order of a few microns, Faraday waves develop at the thinnest region (preceding the rim), which causes the generation of tiny-sized droplets that are ejected perpendicular to the sheet. The corresponding hole formation results in a perforated sheet that causes the detachment of the annular rim, which breaks due to Rayleigh–Plateau (RP) instability. The radial ligaments generated in the sheet, possibly due to Rayleigh–Taylor (RT) instability, break into droplets of different sizes. The secondary breakup exhibits Weber number dependency and includes umbrella, bag, bubble or multi-stage types, ultimately resulting in the complete atomisation of the droplets. Both the primary and the secondary breakup of the droplet involve interfacial instabilities such as Faraday, Kelvin–Helmholtz, RT and RP and are well supported by visual evidence.

Waves are formed on the surface of a sessile drop driven through substrate vibrations oriented at a slanting angle from the normal. A mathematical model is derived, which leads to an infinite system of coupled Mathieu equations governing the wave dynamics that are solved using Floquet theory. The spatial structure of the waves is described by the mode number pair $[\ell,m]$, where $\ell$ and $m$ are the polar and azimuthal mode numbers, respectively. Limiting cases corresponding to horizontal and vertical vibrations are discussed with predictions agreeing well with prior literature. We focus our results on three drop motions – (1) harmonic $[1,1]$ rocking mode, (2) harmonic $[2,0]$ pumping mode, and (3) subharmonic rocking $[1,1]$ mode – as they depend upon the slanting angle, static contact angle, and contact-line conditions, which we assume to be either pinned or freely moving with fixed contact angle. New theoretical predictions are tested through experiments over a range of parameters, showing good agreement.

Existing theoretical analyses of Faraday waves in Hele-Shaw cells rely on the Darcy approximation and assume a parabolic flow profile in the narrow direction. However, Darcy's model is known to be inaccurate when convective or unsteady inertial effects are important. In this work, we propose a gap-averaged Floquet theory accounting for inertial effects induced by the unsteady terms in the Navier–Stokes equations, a scenario that corresponds to a pulsatile flow where the fluid motion reduces to a two-dimensional oscillating Poiseuille flow, similarly to the Womersley flow in arteries. When gap-averaging the linearised Navier–Stokes equation, this results in a modified damping coefficient, which is a function of the ratio between the Stokes boundary layer thickness and the cell's gap, and whose complex value depends on the frequency of the wave response specific to each unstable parametric region. We first revisit the standard case of horizontally infinite rectangular Hele-Shaw cells by also accounting for a dynamic contact angle model. A comparison with existing experiments shows the predictive improvement brought by the present theory and points out how the standard gap-averaged model often underestimates the Faraday threshold. The analysis is then extended to the less conventional case of thin annuli. A series of dedicated experiments for this configuration highlights how Darcy's thin-gap approximation overlooks a frequency detuning that is essential to correctly predict the locations of the Faraday tongues in the frequency–amplitude parameter plane. These findings are well rationalised and captured by the present model.

This paper investigates the steady-state pattern evolution of symmetric Faraday waves excited in a brimful cylindrical container when driving parameters much exceed critical thresholds. In such liquid systems, parametric surface responses are typically considered as the resonant superposition of unstable standing waves. A modified free-surface synthetic Schlieren method is employed to obtain full three-dimensional spatial reconstructions of instantaneous surface patterns. Multi-azimuth structures and localized travelling waves during the small-elevation phases of the oscillation cycle give rise to modal decomposition in the form of $\nu$-basis modes. Two-step surface-fitting results provide insight into the spatiotemporal characteristics of dominant wave components and corresponding harmonics in the experimental observations. Arithmetic combination of modal indices and uniform frequency distributions reveal the nonlinear mechanisms behind pattern formation and the primary pathways of energy transfer. Taking the hypothetical surface manifestation of multiple azimuths as the modal solutions, a linear stability analysis of the inviscid system is utilised to calculate fundamental resonance tongues (FRTs) with non-overlapping bottoms, which correspond to subharmonic or harmonic $\nu$-basis modes induced by surface instability at the air–liquid interface. Close relationships between experimental observations and corresponding FRTs provide qualitative verification of dominant modes identified using surface-fitting results. This supports the validity and rationality of the applied $\nu$-basis modes.

The flow inside a rotating annulus tilted with respect to gravity is characterized experimentally and theoretically. As in the case of a tilted rotating cylinder the flow is forced by the free surface, maintained flat by gravity. It leads to resonances of global inertial modes (Kelvin modes) when the height of fluid is a multiple of half the wavelength of the mode. The divergence of the mode is saturated by viscous effects at the resonance. The maximum amplitude scales as the Ekman number to the power $-1/2$ when surface Ekman pumping is dominant, and to the power $-1$ when volumic damping is dominant. An analytical prediction is given with no fitting parameter, in excellent agreement with experimental results. At lower Ekman numbers, the flow destabilizes with respect to a triadic resonance instability, as already observed by Xu & Harlander (Phys. Rev. Fluids, 2020). We provide here a linear stability analysis leading to the viscous threshold of the instability for small tilt angles. For large tilt angles, a centrifugal instability is observed due to the acceleration of the flow by the inner cylinder. Finally, the features of the turbulent flow and its mixing efficiency are characterized experimentally. We underline the potential interest of this configuration for bioreactors.

We investigate the influence of rotation on the onset and saturation of the Faraday instability in a vertically oscillating two-layer miscible fluid using a theoretical model and direct numerical simulations (DNS). Our analytical approach utilizes Floquet analysis to solve a set of the Mathieu equations obtained from the linear stability analysis. The solution of the Mathieu equations comprises stable and harmonic, and subharmonic unstable regions in a three-dimensional stability diagram. We find that the Coriolis force delays the onset of the subharmonic instability responsible for the growth of the mixing zone size at lower forcing amplitudes. However, at higher forcing amplitudes, the flow is energetic enough to mitigate the instability delaying effect of rotation, and the evolution of the mixing zone size is similar in both rotating and non-rotating environments. These results are corroborated by DNS at different Coriolis frequencies and forcing amplitudes. We also observe that for $(\,f/\omega )^2<0.25$, where $f$ is the Coriolis frequency, and $\omega$ is the forcing frequency, the instability and the turbulent mixing zone size-$L$ saturates. When $(\,f/\omega )^2\geq 0.25$, the turbulent mixing zone size-$L$ never saturates and continues to grow.

The response to harmonic horizontal oscillations of a stably stratified fluid-filled two-dimensional square container is examined as the forcing amplitude is increased. For the studied forcing frequency, the response flow at very small forcing amplitudes is a synchronous periodic flow with piecewise-constant vorticity in regions delineated by the characteristics emanating from the corners of the container, regularized by viscosity. The second temporal harmonic of the forced response flow resonantly excites an intrinsic mode of the stratified container, whose magnitude grows as the square of the forcing amplitude. Above a critical forcing amplitude, a sequence of pairs of other container modes are excited via triadic resonances with the second-harmonic-driven mode. The flows are computed from the Navier–Stokes–Boussinesq equations and the ensuing dynamics is analysed using Fourier techniques, providing a comprehensive picture of the transition to internal wave turbulence.

We study experimentally the dynamics of a water droplet on a tilted and vertically oscillating rigid fibre. As we vary the frequency and amplitude of the oscillations the droplet transitions between different modes: harmonic pumping, subharmonic pumping, a combination of rocking and pumping modes, and a combination of pumping and swinging modes. We characterize these responses and report how they affect the sliding speed of the droplet along the fibre. The swinging mode is explained by a minimal model making an analogy between the droplet and a forced elastic pendulum.

The two-dimensional stability of vertically sheared inertial oscillations at ocean fronts is explored through a linear stability analysis and nonlinear simulations. Baroclinic effects reduce the minimum frequency of inertia-gravity waves to an extent determined by the balanced Richardson number ${{Ri}}$ of the front. Below a critical value of ${{Ri}}$, which depends on the strength of the inertial shear, the inertial oscillations become unstable to parametric subharmonic instability (PSI) resulting in growing perturbations that oscillate at half the inertial frequency $f$. Since the critical value is always greater than 1, PSI can occur at fronts stable to symmetric instability. Although modest in weak inertial shear, growth rates exceeding $f/2$ can be achieved for inertial shear greater than or equal to the thermal wind shear. Our formulation allows for non-hydrostatic perturbations and can be applied to initially unstratified geostrophic adjustment problems. We find that PSI will almost totally damp the transient oscillations that arise during geostrophic adjustment. The perturbations gain energy at the expense of the inertial oscillations through ageostrophic shear production. The perturbations then themselves become unstable to secondary Kelvin–Helmholtz instabilities creating a pathway by which the inertial oscillations can be dissipated rapidly. In contrast to symmetric and baroclinic instabilities that draw on a front's kinetic or potential energy, PSI acts to increase the energy stored in the balanced front as the convergence and divergence of the eddy-momentum fluxes set up a secondary circulation in the sense to stand up the front.

Water stuck in the ear is a common problem during showering, swimming or other water activities. Having water trapped in the ear canal for a long time can lead to ear infections and possibly result in hearing loss. A common strategy for emptying water from the ear canal is to shake the head, where high acceleration helps remove the water. In this present study, we rationalize the underlying mechanism of water ejection/removal from the ear canal by performing experiments and developing a stability theory. From the experiments, we measure the critical acceleration to remove the trapped water inside different sizes of canals. Our theoretical model, modified from the Rayleigh–Taylor instability, can explain the critical acceleration observed in experiments, which strongly depends on the radius of the ear canal. The resulting critical acceleration tends to increase, especially in smaller ear canals, which indicates that shaking heads for water removal can be more laborious and potentially threatening to children due to their small size of the ear canal compared with adults.