Dynamic mode decomposition
Physical systems, such as fluid flow or mechanical vibrations, behave in characteristic patterns, known as modes. In a recirculating flow, for example, one may think of a hierarchy of vortices, a big main vortex driving smaller secondary ones and so on. Most of the motion of such a system can be faithfully described using only a few of those patterns. In a purely mathematical setting similar modes can be extracted form the governing equations using an eigenvalue decomposition. But in many cases the mathematical model is very complicated or not available at all. In an experiment, the mathematical description is not at hand and one has to rely on the measured data only. The dynamic mode decomposition (DMD) is a mathematical method to extract the relevant modes from experimental data, without any recurrence to the governing equations. It can thus be applied to any dynamic phenomenon where appropriate data is available. It is similar, but different from proper orthogonal decomposition which has similar features but lacks dynamical information about the data.
A time-evolving physical situation may be approximated by the action of a linear operator to the instantaneous state vector.
The dynamic mode decomposition strives to approximate the evolution operator from a known sequence of observations, . Thus, we ask the following matrix equation to hold:
Generally, the vectors , and subsequently , are very high-dimensional, and so a strict eigendecomposition of is computationally difficult. However, in DMD it is assumed that the set of does not span the entire vector space (a good assumption, especially if there is spatial structure in the signal). Thus, after a given time , where is much less than the dimensionality of the system, one can write as a linear combination of the previous vectors, i.e., . In matrix form, we then have:
where S is the companion matrix
The eigenvalues of then approximate some of the eigenvalues of . However, since the is small (with dimensions ) as compared to , the eigenvalues and eigenvectors of can be computed with ease.
||This section possibly contains original research. (March 2012)|
Trailing edge of a profile
The wake of an obstacle in the flow may develop a Kármán vortex street. The Fig.1 shows the shedding of a vortex behind the trailing edge of a profile. The DMD-Analysis was applied to 90 sequential Entropy fields (animated gif (1.9MB) ) and yield an approximated Eigenvalue-Spectrum as depicted below. The analysis was applied to the numerical results, without referring to the governing equations. The profile is seen in white. The white arcs are the processor boundaries, since the computation was performed on a parallel computer using different computational blocks.
Roughly a third of the spectrum was highly damped (large, negative ) and is not shown. The dominant shedding mode is shown in the following pictures. The image to the left is the real part, the image to the right, the imaginary part of the Eigenvector.
Again, the Entropy-Eigenvector is shown in this pictures. The acoustic contents of the same mode is seen in the bottom half of the next plot. The top half corresponds to the entropy mode as above.
Synthetic example of a traveling pattern
The DMD analysis assumes a pattern of the form where is any of the independent variables of the problem, but has to be selected in advance. Take for example the pattern
With the time as the preselected exponential factor.
A sample is given in the following figure with , and The left picture shows the pattern without, the right with noise added. The amplitude of the random noise is the same as that of the pattern.
A DMD analysis is performed with 21 synthetically generated fields using a time interval , limiting the analysis to .
The spectrum is symmetric and shows 3 almost undamped modes (small negative Real part), whereas the other modes are heavily damped. Their numerical values are respectively. The real one corresponds to the mean of the field, whereas corresponds to the imposed pattern with . Yielding a relative error of -1/1000. Increasing the noise to 10 times the signal value yields about the same error. The real and imaginary part of one of the latter two Eigenmodes is depicted in the following figure.
Several other decompositions of experimental data exist. If the governing equations are available, an eigenvalue decomposition might be feasible.
- Schmid, Peter J. "Dynamic mode decomposition of numerical and experimental data." Journal of Fluid Mechanics 656.1 (2010): 5-28.
- Hasselmann, K., 1988. POPs and PIPs. The reduction of complex dynamical systems using principal oscillation and interaction patterns. J. Geophys. Res., 93(D9): 10975–10988.