Vertical thermal convection exhibits weak turbulence and spatio-temporally chaotic behaviour. For this configuration, we report seven new equilibria and 26 new periodic orbits. These orbits, together with four previously studied in Zheng et al. (J. Fluid Mech., 2024b, vol. 1000, p. A29) bring the number of periodic-orbit branches computed so far to 30, all solutions to the fully nonlinear three-dimensional Navier–Stokes equations. These new and unstable invariant solutions capture intricate spatio-temporal flow patterns including straight, oblique, wavy, skewed and distorted convection rolls, as well as bursts and defects. These interesting and important fluid mechanical processes in a small flow unit are shown to also appear locally and instantaneously in a chaotic simulation in a large domain. Most of the solution branches show rich spatial and/or spatio-temporal symmetries. The bifurcation-theoretic organisation of these solutions is discussed; the bifurcation scenarios include Hopf, pitchfork, saddle-node, period-doubling, period-halving, global homoclinic and heteroclinic bifurcations, as well as isolas. Furthermore, these orbits are shown to be able to reconstruct statistically the core part of the attractor, so that these results may contribute to a quantitative description of transitional fluid turbulence using periodic orbit theory.

We numerically investigate the cellular detonation dynamics in ethylene/oxygen/ozone/nitrogen mixtures considering detailed chemical kinetics. The aim is to elucidate emergent detonation structures and reveal the transition mechanism from single- to double-cellular structures. Ozone is used to induce two-stage reactions within the mixture. Through systematic initiation strength analysis, we demonstrate two distinct propagation regimes: (i) under strong initiation, a stable double-cellular detonation is established; (ii) weak initiation triggers a multi-stage evolutionary process, beginning with a low-speed single-cellular detonation in the initiation zone. During the initial weak stage, the detonation propagates at a quasi-steady velocity with uniform cellular patterning. The subsequent transition phase features spontaneous acceleration accompanied by structural bifurcation into double cells, ultimately stabilising in a normal stage with sustained double-cellular structures. Further analysis reveals that the weak-stage dynamics is governed exclusively by first-stage chemical reactions, resulting in a single-cellular structure propagating at a velocity much lower than the Chapman–Jouguet speed. In contrast, the double-cellular structure observed at the normal stage results from the two-stage exothermic reactions. Thermodynamic perturbations arising from cellular instability and fluid dynamic instability are identified as critical drivers for the transition from single- to double-cellular detonation. Besides, conditions for the formation of double-cellular detonation are explored, and two qualitative requirements are summarised: the reactions of the two stages must proceed as independently as possible, and both heat releases from the two stages must be high enough to sustain the triple-shock configurations.

Flag flutter frequently features a marked difference between the onset speed of flutter and the speed below which flutter stops. The hysteresis tends to be especially large in experiments as opposed to simulations. This phenomenon has been ascribed to inherent imperfections of flatness in experimental samples, which are thought to inhibit the onset of flutter but have a lesser effect once a flag is already fluttering. In this work, we present an experimental confirmation for this explanation through motion tracking. We also visualize the wake to assess the potential contribution of discrete vortex shedding to hysteresis. We then mould our understanding of the mechanism of bistability and additional observations on flag flutter into a novel, observation-based, semiempirical model for flag flutter in the form of a single ordinary differential equation. Despite its simplicity, the model successfully reproduces key features of the physical system such as bistability, sudden transitions between non-fluttering and fluttering states, amplitude growth and frequency growth.

Finite-amplitude spiral vortex flows are obtained numerically for the Taylor–Couette system in the narrow limit of the gap between two concentric rotating cylinders. These spiral vortex flows bifurcate from circular Couette flow before axisymmetric Taylor vortex flow sets in when the ratio $\mu$ of the angular velocities of the outer to the inner cylinder is less than −0.78, consistent with the results of linear stability analysis by Krueger et al. (J. Fluid Mech., vol. 24, 1966, pp. 521–538), while the boundary of existence of spiral vortex flows is determined not by the linear critical point, but by the saddle-node point of the subcritical spiral vortex flow branch for $\mu \lessapprox -0.75$, when the axial wavenumber $\beta =2.0$. It is found that the nonlinear spiral vortex flows exhibit the mean flow in the axial direction as well as in the azimuthal direction, and that the profiles of both mean-flow components are asymmetric about the centre plane between the gap.

The paper discusses the stochastic dynamics of the vortex shedding process in the presence of external harmonic excitation and coloured multiplicative noise. The situation is encountered in a turbulent practical combustor experiencing combustion instability. Acoustic feedback and turbulent flow are imitated by the harmonic and stochastic excitations, respectively. The Ornstein–Uhlenbeck process is used to generate the noise. A low-order model for vortex shedding is used. The Fokker–Planck framework is used to obtain the evolution of the probability density function of the shedding time period. Stochastic lock-in and resonance characteristics are studied for various parameters associated with the harmonic (amplitude, frequency) and noise (amplitude, correlation time, multiplicative noise factor) excitations. We observed that: (i) the stochastic lock-in (s-lock-in) boundary strongly depends on the noise correlation time; (ii) the parameter sites for s-lock-in can be approximately identified from the noise-induced shedding statistics; and (iii) stochastic resonance is significant for some intermediate correlation times. The effects of the above-mentioned observations are discussed in the context of combustion instability.

In this paper we propose a novel control strategy for modulating nonlinear flapping and symmetry-breaking (S-B) bifurcations of a piezoelectric metamaterial beam behind a circular cylinder subjected to viscous flow. The beam incorporates distributed piezoelectric meta-cells connected via unidirectional circuits to enable self-sensing and adaptive control. A strongly coupled nonlinear fluid-structure-electro-control model within an arbitrary Lagrangian–Eulerian framework is developed for predicting the flapping dynamics of the large deformable piezoelectric metamaterial beam. The system exhibits multiple flow-induced modes, including limit-cycle oscillations, subharmonic responses and S-B deflections. These dynamic regimes arise from nonlinear bifurcations of the system, namely the period-doubling and spontaneous S-B bifurcations. Flapping control and wake topology transition of the system is achieved by suppressing the periodic-doubling bifurcation based on the vibration rebound effect through a self-sensing and adaptive-actuation mechanism of the beam. Floquet stability analysis confirms the effectiveness of control in delaying instability onset and suppressing chaotic transitions. Symmetry modulation of the beam is achieved via the localised perturbations induced from the piezoelectric meta-cells, which reshape the stability of the system. The transition from S-B mode to symmetry-recovery mode reflects a shift from a flow-separation-dominated to vibration-dominated vortex shedding pattern. This symmetry transition reorganises the energy exchange pathways between the flow and the beam. Quantitative analyses of the wake recovery and the energy harvesting efficiency confirm enhanced flow energy conversion under control. These results establish a framework for bifurcation control of slender structures in viscous flow, providing potential applications for underwater energy harvesting and flexible propulsion in unsteady environments.

Equilibrium, travelling-wave and periodic-orbit solutions of the Navier–Stokes equations provide a promising avenue for investigating the structure, dynamics and statistics of transitional flows. Many such invariant solutions have been computed for wall-bounded shear flows, including plane Couette, plane Poiseuille and pipe flow. However, the organisation of invariant solutions is not well understood. In this paper we focus on the role of symmetries in the organisation and computation of invariant solutions of plane Poiseuille flow. We show that enforcing symmetries while computing invariant solutions increases the efficiency of the numerical methods, and that redundancies between search spaces can be eliminated by consideration of equivalence relations between symmetry subgroups. We determine all symmetry subgroups of plane Poiseuille flow in a doubly periodic domain up to translations by half the periodic lengths and classify the subgroups into equivalence classes, each of which represents a physically distinct set of symmetries and an associated set of physically distinct invariant solutions. We calculate fifteen new travelling waves of plane Poiseuille flow in seven distinct symmetry groups and discuss their relevance to the dynamics of transitional turbulence. We present a few examples of subgroups with fractional shifts other than half the periodic lengths and one travelling-wave solution whose symmetry involves shifts by one third of the periodic lengths. We conclude with a discussion and some open questions about the role of symmetry in the behaviour of shear flows.

A literature review suggests that the flows past simply connected bodies with aspect ratio close to unity and symmetries aligned with the flow follow a consistent sequence of regimes (steady, periodic, quasiperiodic) as the Reynolds number increases. However, evidence is fragmented, and studies are rarely conducted using comparable numerical or experimental set-ups. This paper investigates the wake dynamics of two canonical bluff bodies with distinct symmetries: a cube (discrete) and a sphere (continuous). Employing three-dimensional (3-D) global linear stability analysis and nonlinear simulations within a unified numerical framework, we identify the bifurcation sequence driving these regime transitions. The sequence: a pitchfork bifurcation breaks spatial symmetry; a Hopf bifurcation introduces temporal periodicity ($St_1$); a Neimark–Sacker bifurcation destabilises the periodic orbit, leading to quasiperiodic dynamics with two incommensurate frequencies ($St_1, St_2$). A Newton–Krylov method computes the unstable steady and periodic base flows without imposing symmetry constraints. Linear stability reveals similarities between the cube and sphere in the spatial structure of the leading eigenvectors and in the eigenvalue trajectories approaching instability. This study provides the first confirmation of a Neimark–Sacker bifurcation to quasiperiodicity in these 3-D wakes, using Floquet stability analysis of computed unstable periodic orbits and their Floquet modes. The quasiperiodic regime is described in space and time by the Floquet modes’ effects on the base flow and a spectrum dominated by the two incommensurate frequencies and tones arising from nonlinear interactions. Although demonstrated for a cube and a sphere, this bifurcation sequence, leading from steady state to quasiperiodic dynamics, suggests broader applicability beyond these geometries.

The non-uniform evaporation rate at the liquid–gas interface of binary droplets induces solutal Marangoni flows. In glycerol–water mixtures (positive Marangoni number, where the more volatile fluid has higher surface tension), these flows stabilise into steady patterns. Conversely, in water–ethanol mixtures (negative Marangoni number, where the less volatile fluid has higher surface tension), Marangoni instabilities emerge, producing seemingly chaotic flows. This behaviour arises from the opposing signs of the Marangoni number. Perturbations locally reducing surface tension at the interface drive Marangoni flows away from the perturbed region. Continuity of the fluid enforces a return flow, drawing fluid from the bulk towards the interface. In mixtures with a negative Marangoni number, preferential evaporation of the lower-surface-tension component leads to a higher concentration of the higher-surface-tension component at the interface as compared with the bulk. The return flow therefore creates a positive feedback loop, further reducing surface tension in the perturbed region and enhancing the instability. This study investigates bistable quasi-stationary solutions in evaporating binary droplets with negative Marangoni numbers (e.g. water–ethanol) and examines symmetry breaking across a range of Marangoni numbers and contact angles. Bistable domains exhibit hysteresis. Remarkably, flat droplets (small contact angles) show instabilities at much lower critical Marangoni numbers than droplets with larger contact angles. Our numerical simulations reveal that interactions between droplet height profiles and non-uniform evaporation rates trigger azimuthal Marangoni instabilities in flat droplets. This geometrically confined instability can even destabilise mixtures with positive Marangoni numbers, particularly for concave liquid–gas interfaces, as in wells. Finally, through a Lyapunov exponent analysis, we confirm the chaotic nature of flows in droplets with a negative Marangoni number. We emphasise that the numerical models are intentionally simplified to isolate and clarify the underlying mechanisms, rather than to quantitatively predict specific experimental outcomes; in particular, the model becomes increasingly limited in regimes of rapid evaporation.

Flame–flame interactions in continuous combustion systems can induce a range of nonlinear dynamical behaviours, particularly in the thermoacoustic context. This study examines the mutual coupling and synchronisation dynamics of two thermoacoustic oscillators in a model gas-turbine combustor operating within a stochastic environment and subjected to external sinusoidal forcing. Experimental observations from two flames in an annular combustor reveal the emergence of dissimilar limit cycles, indicating localised lock-in of thermoacoustic oscillators. To interpret these dynamics, we introduce a coupled stochastic oscillator model with sinusoidal forcing terms, which highlights the critical role of individual synchronisation in enabling local lock-in. Furthermore, through stochastic system identification using this phenomenological low-order model, we mathematically demonstrate that a transition towards self-sustained oscillations can be driven solely by enhanced mutual coupling under external forcing. This combined experimental and modelling effort offers a novel framework for characterising complex coupled flame dynamics in practical combustion systems.

A rotating detonation combustor exhibits corotating $N$-wave modes with $N$ detonation waves propagating in the same direction. These modes and their responses to ignition conditions and disturbances were studied using a surrogate model. Through numerical continuation, a mode curve (MC) is obtained, depicting the relationship between the wave speed of the one-wave mode and a defined baseline of the combustor circumference ($L_{{base}}$) under fixed equation parameters, limited by deflagration and flow choking. The modes’ existence is confirmed by the equivalence between a one-wave mode within a combustor with circumference $L_{{base}}$/$N$ on the MC and an $N$-wave mode in an $L_{{base}}$ combustor. The stability, measured by the real part of the eigenvalue from linear stability analysis (LSA), revealed the dynamic properties. When multiple stable modes exist under the same parameters, ignition conditions with a spatial period of $L_{{base}}$/$N$ are more likely to form $N$-wave modes. An unstable evolution in formed modes, occurs in the dynamics from stable to unstable modes through saddle-node bifurcation and Hopf bifurcation induced by parameter perturbations and from unstable to stable modes induced by state disturbances. Eigenmodes from LSA reveal mechanisms of the unstable evolution, including the effect of secondary deflagration in the unstable one-wave mode and competitive interaction between detonation waves in the unstable multiwave mode, crucial for the combustor to mode transition.

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.

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.

When an evaporating water droplet is deposited on a thermally conductive substrate, the minimum temperature will be at the apex due to evaporative cooling. Consequently, density and surface tension gradients emerge within the droplet and at the droplet–gas interface, giving rise to competing flows from, respectively, the apex towards the contact line (thermal-buoyancy-driven flow) and the other way around (thermal Marangoni flow). In small droplets with diameter below the capillary length, the thermal Marangoni effects are expected to dominate over thermal buoyancy (‘thermal Rayleigh’) effects. However, contrary to these theoretical predictions, our experiments show mostly a dominant circulation from the apex towards the contact line, indicating a prevailing of thermal Rayleigh convection. Furthermore, our experiments often show an unexpected asymmetric flow that persisted for several minutes. We hypothesise that a tiny amount of contaminants, commonly encountered in experiments with water/air interfaces, act as surfactants and counteract the thermal surface tension gradients at the interface and thereby promote the dominance of Rayleigh convection. Our finite element numerical simulations demonstrate that under our specified experimental conditions, a mere 0.5 % reduction in the static surface tension caused by surfactants leads to a reversal in the flow direction, compared to the theoretical prediction without contaminants. Additionally, we investigate the linear stability of the axisymmetric solutions, revealing that the presence of surfactants also affects the axial symmetry of the flow.

In this work, a systematic study is carried out concerning the dynamic behaviour of finite-size spheroidal particles in non-isothermal shear flows between parallel plates. The simulations rely on a hybrid method combining the lattice Boltzmann method with a finite-difference solver. Fluid–particle and heat–particle interactions are accounted for by using the immersed boundary method. The effect of particle Reynolds number ($\textit{Re}_p=1{-}90$), Grashof number (${Gr}=0{-}200$), initial position and initial orientation of the particle are thoroughly examined. For the isothermal prolate particle, we observed that above a certain Reynolds number, the particle undergoes a pitchfork bifurcation; at an even higher Reynolds number, it returns to the centre position. In contrast, the hot particle behaves differently, with no pitchfork bifurcation. Instead, the Reynolds and Grashof numbers can induce oscillatory tumbling or log-rolling motions in either the lower or upper half of the channel. Heat transfer also plays an important role: at low Grashof numbers, the particle settles near the lower wall, while increasing the Grashof number shifts it towards the upper side. Moreover, the presence of thermal convection increases the rotational speed of the particle. Surprisingly, beyond the first critical Reynolds number, the equilibrium position of the thermal particle shifts closer to the centreline compared with that of a neutrally buoyant isothermal particle. Moreover, higher Grashof numbers can cause the particle to transition from tumbling to log-rolling or even a no-rotation mode. The initial orientation has a stronger influence at low Grashof numbers, while the initial position shows no strong effect in non-isothermal cases.

A nonlinear Schrödinger equation for pure capillary waves propagating at the free surface of a vertically sheared current has been used to study the stability and bifurcation of capillary Stokes waves on arbitrary depth. A linear stability analysis of weakly nonlinear capillary Stokes waves on arbitrary depth has shown that (i) the growth rate of modulational instability increases as the vorticity decreases whatever the dispersive parameter $kh$, where $k$ is the carrier wavenumber and $h$ the depth; (ii) the growth rate is significantly amplified for shallow water depths; and (iii) the instability bandwidth widens as the vorticity decreases. Particular attention has been paid to damping due to viscosity and forcing effects on modulational instability. In addition, a linear stability analysis to transverse perturbations in deep water has been carried out, demonstrating that the dominant modulational instability is two-dimensional whatever the vorticity. Near the minimum of linear phase velocity in deep water, we have shown that generalised capillary solitary waves bifurcate from linear capillary Stokes waves when the vorticity is positive. Moreover, we have shown that the envelope of pure capillary waves in deep water is unstable to transverse perturbations. Consequently, deep-water generalised capillary solitary waves are expected to be unstable to transverse perturbations.

This paper investigates the linear and nonlinear dynamics of two-dimensional penetrative convection subjected to radiative volumetric thermal forcing, focusing on ice-covered freshwater systems. Linear stability analysis reveals how critical wavenumbers $k_c$ and Rayleigh numbers $Ra_c$ are influenced by the attenuation lengths and incoming heat flux. In this configuration, the system easily becomes unstable with a small $Ra_c$, which is two decades smaller than that of the classical Rayleigh–Bénard convection problem, with typically $O(10)$. Weakly nonlinear analysis figures out that this configuration is supercritical, contrasting with the subcritical case by Veronis (Astrophys. J., vol. 137, 1963, 641–663). Numerical bifurcation solutions are performed from the critical points and over several decades, up to $Ra \sim O(10^6)$. This paper found that the system exhibits multiple steady solutions, and under certain specific conditions, a staircase temperature profile emerges. Meanwhile, we further discuss the influence of incoming heat flux and the Prandtl number $Pr$ on the primary bifurcation. Direct numerical simulations are also carried out, showing that heat is transported more efficiently via unsteady convection.

The critical points of vorticity in a two-dimensional viscous flow are essential for identifying coherent structures in the vorticity field. Their bifurcations as time progresses can be associated with the creation, destruction or merging of vortices, and we analyse these processes using the equation of motion for these points. The equation decomposes the velocity of a critical point into advection with the fluid and a drift proportional to viscosity. Conditions for the drift to be small or vanish are derived, and the analysis is extended to cover bifurcations. We analyse the dynamics of vorticity extrema in numerical simulations of merging of two identical vortices at Reynolds numbers ranging from 5 to 1500 in the light of the theory. We show that different phases of the merging process can be identified on the basis of the balance between advection and drift of the critical points, and identify two types of merging, one for low and one for high values of the Reynolds number. In addition to local maxima of positive vorticity and minima of negative vorticity, which can be considered centres of vortices, minima of positive vorticity and maxima of negative vorticity can also exist. We find that such anti-vortices occur in the merging process at high Reynolds numbers, and discuss their dynamics.

Dean’s approximation for curved pipe flow, valid under loose coiling and high Reynolds numbers, is extended to study three-dimensional travelling waves. Two distinct types of solutions bifurcate from the Dean’s classic two-vortex solution. The first type arises through a supercritical bifurcation from inviscid linear instability, and the corresponding self-consistent asymptotic structure aligns with the vortex–wave interaction theory. The second type emerges from a subcritical bifurcation by curvature-induced instabilities and satisfies the boundary region equations. A connection to the zero-curvature limit was not found. However, by continuing from known self-sustained exact coherent structures in the straight pipe flow problem, another family of three-dimensional travelling waves can be shown to exist across all Dean numbers. The self-sustained solutions also possess the two high-Reynolds-number limits. While the vortex–wave interaction type of solutions can be computed at large Dean numbers, their branch remains unconnected to the Dean vortex solution branch.

The nonlinear stability of two-dimensional (2-D) plane Couette flow subject to a constant throughflow is analysed at finite and asymptotically large Reynolds numbers $\textit {Re}$. The speed of this throughflow is quantified by the non-dimensional throughflow number $\eta$. The base flow exhibits a linear instability provided $\eta \gtrsim 3.35$, with multi-deck upper and lower branch structures developing in the limit $1\ll \eta \ll \mathit {O}(\textit {Re})$. This instability provides a springboard for the computation of nonlinear travelling waves which bifurcate subcritically from the linear neutral curve, allowing us to map out a neutral surface at different values of $\eta$. Using strongly nonlinear critical layer theory, we investigate the waves that bifurcate from the upper branch at asymptotically large $\textit {Re}$. This asymptotic structure exists provided the throughflow number is larger than the critical value of $\eta _c\approx 1.20$ and is shown to give quantitatively similar results to the numerical solutions at Reynolds numbers of $\mathit {O}(10^5)$.