The emergence of large-scale spatial modulations of turbulent channel flow, as the Reynolds number is decreased, is addressed numerically using the framework of linear stability analysis. Such modulations are known as the precursors of laminar–turbulent patterns found near the onset of relaminarisation. A synthetic two-dimensional base flow is constructed by adding finite-amplitude streaks to the turbulent mean flow. The streak mode is chosen as the leading resolvent mode from linear response theory. In addition, turbulent fluctuations can be taken into account or not by using a simple Cess eddy viscosity model. The linear stability of the base flow is considered by searching for unstable eigenmodes with wavelengths larger than the base flow streaks. As the streak amplitude is increased in the presence of the turbulent closure, the base flow loses its stability to a large-scale modulation below a critical Reynolds-number value. The structure of the corresponding eigenmode, its critical Reynolds number, its critical angle and its wavelengths are all fully consistent with the onset of turbulent modulations from the literature. The existence of a threshold value of the Reynolds number is directly related to the presence of an eddy viscosity, and is justified using an energy budget. The values of the critical streak amplitudes are discussed in relation with those relevant to turbulent flows.

The Cahn–Hilliard model with reaction terms can lead to situations in which no coarsening is taking place and, in contrast, growth and division of droplets occur which all do not grow larger than a certain size. This phenomenon has been suggested as a model for protocells, and a model based on the modified Cahn–Hilliard equation has been formulated. We introduce this equation and show the existence and uniqueness of solutions. Then, formally matched asymptotic expansions are used to identify a sharp interface limit using a scaling of the reaction term, which becomes singular when the interfacial thickness tends to zero. We compute planar solutions and study their stability under non-planar perturbations. Numerical computations for the suggested model are used to validate the sharp interface asymptotics. In addition, the numerical simulations show that the reaction terms lead to diverse phenomena such as growth and division of droplets in the obtained solutions, as well as the formation of shell-like structures.

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.

Interfaces subjected to strong time-periodic horizontal accelerations exhibit striking patterns known as frozen waves. In this study, we experimentally and numerically investigate the formation of such structures in immiscible fluids under high-frequency forcing. In the inertial regime – characterised by large Reynolds and Weber numbers, where viscous and surface tension effects become negligible – we demonstrate that the amplitude of frozen waves scales proportionally with the square of the forcing velocity. These results are consistent with vibro-equilibria theory and extend the theoretical framework proposed by Gréa & Briard (2019 Phys. Rev. Fluids 4, 064608) to immiscible fluids with large density contrasts. Furthermore, we examine the influence of both Reynolds and Weber numbers, not only in the onset of secondary Faraday instabilities – which drive the transition of frozen wave patterns toward a homogenised turbulent state – but also in selecting the dominant wavelength in the final saturated regime.

Understanding the interplay between thermal, elastic and hydrodynamic effects is crucial for a variety of applications, including the design of soft materials and microfluidic systems. Motivated by these applications, we investigate the emergence of natural convection in a fluid layer that is supported from below by a rigid surface, and covered from above by a thin elastic sheet. The sheet is laterally compressed and is maintained at a constant temperature lower than that of the rigid surface. We show that for very stiff sheets, and below a certain magnitude of the lateral compression, the system behaves as if the fluid were confined between two rigid walls, where the emergent flow exhibits a periodic structure of vortices with a typical length scale proportional to the depth of the fluid, similar to patterns observed in Rayleigh–Bénard convection. However, for more compliant sheets, and above a certain threshold of the lateral compression, a new local minimum appears in the stability diagram, with a corresponding wavenumber that depends solely on the bending modulus of the sheet and the specific weight of the fluid, as in wrinkling instability of thin sheets. The emergent flow field in this region synchronises with the wrinkle pattern. We investigate the exchange of stabilities between these two solutions, and construct a stability diagram of the system.

The constant temperature and constant heat flux thermal boundary conditions, both developing distinct flow patterns, represent limiting cases of ideally conducting and insulating plates in Rayleigh–Bénard convection flows, respectively. This study bridges the gap in between, using a conjugate heat transfer (CHT) set-up and studying finite thermal diffusivity ratios $\kappa _s \! / \! \kappa _f$ to better represent real-life conditions in experiments. A three-dimensional Rayleigh–Bénard convection configuration including two fluid-confining plates is studied via direct numerical simulations given a Prandtl number ${Pr}=1$. The fluid layer of height $H$ and horizontal extension $L$ obeys no-slip boundary conditions at the two solid–fluid interfaces and an aspect ratio of ${\Gamma }=L/H=30$ while the relative thickness of each plate is ${\Gamma _s}=H_s/H=15$. The entire domain is laterally periodic. Here, different $\kappa _s \! / \! \kappa _f$ are investigated for moderate Rayleigh numbers $Ra=\left \{ 10^4, 10^5 \right \}$. We observe a gradual shift of the size of the characteristic flow patterns and their induced heat and mass transfer as $\kappa _s \! / \! \kappa _f$ is varied, suggesting a relation between the recently studied turbulent superstructures and supergranules for constant temperature and constant heat flux boundary conditions, respectively. Performing a linear stability analysis for this CHT configuration confirms these observations theoretically while extending previous studies by investigating the impact of a varying solid plate thickness $\Gamma _s$. Moreover, we study the impact of $\kappa _s \! / \! \kappa _f$ on both the thermal and viscous boundary layers. Given the prevalence of finite $\kappa _s \! / \! \kappa _f$ in nature, this work is a starting point to extend our understanding of pattern formation in geo- and astrophysical convection flows.

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.

A linear stability model based on a phase-field method is established to study the formation of ripples on the ice surface. The pattern on horizontal ice surfaces, e.g. glaciers and frozen lakes, is found to be originating from a gravity-driven instability by studying ice–water–air flows with a range of water and ice thicknesses. Contrary to gravity, surface tension and viscosity act to suppress the instability. The results demonstrate that a larger value of either water thickness or ice thickness corresponds to a longer dominant wavelength of the pattern, and a favourable wavelength of 90 mm is predicted, in agreement with observations from nature. Furthermore, the profiles of the most unstable perturbations are found to be with two peaks at the ice–water and water–air interfaces whose ratio decreases exponentially with the water thickness and wavenumber.

This paper is focused on the existence and uniqueness of nonconstant steady states in a reaction–diffusion–ODE system, which models the predator–prey interaction with Holling-II functional response. Firstly, we aim to study the occurrence of regular stationary solutions through the application of bifurcation theory. Subsequently, by a generalized mountain pass lemma, we successfully demonstrate the existence of steady states with jump discontinuity. Furthermore, the structure of stationary solutions within a one-dimensional domain is investigated and a variety of steady-state solutions are built, which may exhibit monotonicity or symmetry. In the end, we create heterogeneous equilibrium states close to a constant equilibrium state using bifurcation theory and examine their stability.

Parametric oscillations of an interface separating two fluid phases create nonlinear surface waves, called Faraday waves, which organise into simple patterns, such as squares and hexagons, as well as complex structures, such as double hexagonal and superlattice patterns. In this work, we study the influence of surfactant-induced Marangoni stresses on the formation and transition of Faraday-wave patterns. We use a control parameter, $B$, that assesses the relative importance of Marangoni stresses as compared with the surface-wave dynamics. Our results show that the threshold acceleration required to destabilise a surfactant-covered interface through vibration increases with increasing $B$. For a surfactant-free interface, a square-wave pattern is observed. As $B$ is incremented, we report transitions from squares to asymmetric squares, weakly wavy stripes and ultimately to ridges and hills. These hills are a consequence of the bidirectional Marangoni stresses at the neck of the ridges. The mechanisms underlying the pattern transitions and the formation of exotic ridges and hills are discussed.

We report pattern formation in an otherwise non-uniform and unsteady flow arising in high-speed liquid entrainment conditions on the outer wall of a wide rotating drum. We show that the coating flow in this rotary dragout undergoes axial modulations to form an array of roughly vertical thin liquid sheets which slowly drift from the middle of the drum towards its sidewalls. Thus, the number of sheets fluctuates in time such that the most probable rib spacing varies ever so slightly with the speed, and a little less weakly with the viscosity. We propose that these axial patterns are generated due to a primary instability driven by an adverse pressure gradient in the meniscus region of the rotary drag-out flow, similar to the directional Saffman–Taylor instability, as is wellknown for ribbing in film-splitting flows. Rib spacing based on this mechanistic model turns out to be proportional to the capillary length, wherein the scaling factor can be determined based on existing models for film entrainment at both low and large capillary numbers. In addition, we performed direct numerical simulations, which reproduce the experimental phenomenology and the associated wavelength. We further include two numerical cases wherein either the liquid density or the liquid surface tension is quadrupled while keeping all other parameters identical with experiments. The rib spacings of these cases are in agreement with the predictions of our model.

In this paper, we report the spatiotemporal dynamics of an intraguild predation (IGP)-type predator–prey model incorporating harvesting and prey-taxis. We first discuss the local and global existence of the classical solutions in N-dimensional space. It is found that the model has a global classical solution when controlling the prey-taxis coefficient in a certain range. Thereafter, we focus on the existence of the steady-state bifurcation. Moreover, we theoretically investigate the properties of the bifurcating solution near the steady-state bifurcation critical threshold. As a consequence, the spatial pattern formation of this model can be theoretically confirmed. Importantly, by means of rigorous theoretical derivation, we provide discriminant criteria on the stability of the bifurcating solution. Finally, the complicated patterns are numerically displayed. It is demonstrated that the harvesting and prey-taxis significantly affect the pattern formation of this IGP-type predator–prey model. Our main results of this paper reveal that (i) The repulsive prey-taxis could destabilize the spatial homogeneity, while the attractive prey-taxis effect and self-diffusion will stabilize the spatial homogeneity of this model. (ii) Numerical results suggest that over-harvesting for prey or predators is not advisable, it can lead to an ecological imbalance due to a significant reduction in population numbers. However, harvesting within a certain range is a feasible approach.

‘Print stone’ is an iron-banded siltstone from the Pilbara Province of Western Australia that bears partial resemblance to iconic East Kimberley ‘zebra rock’ in both pattern morphology and mineralogical composition. Using a combination of mineralogy and elemental geochemistry, this investigation examines the mechanisms underlying the formation of periodic iron-oxide banding in print stone. We demonstrate that print stone patterns probably arose from the periodic deposition of hydrothermal pyrite during the early Palaeoproterozoic, as evidenced by the distinctive cuboid morphology of the hematite pigment, the deposition of iron oxides along fluid-transport pathways, the presence of extensive hydrothermal pyrite elsewhere in the formation, and the presence of a positive europium anomaly. Through spatial analysis of the iron-oxide banding, we further show that print stone adheres to the Liesegang spacing law with a spacing coefficient of 0.018. This suggests that the periodic deposition of pyrite in print stone arose due to the Liesegang phenomenon, which was probably triggered by the infiltration of near-neutral, sulfidic hydrothermal fluids into a ferruginous, feldspathic shale. Altogether, the present findings demonstrate the opportunity for iron-oxide Liesegang bands to develop in hydrothermal systems, providing additional insight into the mechanisms underlying the formation of the East Kimberley zebra rock and other banded geological material.

In ‘The chemical basis of morphogenesis’ (1952), Alan Turing introduced an idea that revolutionised our thinking about pattern formation. He proposed that diffusion could lead to the spontaneous formation of regular patterns. Here, we discuss the impact of Turing’s idea on plant science using three well-established examples at different scales: ROP patterning inside single cells, epidermal patterning across several cells and whole vegetation patterns. Also at intermediate levels, e.g., organ spacing, plants look surprisingly regular. But not all regular patterns are Turing patterns, careful observation and prediction of the patterning process—not just the final pattern—is critical to distinguish between mechanisms.

Liquid flowing down a fibre readily destabilises into a train of beads, commonly called a bead-on-fibre pattern. Bead formation results from capillary-driven instability and gives rise to patterns with constant velocity and time-invariant bead frequency $f$ whenever the instability is absolute. In this study, we develop a scaling law for $f$ that relates the Strouhal number $St$ and capillary number $Ca$ for Ostwaldian power-law liquids with Newtonian liquids recovered as a limiting case. We validate our proposed scaling law by comparing it with prior experimental data and new experimental data using xanthan gum solutions to produce a low capillary number $Ca \leq 10^{-2}$ regime. The experimental data encompasses both Ostwaldian and Newtonian flow, as well as symmetric and asymmetric patterns, and we find the data collapses along the predicted trend across seven orders of magnitude in $Ca$. Our proposed scaling law is a powerful tool for studying and applying bead-on-fibre flows where $f$ is critical, such as heat and mass transfer systems.