If the square matrix A is real, its eigenvalues are either real or come in complex conjugate pairs. c In such case, A is called diagonalizable or non-defective. If A is defective, we have A = V JV −1 with J being the canonical Jordan form [31,33].

Advances in experimental techniques and the ever-increasing fidelity of numerical simulations have led to an abundance of data describing fluid flows. This review discusses a range of techniques for analyzing such data, with the aim of extracting simplified models that capture the essential features of these flows, in order to gain insight into the flow physics, and potentially identify mechanisms for controlling these flows. We review well-developed techniques, such as proper orthogonal decomposition and Galerkin projection, and discuss more recent techniques developed for linear systems, such as balanced truncation and dynamic mode decomposition (DMD). We then discuss some of the methods available for nonlinear systems, with particular attention to the Koopman operator, an infinite-dimensional linear operator that completely characterizes the dynamics of a nonlinear system and provides an extension of DMD to nonlinear systems.

Dynamic mode decomposition (DMD) provides a practical means of extracting insightful dynamical information from fluids datasets. Like any data processing technique, DMD's usefulness is limited by its ability to extract real and accurate dynamical features from noise-corrupted data. Here we show analytically that DMD is biased to sensor noise, and quantify how this bias depends on the size and noise level of the data. We present three modifications to DMD that can be used to remove this bias: (i) a direct correction of the identified bias using known noise properties, (ii) combining the results of performing DMD forwards and backwards in time, and (iii) a total least-squares-inspired algorithm. We discuss the relative merits of each algorithm, and demonstrate the performance of these modifications on a range of synthetic, numerical, and experimental datasets. We further compare our modified DMD algorithms with other variants proposed in recent literature.

No abstract

No abstract

Unsteady aerodynamic models are necessary to accurately simulate forces and develop feedback controllers for wings in agile motion; however, these models are often high dimensional or incompatible with modern control techniques. Recently, reduced-order unsteady aerodynamic models have been developed for a pitching and plunging airfoil by linearizing the discretized Navier-Stokes equation with lift-force output. In this work, we extend these reduced-order models to include multiple inputs (pitch, plunge, and surge) and explicit parameterization by the pitch-axis location, inspired by Theodorsen's model. Next, we investigate the naïve application of system identification techniques to input-output data and the resulting pitfalls, such as unstable or inaccurate models. Finally, robust feedback controllers are constructed based on these low-dimensional state-space models for simulations of a rigid flat plate at Reynolds number 100. Various controllers are implemented for models linearized at base angles of attack α0 = 0 • , α0 = 10 • , and α0 = 20 • . The resulting control laws are able to track an aggressive reference lift trajectory while attenuating sensor noise and compensating for strong nonlinearities.

We seek to quantify non-normality of the most amplified resolvent modes and predict their features based on the characteristics of the base or mean velocity profile. A 2-by-2 model linear Navier-Stokes (LNS) operator illustrates how non-normality from mean shear distributes perturbation energy in different velocity components of the forcing and response modes. The inverse of their inner product, which is unity for a purely normal mechanism, is proposed as a measure to quantify non-normality. When the operator is normal but not self-adjoint, a phase difference between the forcing and response modes is indicated. In flows where there is downstream spatial dependence of the base/mean, mean flow advection separates the spatial support of forcing and response modes which impacts the inner product and leads to non-normal amplification. Success of mean flow (linear) stability analysis for a particular frequency depends on the normality of amplification. If the amplification is normal, the resolvent operator can be rewritten in its dyadic representation to reveal that the adjoint and forward stability modes are proportional to the forcing and response resolvent modes at that frequency. If the amplification is non-normal, then resolvent analysis is required to understand the origin of observed flow structures. Eigenspectra and pseudospectra are used to characterize these phenomena. Two test cases are studied: low Reynolds number cylinder flow and turbulent channel flow. The first deals mainly with normal mechanisms hence the success of both classical and mean stability analysis with respect to predicting the critical Reynolds number and global frequency of the saturated flow. Quantification of non-normality using the inverse inner product of the leading forcing and response modes agrees well with the product of the resolvent norm and distance between the imaginary axis and least stable eigenvalue. An approximate wavemaker computed from the resolvent modes scales with the length of the mean recirculation bubble. In the case of turbulent channel flow, structures result from both normal and non-normal mechanisms, suggesting that both are necessary elements to sustain turbulence. Mean shear is exploited most efficiently by stationary disturbances while bounds on the pseudospectra illustrate how non-normality is responsible for the most amplified disturbances at spatial wavenumbers and temporal frequencies corresponding to well-known turbulent structures. Some implications for flow control are discussed. * ssymon@caltech.edu arXiv:1712.05473v1 [physics.flu-dyn]

High-fidelity blood flow modelling is crucial for enhancing our understanding of cardiovascular disease. Despite significant advances in computational and experimental characterization of blood flow, the knowledge that we can acquire from such investigations remains limited by the presence of uncertainty in parameters, low resolution, and measurement noise. Additionally, extracting useful information from these datasets is challenging. Data-driven modelling techniques have the potential to overcome these challenges and transform cardiovascular flow modelling. Here, we review several data-driven modelling techniques, highlight the common ideas and principles that emerge across numerous such techniques, and provide illustrative examples of how they could be used in the context of cardiovascular fluid mechanics. In particular, we discuss principal component analysis (PCA), robust PCA, compressed sensing, the Kalman filter for data assimilation, low-rank data recovery, and several additional methods for reduced-order modelling of cardiovascular flows, including the dynamic mode decomposition and the sparse identification of nonlinear dynamics. All techniques are presented in the context of cardiovascular flows with simple examples. These data-driven modelling techniques have the potential to transform computational and experimental cardiovascular research, and we discuss challenges and opportunities in applying these techniques in the field, looking ultimately towards data-driven patient-specific blood flow modelling.