Spectral proper orthogonal decomposition of rapid snapshot pairs sampled at sub-Nyquist intervals

TL;DR

This paper introduces a method to compute spectral proper orthogonal decomposition (SPOD) modes from pairwise data sampled at sub-Nyquist intervals, enabling efficient analysis of turbulent flows with less data.

Contribution

The authors develop a novel approach combining dynamic mode decomposition (DMD) with SPOD to analyze pairwise data with large gaps, reducing data requirements significantly.

Findings

Accurately de-aliases the SPOD spectrum in turbulent jet data.

Estimates mode shapes at frequencies up to St = 1.0 with over 90% less data.

Validates the method on numerical and experimental flow data.

Abstract

Modal decomposition methods are important for characterizing the low-dimensional dynamics of complex systems, including turbulent flows. Different methods have varying data requirements and produce modes with different properties. Spectral proper orthogonal decomposition (SPOD) produces orthogonal, energy-ranked spatial modes at discrete temporal frequencies for statistically stationary flows. However, SPOD requires long stretches of sequential, uniformly sampled, time-resolved data. These data requirements limit SPOD's use in experimental settings where the maximum capture rate of a camera is often slower than the Nyquist sampling rate required to resolve the highest turbulent frequencies. However, if two PIV systems operate in tandem, pairs of data can be acquired that are arbitrarily close in time. The dynamic mode decomposition (DMD) uses this pairwise data to resolve frequencies up to the Nyquist frequency associated with the small time step within a pair. However, these modes do not form a basis and have no set ranking. The present work attempts to compute SPOD modes from pairwise data with a small time step but with large gaps between pairs. We use DMD on pairwise data to estimate segment-wise, uniformly sampled series that can then be used to estimate the SPOD modes, intending to resolve frequencies between the gap and pair Nyquist limits. The method is tested on numerically obtained data of the linearized complex Ginzburg-Landau equation, as well as a Mach 0.4 isothermal turbulent jet. For the jet, pairwise SPOD can accurately de-alias the SPOD spectrum and estimate mode shapes at frequencies up to St = 1.0, while using over 90% less data.