The primary bifurcation of the flow past three-dimensional axisymmetric bodies is investigated. We show that the azimuthal vorticity generated at the body surface is at the root of the instability, and that the mechanism proposed by Magnaudet & Mougin (2007, J. Fluid Mech., vol. 572, 311–337) in the context of spheroidal bubbles extends to axisymmetric bodies with a no-slip surface. The instability arises in a thin region of the flow in the near wake, and is associated with the occurrence of strong vorticity gradients. We propose a simple yet effective scaling law for the prediction of the instability, based on a measure of the near-wake vorticity and of the radial extent of the separation bubble. At criticality, the resulting Reynolds number collapses approximately to a constant value for bodies with different geometries and aspect ratios, with a relative variation that is one order of magnitude smaller than that of the standard Reynolds number based on the free-stream velocity and body diameter. The new scaling can be useful to assess whether the steady flow past axisymmetric bodies is globally unstable, without the need for an additional stability analysis.

Yaw control can effectively enhance wind farm power output, but the vorticity distribution and coherent structures in yawed turbine wakes remain poorly understood. We propose a physical model capable of accurately predicting tip vortex dynamics from their generation to destabilisation. This model integrates a point vortex framework with advanced blade element momentum theory and vortex cylinder theory for yawed turbines. Comparisons with large eddy simulations demonstrate that the model effectively predicts the vorticity distribution of tip vortices and the wake profile of yawed turbines. Finally, we employ sparsity-promoting dynamic mode decomposition to analyse the dynamics of the far wake. Our analysis reveals four primary mode types: (i) the averaged mode; (ii) shear modes; (iii) harmonic modes; and (iv) merging modes. Under yawed conditions, these modes become asymmetric, leading to interactions between the tip and root vortex modes. This direct interaction plays a critical role during the formation process of the counter-rotating vortex pair observed in yawed wakes.

This study investigated the cylindrically divergent Rayleigh–Taylor instability (RTI) on a liquid–gas interface and its dependence on initial conditions. A novel hydrophobic technique was developed to generate a two-dimensional water–air interface with controlled initial conditions. The experimental configuration utilised high-pressure air injection to produce uniform circumferential acceleration. Amplitude measurements over time revealed that the cylindrical RTI growth depends strongly on the azimuthal wavenumber. Experimental results demonstrated that surface tension significantly suppresses the liquid–gas cylindrical RTI, even inducing a freeze-out and oscillatory perturbation growth – a phenomenon observed for the first time. Spectrum analysis of the interface contours demonstrated that the cylindrical RTI evolves in a weakly nonlinear regime. Linear and weakly nonlinear models were derived to accurately predict the time-varying interface amplitudes and high-order modes. The linear model was further used to determine conditions for unstable, freeze-out and oscillatory solutions of the cylindrically divergent RTI. These findings offer valuable insights into manipulating hydrodynamic instabilities in contracting/expanding geometries using surface tension.

The linear stability of a thermally stratified fluid layer between horizontal walls, where non-Brownian thermal particles are injected continuously at one boundary and extracted at the other – a system known as particulate Rayleigh–Bénard (pRB) – is studied. For a fixed volumetric particle flux and minimal thermal coupling, reducing the injection velocity stabilises the system when heavy particles are introduced from above, but destabilises it when light particles are injected from below. For very light particles (bubbles), low injection velocities can shift the onset of convection to negative Rayleigh numbers, i.e. heating from above. Particles accumulate non-uniformly near the extraction wall and in regions of strong vertical flow, aligning with either wall-impinging or wall-detaching zones depending on whether injection is at sub- or super-terminal velocity. The increase of the volumetric particle flux always enhances these effects.

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.

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 paper explores the construction of quadratic Lyapunov functions for establishing the conditional stability of shear flows described by truncated ordinary differential equations, addressing the limitations of traditional methods like the Reynolds–Orr equation and linear stability analysis. The Reynolds–Orr equation, while effective for predicting unconditional stability thresholds in shear flows due to the non-contribution of nonlinear terms, often underestimates critical Reynolds numbers. Linear stability analysis, conversely, can yield impractically high limits due to subcritical transitions. Quadratic Lyapunov functions offer a promising alternative, capable of proving conditional stability, albeit with challenges in their construction. Typically, sum-of-squares programs are employed for this purpose, but these can result in sizable optimisation problems as system complexity increases. This study introduces a novel approach using linear transformations described by matrices to define quadratic Lyapunov functions, validated through nonlinear optimisation techniques. This method proves particularly advantageous for large systems by leveraging analytical gradients in the optimisation process. Two construction methods are proposed: one based on general optimisation of transformation matrix coefficients, and another focusing solely on the system’s linear aspects for more efficient Lyapunov function construction. These approaches are tested on low-order models of subcritical transition and a two-dimensional Poiseuille flow model with degrees of freedom nearing 1000, demonstrating their effectiveness and efficiency compared with sum-of-squares programs.

The evolution of a Lamb–Oseen vortex is studied in a stratified rotating fluid under the complete Coriolis force. In a companion paper, it was shown that the non-traditional Coriolis force generates a vertical velocity field and a vertical vorticity anomaly at a critical radius when the Froude number is larger than unity. Below a critical non-traditional Rossby number $\widetilde {Ro}$, a two-dimensional shear instability was next triggered by the vorticity anomaly. Here, we test the robustness of this two-dimensional evolution against small three-dimensional perturbations. Direct numerical simulations (DNS) show that the two-dimensional shear instability then develops only in an intermediate range of non-traditional Rossby numbers for a fixed Reynolds number $Re$. For lower $\widetilde {Ro}$, the instability is three-dimensional. Stability analyses of the flows in the DNS prior to the instability onset fully confirm the existence of these two competing instabilities. In addition, stability analyses of the local theoretical flows at leading order in the critical layer demonstrate that the three-dimensional instability is due to the shear of the vertical velocity. For a given Froude number, its growth rate scales as $Re^{2/3}/\widetilde {Ro}$, whereas the growth rate of the two-dimensional instability depends on $Re/\widetilde {Ro}^2$, provided that the critical layer is smoothed by viscous effects. However, the growth rate of the three-dimensional instability obtained from such local stability analyses agrees quantitatively with those of the DNS flows only if second-order effects due to the traditional Coriolis force and the buoyancy force are taken into account. These effects tend to damp the three-dimensional instability.

Stall cells are transverse cellular patterns that often appear on the suction side of airfoils near stalling conditions. Wind-tunnel experiments on a NACA4412 airfoil at Reynolds number ${Re}=3.5 \times 10^5$ show that they appear for angles of attack larger than $\alpha = 11.5^{\circ }\ (\pm 0.5^{\circ })$. Their onset is further investigated based on global stability analyses of turbulent mean flows computed with the Reynolds-averaged Navier–Stokes (RANS) equations. Using the classical Spalart–Allmaras turbulence model and following Plante et al. (J. Fluid Mech., vol. 908, 2021, A16), we first show that a three-dimensional stationary mode becomes unstable for a critical angle of attack $\alpha = 15.5^{\circ }$ which is much larger than in the experiments. A data-consistent RANS model is then proposed to reinvestigate the onset of these stall cells. Through an adjoint-based data-assimilation approach, several corrections in the turbulence model equation are identified to minimize the differences between assimilated and reference mean-velocity fields, the latter reference field being extracted from direct numerical simulations. Linear stability analysis around the assimilated mean flow obtained with the best correction is performed first using a perturbed eddy-viscosity approach which requires the linearization of both RANS and turbulence model equations. The three-dimensional stationary mode becomes unstable for angle $\alpha = 11^{\circ }$ which is in significantly better agreement with the experimental results. The interest of this perturbed eddy-viscosity approach is demonstrated by comparing with results of two frozen eddy-viscosity approaches that neglect the perturbation of the eddy viscosity. Both approaches predict the primary destabilization of a higher-wavenumber mode which is not experimentally observed. Uncertainties in the stability results are quantified through a sensitivity analysis of the stall cell mode's eigenvalue with respect to residual mean-flow velocity errors. The impact of the correction field on the results of stability analysis is finally assessed.

Understanding interfacial instability in a coflow system has relevance in the effective manipulation of small objects in microfluidic applications. We experimentally elucidate interfacial instability in stratified coflow systems of Newtonian and viscoelastic fluid streams in microfluidic confinements. By performing a linear stability analysis, we derive equations that describe the complex wave speed and the dispersion relationship between wavenumber and angular frequency, thus categorizing the behaviour of the systems into two main regimes: stable (with a flat interface) and unstable (with either a wavy interface or droplet formation). We characterize the regimes in terms of the capillary numbers of the phases in a comprehensive regime plot. We decipher the dependence of interfacial instability on fluidic parameters by decoupling the physics into viscous and elastic components. Remarkably, our findings reveal that elastic stratification can both stabilize and destabilize the flow, depending on the fluid and flow parameters. We also examine droplet formation, which is important for microfluidic applications. Our findings suggest that adjusting the viscous and elastic properties of the fluids can control the transition between wavy and droplet-forming unstable regimes. Our investigation uncovers the physics behind the instability involved in interfacial flows of Newtonian and viscoelastic fluids in general, and the unexplored behaviour of interfacial waves in stratified liquid systems. The present study can lead to a better understanding of the manipulation of small objects and production of droplets in microfluidic coflow systems.

We study the existence and strong instability of standing wave solutions for the fractional nonlinear Schrödinger equation

\begin{equation*}\begin{cases}\mathrm{i} \psi_{t}=\left(-\Delta\right)^s \psi-\left(|\psi|^{p-2}\psi+\mu|\psi|^{q-2}\psi\right), & \quad(t,x)\in(0,\infty)\times\mathbb{R}^N,\\\psi(0,x)=\psi_0(x),&\quad x\in\mathbb{R}^N,\end{cases}\end{equation*}

where $N\geq2$, $0 \lt s \lt 1$, $2 \lt q \lt p \lt 2_s^*=2N/(N-2s)$, and $\mu\in\mathbb{R}$. The primary challenge lies in the inhomogeneity of the nonlinearity.We deal with the following three cases: (i) for $2 \lt q \lt p \lt 2+4s/N$ and µ < 0, there exists a threshold mass a0 for the existence of the least energy normalized solution; (ii) for $2+4s/N \lt q \lt p \lt 2_s^*$ and µ > 0, we reveal the existence of the ground state solution, explore the strong instability of standing waves, and provide a blow-up criterion; (iii) for $2 \lt q\leq2+4s/N \lt p \lt 2_s^*$ and µ < 0, the strong instability of standing wave solutions is demonstrated. These findings are illuminated through variational characterizations, the profile decomposition, and the virial estimate.

The resolvent analysis reveals the worst-case disturbances and the most amplified response in a fluid flow that can develop around a stationary base state. The recent work by Padovan et al. (J. Fluid Mech., vol. 900, 2020, A14) extended the classical resolvent analysis to the harmonic resolvent analysis framework by incorporating the time-varying nature of the base flow. The harmonic resolvent analysis can capture the triadic interactions between perturbations at two different frequencies through a base flow at a particular frequency. The singular values of the harmonic resolvent operator act as a gain between the spatiotemporal forcing and the response provided by the singular vectors. In the current study, we formulate the harmonic resolvent analysis framework for compressible flows based on the linearized Navier–Stokes equation (i.e. operator-based formulation). We validate our approach by applying the technique to the low-Mach-number flow past an airfoil. We further illustrate the application of this method to compressible cavity flows at Mach numbers of 0.6 and 0.8 with a length-to-depth ratio of $2$. For the cavity flow at a Mach number of 0.6, the harmonic resolvent analysis reveals that the nonlinear cross-frequency interactions dominate the amplification of perturbations at frequencies that are harmonics of the leading Rossiter mode in the nonlinear flow. The findings demonstrate a physically consistent representation of an energy transfer from slow-evolving modes toward fast-evolving modes in the flow through cross-frequency interactions. For the cavity flow at a Mach number of 0.8, the analysis also sheds light on the nature of cross-frequency interaction in a cavity flow with two coexisting resonances.

Zonal flows are mean flows in the east–west direction, which are ubiquitous on planets, and can be formed through ‘zonostrophic instability’: within turbulence or random waves, a weak large-scale zonal flow can grow exponentially to become prominent. In this paper, we study the statistical behaviour of the zonostrophic instability and the effect of magnetic fields. We use a stochastic white noise forcing to drive random waves, and study the growth of a mean flow in this random system. The dispersion relation for the growth rate of the expectation of the mean flow is derived, and properties of the instability are discussed. In the limits of weak and strong magnetic diffusivity, the dispersion relation reduces to manageable expressions, which provide clear insights into the effect of the magnetic field and scaling laws for the threshold of instability. The magnetic field mainly plays a stabilising role and thus impedes the formation of the zonal flow, but under certain conditions it can also have destabilising effects. Numerical simulation of the stochastic flow is performed to confirm the theory. Results indicate that the magnetic field can significantly increase the randomness of the zonal flow. It is found that the zonal flow of an individual realisation may behave very differently from the expectation. For weak magnetic diffusivity and moderate magnetic field strengths, this leads to considerable variation of the outcome, that is whether zonostrophic instability takes place or not in individual realisations.

We examine the mechanisms responsible for the onset of the three-dimensional mode B instability in the wake behind a circular cylinder. We show that it is possible to explicitly account for the stabilising effect of spanwise viscous diffusion and then demonstrate that the remaining mechanisms involved in this short-wavelength instability are preserved in the limit of zero wavelength. Using the resulting simplified equations, we show that perturbations in different fluid particles interact only through the in-plane viscous diffusion which turns out to have a destabilising effect. We also show that in the presence of viscous diffusion, the closed trajectories which had been conjectured to play a crucial role in the onset of the mode B instability are not actually a prerequisite for the growth of mode B type perturbations. We combine these observations to identify the three essential ingredients for the development of the mode B instability: (i) the amplification of perturbations in the braid regions due to the stretching mechanism; and the spreading of perturbations through (ii) viscous diffusion, and (iii) cross-flow advection which transports fluid between the two braid regions on either side of the cylinder. Finally, we develop a simple criterion that allows the prediction of the regions where three-dimensional short-wavelength perturbations are amplified by the stretching mechanism. The approach used in our study is general and has the potential to give insights into the onset of three-dimensionality via short-wavelength instabilities in other flows.

In this paper it is shown that a modal detuned instability of periodic near-wall streaks originates a large-scale structure in the bulk of the turbulent channel flow. The effect of incoherent turbulent fluctuations is included in the linear operator by means of an eddy viscosity. The base flow is an array of periodic two-dimensional streaks, extracted from numerical simulations in small domains, superposed to the turbulent mean profile. The stability problem for a large number of periodic units is efficiently solved using the block-circulant matrix method proposed by Schmid et al. (Phys. Rev. Fluids, vol. 2, 2017, 113902). For friction Reynolds numbers equal or higher than $590$, it is shown that an unstable branch is present in the eigenspectra. The most unstable eigenmodes display large-scale modulations whose characteristic wavelengths are compatible with the large-scale end of the premultiplied velocity fluctuation spectra reported in previous computational studies. The wall-normal location of the large-wavelength near-wall peak in the spanwise spectrum of the eigenmode exhibits a power-law dependence on the friction Reynolds number, similarly to that found in experiments of pipes and boundary layers. Lastly, the shape of the eigenmode in the streamwise-wall-normal plane is reminiscent of the superstructures reported in the recent experiments of Deshpande et al. (J. Fluid Mech., vol. 969, 2023, A10). Therefore, there is evidence that such large-wavelength instabilities generate large-scale motions in wall-bounded turbulent flows.

The asymptotic analysis of steady azimuthally invariant electromagnetically driven flows occurring in a shallow annular layer of electrolyte undertaken in Part 1 of this study (McCloughan & Suslov, J. Fluid Mech., vol. 980, 2024, A59) predicted the existence of a two-tori flow state that has not been detected previously. In Part 2 of the study we confirm its existence by numerical time integration of the governing equations. We observe a hysteresis, where the type of solution obtained for the same set of governing parameters depends on the choice of the initial conditions and the way the governing parameters change, which is fully consistent with the analytic results of Part 1. Subsequently, we perform a linear stability analysis of the newly obtained steady state and deduce that the experimentally observed anti-cyclonic free-surface vortices appear on its background as a result of a centrifugal (Rayleigh-type) instability of the interface separating two counter-rotating toroidal structures that form the newly found flow solution. The quantitative characteristics of such instability structures are determined. It is shown that such structures can only exist in sufficiently thin layers with the depth not exceeding a certain critical value.

Resistive tearing instabilities are common in fluids that are highly electrically conductive and carry strong currents. We determine the effect of stable stratification on the tearing instability under the Boussinesq approximation. Our results generalise previous work that considered only specific parameter regimes, and we show that the length scale of the fastest-growing mode depends non-monotonically on the stratification strength. We confirm our analytical results by solving the linearised equations numerically, and we discuss whether the instability could operate in the solar tachocline.

This work is devoted to a theoretical and numerical study of the dynamics of a two-phase system vapour bubble in equilibrium with its liquid phase under translational vibrations in the absence of gravity. The bubble is initially located in the container centre. The liquid and vapour phases are considered as viscous and incompressible. Analysis focuses on the vibrational conditions used in experiments with the two-phase system SF$_6$ in the MIR space station and with the two-phase system para-Hydrogen (p-H$_2$) under magnetic compensation of Earth's gravity. These conditions correspond to small-amplitude high-frequency vibrations. Under vibrations, additionally to the forced oscillations, an average displacement of the bubble to the wall is observed due to an average vibrational attraction force related to the Bernoulli effect. Vibrational conditions for SF$_6$ correspond to much smaller average vibrational force (weak vibrations) than for p-H$_2$ (strong vibrations). For weak vibrations, the role of the initial vibration phase is crucial. The difference in the behaviour at different initial phases is explained using a simple mechanical model. For strong vibrations, the average displacement to the wall stops when the bubble reaches a quasi-equilibrium position where the resulting average force is zero. At large vibration velocity amplitudes this position is near the wall where the bubble performs only forced oscillations. At moderate vibration velocity amplitudes the bubble average displacement stops at a finite distance from the wall, then large-scale damped oscillations around this position accompanied by forced oscillations are observed. Bubble shape oscillations and the parametric resonance of forced oscillations are also studied.