Skip to main content Accessibility help

Frequency–wavenumber spectral analysis of spatio-temporal flows

  • Christopher J. Geoga (a1), Charlotte L. Haley (a1), Andrew R. Siegel (a1) (a2) and Mihai Anitescu (a1) (a3)


We propose a fully spatio-temporal approach for identifying spatially varying modes of oscillation in fluid dynamics simulation output by means of multitaper frequency–wavenumber spectral analysis. One-dimensional spectrum estimation has proven to be a valuable tool in the analysis of turbulence data applied spatially to determine the rate of energy transport between spatial scales, or temporally to determine frequencies of oscillatory flows. It also allows for the quantitative comparison of flow characteristics between two scenarios using a standard basis. It has the limitation, however, that it neglects coupling between spatial and temporal structures. Two-dimensional frequency–wavenumber spectral analysis allows one to decompose waveforms into standing or travelling variety. The extended higher-dimensional multitaper method proposed here is shown to have improved statistical properties over conventional non-parametric spectral estimators, and is accompanied by confidence intervals which estimate their uncertainty. Multitaper frequency–wavenumber analysis is applied to a canonical benchmark problem, namely, a direct numerical simulation of von Kármán vortex shedding off a square wall-mounted cylinder with two inflow scenarios with matching momentum-thickness Reynolds numbers $Re_{\unicode[STIX]{x1D703}}\approx 1000$ at the obstacle. Frequency–wavenumber analysis of a two-dimensional section of these data reveals that although both the laminar and turbulent inflow scenarios show a turbulent $-5/3$ cascade in wavenumber ( $\unicode[STIX]{x1D708}$ ) and frequency ( $f$ ), the flow characteristics differ in that there is a significantly more prominent discrete harmonic oscillation near $(f,\unicode[STIX]{x1D708})=(0.2,0.21)$ in wavenumber and frequency in the laminar inflow scenario than the turbulent scenario. This frequency–wavenumber pair corresponds to a travelling wave with velocity near one near the centre path of the vortex street.


Corresponding author

Email address for correspondence:


Hide All
Berkooz, G., Holmes, P. & Lumley, J. L. 1993 The proper orthogonal decomposition in the analysis of turbulent flows. Annu. Rev. Fluid Mech. 25, 539575.
Cressie, N. 1981 Transformations and the jackknife. J. Roy. Statist. Soc. B 43, 177182.
Cressie, N. & Wikle, C. K. 2011 Statistics for Spatio-Temporal Data. Wiley.
Efron, B. & Stein, C. 1981 The jackknife estimate of variance. Ann. Stat. 9, 586596.
Grünbaum, F. A. 1981 Toeplitz matrices commuting with tridiagonal matrices. Linear Algebr. Applics. 40, 2536.
Hanssen, A. 1997 Multidimensional multitaper spectral estimation. Signal Process. 58, 327332.
Hayashi, Y. 1979 Space–time spectral analysis of rotary vector series. J. Atmos. Sci. 36, 757766.
Holmes, P., Lumley, J. L. & Berkooz, G. 1998 Turbulence, Coherent Structures, Dynamical Systems and Symmetry. Cambridge University Press.
Kirby, J. F. 2014 Estimation of the effective elastic thickness of the lithosphere using inverse spectral methods: the state of the art. Tectonophysics 631, 87116.
Mallows, C. L. 1967 Linear processes are nearly Gaussian. J. Appl. Probab. 4, 313329.
Percival, D. B. & Walden, A. T. 1993 Spectral Analysis for Physical Applications: Multitaper and Conventional Univariate Techniques. Cambridge University Press.
Pratt, R. W. 1976 The interpretation of space–time spectral quantities. J. Atmos. Sci. 33 (6), 10601066.
Rowley, C. W. 2005 Model reduction for fluids, using balanced proper orthogonal decomposition. Intl J. Bifurcation Chaos 15, 9971013.
Rowley, C. W., Mezić, I., Bagheri, S., Schlatter, P. & Henningson, D. S. 2009 Spectral analysis of nonlinear flows. J. Fluid Mech. 641, 115127.
Schmidt, O. T., Towne, A., Colonius, T., Cavalieri, A. V. G., Jordan, P. & Brès, G. A. 2017 Wavepackets and trapped acoustic modes in a turbulent jet: coherent structure eduction and global stability. J. Fluid Mech. 825, 11531181.
Simons, F. J. & Wang, D. V. 2011 Spatiospectral concentration in the Cartesian plane. Intl J. Geomath. 2 (1), 136.
Slepian, D. 1964 Prolate spheroidal wave functions, Fourier analysis and uncertainty IV. Bell System Tech. J. 43, 30093057.
Slepian, D. 1978 Prolate spheroidal wave functions, Fourier analysis, and uncertainty V: the discrete case. Bell System Tech. J. 57 (5), 13711429.
von Storch, H. & Zwiers, F. W. 1999 Statistical Analysis in Climate Research. Cambridge University Press.
Thomson, D. J. 1982 Spectrum estimation and harmonic analysis. Proc. IEEE 70 (9), 10551096.
Thomson, D. J. 1990 Quadratic-inverse spectrum estimates: applications to paleoclimatology. Phil. Trans. R. Soc. Lond. A 332, 539597.
Thomson, D. J. 1994 An overview of multiple-window and quadratic-inverse spectrum estimation methods. Proc. ICASSP VI, 185194; invited plenary lecture.
Thomson, D. J. & Chave, A. D. 1991 Jackknifed error estimates for spectra, coherences, and transfer functions. In Advances in Spectrum Analysis and Array Processing (ed. Haykin, S.), vol. 1, chap. 2, pp. 58113. Prentice-Hall.
Vinuesa, R., Schlatter, P., Malm, J., Mavriplis, C. & Henningson, D. S. 2015 Direct numerical simulation of the flow around a wall-mounted square cylinder under various inflow conditions. J. Turbul. 16, 555587.
Wilks, D. S. 2011 Statistical Methods in the Atmospheric Sciences. Academic Press.
MathJax is a JavaScript display engine for mathematics. For more information see

JFM classification

Frequency–wavenumber spectral analysis of spatio-temporal flows

  • Christopher J. Geoga (a1), Charlotte L. Haley (a1), Andrew R. Siegel (a1) (a2) and Mihai Anitescu (a1) (a3)


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed