In a combined experimental and numerical effort, we investigate the generation and reduction of airfoil tonal noise. The means of noise control are streak generators in the form of cylindrical roughness elements. These elements are placed periodically along the span of the airfoil at the mid-chord streamwise position. Experiments are performed for a wide range of Reynolds numbers and angles of attack in a companion work (Alva et al., AIAA Aviation Forum, 2023). In the present work, we concentrate on numerical investigations for a further investigation of selected cases. We have performed wall-resolved large-eddy simulations for a NACA 0012 airfoil at zero angle of attack and Mach 0.3. Two Reynolds numbers (${0.8\times 10^{5}}$ and ${1.0 \times 10^{5}}$) have been investigated, showing acoustic results consistent with experiments at the same Reynolds but lower Mach numbers. Roughness elements attenuate tones in the acoustic field and, for the higher Reynolds number, suppress them. Through Fourier decomposition and spectral proper orthogonal decomposition analysis of streamwise velocity data, dominating structures have been identified. Further, the coupling between the structures generated by the surface roughness and the instability modes (Kelvin–Helmholtz) of the shear layer has been identified through stability analysis, suggesting stabilisation mechanisms by which the sound generation by the airfoil is reduced by the roughness elements.

This study investigates the transport of particles in turbulent channel flow with friction Reynolds number $Re_\tau = 1000$ by direct numerical simulation. We focus on how large-scale flow structures, namely the $Qs$ structures (Lozano-Durán et al. 2012, J. Fluid Mech., vol. 694, pp. 100–130), affect the wall-normal transport of particles. Despite occupying less than $10\,\%$ of the physical domain, our results highlight the critical role played by $Qs$ structures in the particle transport, namely that the particle number and momentum flux inside the $Qs$ structures are substantially higher than outside. The fraction of particle wall-normal momentum flux inside $Qs$ structures is considerably larger than their volume fraction, suggesting highly efficient transport inside the $Qs$ structures. This prominent role played by $Qs$ structures in the transport of inertial particles is more effective by diminishing the inertia of particles. Notably, the long-distance transport of particles in the wall-normal direction is driven primarily by the continuous effect of $Qs$ structures. In summary, our findings advance the understanding of the effects of $Qs$ structures on particle transport, and demonstrate their significant role in the process.

Volcanic fissure eruptions typically start with the opening of a linear fissure that erupts along its entire length, following which, activity localises to one or more isolated vents within a few hours or days. Localisation is important because it influences the spatiotemporal evolution of the hazard posed by the eruption. Previous work has proposed that localisation can arise through a thermoviscous fingering instability driven by the strongly temperature dependent viscosity of the rising magma. Here, we explore how thermoviscous localisation is influenced by the irregular geometry of natural volcanic fissures. We model the pressure-driven flow of a viscous fluid with temperature-dependent viscosity through a narrow fissure with either sinusoidal or randomised deviations from a uniform width. We identify steady states, determine their stability and quantify the degree of flow enhancement associated with localised flow. We find that, even for relatively modest variations of the fissure width (${\lt } 10$ %), the non-planar geometry supports strongly localised steady states, in which the wider parts of the fissure host faster, hotter flow, and the narrower parts of the fissure host slower, cooler flow. This geometrically driven localisation differs from the spontaneous thermoviscous fingering observed in planar geometries and can strongly impact the localisation process. We delineate the regions of parameter space under which geometrically driven localisation is significant, showing that it is a viable mechanism for the observed localisation under conditions typical of basaltic eruptions, and that it has the potential to dominate the effects of spontaneous thermoviscous fingering in these cases.

Ocean submesoscales, flows with characteristic size $10\,\text{m}{-}10\,\text{km}$, are transitional between the larger, rotationally constrained mesoscale and three-dimensional turbulence. In this paper, we present simulations of a submesoscale ocean filament. In our case, the filament is strongly sheared in both vertical and cross-filament directions, and is unstable. Instability indeed dominates the early behaviour with a fast extraction of kinetic energy from the vertically sheared thermal wind. However, the instability that emerges does not exhibit characteristics that match the perhaps expected symmetric or Kelvin–Helmholtz instabilities, and appears to be non-normal in nature. The prominence of the transient response depends on the initial noise, and for large initial noise amplitudes, saturates before symmetric instability normal modes are able to develop. The action of the instability is sufficiently rapid – with energy extraction from the mean flow emerging and peaking within the first inertial period ($\sim\! 18\ \text{h}$) – that the filament does not respond in a geostrophically balanced sense. Instead, at all initial noise levels, it later exhibits vertically sheared near-inertial oscillations with higher amplitude as the initial minimum Richardson number decreases. Horizontal gradients strengthen only briefly as the fronts restratify. These unstable filaments can be generated by strong mixing events at pre-existing stable structures; we also caution against inadvertently triggering this response in idealised studies that start in a very unstable state.

Direct numerical simulations have been conducted to explore the coupling effect of the thermoelectric effect and vertical convection (VC) in a square cavity composed of liquid lithium and stainless steel under different Hartmann numbers at $Ra=10^5$. By leveraging thermoelectric phenomena, an innovative approach is proposed to actively modulate heat transfer efficiency. The core concept lies in modulating the intensity of large-scale circulation (LSC) in VC systems through the torque generated by the interaction between thermoelectric currents and magnetic fields via Lorentz forces. The findings reveal that when the torque aligns with the direction of LSC induced by pure buoyancy, both momentum and heat transfer are enhanced. However, due to the magnetic damping itself, this enhancement is not sustained indefinitely, resulting in a trend of initial increase followed by decline in both momentum and heat transfer efficiency. Conversely, when the magnetic field direction is reversed, causing the Lorentz force torque to oppose the buoyancy-driven circulation, both momentum and heat transfer efficiency diminish until the flow reverses. By varying the magnetic field intensity, three distinct flow regimes are identified: the buoyancy-dominated regime, the thermoelectric-dominated regime and the magnetic-damping-dominated regime. The transition between the buoyancy-dominated regime and thermoelectric-dominated regime – specifically, the onset of flow reversal – is analysed through a boundary-layer–bulk–boundary-layer coupling model. This model enables precise prediction of the critical $Ha$ based on the torque balance between buoyancy forces and thermoelectrically induced Lorentz forces, and demonstrates close agreement with numerical simulations.

We present the results of a theoretical investigation of orbital stability in pilot-wave hydrodynamics, wherein a droplet bounces and self-propels across the surface of a vertically vibrating liquid bath. A critical notion in pilot-wave hydrodynamics is that the bath plays the role of the system memory, recording the history of the droplet in its wave field. Quantised orbital motion may arise when the droplet is confined by either an axisymmetric potential or the Coriolis force induced by system rotation. We here elucidate the dependence of the stability of circular orbits on both the form of the confining force and the system memory. We first provide physical insight by distinguishing between potential- and wave-driven instabilities. We demonstrate that the former are a generic feature of classical orbital dynamics at constant speed, while the latter are peculiar to pilot-wave systems. The wave-driven instabilities are marked by radial perturbations that either grow monotonically or oscillate at an integer multiple of the orbital frequency, in which case they are said to be resonant. Conversely, for potential-driven wobbling, the instability frequency may be resonant or non-resonant according to the form of the applied potential. Asymptotic analysis rationalises the different stability characteristics for linear-spring and Coriolis forces, the two cases that have been explored experimentally. Our results are generalised to consider other potentials of interest in pilot-wave hydrodynamics, and elucidate the distinct roles of wave- and potential-driven instabilities. Our study highlights the limitations of prior heuristic arguments for predicting the onset of orbital instability.

The present work experimentally investigates the interaction of a buoyant (rigid) spherical particle with a single translating (water) vortex ring, focusing on the effects of particle-to-vortex core size ratio ($D_p/D_{c,o}$) on both the particle dynamics and ring dynamics ($D_p$ = particle diameter, $D_{c,o}$ = vortex core diameter). These interactions are studied for $D_p/D_{c,o}$ = 0.6–1.7, over ring Reynolds numbers ($Re={\varGamma }/{\nu }$; $\varGamma$ = ring circulation) of 6000–67 300. As the buoyant particle comes close to the ring, it gets captured into the low-pressure vortex core, and the interaction begins. The particle within the core undergoes radial oscillation, spins and translates along the ring’s azimuthal axis. As $D_p/D_{c,o}$ increases, the particle undergoes higher-amplitude radial oscillation and a relatively shorter azimuthal translation. The differences in the particle size and its motion within the ring lead to large differences in the ring’s dynamics. A larger particle is seen to lead to a higher ring disruption, substantially reducing the ring’s convection speed and azimuthal enstrophy, which are seen to scale as $(D_p/D_{c,o})^{2.3}Re^{-0.37}$ and $(D_p/D_{c,o})^{1.3}Re^{-0.25}$, respectively. The ring disruption is significant above $D_p/D_{c,o}\approx$ 1.0, beyond which the ring fragments, with up to 60 % drop in convection speed and 90 % drop in enstrophy, at low $Re$, as compared with the base ring. These results for the rigid particle size effects on the vortex ring dynamics are more dramatic than for a deforming bubble. Our results could help to better understand and model buoyant particle (and bubble) interactions with coherent structures in turbulence.

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.

This paper presents a theoretical investigation of vortex modes in acoustofluidic cylindrical resonators with rigid boundaries and viscous fluids. By solving the Helmholtz equation for linear pressure, incorporating boundary conditions that account for no-slip surfaces and vortex and non-vortex excitation at the base, we analyse both single- and dual-eigenfunction modes near system resonance. The results demonstrate that single-vortex modes generate spin angular momentum exclusively along the axial direction, while dual modes introduce a transverse spin component due to the nonlinear interaction between axial and transverse ultrasonic waves, even in the absence of vortex excitation. We find that nonlinear acoustic fields, including energy density, radiation force potential and spin, scale with the square of the shear wave number, defined as the ratio of the cavity radius to the thickness of the viscous boundary layer. Theoretical predictions align closely with finite element simulations based on a model for an acoustofluidic cavity with adiabatic and rigid walls. These findings hold particular significance for acoustofluidic systems, offering potential applications in the precise control of cells and microparticles.

The existing studies on vortex rings have concentrated on non-zero circulation. However, the cases of zero circulation may also be significantly noteworthy on both theoretical and practical grounds. As the first attempt on this subject, in this paper a family of viscous laminar vortex rings with zero circulation and a moderate ratio of core radius to ring radius is studied using numerical simulations of the incompressible Navier–Stokes equations. This unusual zero circulation is achieved by assigning a special layered vorticity distribution with alternate signs to the vortex core. At the initial moment, the ring is axisymmetric, swirl-free and of a circular cross-section. It is found that the axial symmetry and the non-swirl nature of the vortex ring are preserved during the evolution, and the vortex ring endures a transition from the initial layered structure to a shell structure, then degenerates to an ordinary vortex ring with non-zero circulation at last. Significant vorticity cancellation is observed due to the interactions among the layered structures. A new Reynolds number, based on the absolute value of vorticity, is applied to the zero-circulation vortex rings in the present work. For such vortex rings, cases of both zero and non-zero vortical impulse can happen, unlike the ordinary ones with only non-zero vortical impulse. Additionally, it is found that the vortical impulse can be irrelevant to the ring diameter. The study may shed light on modelling certain real flows characterised by distinct vortex structures or configurations.

Generation of steady streaming vortices is usually accomplished by mechanically vibrating bodies, as occurs in several microfluidic applications for mixing, as well as for transport and handling of microparticles. Here, we propose the generation of streaming from the harmonic electromagnetic forcing of a free-moving circular magnet held afloat on a shallow electrolytic layer, and show that the sense of rotation of steady vortices is the opposite to that of the classical streaming flow. Reverse streaming is attributed to the coupling between the fluid and the free-moving body. Analytical solutions offer a physical rationale for the observed flow dynamics, while numerical simulation reproduces the experimental observations satisfactorily.

Non-spherical bubble collapses near solid boundaries, generating water hammer pressures and shock waves, were recognized as key mechanisms for cavitation erosion. However, there is no agreement on local erosion patterns, and cavitation erosion damage lacks quantitative analysis. In our experiments, five distinct local erosion patterns were identified on aluminium sample surfaces, resulting from the collapse of laser-induced cavitation bubbles at moderate stand-off distances of $0.4\leqslant \gamma \leqslant 2.2$, namely bipolar, monopolar, annular, solar-halo and central. Among them, the bipolar and monopolar patterns exhibit the most severe cavitation erosion when the toroidal bubbles undergo asymmetrical collapse along the circumferential direction during the second cycle. Shadowgraphy visualization revealed that asymmetrical collapse caused shockwave focusing through head-on collision and oblique superposition of wavefronts. This led to the variations in toroidal bubble radii and the positions of maximum erosion depth not matching at certain stand-off distances. Both initial plasma asymmetry and bubble–wall stand-off distance were critical in determining circumferential asymmetrical collapse behaviours. At large initial aspect ratios, the elliptical jet tips form during the contraction process, resulting in the toroidal bubble collapsing from regions with smaller curvature radii, ultimately converging to the colliding point along the circumferential direction. Our three-dimensional simulations using OpenFOAM successfully reproduce the key features of circumferentially asymmetrical bubble collapse. This study provides new insights into the non-spherical near-wall bubble collapse dynamics and provides a foundation for developing predictive models for cavitation erosion.

Periodic travelling waves at the free surface of an incompressible inviscid fluid in two dimensions under gravity are numerically computed for an arbitrary vorticity distribution. The fluid domain over one period is conformally mapped from a fixed rectangular one, where the governing equations along with the conformal mapping are solved using a finite-difference scheme. This approach accommodates internal stagnation points, critical layers and overhanging profiles, thereby overcoming limitations of previous studies. The numerical method is validated through comparisons with known solutions for zero and constant vorticity. Novel solutions are presented for affine vorticity functions and a two-layer constant-vorticity scenario.

We consider the conceptual two-layered oscillating tank of Inoue & Smyth (2009 J. Phys. Oceanogr. vol. 39, no. 5, pp. 1150–1166), which mimics the time-periodic parallel shear flow generated by low-frequency (e.g. semi-diurnal tides) and small-angle oscillations of the density interface. Such self-induced shear of an oscillating pycnocline may provide an alternate pathway to pycnocline turbulence and diapycnal mixing in addition to the turbulence and mixing driven by wind-induced shear of the surface mixed layer. We theoretically investigate shear instabilities arising in the inviscid two-layered oscillating tank configuration and show that the equation governing the evolution of linear perturbations on the density interface is a Schrödinger-type ordinary differential equation with a periodic potential. The necessary and sufficient stability condition is governed by a non-dimensional parameter $\beta$ resembling the inverse Richardson number; for two layers of equal thickness, instability arises when $\beta \,{\gt}\,1/4$. When this condition is satisfied, the flow is initially stable but finally tunnels into the unstable region after reaching the time marking the turning point. Once unstable, perturbations grow exponentially and reveal characteristics of Kelvin–Helmholtz (KH) instability. The modified Airy function method, which is an improved variant of the Wentzel–Kramers–Brillouin theory, is implemented to obtain a uniformly valid, composite approximate solution to the interface evolution. Next, we analyse the fully nonlinear stages of interface evolution by modifying the circulation evolution equation in the standard vortex blob method, which reveals that the interface rolls up into KH billows. Finally, we undertake real case studies of Lake Geneva and Chesapeake Bay to provide a physical perspective.

This paper numerically investigates the heat transport and bifurcation of natural convection in a differentially heated cavity filled with entangled polymer solution combined with the boundary layer and kinetic energy budget analysis. The polymers are described by the Rolie-Poly model, which effectively captures the rheological response of entangled polymers. The results indicate that the competition between its shear-thinning and elasticity dominates the flow structures and heat transfer rate. The addition of polymers tends to enhance the heat transfer as the polymer viscosity ratio ($\beta$) decreases or the relaxation time ratio ($\xi$) increases. The amount of heat transfer enhancement (HTE) behaves non-monotonically, which first increases significantly and then remains almost constant or decreases slightly with the Weissenberg number ($Wi$). The critical $Wi$ gradually increases with the increasing $\xi$, where the maximum HTE reaches approximately $64.9\,\%$ at $\beta = 0.1$. It is interesting that even at low Rayleigh numbers, the flow transitions from laminar to periodic flows in scenarios with strong elasticity. The bifurcation is subcritical and exhibits a typical hysteresis loop. Then, the bifurcation routes driven by inertia and elasticity are examined by direct numerical simulations. These results are illustrated by time histories, Fourier spectra analysis and spatial structures observed at varying time intervals. The kinetic energy budget indicates that the stretch of the polymers leads to great energy exchange between polymers and flow structures, which plays a crucial role in the hysteresis phenomenon. This dynamic behaviour contributes to the strongly self-sustained and self-enhancing processes in the flow.

An analytical formulation is provided that describes the first two natural modes of the fluid–structure interaction of an incompressible current with a pitching and heaving flexible plate. The objective is twofold: first, to present a general derivation of analytical expressions for the lift, moment and the flexural moments exerted by an inviscid flow on a pitching and heaving plate whose deformation is general enough that the coupling of the flexural moments with the structural equations allows solving analytically the first two natural modes of the system; second, to analyse the propulsion performance of the foil when actuated near the first two natural frequencies. For the second purpose, one also needs the thrust force generated through the motion and the general deformation of the foil considered, which is analytically derived using the linearized vortex impulse theory, extending and systematizing previous works. The analytical expressions, once viscous effects are taken into consideration through nonlinear transverse damping and offset drag coefficients, are compared with small-amplitude available experimental data, discussing their limitations. It is found that low stiffness pitching and heaving are quite different, with a pitching flexible foil only generating thrust near the second resonant frequency, whereas heaving always generates thrust, with the maximum slightly below the second natural frequency. Maximum thrust for large stiffness pitching is around the first natural frequency. The maximum efficiency occurs at frequencies close to the first natural mode if the foil is sufficiently rigid, but it is not related to the natural frequencies as the rigidity decreases.

Rayleigh–Taylor (RT) stability occurs when a single-mode light/heavy interface is accelerated by rarefaction waves, exhibiting a sustained oscillation in perturbation amplitude. If the perturbation is accelerated again by a shock propagating in the same direction as the rarefaction waves, the interface evolution will shift from RT stability to Richtmyer–Meshkov (RM) instability. Depending upon the interface state when the shock arrives, the perturbation growth can be actively manipulated through controlling the magnitudes of vorticity deposited by rarefaction and shock waves. The present work first theoretically analyses the 12 different growth possibilities of a light/heavy interface accelerated by co-directional rarefaction and shock waves. A theoretical model is established by combining the RT growth rate with the RM growth rate, providing the conditions for the different possibilities of the perturbation growth. Based on the model, extensive experiments are designed and conducted in the specially designed rarefaction-shock tube. By precisely controlling the shock arrival time at the interface, the different growth possibilities, including promotion, reduction and freeze-out, are realised in experiments. This work verifies the feasibility of manipulating the light/heavy perturbation via co-directional rarefaction and shock waves, which sheds light on control of hydrodynamic instabilities in practical applications.

This paper presents an experimental and analytical investigation into the use of trailing edge slits for the reduction of aerofoil trailing edge noise. The noise reduction mechanism is shown to be fundamentally different from conventional trailing edge serrations, relying on destructive interference from highly compact and coherent sources generated at either ends of the slit. This novel approach is the first to exploit the coherence intrinsic to the boundary layer turbulence. Furthermore, the study demonstrates that trailing edge slits not only achieve superior noise reductions compared with sawtooth serrations of the same amplitude at certain conditions, but also offer frequency-tuning capability for noise reduction. Noise reduction is driven by the destructive interference between acoustic sources at the root and tip of the slit, which radiate with a phase difference determined by the difference in times taken for the boundary layer flow to convect between the root and tip. Maximum noise reductions occur at frequencies where the phase difference between these sources is $180^\circ$. The paper also presents a detailed parametric study into the variation in noise reductions due to the slit length, slit wavelength and slit root width. Additionally, a simple two-source analytic model is proposed to explain the observed results. Wind tunnel measurements of the unsteady flow field around the trailing edge slits are also presented, providing insights into the underlying flow physics.

The Richtmyer–Meshkov instability at gas interfaces with controllable initial perturbation spectra under reshock conditions is investigated both experimentally and theoretically. A soap-film method is adopted to generate well-defined single-, dual- and triple-mode air/SF$_6$ interfaces. By inserting an acrylic block into the test section, a reflected shock with controllable reshock timing is created. The results reveal a complex relationship between the post-reshock perturbation growth rate and the pre-reshock interface morphology. For single-mode interfaces, the post-reshock growth rate exhibits a strong dependence on pre-reshock conditions. In contrast, for multi-mode interfaces, this dependence weakens significantly due to mode-coupling effects. It is found that, following reshock, each fundamental mode develops independently and later is significantly influenced by mode-coupling effects. Based on this finding, we propose an empirical model that matches the initial linear growth rate and the asymptotic growth rate, accurately predicting the evolution of fundamental modes from early to late stages across all three configurations. Furthermore, a theoretical formula is derived, linking the empirical coefficient in the model of Charakhch’An (2020 J. Appl. Mech. Tech. Phys. vol. 41, no. 1, pp. 23–31) to the initial perturbation. This provides a unified framework to explain the varying dependence of post-reshock growth rates on pre-reshock morphology observed in previous experiments.

A new model is presented for the decay of plane shock waves in equilibrium flows with an arbitrary equation of state. A fundamental challenge for the accurate prediction of shock propagation using analytical modelling is to account for the coupling between a shock’s motion and the post-shock flow. Our model accomplishes this by neglecting only higher-order perturbations to the second velocity gradient, $u_{xx}$, in the incident simple wave. The second velocity gradient is generally small and exactly zero for centred expansion waves in a perfect gas, so neglecting its effect on the shock motion provides an accurate closure criterion for a shock-change equation. This second-order shock-change equation is derived for a general equation of state. The model is tested by comparison with numerical simulations for three problems: decay by centred waves in a perfect gas, decay by centred waves in equilibrium air and decay by the simple wave generated from the constant deceleration of piston in a perfect gas. The model is shown to be exceptionally accurate for a wide range of conditions, including small $\gamma$ and large shock Mach numbers. For a Mach 15 shock in equilibrium air, model errors are less than 2 % in the first 60 % of the shock’s decay. The analytical results possess a simple formulation but are applicable to fluids with a general equation of state, enabling new insight into this fundamental problem in shock wave physics.