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.

Experimental studies of natural convection in yield stress fluids have revealed transient behaviours that contradict predictions from viscoplastic models. For example, at a sufficiently large yield stress, these models predict complete motionlessness; below a critical value, yielding and motion onset can be delayed in viscoplastic models. In both cases, however, experiments observe immediate motion onset. We present numerical simulations of the transient natural convection of elastoviscoplastic (EVP) fluids in a square cavity with differentially heated side walls, exploring the role of elasticity in reconciling theoretical predictions with experimental observations. We consider motion onset in EVP fluids under two initial temperature distributions: (i) a linear distribution characteristic of steady pure conduction, and (ii) a uniform distribution representative of experimental conditions. The Saramito EVP model exhibits an asymptotic behaviour similar to the Kelvin-Voigt model as $t\to 0^+$, where material behaviour is primarily governed by elasticity and solvent viscosity. The distinction between motion onset and yielding, a hallmark of EVP models, is the key feature that bridges theoretical predictions with experimental observations. While motion onset is consistently immediate (as seen in experiments), yielding occurs with a delay (as predicted by viscoplastic models). Scaling analysis suggests that this delay varies logarithmically with the yield stress and is inversely proportional to the elastic modulus. The intensity of the initial pre-yield motion increases with higher yield stress and lower elastic modulus. The observed dynamics resemble those of under- and partially over-damped systems, with a power-law fit providing an excellent match for the variation of oscillation frequency with the elastic modulus.