航空学报 > 2014, Vol. 35 Issue (5): 1181-1192   doi: 10.7527/S1000-6893.2013.0430

基于高阶耗散紧致格式的GMRES方法收敛特性研究

燕振国1, 刘化勇1, 毛枚良1,2, 邓小刚3, 朱华君1   

  1. 1. 中国空气动力研究与发展中心 空气动力学国家重点实验室, 四川 绵阳 621000;
    2. 中国空气动力研究与发展中心 计算空气动力研究所, 四川 绵阳 621000;
    3. 国防科学技术大学, 湖南 长沙 410073
  • 收稿日期:2013-06-25 修回日期:2013-10-12 出版日期:2014-05-25 发布日期:2013-11-16
  • 通讯作者: 毛枚良,Tel.:0816-2463102 E-mail:mml219@163.com E-mail:mml219@163.com
  • 作者简介:燕振国男,硕士研究生。主要研究方向:计算流体力学,高阶精度格式,隐式时间推进方法。Tel:0816-7067915E-mail:zhgyan@outlook.com;刘化勇男,博士。主要研究方向:计算流体力学,高阶精度格式,高超声速空气动力学。Tel:0816-7067915E-mail:hyliu@skla.com.cn;毛枚良男,博士,研究员。主要研究方向:计算流体力学,高阶精度格式,高超声速空气动力学。Tel:0816-2463102E-mail:mml219@163.com
  • 基金资助:

    国家自然科学基金(11072259);国家“973”计划(2009CB723801)

Convergence Property Investigation of GMRES Method Based on High-order Dissipative Compact Scheme

YAN Zhenguo1, LIU Huayong1, MAO Meiliang1,2, DENG Xiaogang3, ZHU Huajun1   

  1. 1. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, China;
    2. Institute of Computational Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, China;
    3. National University of Defense Technology, Changsha 410073, China
  • Received:2013-06-25 Revised:2013-10-12 Online:2014-05-25 Published:2013-11-16
  • Supported by:

    National Natural Science Foundation of China (11072259); National Basic Research Program of China (2009CB723801)

摘要:

计算效率较低是当前限制高阶精度计算方法应用的重要因素。为了提高高阶精度混合型耗散紧致格式(HDCS)的计算效率,发展了适合多块对接网格的广义最小残值(GMRES)方法,并利用GMRES方法开展了HDCS格式的加速收敛研究。首先研究了GMRES的预处理方法、CFL数和内层迭代步数对HDCS数值模拟收敛特性的影响,计算结果显示:点松弛方法是一种高效的预处理方法;CFL数对计算收敛速度影响较大;GMRES方法存在最优的内层迭代步数。利用GMRES方法完成了NACA 0012翼型绕流、NLR 7301翼型绕流和DLR-F4翼身组合体绕流的数值模拟,并与其他隐式时间推进方法进行了对比,GMRES方法计算更加稳定,并且计算效率相对LU-SGS(Lower-Upper Symmetric Gauss-Seidel)方法可以提高5倍以上。研究结果表明,本文发展的GMRES方法在多块对接网格中具有良好的计算稳定性,计算结果的残差可以收敛到更低的量级,并且可以较大幅度地提高高阶精度数值模拟的计算效率。

关键词: 隐式时间推进方法, 数值方法, 高阶精度格式, HDCS, GMRES, 收敛特性, 计算效率

Abstract:

Low computational efficiency is an important factor constraining the application of high-order numerical methods. To improve the computational efficiency of hybrid cell-edge and cell-node dissipative compact scheme (HDCS), a generalized minimum residual (GMRES) algorithm suitable for multi-block structured grids is developed to accelerate simulations. The influence of GMRES's precondition methods, CFL number and sub-iteration number on convergence property of HDCS high-order simulations is investigated. It is shown that the point relaxation method is an efficient precondition method, that the CFL number can greatly affect the computational efficiency, and that GMRES has an optimal sub-iteration number. GMRES is applied to simulations of NACA 0012 airfoil, NLR 7301 airfoil and DLR-F4 wing/body configuration, and is compared with other implicit time integration methods. By using GMRES, the computation becomes more stable, and the computational efficiency can be improved by more than 5 times when compared with the LU-SGS(Lower-Upper Symmetric Gauss-Seidel) method. The results indicate that the GMRES method developed in this paper has good stability in multi-block structured grids, the residual can converge to lower levels, and GMRES can greatly improve the computational efficiency of high-order simulations.

Key words: implicit time integration method, numerical methods, high-order scheme, HDCS, GMRES, convergence property, computational efficiency

中图分类号: