EI、Scopus 收录
中文核心期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Gappy POD方法重构湍流数据的研究

李天一 BuzzicottiMichele BiferaleLuca 万敏平 陈十一

李天一, Buzzicotti Michele, Biferale Luca, 万敏平, 陈十一. Gappy POD方法重构湍流数据的研究. 力学学报, 2021, 53(10): 2703-2711 doi: 10.6052/0459-1879-21-464
引用本文: 李天一, Buzzicotti Michele, Biferale Luca, 万敏平, 陈十一. Gappy POD方法重构湍流数据的研究. 力学学报, 2021, 53(10): 2703-2711 doi: 10.6052/0459-1879-21-464
Li Tianyi, Buzzicotti Michele, Biferale Luca, Wan Minping, Chen Shiyi. Reconstruction of turbulent data with gappy POD method. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(10): 2703-2711 doi: 10.6052/0459-1879-21-464
Citation: Li Tianyi, Buzzicotti Michele, Biferale Luca, Wan Minping, Chen Shiyi. Reconstruction of turbulent data with gappy POD method. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(10): 2703-2711 doi: 10.6052/0459-1879-21-464

Gappy POD方法重构湍流数据的研究

doi: 10.6052/0459-1879-21-464
基金项目: 国家数值风洞工程(NNW2019ZT1-A01)和深圳市科技计划(KQTD20180411143441009)资助项目
详细信息
    作者简介:

    万敏平, 教授, 主要研究方向: 湍流数值模拟和理论, 磁流体和太阳风. E-mail: wanmp@sustech.edu.cn

  • 中图分类号: O357.5

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基函数在已知点上的值组成的矩阵的非列满秩性, 它主要取决于破损区域的面积大小和几何形状. 破损区域的面积越大, 或者破损面积相同时, 破损区域内信息所包含的相关性越大, 第二部分的重构误差越大.

     

  • 图  1  POD模态的特征值曲线

    Figure  1.  Eigenvalues of the POD modes

    2  (a)-(e) 破损面积相同但破损区域几何不同的部分缺失流场, (f) 真实的完整流场

    2.  (a)-(e) Damaged flow fields with gaps of the same area but different geometries. (f) The complete flow field

    图  2  (a)-(e) 破损面积相同但破损区域几何不同的部分缺失流场, (f) 真实的完整流场 (续)

    Figure  2.  (a)-(e) Damaged flow fields with gaps of the same area but different geometries. (f) The complete flow field (continued)

    图  3  破损区域中的均方重构误差$ {\rm{MS}}{{\rm{E}}_{{\rm{gap}}}} $关于${{{N_{{\rm{comp}}}}} \mathord{\left/ {\vphantom {{{N_{{\rm{comp}}}}} {{N_{{\rm{flow}}}}}}} \right. } {{N_{{\rm{flow}}}}}}$的变化曲线

    Figure  3.  Normalized mean square error in the gap, $ {{MS}}{{{E}}_{{\rm{gap}}}} $, as a function of ${{{N_{{\rm{comp}}}}} \mathord{\left/ {\vphantom {{{N_{{\rm{comp}}}}} {{N_{{\rm{flow}}}}}}} \right. } {{N_{{\rm{flow}}}}}}$

    图  4  $\boldsymbol{\tilde X}$的最小奇异值$ {\sigma _{\min }} $关于${{N_{{\rm{comp}}}}}$/$ {\vphantom {{{N_{{\rm{comp}}}}} {{N_{{\rm{flow}}}}}}} {{N_{{\rm{flow}}}}} $的变化曲线

    Figure  4.  The minimum singular value of $ \boldsymbol{\tilde X} $, $ {\sigma _{\min }} $, as a function of ${{N_{{\rm{comp}}}}}$/${ {\vphantom {{{N_{{\rm{comp}}}}} {{N_{{\rm{flow}}}}}}} } {{N_{{\rm{flow}}}}} $

    图  5  对于不同几何形状的破损, 破损区域中的均方重构误差$ {{MS}}{{{E}}_{{\rm{gap}}}} $关于破损尺寸的变化曲线

    Figure  5.  Normalized mean square error in the gap, $ {{MS}}{{{E}}_{{\rm{gap}}}} $, as a function of the gap size for different gap geometries

    图  6  对于不同几何形状的破损, ${N_{{\rm{flow}}}} - r\left( {\boldsymbol{\tilde X}} \right)$关于破损尺寸的变化曲线

    Figure  6.  ${N_{{\rm{flow}}}} - r\left( {\boldsymbol{\tilde X}} \right)$ as a function of the gap size for different gap geometries

    图  7  相同破损面积, 不同破损几何的gappy POD重构结果(对应不同行). 第1列: 破损的流场; 第2列: 重构的流场; 第3列: 原始流场; 第4列: 第一列中红色参考线上重构结果(虚线)与真实结果(实线)的分布; 第5列: 平均重构误差(红色), $ {\varDelta _u}\left( {{x_1}} \right) $, 及其在测试集上的平均(黑色), $ \left\langle {{\varDelta _u}\left( {{x_1}} \right)} \right\rangle $. ${x_1}$${x_2}$分别代表水平和竖直方向

    Figure  7.  Gappy POD reconstruction results for the same gap area and different gap geometries (one for each row). 1st column: damaged image in input. 2nd column: image generated in output. 3rd column: ground truth. 4th column: generated (dashed) and ground truth (solid) profiles along the vertical line shown in the 1st column. 5th column: mean reconstruction error, $ {\varDelta _u}\left( {{x_1}} \right) $, for each image (red line) and the average error, $ \left\langle {{\varDelta _u}\left( {{x_1}} \right)} \right\rangle $, over the test images (black curve). Note that ${x_1}$ and ${x_2}$ denotes the horizontal and vertical directions, respectively

  • [1] Holmes P, Lumley JL, Berkooz G, et al. Turbulence, Coherent Structures, Dynamical Systems and Symmetry. Cambridge: Cambridge University Press, 2012
    [2] Moin P, Moser RD. Characteristic-eddy decomposition of turbulence in a channel. Journal of Fluid Mechanics, 1989, 200: 471-509 doi: 10.1017/S0022112089000741
    [3] De Cillis G, Cherubini S, Semeraro O, et al. POD analysis of the recovery process in wind turbine wakes. Journal of Physics: Conference Series, 2020, 1618(6): 062016
    [4] De Cillis G, Cherubini S, Semeraro O, et al. POD-based analysis of a wind turbine wake under the influence of tower and nacelle. Wind Energy, 2021, 24(6): 609-633 doi: 10.1002/we.2592
    [5] Lacassagne T, Simoëns S, El Hajem M, et al. POD analysis of oscillating grid turbulence in water and shear thinning polymer solution. AIChE Journal, 2021, 67(1): e17044
    [6] 陈鑫, 刘莉, 岳振江. 基于本征正交分解和代理模型 的高超声速气动热模型降阶研究. 航空学报, 2015, 36(2): 462-472 (Chen Xin, Liu Li, Yue Zhenjiang. Reduced order aerothermodynamic modeling research for hypersonic vehicles based on proper orthogonal decomposition and surrogate method. Acta Aeronautica et Astronautica Sinica, 2015, 36(2): 462-472 (in Chinese)
    [7] 李卓越, 白宏磊, 高南. 基于 POD 的受激励后向台阶分离流低阶模型. 气体物理, 2017, 2(3): 33-43 (Li Zhuoyue, Bai Honglei, Gao Nan. Reduced order models of forced flow over a backward facing step based on POD method. Physics of Gases, 2017, 2(3): 33-43 (in Chinese)
    [8] Abbaszadeh M, Dehghan M. Reduced order modeling of time-dependent incompressible Navier–Stokes equation with variable density based on a local radial basis functions-finite difference (LRBF-FD) technique and the POD/DEIM method. Computer Methods in Applied Mechanics and Engineering, 2020, 364: 112914 doi: 10.1016/j.cma.2020.112914
    [9] Kean K, Schneier M. Error analysis of supremizer pressure recovery for POD based reduced-order models of the time-dependent navier-stokes equations. SIAM Journal on Numerical Analysis, 2020, 58(4): 2235-2264 doi: 10.1137/19M128702X
    [10] Girfoglio M, Quaini A, Rozza G. A POD-Galerkin reduced order model for a LES filtering approach. Journal of Computational Physics, 2021, 436: 110260 doi: 10.1016/j.jcp.2021.110260
    [11] Everson R, Sirovich L. Karhunen–Loeve procedure for gappy data. JOSA A, 1995, 12(8): 1657-1664 doi: 10.1364/JOSAA.12.001657
    [12] Willcox K. Unsteady flow sensing and estimation via the gappy proper orthogonal decomposition. Computers & Fluids, 2006, 35(2): 208-226
    [13] Yu B, Tong Y, Hu P, et al. A novel inversion approach for identifying the shape of cavity by combining gappy POD with direct inversion scheme. International Journal of Heat and Mass Transfer, 2020, 150: 119365 doi: 10.1016/j.ijheatmasstransfer.2020.119365
    [14] Dechanubeksa C, Chaturantabut S. An Application of a modified gappy proper orthogonal decomposition on complexity reduction of allen-cahn equation. Algorithms, 2020, 13(6): 148 doi: 10.3390/a13060148
    [15] 李静, 张伟伟. 基于 Gappy POD 的流场数据填补方法. 气体物理, 2020, 5(4): 1-10 (Li Jing, Zhang Weiwei. Gappy proper orthogonal decomposition for flow data reconstruction. Physics of Gases, 2020, 5(4): 1-10 (in Chinese)
    [16] Venturi D, Karniadakis GE. Gappy data and reconstruction procedures for flow past a cylinder. Journal of Fluid Mechanics, 2004, 519: 315-336 doi: 10.1017/S0022112004001338
    [17] Gunes H, Sirisup S, Karniadakis GE. Gappy data: to Krig or not to Krig? Journal of Computational Physics, 2006, 212(1): 358-382 doi: 10.1016/j.jcp.2005.06.023
    [18] Cambon C, Mansour NN, Godeferd FS. Energy transfer in rotating turbulence. Journal of Fluid Mechanics, 1997, 337: 303-332 doi: 10.1017/S002211209700493X
    [19] Sharma MK, Verma MK, Chakraborty S. On the energy spectrum of rapidly rotating forced turbulence. Physics of Fluids, 2018, 30(11): 115102 doi: 10.1063/1.5051444
    [20] Alexakis A, Biferale L. Cascades and transitions in turbulent flows. Physics Reports, 2018, 767: 1-101
    [21] Buzzicotti M, Aluie H, Biferale L, et al. Energy transfer in turbulence under rotation. Physical Review Fluids, 2018, 3(3): 034802 doi: 10.1103/PhysRevFluids.3.034802
    [22] Li T, Wan M, Wang J, et al. Flow structures and kinetic-potential exchange in forced rotating stratified turbulence. Physical Review Fluids, 2020, 5(1): 014802 doi: 10.1103/PhysRevFluids.5.014802
    [23] van Kan A, Alexakis A. Critical transition in fast-rotating turbulence within highly elongated domains. Journal of Fluid Mechanics, 2020, 899: A33
    [24] Di Leoni PC, Alexakis A, Biferale L, et al. Phase transitions and flux-loop metastable states in rotating turbulence. Physical Review Fluids, 2020, 5(10): 104603 doi: 10.1103/PhysRevFluids.5.104603
    [25] Pedlosky J. Geophysical Fluid Dynamics. New York: Springer, 1987
    [26] McWilliams JC. Fundamentals of Geophysical Fluid Dynamics. Cambridge: Cambridge University Press, 2006
    [27] Yang CF, Chi WC, van Haren H. Deep-sea turbulence evolution observed by multiple closely spaced instruments. Scientific Reports, 2021, 11(1): 1-12 doi: 10.1038/s41598-020-79139-8
    [28] 刘式适, 刘式达. 大气动力学 (第二版). 北京: 北京大学出版社, 2011

    (Liu Shishi, Liu Shida. Atmospheric Dynamics (second edition). Beijing: Peking University Press, 2011 (in Chinese))
    [29] 任福安. 海洋动力学. 大连: 大连海事大学出版社, 2001

    (Ren Fuan. Ocean Dynamics. Dalian: Dalian Maritime University Press, 2001 (in Chiniese))
    [30] Biferale L, Bonaccorso F, Buzzicotti M, et al. TURB-Rot. A large database of 3D and 2D snapshots from turbulent rotating flows. ArXiv Preprint ArXiv: 2006.07469, 2020
    [31] Canuto C, Hussaini MY, Quarteroni A, et al. Spectral Methods in Fluid Dynamics. Springer Science & Business Media, 2012
  • 加载中
图(8)
计量
  • 文章访问数:  981
  • HTML全文浏览量:  274
  • PDF下载量:  266
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-09-10
  • 录用日期:  2021-10-07
  • 网络出版日期:  2021-10-08
  • 刊出日期:  2021-10-26

目录

    /

    返回文章
    返回