Using direct numerical simulations, we systematically investigate the inner-layer turbulence of a turbulent vertical buoyancy layer (a model for a vertical natural convection boundary layer) at a constant Prandtl number of $0.71$. Near-wall streaky structures of streamwise velocity fluctuations, synonymous with the buffer layer streaks of canonical wall turbulence, are not evident at low and moderate Reynolds numbers (${\textit{Re}}$) but manifest at high ${\textit{Re}}$. At low ${\textit{Re}}$, the turbulent production in the near-wall region is negligible; however, this increases with increasing ${\textit{Re}}$. By using domains truncated in the streamwise, spanwise and wall-normal directions, we demonstrate that the turbulence production in the near-wall region at moderate and high ${\textit{Re}}$ is largely independent of large-scale motions and outer-layer turbulence. On a fundamental level, the near-wall turbulence production is autonomous and self-sustaining, and a well-developed bulk is not needed to drive the inner-layer turbulence. Near-wall streaks are also not essential for this autonomous process. The type of thermal boundary condition only marginally influences the velocity fluctuations, revealing that the turbulence dynamics are primarily governed by the mean-shear induced by the buoyancy field and not by the thermal fluctuations, despite the current flow being solely driven by buoyancy. In the inner layer, the spanwise wavelength of the eddies responsible for positive shear production is remarkably similar to that of canonical wall turbulence at moderate and high ${\textit{Re}}$ (irrespective of near-wall streaks). Based on these findings, we propose a mechanistic model that unifies the near-wall shear production of vertical buoyancy layers and canonical wall turbulence.

We derive a mathematical model for the overflow fusion glass manufacturing process. In the limit of zero wedge angle, the model leads to a canonical fluid mechanics problem in which, under the effects of gravity and surface tension, a free-surface viscous flow transitions from lubrication flow to extensional flow. We explore the leading-order behaviour of this problem in the limit of small capillary number, and find that there are four distinct regions where different physical effects control the flow. We obtain leading-order governing equations, and determine the solution in each region using asymptotic matching. The downstream behaviour reveals appropriate far-field conditions to impose on the full problem, resulting in a simple governing equation for the film thickness that holds at leading order across the entire domain.

Induced diffusion (ID), an important mechanism of spectral energy transfer due to interacting internal gravity waves (IGWs), plays a significant role in driving turbulent dissipation in the ocean interior. In this study, we revisit the ID mechanism to elucidate its directionality and role in ocean mixing under varying IGW spectral forms, with particular attention to deviations from the standard Garrett–Munk spectrum. The original interpretation of ID as an action diffusion process, as proposed by McComas et al., suggests that ID is inherently bidirectional, with its direction governed by the vertical-wavenumber spectral slope $\sigma$ of the IGW action spectrum, $n \propto m^\sigma$. However, through the direct evaluation of the wave kinetic equation, we reveal a more complete depiction of ID, comprising both a diffusive and a scale-separated transfer rooted in the energy conservation within wave triads. Although the action diffusion may reverse direction depending on the sign of $\sigma$ (i.e. red or blue spectra), the net transfer by ID consistently leads to a forward energy cascade at the dissipation scale, contributing positively to turbulent dissipation. This supports the viewpoint of ID as a dissipative mechanism in physical oceanography. This study presents a physically grounded overview of ID, and offers insights into the specific types of wave–wave interactions responsible for turbulent dissipation.

We analyse the long-time dynamics of trajectories within the stability boundary between laminar and turbulent square duct flow. If not constrained to a symmetric subspace, the edge trajectories exhibit a chaotic dynamics characterised by a sequence of alternating quiescent phases and intense bursting episodes. The dynamics reflects the different stages of the well-known near-wall streak–vortex interaction. Most of the time, the edge states feature a single streak with a number of flanking vortices attached to one of the four surrounding walls. The initially straight streak undergoes a linear instability and eventually breaks in an intense bursting event. At the same time, the downstream vortices give rise to a new low-speed streak at one of the neighbouring walls, thereby causing the turbulent activity to ‘switch’ from one wall to the other. If the edge dynamics is restricted to a single or twofold mirror-symmetric subspace, the bursting and wall-switching episodes become self-recurrent in time, representing the first periodic orbits found in square duct flow. In contrast to the chaotic edge states in the non-symmetric case, the imposed symmetries enforce analogue bursting cycles to simultaneously appear at two parallel opposing walls in a mirror-symmetric configuration. Both the localisation of turbulent activity to one or two walls and the wall-switching dynamics are shown to be common phenomena in marginally turbulent duct flows. We argue that such episodes represent transient visits of marginally turbulent trajectories to some of the edge states detected here.

Symmetry-based analyses of multiscale velocity gradients highlight that strain self-amplification (SS) and vortex stretching (VS) drive forward energy transfer in turbulent flows. By contrast, a strain–vorticity covariance mechanism produces backscatter that contributes to the bottleneck effect in the subinertial range of the energy cascade. We extend these analyses by using a normality-based decomposition of filtered velocity gradients in forced isotropic turbulence to distinguish contributions from normal straining, pure shearing and rigid rotation at a given scale. Our analysis of direct numerical simulation (DNS) data illuminates the importance of shear layers in the inertial range and (especially) the subinertial range of the cascade. Shear layers contribute significantly to SS and VS and play a dominant role in the backscatter mechanism responsible for the bottleneck effect. Our concurrent analysis of large-eddy simulation (LES) data characterizes how different closure models affect the flow structure and energy transfer throughout the resolved scales. We thoroughly demonstrate that the multiscale flow features produced by a mixed model closely resemble those in a filtered DNS, whereas the features produced by an eddy viscosity model resemble those in an unfiltered DNS at a lower Reynolds number. This analysis helps explain how small-scale shear layers, whose imprint is mitigated upon filtering, amplify the artificial bottleneck effect produced by the eddy viscosity model in the inertial range of the cascade. Altogether, the present results provide a refined interpretation of the flow structures and mechanisms underlying the energy cascade and insight for designing and evaluating LES closure models.

Transonic buffet presents time-dependent aerodynamic characteristics associated with shock, turbulent boundary layer and their interactions. Despite strong nonlinearities and a large degree of freedom, there exists a dominant dynamic pattern of a buffet cycle, suggesting the low dimensionality of transonic buffet phenomena. This study seeks a low-dimensional representation of transonic airfoil buffet at a high Reynolds number with machine learning. Wall-modelled large-eddy simulations of flow over the OAT15A supercritical airfoil at two Mach numbers, $M_\infty = 0.715$ and 0.730, respectively producing non-buffet and buffet conditions, at a chord-based Reynolds number of ${Re} = 3\times 10^6$ are performed to generate the present datasets. We find that the low-dimensional nature of transonic airfoil buffet can be extracted as a sole three-dimensional latent representation through lift-augmented autoencoder compression. The current low-order representation not only describes the shock movement but also captures the moment when the separation occurs near the trailing edge in a low-order manner. We further show that it is possible to perform sensor-based reconstruction through the present low-dimensional expression while identifying the sensitivity with respect to aerodynamic responses. The present model trained at ${Re} = 3\times 10^6$ is lastly evaluated at the level of a real aircraft operation of ${Re} = 3\times 10^7$, exhibiting that the phase dynamics of lift is reasonably estimated from sparse sensors. The current study may provide a foundation towards data-driven real-time analysis of transonic buffet conditions under aircraft operation.

We explore the mechanisms and regimes of mixing in yield-stress fluids by simulating the stirring of an infinite, two-dimensional domain filled with a Bingham fluid. A cylindrical stirrer moves along a circular path at constant speed, with the path radius fixed at twice the stirrer diameter; the domain is initially quiescent and marked by a passive dye in the lower half. We first examine the mixing process in Newtonian fluids, identifying three key mechanisms: interface stretching and folding around the stirrer’s path, diffusion across streamlines and dye advection and interface stretching due to vortex shedding. Introducing yield stress leads to notable mixing localisation, manifesting through three mechanisms: advection of vortices within a finite distance of the stirrer, vortex entrapment near the stirrer and complete suppression of vortex shedding at high yield stresses. Based on these mechanisms, we classify three distinct mixing regimes: (i) regime SE, where shed vortices escape the central region, (ii) regime ST, where shed vortices remain trapped near the stirrer and (iii) regime NS, where no vortex shedding occurs. These regimes are quantitatively distinguished through spectral analysis of energy oscillations, revealing transitions and the critical Bingham and Reynolds numbers. The transitions are captured through effective Reynolds numbers, supporting the hypothesis that mixing regime transitions in yield-stress fluids share fundamental characteristics with bluff-body flow dynamics. The findings provide a mechanistic framework for understanding and predicting mixing behaviours in yield-stress fluids, suggesting that the localisation mechanisms and mixing regimes observed here are archetypal for stirred-tank applications.

We analyse the dynamics of a weakly elastic spherical particle translating parallel to a rigid wall in a quiescent Newtonian fluid in the Stokes limit. The particle motion is constrained parallel to the wall by applying a point force and a point torque at the centre of its undeformed shape. The particle is modelled using the Navier elasticity equations. The series solutions to the Navier and the Stokes equations are used to obtain the displacement and velocity fields in the solid and fluid, respectively. The point force and the point torque are calculated as series in small parameters $\alpha$ and $1/H$, using the domain perturbation method and the method of reflections. Here, $\alpha$ is the measure of elastic strain induced in the particle resulting from the fluid’s viscous stress and $H$ is the non-dimensional gap width, defined as the ratio of the distance of the particle centre from the wall to its radius. The results are presented up to $\textit {O}(1/H^3)$ and $\textit {O}(1/H^2)$, assuming $\alpha \sim 1/H$, for cases where gravity is aligned and non-aligned with the particle velocity, respectively. The deformed shape of the particle is determined by the force distribution acting on it. The hydrodynamic lift due to elastic effects (acting away from the wall) appears at $\textit {O}(\alpha /H^2)$ in the former case. In an unbounded domain, the elastic effects in the latter case generate a hydrodynamic torque at O($\alpha$) and a drag at O($\alpha ^2$). Conversely, in the former case, the torque is zero, while the drag still appears at O($\alpha ^2$).

Nearly fifty years ago, Roberts (1978) postulated that the Earth’s magnetic field, which is generated by turbulent motions of liquid metal in its outer core, likely results from a subcritical dynamo instability characterised by a dominant balance between Coriolis, pressure and Lorentz forces (requiring a finite-amplitude magnetic field). Here, we numerically explore subcritical convective dynamo action in a spherical shell, using techniques from optimal control and dynamical systems theory to uncover the nonlinear dynamics of magnetic field generation. Through nonlinear optimisation, via direct-adjoint looping, we identify the minimal seed – the smallest magnetic field that attracts to a nonlinear dynamo solution. Additionally, using the Newton-hookstep algorithm, we converge stable and unstable travelling wave solutions to the governing equations. By combining these two techniques, complex nonlinear pathways between attracting states are revealed, providing insight into a potential subcritical origin of the geodynamo. This paper showcases these methods on the widely studied benchmark of Christensen et al. (2001, Phys. Earth Planet. Inter., vol. 128, pp. 25–34), laying the foundations for future studies in more extreme and realistic parameter regimes. We show that the minimal seed reaches a nonlinear dynamo solution by first approaching an unstable travelling wave solution, which acts as an edge state separating a hydrodynamic solution from a magnetohydrodynamic one. Furthermore, by carefully examining the choice of cost functional, we establish a robust optimisation procedure that can systematically locate dynamo solutions on short time horizons with no prior knowledge of its structure.

We present a framework for parametric proper orthogonal decomposition (POD)-Galerkin reduced-order modelling (ROM) of fluid flows that accommodates variations in flow parameters and control inputs. As an initial step, to explore how the locally optimal POD modes vary with parameter changes, we demonstrate a sensitivity analysis of POD modes and their spanned subspace, respectively rooted in Stiefel and Grassmann manifolds. The sensitivity analysis, by defining distance between POD modes for different parameters, is applied to the flow around a rotating cylinder with varying Reynolds numbers and rotation rates. The sensitivity of the subspace spanned by POD modes to parameter changes is represented by a tangent vector on the Grassmann manifold. For the cylinder case, the inverse of the subspace sensitivity on the Grassmann manifold is proportional to the Roshko number, highlighting the connection between geometric properties and flow physics. Furthermore, the Reynolds number at which the subspace sensitivity approaches infinity corresponds to the lower bound at which the characteristic frequency of the Kármán vortex street exists (Noack & Eckelmann, J. Fluid Mech., 1994, vol. 270, pp. 297–330). From the Stiefel manifold perspective, sensitivity modes are derived to represent the flow field sensitivity, comprising the sensitivities of the POD modes and expansion coefficients. The temporal evolution of the flow field sensitivity is represented by superposing the sensitivity modes. Lastly, we devise a parametric POD-Galerkin ROM based on subspace interpolation on the Grassmann manifold. The reconstruction error of the ROM is intimately linked to the subspace-estimation error, which is in turn closely related to subspace sensitivity.