[论文解读] Optimal frequency resolution for spectral proper orthogonal decomposition

Liam Heidt, Tim Colonius|arXiv (Cornell University)|Feb 24, 2024

Image and Signal Denoising Methods被引用 4

一句话总结

本文表明，频谱本征正交分解（SPOD）中的频率分辨率对模态精度具有决定性影响，先前使用的固定分辨率在物理上重要的频率处引入了显著偏差。作者提出一种基于物理的自适应频率分辨率算法，根据模态的空间变化动态调整分辨率，与标准的恒定分辨率方法相比，显著降低了SPOD模态的偏差，尤其在湍流喷流中表现更优。

ABSTRACT

We demonstrate that accurate computation of the spectral proper orthogonal decomposition (SPOD) critically depends on the choice of frequency resolution. Using both artificially generated data and large-eddy simulation data of a turbulent subsonic jet, we show that the optimal choice depends on how rapidly the SPOD modes change in space at adjacent frequencies. Previously employed values are found to be too high, resulting in unnecessarily biased results at physically important frequencies. A physics-informed adaptive frequency-resolution SPOD algorithm is developed that provides substantially less biased SPOD modes than the standard constant resolution method.

研究动机与目标

• 识别固定频率分辨率SPOD在准确捕捉统计平稳流中相干结构方面的局限性。
• 研究SPOD模态在相邻频率之间的空间变化如何影响最优频率分辨率。
• 开发一种基于物理的自适应频率分辨率算法，以减少SPOD模态估计中的偏差。
• 评估广泛使用的SPOD与响应函数模态之间对齐度量的统计收敛性和敏感性。
• 通过人工数据和湍流亚音速喷流的大涡模拟，证明自适应方法优于恒定分辨率SPOD。

提出的方法

• 作者使用具有已知SPOD谱和模态的人工数据，系统测试频率分辨率Δf对模态偏差的影响。
• 采用可变块大小和窗函数的Welch周期图法，估计每个频率箱的互功率谱密度（CSD）。
• 开发一种基于物理的自适应频率分辨率算法，根据相邻频率间SPOD模态变化率调整Δf，并利用收敛准则指导分辨率调整。
• 该方法采用改进的自适应过程，在低Strouhal数区域降低分辨率，高频频段则提高分辨率，依据模态相似性和方差减少原则。
• 使用快照法标准SPOD计算求解CSD矩阵的特征值问题，从特征向量中提取按能量排序的模态。
• 利用SPOD与响应函数模态之间的对齐度量来量化模态相似性，并通过不同数据长度评估其统计收敛性。

实验结果

研究问题

• RQ1 频率分辨率Δf的选择如何影响SPOD模态估计中的偏差，特别是在模态在相邻频率间快速变化时？
• RQ2 当底层模态在整个频带内表现出显著的空间变化时，SPOD的最优频率分辨率是什么？
• RQ3 与标准恒定分辨率方法相比，基于物理的自适应频率分辨率算法能否减少SPOD模态的偏差？
• RQ4 SPOD与响应函数模态之间的对齐度量对SPOD估计中的偏差和统计收敛性有多敏感？
• RQ5 增加数据长度在多大程度上能改善SPOD与响应函数模态之间对齐度量的收敛性和可靠性？

主要发现

• 先前使用的固定频率分辨率过高，在物理上重要的频率处引入了显著偏差，尤其是在模态形状快速变化的区域。
• 自适应频率分辨率算法显著降低了偏差，改进的自适应方法在所有频率和数据长度下均实现了SPOD与响应函数模态之间极佳的对齐。
• SPOD与响应函数模态之间的对齐度量对偏差和统计收敛性高度敏感，当数据长度从5,000增加到20,000个快照时，对齐度量最高可提升≈0.2。
• 最优Δf随数据长度增加而减小，与方差降低一致，实践者应在模态形状保持稳定且噪声减少时提高分辨率。
• 自适应方法在捕捉低频和中频Strouhal数区域的模态结构方面优于固定分辨率SPOD，这些区域的模态变化最为显著。
• 本研究表明，对齐度量需要大量数据才能收敛，仅依赖对齐度量得出的结论应谨慎对待，因其方差高且收敛缓慢。