Gappy POD方法重构湍流数据的研究
RECONSTRUCTION OF TURBULENT DATA WITH GAPPY POD METHOD
-
摘要:Gappy POD 是一种基于本征正交分解(proper orthogonal decomposition, POD)的数据重构方法. 本文研究了gappy POD在湍流数据重构中的应用, 主要关注了以下两个因素的影响: 第一, 数据本身的复杂程度, 即构成流场的POD模态数量; 第二, 破损区域的面积大小和几何形状. 考虑到上述因素, 本文重新严格地表述了gappy POD的重构过程, 并推导出gappy POD重构误差的公式. 论文选取旋转湍流数据为案例进行了gappy POD重构的研究, 并解释了构成gappy POD重构误差的两个部分. 第一部分来自流场POD展开的截断误差, 该截断误差会被POD基函数在已知点上的值组成的矩阵的最小特征值放大. 这部分误差主要取决于流场的复杂程度, 当流场复杂程度较低时, 相应误差随采用的POD模态数目增大而减小. 当流场复杂程度较高时, 很小的POD截断误差也会导致很大的重构误差, 此时需要采用流场所有的POD模态进行重构以消除截断误差. 重构误差的第二部分来自POD基函数在已知点上的值组成的矩阵的非列满秩性, 它主要取决于破损区域的面积大小和几何形状. 破损区域的面积越大, 或者破损面积相同时, 破损区域内信息所包含的相关性越大, 第二部分的重构误差越大.Abstract:Gappy POD is a method of data reconstruction based on the proper orthogonal decomposition (POD). We study the applicability of gappy POD to the reconstruction of fluid turbulence configurations and focus mainly on two factors. The first factor is the complexity of the data, which mostly depends on the number of POD modes with non-zero eigenvalues. The second factor is the area and the geometry of the gap. By taking these factors into account, we reformulate the gappy POD reconstruction and derive a formula to compute the reconstruction error. Rotating turbulence data is used as a case study of gappy POD reconstruction, where the reconstruction error can be separated into two parts. The first contribution to the reconstruction error is from the truncation error during the POD expansion and it is amplified by the smallest eigenvalue of the matrix, which consists of POD modes at known indexes. This error mainly depends on the flow complexity, e.g. for flow of moderate complexity, this error decreases with the increase in number of POD modes employed during the reconstruction process. For flow of large complexity, a small POD truncation error can be detrimental and contribute signification to the reconstruction error. Therefore, all POD modes should be considered when utilizing Gappy POD reconstruction to eliminate the truncation error, especially for the turbulent flow field. The second part of the reconstruction error appears when the matrix composed of POD modes at the known points is not of full column rank. This part of error depends on the area and the geometry of the gap. The gap area determines the amount of the lost information. For the same gap area, the gap geometry determines the correlation of the lost information. Gappy POD reconstruction works well when both the amount and the correlation of the lost information are small.