航空学报 > 2017, Vol. 38 Issue (7): 120894-120894   doi: 10.7527/S1000-6893.2016.120894

时间谱方法中的高效GMRES算法

贡伊明1, 刘战合2, 刘溢浪1, 张伟伟1   

  1. 1. 西北工业大学 航空学院, 西安 710072;
    2. 郑州航空工业管理学院 航空工程学院, 郑州 450046
  • 收稿日期:2016-10-27 修回日期:2016-12-04 出版日期:2017-07-15 发布日期:2016-12-26
  • 通讯作者: 张伟伟,E-mail:aeroelastic@nwpu.edu.cn E-mail:aeroelastic@nwpu.edu.cn
  • 基金资助:

    国家自然科学基金优秀青年基金(11622220);高等学校创新引智计划(B17037)

Efficient GMRES algorithm in time spectral method

GONG Yiming1, LIU Zhanhe2, LIU Yilang1, ZHANG Weiwei1   

  1. 1. School of Aeronautics, Northwestern Polytechnical University, Xi'an 710072, China;
    2. School of Aeronautic Engineering, Zhengzhou University of Aeronautics, Zhengzhou 450046, China
  • Received:2016-10-27 Revised:2016-12-04 Online:2017-07-15 Published:2016-12-26
  • Supported by:

    National Natural Science Foundation of China for Excellent Young Scholars (11622220);Programme of Introducing Talents of Discipline to Universities (B17037)

摘要:

研究了时间谱方法求解周期性非定常流场的计算效率,并对时间谱方法应用于周期性非定常流动的隐式求解方法进行探讨。当采样点数增加或减缩频率增大时,时间谱方法对应的雅可比矩阵对角占优性质迅速恶化,导致很多传统的迭代方法失效。为了解决上述问题,论文采用带预处理的广义极小残差(GMRES)算法来提高雅可比系数矩阵的计算收敛性。使用时间谱方法对NACA0012翼型强迫振荡算例进行计算,并与时域差分方法的计算效率和精度进行对比。研究表明在保证计算精度的同时,时间谱方法普遍可将计算效率提高一个量级左右。对于跨声速周期性流动,广义极小残差算法不论是稳定性还是收敛性都优于对称SGS迭代算法。

关键词: 时间谱方法, 广义极小值残差(GMRES)算法, 周期性非定常流动, 预处理, 计算效率

Abstract:

In this paper, the computational efficiency of the time-spectral method for solving the periodic unsteady flow field is studied, and the implicit method of time spectral method for solving the periodic unsteady flow is discussed. When the number of sampling points increases or the reduced frequency magnifies, the diagonal dominant property of the Jacobian matrix corresponding to the time spectral method deteriorates rapidly, resulting in the failure of many traditional iterative methods. In order to solve the problems above, the generalized minimum residual (GMRES) algorithm with preprocessing is used to improve the computational convergence of the Jacobian matrix. The time spectral method is used to compute the NACA0012 airfoil forced oscillation, and the computational efficiency and accuracy is compared with that of the time-domain difference method. The results show that the time spectral method can generally improve the computational efficiency an order of magnitude with saturated computational accuracy. For the transonic periodic flow, the GMRES algorithm is superior to SGS iterative algorithm both in stability and computational convergence.

Key words: time spectral method, generalized minimum residual (GMRES) algorithm, periodic unsteady flow, preprocess, computational efficiency

中图分类号: