A hierarchy of low-dimensional Galerkin models is proposed for the viscous, incompressible flow around a circular cylinder building on the pioneering works of Stuart (1958), Deane et al. (1991), and Ma & Karniadakis (2002). The empirical Galerkin model is based on an eight-dimensional Karhunen–Loève decomposition of a numerical simulation and incorporates a new ‘shift-mode’ representing the mean-field correction. The inclusion of the shift-mode significantly improves the resolution of the transient dynamics from the onset of vortex shedding to the periodic von Kármán vortex street. In addition, the Reynolds-number dependence of the flow can be described with good accuracy. The inclusion of stability eigenmodes further enhances the accuracy of fluctuation dynamics. Mathematical and physical system reduction approaches lead to invariant-manifold and to mean-field models, respectively. The corresponding two-dimensional dynamical systems are further reduced to the Landau amplitude equation.

Closed-loop turbulence control is a critical enabler of aerodynamic drag reduction, lift increase, mixing enhancement, and noise reduction. Current and future applications have epic proportion: cars, trucks, trains, airplanes, wind turbines, medical devices, combustion, chemical reactors, just to name a few. Methods to adaptively adjust open-loop parameters are continually improving toward shorter response times. However, control design for in-time response is challenged by strong nonlinearity, high-dimensionality, and time-delays. Recent advances in the field of model identification and system reduction, coupled with advances in control theory (robust, adaptive, and nonlinear) are driving significant progress in adaptive and in-time closed-loop control of fluid turbulence. In this review, we provide an overview of critical theoretical developments, highlighted by compelling experimental success stories. We also point to challenging open problems and propose potentially disruptive technologies of machine learning and compressive sensing.

We propose a novel cluster-based reduced-order modelling (CROM) strategy for unsteady flows. CROM combines the cluster analysis pioneered in Gunzburger’s group (Burkardt, Gunzburger & Lee,Comput. Meth. Appl. Mech. Engng, vol. 196, 2006a, pp. 337–355) and transition matrix models introduced in fluid dynamics in Eckhardt’s group (Schneider, Eckhardt & Vollmer,Phys. Rev. E, vol. 75, 2007, art. 066313). CROM constitutes a potential alternative to POD models and generalises the Ulam–Galerkin method classically used in dynamical systems to determine a finite-rank approximation of the Perron–Frobenius operator. The proposed strategy processes a time-resolved sequence of flow snapshots in two steps. First, the snapshot data are clustered into a small number of representative states, called centroids, in the state space. These centroids partition the state space in complementary non-overlapping regions (centroidal Voronoi cells). Departing from the standard algorithm, the probabilities of the clusters are determined, and the states are sorted by analysis of the transition matrix. Second, the transitions between the states are dynamically modelled using a Markov process. Physical mechanisms are then distilled by a refined analysis of the Markov process, e.g. using finite-time Lyapunov exponent (FTLE) and entropic methods. This CROM framework is applied to the Lorenz attractor (as illustrative example), to velocity fields of the spatially evolving incompressible mixing layer and the three-dimensional turbulent wake of a bluff body. For these examples, CROM is shown to identify non-trivial quasi-attractors and transition processes in an unsupervised manner. CROM has numerous potential applications for the systematic identification of physical mechanisms of complex dynamics, for comparison of flow evolution models, for the identification of precursors to desirable and undesirable events, and for flow control applications exploiting nonlinear actuation dynamics.

Low-dimensional empirical Galerkin models are developed for spatially evolving laminar and transitional shear layers, based on a Karhunen–Loève decomposition of incompressible two- and three-dimensional Navier–Stokes simulations. It is shown that the key to an accurate Galerkin model is a novel analytical pressure-term representation. The effect of the pressure term is elucidated by a modal energy-flow analysis in a mixing layer, which generalizes the framework developed by Rempfer (1991). In convectively unstable shear layers, it is shown in particular that neglecting small energy terms leads to large amplitude errors in the Galerkin model. The effect of the pressure term and small neglected energy flows is very important for a two-dimensional mixing layer, is less pronounced for the three-dimensional analogue, and can be considered as small in an absolutely unstable wake flow.

The spatio-temporal evolution of a turbulent swirling jet undergoing vortex breakdown has been investigated. Experiments suggest the existence of a self-excited global mode having a single dominant frequency. This oscillatory mode is shown to be absolutely unstable and leads to a rotating counter-winding helical structure that is located at the periphery of the recirculation zone. The resulting time-periodic 3D velocity field is predicted theoretically as being the most unstable mode determined by parabolized stability analysis employing the mean flow data from experiments. The 3D oscillatory flow is constructed from uncorrelated 2D snapshots of particle image velocimetry data, using proper orthogonal decomposition, a phase-averaging technique and an azimuthal symmetry associated with helical structures. Stability-derived modes and empirically derived modes correspond remarkably well, yielding prototypical coherent structures that dominate the investigated flow region. The proposed method of constructing 3D time-periodic velocity fields from uncorrelated 2D data is applicable to a large class of turbulent shear flows.

We propose a general dynamic reduced-order modelling framework for typical experimental data: time-resolved sensor data and optional non-time-resolved particle image velocimetry (PIV) snapshots. This framework can be decomposed into four building blocks. First, the sensor signals are lifted to a dynamic feature space without false neighbours. Second, we identify a sparse human-interpretable nonlinear dynamical system for the feature state based on the sparse identification of nonlinear dynamics (SINDy). Third, if PIV snapshots are available, a local linear mapping from the feature state to the velocity field is performed to reconstruct the full state of the system. Fourth, a generalized feature-based modal decomposition identifies coherent structures that are most dynamically correlated with the linear and nonlinear interaction terms in the sparse model, adding interpretability. Steps 1 and 2 define a black-box model. Optional steps 3 and 4 lift the black-box dynamics to a grey-box model in terms of the identified coherent structures, if non-time-resolved full-state data are available. This grey-box modelling strategy is successfully applied to the transient and post-transient laminar cylinder wake, and compares favourably with a proper orthogonal decomposition model. We foresee numerous applications of this highly flexible modelling strategy, including estimation, prediction and control. Moreover, the feature space may be based on intrinsic coordinates, which are unaffected by a key challenge of modal expansion: the slow change of low-dimensional coherent structures with changing geometry and varying parameters.