In a vertical channel driven by an imposed horizontal temperature gradient, numerical simulations (Gao et al., Phys. Rev. E, vol. 88, 2013, 023010; Phys. Rev. E, vol. 91, 2015, 013006; Phys. Rev. E, vol. 97, 2018, 053107) have previously shown steady, time-periodic and chaotic dynamics. We explore the observed dynamics by constructing invariant solutions of the three-dimensional Oberbeck–Boussinesq equations, characterizing the stability of these equilibria and periodic orbits, and following the bifurcation structure of the solution branches under parametric continuation in Rayleigh number. We find that in a narrow vertically periodic domain of aspect ratio 10, the flow is dominated by the competition between three and four co-rotating rolls. We demonstrate that branches of three- and four-roll equilibria are connected and can be understood in terms of their discrete symmetries. Specifically, the $D_4$ symmetry of the four-roll branch dictates the existence of qualitatively different intermediate branches that themselves connect to the three-roll branch in a transcritical bifurcation due to $D_3$ symmetry. The physical appearance, disappearance, merging and splitting of rolls along the connecting branch provide a physical and phenomenological illustration of the equivariant theory of $D_3$–$D_4$ mode interaction. We observe other manifestations of the competition between three and four rolls, in which the symmetry in time or in the transverse direction is broken, leading to limit cycles or wavy rolls, respectively. Our work highlights the interest of combining numerical simulations, bifurcation theory and group theory, in order to understand the transitions between and origin of flow patterns.

Vertical thermal convection is a non-equilibrium system in which both buoyancy and shear forces play a role in driving the convective flow. Beyond the onset of convection, the driven dissipative system exhibits chaotic dynamics and turbulence. In a three-dimensional domain extended in both the vertical and the transverse dimensions, Gao et al. (Phys. Rev. E, vol. 97, 2018, 053107) have observed a variety of convection patterns which are not described by linear stability analysis. We investigate the fully nonlinear dynamics of vertical convection using a dynamical-systems approach based on the Oberbeck–Boussinesq equations. We compute the invariant solutions of these equations and the bifurcations that are responsible for the creation and termination of various branches. We map out a sequence of local bifurcations from the laminar base state, including simultaneous bifurcations involving patterned steady states with different symmetries. This atypical phenomenon of multiple branches simultaneously bifurcating from a single parent branch is explained by the role of $D_4$ symmetry. In addition, two global bifurcations are identified: first, a homoclinic cycle from modulated transverse rolls and second, a heteroclinic cycle linking two symmetry-related diamond-roll patterns. These are confirmed by phase space projections as well as the functional form of the divergence of the period close to the bifurcation points. The heteroclinic orbit is shown to be robust and to result from a 1:2 mode interaction. The intricacy of this bifurcation diagram highlights the essential role played by dynamical systems theory and computation in hydrodynamic configurations.

The nonlinear evolution of free-stream vortical disturbances entrained in the entrance region of a circular pipe is investigated using asymptotic and numerical methods. Attention is focused on the low-frequency disturbances that induce streamwise elongated structures. A pair of vortical modes with opposite azimuthal wavenumbers is used to model the free-stream disturbances. Their amplitude is assumed to be intense enough for nonlinear interactions to occur inside the pipe. The formation and evolution of the perturbation flow are described by the nonlinear unsteady boundary-region equations in the cylindrical coordinate system, derived and solved herein for the first time. Matched asymptotic expansions are employed to construct appropriate initial conditions and the initial–boundary value problem is solved numerically by a marching procedure in the streamwise direction. Numerical results show the stabilising effect of nonlinearity on the intense algebraic growth of the disturbances and an increase of the wall-shear stress due to the nonlinear interactions. A parametric study is carried out to evince the effect of the Reynolds number, the streamwise and azimuthal wavelengths, and the radial length scale of the inlet disturbance on the nonlinear flow evolution. Elongated pipe-entrance nonlinear structures (EPENS) occupying the whole pipe cross-section are discovered. EPENS with $h$-fold rotational symmetry comprise $h$ high-speed streaks positioned near the wall, and $h$ low-speed streaks centred around the pipe core. These distinct structures display a striking resemblance to nonlinear travelling waves found numerically and observed experimentally in fully developed pipe flow. Good agreement of our mean-flow and root mean square data with experimental measurements is obtained.

Liquid bridges are formed when a flowing liquid interacts with multiple parallel fibres, as relevant to heat and mass transfer applications that utilize flow down fibre arrays. We perform a comprehensive experimental study of flowing liquid bridges between two vertical fibres whose spacing is controlled dynamically in our experimental apparatus. The bridge patterns exhibit a regular periodic spacing typical of absolute instability for low flow rates, but become spatially inhomogeneous above a critical flow rate where the base flow is convectively unstable. The shapes of individual bridges and their associated dynamics are measured, as they depend upon the liquid properties, and fibre geometry/spacing. The bridge length scales similarly to static bridges between parallel fibres. The bridge dynamics exhibits a dependence on viscosity and scale with the impedance. A simple energy balance is used to derive a scaling relationship for the bridge velocity that captures the general trend of our experimental data. Finally, we demonstrate that these scalings similarly apply when the fibres are dynamically separated or brought together.

The first experimental results on pattern transitions in the co-rotation regime (i.e. the rotation ratio $\varOmega = \omega _o/\omega _i > 0$, where $\omega _i$ and $\omega _o$ are the angular speeds of the inner and outer cylinders, respectively) of the Taylor–Couette flow (TCF) are reported for a neutrally buoyant suspension of non-colloidal particles, up to a particle volume fraction of $\phi = 0.3$. While the stationary Taylor vortex flow (TVF) is the primary bifurcating state in dilute suspensions ($\phi \leq ~0.05$), the non-axisymmetric oscillatory states, such as the spiral vortex flow (SVF) and the ribbon (RIB), appear as primary bifurcations with increasing particle loading, with an overall de-stabilization of the primary bifurcating states (TVF/SVF/RIB) being found with increasing $\phi$ for all $\varOmega \geq ~0$. At small co-rotations ($\varOmega \sim 0$), the particles play the dual role of stabilization ($\phi < 0.1$) and destabilization ($\phi \geq ~0.1$) on the secondary/tertiary oscillatory states. The distinctive features of the ‘particle-induced’ spiral vortices are identified and contrasted with those of the ‘fluid-induced’ spirals that operate in the counter-rotation regime.

Symmetry-breaking bifurcations, where a flow state with a certain symmetry undergoes a transition to a state with a different symmetry, are ubiquitous in fluid mechanics. Much can be understood about the nature of these transitions from symmetry alone, using the theory of groups and their representations. Here, we show how the extensive databases on groups in crystallography can be exploited to yield insights into fluid dynamical problems. In particular, we demonstrate the application of the crystallographic layer groups to problems in fluid layers, using thermal convection as an example. Crystallographic notation provides a concise and unambiguous description of the symmetries involved, and we advocate its broader use by the fluid dynamics community.

Although the roll-streak (R-S) is fundamentally involved in the dynamics of wall turbulence, the physical mechanism responsible for its formation and maintenance remains controversial. In this work we investigate the dynamics maintaining the R-S in turbulent Poiseuille flow at $R=1650$. Spanwise collocation is used to remove spanwise displacement of the streaks and associated flow components, which isolates the streamwise-mean flow R-S component and the second-order statistics of the streamwise-varying fluctuations that are collocated with the R-S. This partition of the dynamics into streamwise-mean and fluctuation components facilitates exploiting insights gained from the analytic characterization of turbulence in the second-order statistical state dynamics (SSD), referred to as S3T, and its closely associated restricted nonlinear dynamics (RNL) approximation. Symmetry of the statistics about the streak centreline permits separation of the fluctuations into sinuous and varicose components. The Reynolds stress forcing induced by the sinuous and varicose fluctuations acting on the R-S is shown to reinforce low- and high-speed streaks, respectively. This targeted reinforcement of streaks by the Reynolds stresses occurs continuously as the fluctuation field is strained by the streamwise-mean streak and not intermittently as would be associated with streak-breakdown events. The Reynolds stresses maintaining the streamwise-mean roll arise primarily from the dominant proper orthogonal decomposition (POD) modes of the fluctuations, which can be identified with the time average structure of optimal perturbations growing on the streak. These results are consistent with a universal process of R-S growth and maintenance in turbulent shear flow arising from roll forcing generated by straining turbulent fluctuations, which was identified using the S3T SSD.

Snap-through is a buckling instability that allows slender objects, including those in plant and biological systems, to generate rapid motion that would be impossible if they were to use their internal forces exclusively. In microfluidic devices, such as micromechanical switches and pumps, this phenomenon has practical applications for manipulating fluids at small scales. The onset of this elastic instability often drives the surrounding fluid into motion – a process known as snap-induced flow. To analyse the complex dynamics resulting from the interaction between a sheet and a fluid, we develop a prototypical model of a thin sheet that is compressed between the two sides of a closed channel filled with an inviscid fluid. At first, the sheet bends towards the upstream direction and the system is at rest. However, once the pressure difference in the channel exceeds a critical value, the sheet snaps to the opposite side and drives the fluid dynamics. We formulate an analytical model that combines the elasticity of thin sheets with the hydrodynamics of inviscid fluids to explore how external pressure differences, material properties and geometric factors influence the system's behaviour. To analyse the early stages of the evolution, we perform a linear stability analysis and obtain the growth rate and the critical pressure difference for the onset of the instability. A weakly nonlinear analysis suggests that the system can exhibit a pressure spike in the vicinity of the inverted configuration.

A dense particle suspension under shear may lose its uniform state to large local density and stress fluctuations, which challenge the mean-field description of the system. Here, we explore the novel dynamics of a non-Brownian suspension under orbital oscillations, where localized density waves along the flow direction appear beyond an excitation frequency threshold and self-organize into a hexagonal pattern across the system. The spontaneous occurrence of the inhomogeneity pattern arises from a coupling between particle advection and the shear-thickening nature of the suspension. Through linear stability analysis, we show that they overcome the stabilizing effects of particle pressure at sufficient particle volume fraction and oscillation frequency. In addition, the long-standing density waves degenerate into random fluctuations when replacing the free surface with rigid confinement. It indicates that the shear-thickened state is intrinsically heterogeneous, and the boundary conditions are crucial for developing local disturbance.

We study structure formation in two-dimensional turbulence driven by an external force, interpolating between linear instability forcing and random stirring, subject to nonlinear damping. Using extensive direct numerical simulations, we uncover a rich parameter space featuring four distinct branches of stationary solutions: large-scale vortices, hybrid states with embedded shielded vortices (SVs) of either sign, and two states composed of many similar SVs. Of the latter, the first is a dense vortex gas where all SVs have the same sign and diffuse across the domain. The second is a hexagonal vortex crystal forming from this gas when the linear instability is sufficiently weak. These solutions coexist stably over a wide parameter range. The late-time evolution of the system from small-amplitude initial conditions is nearly self-similar, involving three phases: initial inverse cascade, random nucleation of SVs from turbulence and, once a critical number of vortices is reached, a phase of explosive nucleation of SVs, leading to a statistically stationary state. The vortex gas is continued in the forcing parameter, revealing a sharp transition towards the crystal state as the forcing strength decreases. This transition is analysed in terms of the diffusivity of individual vortices using ideas from statistical physics. The crystal can also decay via an inverse cascade resulting from the breakdown of shielding or insufficient nonlinear damping acting on SVs. Our study highlights the importance of the forcing details in two-dimensional turbulence and reveals the presence of non-trivial SV states in this system, specifically the emergence and melting of a vortex crystal.

In this paper, we analyse Turing instability and bifurcations in a host–parasitoid model with nonlocal effect. For a ordinary differential equation model, we provide some preliminary analysis on Hopf bifurcation. For a reaction–diffusion model with local intraspecific prey competition, we first explore the Turing instability of spatially homogeneous steady states. Next, we show that the model can undergo Hopf bifurcation and Turing–Hopf bifurcation, and find that a pair of spatially nonhomogeneous periodic solutions is stable for a (8,0)-mode Turing–Hopf bifurcation and unstable for a (3,0)-mode Turing–Hopf bifurcation. For a reaction–diffusion model with nonlocal intraspecific prey competition, we study the existence of the Hopf bifurcation, double-Hopf bifurcation, Turing bifurcation, and Turing–Hopf bifurcation successively, and find that a spatially nonhomogeneous quasi-periodic solution is unstable for a (0,1)-mode double-Hopf bifurcation. Our results indicate that the model exhibits complex pattern formations, including transient states, monostability, bistability, and tristability. Finally, numerical simulations are provided to illustrate complex dynamics and verify our theoretical results.

Supergranule aggregation, i.e. the gradual aggregation of convection cells to horizontally extended networks of flow structures, is a unique feature of constant heat flux-driven turbulent convection. In the present study, we address the question if this mechanism of self-organisation of the flow is present for any fluid. Therefore, we analyse three-dimensional Rayleigh–Bénard convection at a fixed Rayleigh number ${Ra} \approx 2.0 \times 10^{5}$ across $4$ orders of Prandtl numbers ${Pr} \in [10^{-2}, 10^{2}]$ by means of direct numerical simulations in horizontally extended periodic domains with aspect ratio $\varGamma = 60$. Our study confirms the omnipresence of the mechanism of supergranule aggregation for the entire range of investigated fluids. Moreover, we analyse the effect of ${Pr}$ on the global heat and momentum transport, and clarify the role of a potential stable stratification in the bulk of the fluid layer. The ubiquity of the investigated mechanism of flow self-organisation underlines its relevance for pattern formation in geophysical and astrophysical convection flows, the latter of which are often driven by prescribed heat fluxes.

In this paper a fully coupled particle-size segregation model for granular flows (Barker et al., J. Fluid Mech., vol. 909, 2021, p. A22) is used to simulate the development of the patterns in a triangular rotating drum. The results are compared with the experimental patterns formed with bidisperse and tridisperse granular mixtures, and with varying compositions and fill heights. In all cases the agreement between the simulations and experiments is remarkably good. The experimental patterns are generated in a narrow gap between transparent front and back sidewalls. These prevent three-dimensional motion, but also impose friction on the flow, making it thinner and faster than it would otherwise be. This promotes segregation, as it simultaneously increases the shear rate and reduces the local pressure. To obtain the correct flow dynamics and segregation, width-averaged sidewall friction is incorporated into the two-dimensional simulations, which are performed in OpenFOAM$^{\circledR}$. The free-surface avalanche forms a boundary layer within which all the segregation occurs. Material in the lower reach of the avalanche is continuously deposited into an underlying solid body of grains, which rotates with the drum, and is eventually re-entrained into the avalanche along its upper reach. The changing geometry of the granular region (as the drum rotates) implies that the avalanche is constantly adjusting its length, position and depth. This generates a complex quasi-periodic flow, which when combined with particle-size segregation generates amazing patterns in the solid rotating granular body after only two drum rotations.

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

An experiment on convective flows induced by the dielectrophoretic force was performed under the microgravity condition provided during a sounding rocket flight. The dielectrophoretic force possesses a non-conservative term that can be seen as resulting from an electric gravity. That gravity can be responsible for an electric Rayleigh–Bénard convection between a hot inner cylinder and a cold outer cylinder when an electric field is applied in the radial direction. Four cells with independent temperature and electric field controls allowed the investigation of eight different values of the electric Rayleigh number relatively close to the onset of the thermo-electric instability. A linear stability analysis is performed to predict the stability threshold and the evolution of the growth rate of the instability. The three-dimensional structure of the flow is captured by simultaneous particle image velocimetry and by shadowgraphy. The amplitude of the instability modes and the time evolution of the flow is analysed, and various methods are proposed to extrapolate the experimental critical value of the electric Rayleigh number for the onset of convection. The measured critical electric Rayleigh number is in agreement with the prediction of the linear stability theory. The comparison of the new experimental results with previous ones from parabolic flight campaigns highlights the importance of long-term microgravity for the achievement of thermal convection at low values of the control parameters.

This paper reports several new classes of unstable recurrent solutions of the two-dimensional Euler equation on a square domain with periodic boundary conditions. These solutions are in many ways analogous to recurrent solutions of the Navier–Stokes equation which are often referred to as exact coherent structures. In particular, we find that recurrent solutions of the Euler equation are dynamically relevant: they faithfully reproduce large-scale flows in simulations of turbulence at very high Reynolds numbers. On the other hand, these solutions have a number of properties which distinguish them from their Navier–Stokes counterparts. First of all, recurrent solutions of the Euler equation come in infinite-dimensional continuous families. Second, solutions of different types are connected, e.g. an equilibrium can be smoothly continued to a travelling wave or a time-periodic state. Third, and most important, they are only weakly unstable and, as a result, fully developed turbulence mimics some of these solutions remarkably frequently and over unexpectedly long temporal intervals.

We develop a frame-invariant theory of material spike formation during flow separation over a no-slip boundary in three-dimensional flows with arbitrary time dependence. Based on the exact evolution of the largest principal curvature on near-wall material surfaces, our theory identifies fixed and moving separation. Our approach is effective over short time intervals and admits an instantaneous limit. As a byproduct, we derive explicit formulas for the evolution of the Weingarten map and the principal curvatures of any surface advected by general three-dimensional flows. The material backbone we identify acts first as a precursor and later as the centrepiece of Lagrangian flow separation. We discover previously undetected spiking points and curves where the separation backbones connect to the boundary and provide wall-based analytical formulas for their locations. We illustrate our results on several steady and unsteady flows.

Our understanding of the material organization of complex fluid flows has benefited recently from mathematical developments in the theory of objective coherent structures. These methods have provided a wealth of approaches that identify transport barriers in three-dimensional (3-D) turbulent flows. Specifically, theoretical advances have been incorporated into numerical algorithms that extract the most influential advective, diffusive and active barriers to transport from data sets in a frame-indifferent fashion. To date, however, there has been very limited investigation into these objectively defined transport barriers in 3-D unsteady flows with complicated spatiotemporal dynamics. Similarly, no systematic comparison of advective, diffusive and active barriers has been carried out in a 3-D flow with both thermally driven and mechanically modified structures. In our study, we utilize simulations of turbulent rotating Rayleigh–Bénard convection to uncover the interplay between advective transport barriers (Lagrangian coherent structures), material barriers to diffusive heat transport, and objective Eulerian barriers to momentum transport. For a range of (inverse) Rossby numbers, we identify each type of barrier and find intriguing relationships between momentum and heat transport that can be related to changes in the relative influence of mechanical and thermal forces. Further connections between bulk behaviours and structure-specific behaviours are also developed.

Thermoacoustic systems are complex systems where the interactions between the hydrodynamic, acoustic and heat release rate fluctuations lead to diverse dynamics such as chaos, intermittency and ordered dynamics. Such complex interactions cause catastrophically high-amplitude acoustic pressure oscillations and the emergence of order in the spatiotemporal dynamics, referred to as thermoacoustic instability. In this work, we use multilayer networks to study the spatial pattern of inter-subsystem interactions between the vorticity dynamics and thermoacoustic power generated due to acoustically coupled combustion in a bluff-body-stabilised turbulent dump combustor. We construct a two-layered network where the layers represent the thermoacoustic power and vorticity fields. The inter-layer links are determined using cross-variable short-window correlations between vorticity and thermoacoustic power fluctuations at any two locations in the flow field. Analysing the topology of inter-layer networks, using network properties such as degree correlations and link-rank distributions, helps us infer the spatial inhomogeneities in inter-subsystem interactions and unravel the fluid mechanical processes involved during different dynamical states. We show that, during chaotic dynamics, interactions between subsystems are non-localised and spread throughout the flow field of the combustor. During the state of thermoacoustic instability (order), we find that intense interactions occur in between regions of coherent vortex shedding and thermoacoustic power generation and we understand that these processes are strongly and locally coupled. Moreover, we discover that such dense inter-layer connections emerge in spatial pockets in the dump plane of the combustor during the state of intermittency much prior to the onset of order.

Low Reynolds number turbulence in wall-bounded shear flows en route to laminar flow takes the form of spatially intermittent turbulent structures. In plane shear flows, these appear as a regular pattern of alternating turbulent and quasi-laminar flow. Both the physical and the spectral energy balance of a turbulent–laminar pattern in plane Couette flow are computed and compared to those of uniform turbulence. In the patterned state, the mean flow is strongly modulated and is fuelled by two mechanisms: primarily, the nonlinear self-interaction of the mean flow (via mean advection), and secondly, the extraction of energy from turbulent fluctuations (via negative spectral production, associated with an energy transfer from small to large scales). Negative production at large scales is also found in the uniformly turbulent state. Important features of the energy budgets are surveyed as a function of $Re$ through the transition between uniform turbulence and turbulent–laminar patterns.