In this work, we study the effect of flow curvature, or angular momentum, on the propagation and trapping characteristics of near-inertial waves (NIWs) in a curved front. The curved front is idealised as a baroclinic vortex in cyclogeostrophic balance. Motivated by ocean observations, we employ a Gaussian base flow, which by construction possesses a shield of oppositely signed vorticity surrounding its core, and we consider both cyclonic and anticyclonic representations of this flow. Following two main assumptions, i.e. that (i) the horizontal wavelength of the NIW is smaller than the length scale of the background flow (the WKBJ approximation), and (ii) the vertical wavelength of the NIW is smaller than the radial distance of interest, we derive the NIW dispersion relation and discuss the group velocity and direction of energy propagation. We show that the curvature can (i) increase the critical depth and horizontal extent of the trapping region, (ii) reduce NIW activity at the centre of the anticyclonic vortex core and enhance it in the cyclonic shield surrounding the core for high curvatures, (iii) lead to NIW trapping in the anticyclonic shield surrounding the cyclonic core, and (iv) increase the available band of NIW frequencies that are trapped. The solutions from the ray-tracing method are supported by numerical solutions of the governing equations linearised about the cyclogeostrophic base state. Finally, these methods are applied to an idealised model of oceanic mesoscale Arctic eddies showing an increase in the critical depth of trapping. Our results – while applied to polar eddies – equally apply at lower latitudes in both oceans and atmospheres, highlighting the potential importance of flow curvature in controlling the propagation of NIW energy.

The aerodynamic sound generated by the oblique collision of two vortex rings is featured by the asymmetric emission associated with the octupole mode, which differs from the symmetric emission associated with the quadrupole mode observed in the coaxial collision of two vortex rings. This distinctive feature of aerodynamic sound is closely related to the tilting and reconnecting of the vortex rings. While previous studies have explored the effects of reconnecting on aerodynamic sound, this study specifically addresses the impact of vortex ring tilting. We propose a novel vortex sound formula to quantitatively assess the role of tilting in aerodynamic sound generation. The proposed formula relates the far-field sound pressure to equivalent circulations and vorticity centroids by referring to Truesdell’s consistency conditions for vorticity moments. The variations of the equivalent circulations and vorticity centroids in the oblique collision of two vortex rings under different configurations are analysed based on the numerical solution of the Navier–Stokes equations in the source region. It is found that the tilting of vortex rings results in a rapid change of the equivalent circulation associated with the vorticity in the collision direction. However, the change caused by titling is almost out of phase with that caused by reconnecting and deforming. The vortex tilting significantly reduces the aerodynamic sound associated with the longitudinal quadrupole and octupole modes, which is opposite to the role of vortex reconnecting that was reported in the oblique collision of vortex rings.

Inspired by the need to theoretically understand the naturally occurring interactions between internal waves and mesoscale phenomena in the ocean, we derive a novel model equation from the primitive rotational Euler equations using the multi-scale asymptotic expansion method. By applying the classic balance $\epsilon =\mu ^2$ between nonlinearity (measured by $\epsilon$) and dispersion (measured by $\mu$), along with the assumption that variations in the transverse direction are of order $\mu$, which is smaller than those in the propagation direction, we arrive at terms from the classic Kadomtsev–Petviashvili equation. However, when incorporating background shear currents in two horizontal dimensions and accounting for Earth’s rotation, we introduce three additional terms that, to the best of the authors’ knowledge, have not been addressed in the previous literature. Theoretical analyses and numerical results indicate that these three terms contribute to a tendency for propagation in the transverse direction and an overall variation in wave amplitudes. The specific effects of these terms can be estimated qualitatively based on the signs of the coefficients for each term and the characteristics of the initial waves. Finally, the potential shortcomings of this proposed equation are illuminated.

Direct numerical simulations are performed to explore the impact of surface roughness on inter-scale energy transfer and interaction in a turbulent open-channel flow over differently arranged rough walls. With friction Reynolds number approximately 540, six distinct configurations of roughness arrangements are examined. The results show that the clustered roughness arrangements yield notable changes in large-scale secondary-flow structures, which manifest in the profiles of dispersive stresses, predominantly near the roughness elements. They are marked by the presence of spanwise alternating high-momentum pathways and low-momentum pathways. From the outer peak in the spanwise energy spectra, the size and intensity of turbulent secondary flows are shown to be related to the spanwise spacing of the roughness heterogeneity. The emergence of turbulent secondary flows serves to suppress the original large-scale structures in the outer region of smooth-wall turbulence, paving the way for the development of new turbulent structures at the second harmonic scale. Furthermore, the spanwise triadic interaction analysis reveals the mutual energy exchange between the secondary harmonic scale and the secondary-flow scale. These findings elucidate the underlying mechanisms behind the attenuation of large-scale structures in the outer region influenced by roughness, offering new insights into the dynamic interplay of scale interactions in rough-wall turbulence.

Synthetic-aperture radar images and mesoscale models show that wind-farm wakes differ from single-turbine wakes. For instance, wind-farm wakes often narrow and do not disperse over long distances, contrasting the broader and more dissipating wakes of individual turbines. In this work, we aim to better understand the mechanisms that govern wind-farm wake behaviour and recovery. Hence we study the wake properties of a $1.6$ GW wind farm operating in conventionally neutral boundary layers with capping-inversion heights $203$, $319$, $507$ and $1001$ m. In shallow boundary layers, we find strong flow decelerations that reduce the Coriolis force magnitude, leading to an anticlockwise wake deflection in the Northern Hemisphere. In deep boundary layers, the vertical turbulent entrainment of momentum adds clockwise-turning flow from aloft into the wake region, leading to a faster recovery rate and a clockwise wake deflection. To estimate the wake properties, we propose a simple function to fit the velocity magnitude profiles along the spanwise direction. In the vertical direction, the wake spreads up to the capping-inversion height, which significantly limits vertical wake development in shallow-boundary-layer cases. In the horizontal direction and for shallow boundary layers, the wake behaves as two distinct mixing layers located at the lateral wake edges, which expand and turn towards their low-velocity side, causing the wake to narrow along the streamwise direction. A detailed analysis of the momentum budget reveals that in deep boundary layers, the wake is predominantly replenished through turbulent vertical entrainment. Conversely, in shallow boundary layers, wakes are mostly replenished by mean flow advection in the spanwise direction.

Bubble–particle collisions in turbulence are key to the froth flotation process that is widely employed industrially to separate hydrophobic from hydrophilic materials. In our previous study (Chan et al., 2023 J. Fluid Mech. 959, A6), we elucidated the collision mechanisms and critically reviewed the collision models in the no-gravity limit. In reality, gravity may play a role since, ultimately, separation is achieved through buoyancy-induced rising of the bubbles. This effect has been included in several collision models, which have remained without a proper validation thus far due to a scarcity of available data. We therefore conduct direct numerical simulations of bubbles and particles in homogeneous isotropic turbulence with various Stokes, Froude and Reynolds numbers, and particle density ratios using the point-particle approximation. Generally, turbulence enhances the collision rate compared with the pure relative settling case by increasing the collision velocity. Surprisingly, however, for certain parameters the collision rate is lower with turbulence compared with without, independent of the history force. This is due to turbulence-induced bubble–particle spatial segregation, which is most prevalent at weak relative gravity and decreases as gravitational effects become more dominant, and reduced bubble slip velocity in turbulence. The existing bubble–particle collision models only qualitatively capture the trends in our numerical data. To improve on this, we extend the model by Dodin & Elperin (2002 Phys. Fluids 14, 2921–2924) to the bubble–particle case and found excellent quantitative agreement for small Stokes numbers when the history force is negligible and segregation is accounted for.

The connection between the drag and vorticity dynamics for viscous flow over a bluff body is explored using the Josephson–Anderson (J–A) relation for classical fluids. The instantaneous rate of work on the fluid, associated with the drag force, is related to the vorticity flux across the streamlines of a background potential flow. The vorticity transport itself is examined by aid of the Huggins vorticity-flux tensor. The analysis is performed for three flows: flow over a sphere at Reynolds numbers $Re=200,3700$, and flow over a prolate spheroid at $Re=3000$ and $20^{\circ }$ incidence. In these flows, the vorticity transport shifts the flow away from and towards the ideal potential flow, with a net balance towards the former effect thus making an appreciable contribution to the drag. The J–A relation is first demonstrated for the flow over a sphere at $Re=200$. The drag power injection is related to the viscous flux of azimuthal vorticity from the wall into the fluid, and the advection of vorticity by the detached shear layer. In the wake, the azimuthal vorticity is advected towards the wake centreline and is annihilated by viscous effects, which contributes a reduction in drag. The analysis of the flow over a sphere at $Re=3700$ is reported for the impulsively started and stationary stages, with emphasis on the effects of unsteady two-dimensional separation and turbulent transport in the transitional wake. The turbulent flux in the wake enhances the transport of mean azimuthal vorticity towards the wake centreline, and is the main driver of the recovery of enthalpy between the rear point of the sphere and far downstream. The rate of work on the fluid by the drag force for a prolate spheroid is mostly due to the transport of vorticity along the separated boundary layers. Primary and secondary separation contribute oppositely to the power injection by the drag force, while the large-scale vortices only re-distribute vorticity without a net contribution. A mechanism for secondary separation is proposed based on the theory of vortex-induced separation.

Hydrodynamic modulation of short ocean surface waves by longer ambient waves significantly influences remote sensing, interpretation of in situ wave measurements and numerical wave forecasting. This paper revisits the wave crest and action conservation laws and derives steady, nonlinear, analytical solutions for the change of short-wave wavenumber, action and gravitational acceleration due to the presence of longer waves. We validate the analytical solutions with numerical solutions of the full crest and action conservation equations. The nonlinear analytical solutions of short-wave wavenumber, amplitude and steepness modulation significantly deviate from the linear analytical solutions of Longuet-Higgins & Stewart (1960 J. Fluid Mech. vol. 8, no. 4, pp. 565–583) and are similar to the nonlinear numerical solutions by Longuet-Higgins (1987 J. Fluid Mech. vol. 177, pp. 293–306) and Zhang & Melville (1990 J. Fluid Mech. vol. 214, pp. 321–346). The short-wave steepness modulation is attributed 5/8 to wavenumber, 1/4 due to wave action and 1/8 due to effective gravity. Examining the homogeneity and stationarity requirements for the conservation of wave action reveals that stationarity is a stronger requirement and is generally not satisfied for very steep long waves. We examine the results of Peureux et al. (2021 J. Geophys. Res.: Oceans vol. 126, no. 1, e2020JC016735) who found through numerical simulations that the short-wave modulation grows unsteadily with each long-wave passage. We show that this unsteady growth only occurs for homogeneous initial conditions as a special case and not generally. The proposed steady solutions are a good approximation of the nonlinear crest-action conservation solutions in long-wave steepness $\lesssim 0.2$. Except for a subset of initial conditions, the solutions to the nonlinearised crest-action conservation equations are mostly steady in the reference frame of the long waves.

In this work, we numerically investigate heat transfer in low-Prandtl-number drop-laden wall-bounded turbulence. These flows are characteristic of nuclear and fusion technologies, where liquid metals – known for their high thermal conductivity – are laden with drops or bubbles of another liquid or pressurised gas. To this end, we consider forced convection turbulence between two differentially heated parallel plates. The carrier phase (i.e. liquid metal) is characterised by a low Prandtl number $Pr_c=0.013$, while for the dispersed phase, we explore a range of Prandtl numbers from $Pr_d=0.013$ (matched case) to $Pr_d=7$ (super-unitary Prandtl number in the dispersed phase). Simulations are conducted at constant friction Reynolds number $Re_\tau =300$, and for each dispersed phase Prandtl number, two volume fractions are examined: $\alpha =5.4\,\%$ and $\alpha =10.6\,\%$. The simulation framework relies on direct numerical simulation of the Navier–Stokes equations, coupled with a phase-field method and the energy equation. Results show that an increase of the dispersed phase Prandtl number reduces heat transfer, leading to a lower Nusselt number for both volume fractions. To explain this behaviour, we analyse how the drops modify the temperature field, and demonstrate that the heat transfer reduction stems from a decreased diffusive heat flux within the dispersed phase. Finally, we propose a phenomenological model to predict the Nusselt number as a function of both the dispersed phase volume fraction and Prandtl number.

Viscoplastic fluids exhibit yield stress, beyond which they flow viscously, while at lower stress levels they behave as solids. Despite their fundamental biological and medical importance, the hydrodynamics of swimming in viscoplastic environments is still evolving. In this study, we investigate the swimming of an ellipsoidal squirmer and the associated tracer diffusion in a Bingham viscoplastic fluid. The results illustrate that neutral squirmers in viscoplastic fluids experience a reduction in swimming speed and an increase in power dissipation as the Bingham number increases, with swimming efficiency peaking at moderate Bingham numbers. As the aspect ratio of a squirmer increases, ellipsoidal squirmers exhibit significantly higher swimming speeds in viscoplastic fluids. The polar and swirling modes can either enhance or reduce swimming speed, depending on the specific scenarios. These outcomes are closely related to the confinement effects induced by the yield surface surrounding the swimmer, highlighting how both swimmer shape and swimming mode can significantly alter the yield surface and, in turn, modify the swimming hydrodynamics. In addition, this study investigates the influence of viscoplasticity on swimmer-induced diffusion in a dilute suspension. The plasticity enforces the velocity far from the swimmer to be zero, thus breaking the assumptions used in Newtonian fluids. The diffusivity reaches its maximum at intermediate aspect ratios and Bingham numbers, and increases with the magnitude of the squirmer’s dipolarity. These findings are important to understand microscale swimming in viscoplastic environments and the suspension properties.

Bubbles entrained by breaking waves rise to the ocean surface, where they cluster before bursting and release droplets into the atmosphere. The ejected drops and dry aerosol particles, left behind after the liquid drop evaporates, affect the radiative balance of the atmosphere and can act as cloud condensation nuclei. The remaining uncertainties surrounding the sea spray emissions function motivate controlled laboratory experiments that directly measure and link collective bursting bubbles and the associated drops and sea salt aerosols. We perform experiments in artificial seawater for a wide range of bubble size distributions, measuring both bulk and surface bubble distributions (measured radii from $30\,\unicode{x03BC} \mathrm{m}$ to $5\,\mathrm{mm}$), together with the associated drop size distribution (salt aerosols and drops of measured radii from $50\,\mathrm{nm}$ to $500\,\unicode{x03BC} \mathrm{m}$) to quantify the link between emitted drops and bursting surface bubbles. We evaluate how well the individual bubble bursting scaling laws describe our data across all scales and demonstrate that the measured drop production by collective bubble bursting can be represented by a single framework integrating individual bursting scaling laws over the various bubble sizes present in our experiments. We show that film drop production by bubbles between $100\,\unicode{x03BC} \mathrm{m}$ and $1\,\mathrm{mm}$ describes the submicron drop production, while jet drop production by bubbles from $30\,\unicode{x03BC} \mathrm{m}$ to $2\,\mathrm{mm}$ describes the production of drops larger than $1\,\unicode{x03BC} \mathrm{m}$. Our work confirms that sea spray emission functions based on individual bursting processes are reasonably accurate as long as the surface bursting bubble size distribution is known.

Sound entering the ear is known not only to transmit signals to the nerve system, but also to generate vortex-like steady streaming in the cochlea. This streaming has been suggested as the primary vehicle for drug delivery in the inner ear (Sumner, Mestel & Reichenbach, 2021, Sci. Rep., vol. 11, 57). An alternative vehicle by pure diffusion alone has also been suggested by Sadreev et al. (2019, Front. Cell. Neurosci., vol. 13, 161). This paper purports to examine both mechanisms analytically, and compare their relative importance, based on the two-dimensional model of Allen (1977, Acoust. Soc. Am., vol. 61, 110–119). First, we reconstruct the fluid mechanics of the Békséy vortices by an asymptotic theory of multiple scales as a complement to the two-dimensional numerical theory of Edom, Obrist & Kleiser (2014, J. Fluid Mech., vol. 753, 254–278). For discerning the difference between Sumner, Mestel & Reichenbach (2021) and Sadreev et al. (2019), we combine sound-induced streaming and molecular diffusion by modeling the drug as a solute of known diffusivity. It will be shown that due to the high frequency of sound, advection is augmented by the Lagrangian velocity, but molecular diffusion still dominates drug transport in the cochlear duct, unlike Taylor dispersion of pollutant by tides in a shallow river.

The Leidenfrost effect occurs when drops are deposited on a highly superheated solid surface, creating a thin vapour film through rapid evaporation that levitates the drops. For drop with a radius exceeding a critical value, a vapour bubble forms and bursts from its bottom centre, a phenomenon known as chimney instability. Despite extensive investigation, the impact of Leidenfrost drop’s rotation on its chimney instability has remained unexplored. This study addresses this gap by providing both numerical and approximate solutions to the theoretical models. We identify two distinct regimes where either gravitational force or centrifugal force is the primary driver of chimney instability. These regimes are characterised by a non-dimensional rotation number, Ro, which represents the ratio between centrifugal force and gravitational force. Our findings reveal clear scaling laws that relate the critical geometrical parameters (radius, volume and height of the drop) for chimney instability to Ro, demonstrating that rotation can induce chimney instability in smaller drops. The scaling laws are elucidated through pressure perturbation analyses under a virtual perturbation to the drop profile at the critical state for chimney instability. Additionally, by varying the evaporation number $Ev$, we demonstrate that while increased superheat reduces the critical radius in the absence of rotation, the scaling laws related to Ro for a rotating drop remain unaffected. Building on these insights, we present a master curve in a simplified form that accurately predicts the critical state for chimney instability under various angular velocities, gravitational accelerations and superheat conditions.

A knowledge gap exists for flows and transport phenomena at the angstrom scale when the Poisson–Nernst–Planck equation based on the concept of the electrical double layer fails. We discovered that streaming conductance becomes pressure-dependent in angstrom channels using latent-track membranes. The streaming current emerges only when the applied pressure exceeds a threshold value, which is inconsistent with the existing knowledge as a constant. With increasing channel size, we found that the pressure-dependent streaming conductance phenomenon weakens and vanishes into a constant streaming conductance regime when the mean channel radius exceeds $\sim$2 nm. The effective surface potential derived from the streaming conductance that divides conduction anomalously increases as the channel narrows. We suspect that the pressure-dependent streaming current is due to the reinforced Coulomb interaction between counterions and deprotonated carboxyl groups at the surface, which is close to the ion channel but different from that of electrified two-dimensional materials. The streaming current emerged due to hydrodynamic friction when the counterions were released from the surface. We approximated the stochastic process of counterion dissociation by a one-dimensional Kramer escape theory framework and defined the Damk$\ddot {\mathrm{o}}$hler number to describe the transition from the nonlinear streaming conductance regime to the linear regime as functions of applied pressure and channel radius and well explained the enhanced effective surface potential in confinement.

Phase change materials (PCMs) hold considerable promise for thermal energy storage applications. However, designing a PCM system to meet a specific performance presents a formidable challenge, given the intricate influence of multiple factors on the performance. To address this challenge, we hereby develop a theoretical framework that elucidates the melting process of PCMs. By integrating stability analysis with theoretical modelling, we derive a transition criterion to demarcate different melting regimes, and subsequently formulate the melting curve that uniquely characterises the performance of an exemplary PCM system. This theoretical melting curve captures the key trends observed in experimental and numerical data across a broad parameter space, establishing a convenient and quantitative relationship between design parameters and system performance. Furthermore, we demonstrate the versatility of the theoretical framework across diverse configurations. Overall, our findings deepen the understanding of thermo-hydrodynamics in melting PCMs, thereby facilitating the evaluation, design and enhancement of PCM systems.

A numerical study is presented on flow-induced vibration of a circular cylinder, under the effect of a downstream stationary cylinder-induced proximity interference. The interference-induced various types of gap-flow regimes and characteristics of vibration and gap-flow rate $Q^*_g$ are presented, by considering various non-dimensional gaps $G^* = 0.1{-}2.5$ and reduced velocities $U^* = 3{-}20$ at a constant Reynolds number $Re = 100$, mass ratio $m^*= 2$ and damping ratio $\zeta = 0.005$. Decreasing $G^*$ or increasing proximity leads to the four gap-flow regimes: bi-directional gap flow at $G^* \geqslant 1.0$, uni-directional non-orthogonal gap flow at $G^* = 1.5{-}1.0$, uni-directional orthogonal gap flow at $G^* \leqslant 0.5$ and uni-directional one-sided gap flow at $G^* \leqslant 0.3$. Further, the respective regimes at larger $U^*$ are associated with proximity-induced modified vortex-induced vibration (PImVIV), proximity-induced galloping (PIG), transitional PImVIV–PIG, and proximity-induced staggered vibration (PISV). Quantitative presentation of maximum gap-flow rate $Q^*_{{g,max}}$, phase $ \phi _g$ (between $Q^*_{g}$ and displacement $y^*$) and phase portraits ($Q^*_{g}$ versus $y^*$) provides clear demarcation between the various gap-flow regimes. Flow mechanisms are presented for the PImVIV, PIG and PISV responses. For the PIG, the mechanism is presented for the first time on generation of galloping instability, asymptotically increasing $A^*$ and existence of optimum gap $G^* = 0.5$ for the maximum amplitude. This work is significant as it provides new insights into the proximity interference-induced gap-flow dynamics between two cylinders, associated flow mechanism for both vibration mitigation and enhancement and promising potential applications for energy harvesting.

We investigate enstrophy variations by collapse of point vortices in an inviscid flow and, in particular, focus on the enstrophy dissipation that is a significant property characterising two-dimensional (2-D) turbulent flows. To reveal the vortex dynamics causing the enstrophy dissipation, we consider the dynamics of point vortices, whose vorticity is concentrated on points and dynamics on the inviscid flow, governed by the point-vortex system. The point-vortex system has self-similar collapsing solutions, which are expected to be a key to understand the enstrophy dissipation, but the collapsing process cannot be described by solutions to the 2-D Euler equations. We thus consider the 2-D filtered-Euler equations, which are a regularised model of the 2-D Euler equations, and their point-vortex solutions. The preceding studies (Gotoda and Sakajo, J. Nonlinear Sci. 2016, vol. 26, pp. 1525–1570, Gotoda and Sakajo, SIAM J. Appl. Math. 2018, vol. 78, 2105–2128) have proven that there exist three point-vortex solutions to the filtered model such that they converge to self-similar collapsing orbits in the three point-vortex system and dissipate the enstrophy at the event of collapse in the zero limit of the filter parameter. In this study, we numerically show that the enstrophy dissipation by the collapse of point vortices could occur for the four and five vortex problems in a filtered model. Moreover, we show the detailed convergence process of the point vortices for gradually decreasing filter parameters, which provides a new insight for the three vortex problem. In addition, numerical computations suggest that the enstrophy dissipation is caused by collapse of separated point vortices with the negative interactive energy.

The electromagnetically driven magnetised spherical Couette flow is studied experimentally, theoretically and numerically in the laminar regime. The working fluid, Galinstan, is contained in the spherical gap between two concentric spheres at rest. The electromagnetic stirring is primarily generated due to the interaction of a direct current, which is injected through two ring-shaped electrodes located at the equatorial zone of each sphere, and a dipolar magnetic field produced by a permanent magnet located inside the inner sphere. The flows were explored experimentally for a Reynolds number ranging from 450 to 2230 and a Hartmann number of 240. Ultrasound Doppler velocimetry and particle image velocimetry were used to characterise the flow. For low Reynolds numbers, given the symmetry of the problem, a one-dimensional analytic solution is obtained in the equatorial plane from the magnetohydrodynamic equations. The analytical solution reproduces the main characteristics of the flow. In addition, a full three-dimensional numerical model is able to reproduce both the analytical solution and the experimental measurements. To the best knowledge of the authors, this is the first time experimental results of the magnetised spherical Couette flow have been reported with electromagnetic forcing using a liquid metal as the working media.

Experimental investigation of the Rayleigh–Taylor instability (RTI) and its dependence on initial conditions has been challenging, primarily due to the difficulty of creating a well-defined gaseous interface. To address this, a novel soap film technique was developed to create a discontinuous two-dimensional SF$_6$air interface with precisely controlled initial conditions. High-order modes were superimposed on a long-wavelength perturbation to study the influence of initial conditions on RTI evolution. Experiments conducted at Atwood numbers ranging from 0.26 to 0.66 revealed that bubble growth shows a weak dependence on both initial conditions and Atwood numbers, whereas spike growth is more influenced by these factors. Spike growth accelerates as the wavenumber of the imposed high-order modes decreases and/or the Atwood number increases. To quantify these effects, a variation on the previously developed potential flow model was applied, capturing the suppression of high-order modes and Atwood number dependence on RTI growth. In turbulent flow, the self-similar factors of bubbles and spikes exhibit minimal sensitivity to initial conditions. However, in relation to the Atwood number, the self-similar factors of bubbles (or spikes) demonstrate negligible (or significant) dependence. Comparisons with literature revealed that two-dimensional flows yield lower self-similar factors than three-dimensional flows. Furthermore, the discontinuity of the initial interface in this study, achieved through the soap film technique, results in faster spike growth compared with previous studies involving a diffusive initial interface. These findings provide critical insights into the nonlinear dynamics of RTI and underscore the importance of well-characterised initial conditions in experimental studies.

This study investigates the interactions between flexural-gravity waves and interfacial waves in a two-layer fluid, focusing on wave blocking. Both liquid layers are of finite depth bounded on top by a viscoelastic thin plate. Both liquids are incompressible and inviscid, and their flows are two-dimensional and potential. Linear wave theory and a linear equation of a thin floating viscoelastic plate of constant thickness are used. We analyse the phenomenon of wave blocking and Kelvin–Helmholtz (KH) instability in a two-layer fluid with a discontinuous background mean flow. A quartic dispersion relation for frequency as a function of wavenumber and other parameters of the problem is derived. Two cases of uniform current and layers moving with different velocities are studied. Wave blocking occurs when roots of the dispersion relation coalesce without accounting for plate viscosity, leading to zero group velocity. Our findings indicate that wave blocking can occur for both flexural-gravity and interfacial waves under various frequency and current speed conditions, provided that plate viscosity is absent. The role of different parameters and the flow velocities of the upper and lower layers are investigated in the occurrence of wave blocking and KH instability. The loci of the roots of the dispersion relation involving plate viscosity depict that no root coalescence occurs irrespective of the values of wavenumber and frequency in the presence of plate viscosity. The amplitude ratio of the interfacial wave elevation to that of floating viscoelastic plate deflection exhibits the dead-water phenomenon as a density ratio approaches unity.