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

doi:10.7527/S1000-6893.2016.120894

流体力学与飞行力学

本期目录 | 过刊浏览 | 高级检索

前一篇 | 后一篇

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

贡伊明¹, 刘战合², 刘溢浪¹, 张伟伟¹

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 Yiming¹, LIU Zhanhe², LIU Yilang¹, ZHANG Weiwei¹

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)

摘要/Abstract

摘要：

研究了时间谱方法求解周期性非定常流场的计算效率，并对时间谱方法应用于周期性非定常流动的隐式求解方法进行探讨。当采样点数增加或减缩频率增大时，时间谱方法对应的雅可比矩阵对角占优性质迅速恶化，导致很多传统的迭代方法失效。为了解决上述问题，论文采用带预处理的广义极小残差（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

中图分类号:

V211.3

贡伊明, 刘战合, 刘溢浪, 张伟伟. 时间谱方法中的高效GMRES算法[J]. 航空学报, 2017, 38(7): 120894-120894.

GONG Yiming, LIU Zhanhe, LIU Yilang, ZHANG Weiwei. Efficient GMRES algorithm in time spectral method[J]. ACTA AERONAUTICAET ASTRONAUTICA SINICA, 2017, 38(7): 120894-120894.

参考文献

[1] HALL K C, CRAWLEY E F. Calculation of unsteady flows in turbomachinery using the linearized Euler equations[J]. AIAA Journal, 1989, 27(6): 777-787.
[2] NING W, HE L. Computation of unsteady flows around oscillation blades using linear and nonlinear harmonic Euler methods[J]. Journal of Turbomachinery, 1998, 120(3): 508-514.
[3] HALL K C, THOMAS J P, CLARK W S. Computation of unsteady nonlinear flows in cascades using a harmonic balance technique[J]. AIAA Journal, 2002, 40(5): 879-886.
[4] MCMULLEN M, JAMESON A, ALONSO J J. Acceleration of convergence to a periodic steady state in turbomachinery flows: AIAA-2001-0152[R]. Reston: AIAA, 2001.
[5] MCMULLEN M, JAMESON A, ALONSO J J. Application of a nonlinear frequency domain solver to the Euler and Navier-Stokes equations: AIAA-2002-0120[R]. Reston: AIAA, 2002.
[6] DAI H H, YUE X K, YUAN J P, et al. A time domain collocation method for studying the aeroelasticity of a two dimensional airfoil with a structural nonlinearity[J]. Journal of Computational Physics, 2014, 270: 214-237.
[7] GOPINATH A K, JAMESON A. Time spectral method for periodic unsteady computations over two- and three- dimensional bodies: AIAA-2005-1220[R]. Reston: AIAA, 2005.
[8] VAN DER WEIDE E, GOPINATH A K, JAMESON A. Turbomachinery applications with the time spectral method: AIAA-2005-4905[R]. Reston: AIAA, 2005.
[9] CHOI S, POTSDAM M, LEE K, et al. Helicopter rotor design using a time-spectral and adjoint-based method: AIAA-2008-5810[R]. Reston: AIAA, 2008.
[10] GOMAR A, BONVY Q, SICOT F, et al. Convergence of Fourier-based time methods for turbomachinery wake passing problems[J]. Journal of Computational Physics, 2014, 278(C): 229-256.
[11] GOPINATH A K, JAMESON A. Application of the time spectral method to periodic unsteady vortex shedding: AIAA-2006-0449[R]. Reston: AIAA, 2006.
[12] 杨小权, 程苏堃, 杨爱明, 等. 基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟[J]. 航空学报, 2013, 34(4): 787-797. YANG X Q, CHENG S K, YANG A M, et al. Time spectral method for numerical simulation of unsteady viscous flow over oscillating airfoil and wing[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(4): 787-797 (in Chinese).
[13] 谢立军, 杨云军, 刘周, 等. 基于时间谱方法的飞行器动导数高效计算技术[J]. 航空学报, 2015, 36(6): 2016-2026. XIE L J, YANG Y J, LIU Z, et al. A high efficient method for computing dynamic derivatives of aircraft based on time spectral method[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(6): 2016-2026 (in Chinese).
[14] 詹磊, 刘锋. 应用傅里叶时间谱方法求解二维跨音速流动问题解的精度研究[J]. 航空工程进展, 2015, 6(4): 395-404. ZHAN L, LIU F.Study on accuracy of the solutions using the Fourier time spectral method for two-dimensional transonic flows[J]. Advances in Aeronautical Science and Engineering, 2015, 6(4): 395-404 (in Chinese).
[15] JAMESON A, SHANKARAN S. An assessment of dual-time stepping, time spectral and artificial compressibility based numerical algorithms for unsteady flow with applications to flapping wings[C]//19th AIAA Computational Fluid Dynamics. Reston: AIAA, 2009.
[16] MAVRIPLIS D J, YANG Z. Time spectral method for periodic and quasi-periodic unsteady computations on unstructured meshes[J]. Mathematical Modelling of Natural Phenomena, 2011, 6(3): 213-236.
[17] SAAD Y, SCHULTZ M H. GMRES: A generalized minimal residual algorithm for solvingnonsymmetric linear systems[J]. SIAM Journal on Scientific and Statistical Computing, 1986, 7(3): 856-869.
[18] JAMESON A, YOON S. Lower-upper implicit schemes with multiple grids for the Euler equations[J]. AIAA Journal, 1987, 25(7): 929-935.
[19] 李春娜, 叶正寅, 王刚. 基于二维非结构网格的GMRES隐式算法[J]. 西北工业大学学报, 2007, 25(5): 630-635. LI C N, YE Z Y, WANG G.GMRES implicit algorithm based on 2D unstructured meshes for solving Euler equations[J]. Journal of Northwestern Polytechnical University, 2007, 25(5): 630-635 (in Chinese).
[20] 刘巍, 张理论, 王勇献, 等. 计算空气动力学并行编程基础[M]. 北京: 国防工业出版社, 2013. LIU W, ZHANG L L, WANG Y X,et al. Computational aerodynamics parallel programming basis[M]. Beijing: National Defence Industry Press, 2013 (in Chinese).
[21] LANDON R H. NACA0012 oscillatory and transient pitching: AGARD-R-702 [R]. Paris: AGARD, 1982.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

主管单位：中国科学技术协会主办单位：中国航空学会北京航空航天大学

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

Efficient GMRES algorithm in time spectral method

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	张宝振, 王汉平, 徐峰, 吴志青. VSV调节机构的仿真提速和精度补偿措施[J]. 航空学报, 2022, 43(9): 226034-226034.
[2]	杨鲤铭, 吴杰, 董昊, 杜银杰. 基于速度空间非结构网格和守恒修正的改进离散速度方法[J]. 航空学报, 2022, 43(12): 627033-627033.
[3]	陈龙, 夏健, 田书玲. 旋翼干扰流场数值模拟中的新型Blend低速预处理方法[J]. 航空学报, 2018, 39(8): 121903-121903.
[4]	杨体浩, 白俊强, 史亚云, 杨一雄. 适用于非周期流固耦合问题的时间谱方法[J]. 航空学报, 2018, 39(5): 121654-121654.
[5]	李欢, 陈江涛, 马明生, 周乃春. 一种基于特征关系式的预处理远场边界条件[J]. 航空学报, 2017, 38(12): 121364-121364.
[6]	谢立军, 杨云军, 刘周, 周伟江. 基于时间谱方法的飞行器动导数高效计算技术[J]. 航空学报, 2015, 36(6): 2016-2026.
[7]	燕振国, 刘化勇, 毛枚良, 马燕凯, 朱华君. 三阶HWCNS的构造及其在高超声速流动中的应用[J]. 航空学报, 2015, 36(5): 1460-1470.
[8]	黄飞, 沈清, 程晓丽, 赵波. 一种DSMC分子仿真下的权因子预定义方法[J]. 航空学报, 2014, 35(8): 2174-2181.
[9]	燕振国, 刘化勇, 毛枚良, 邓小刚, 朱华君. 基于高阶耗散紧致格式的GMRES方法收敛特性研究[J]. 航空学报, 2014, 35(5): 1181-1192.
[10]	杨小权, 程苏堃, 杨爱明, 孙刚. 基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟[J]. 航空学报, 2013, 34(4): 787-797.
[11]	杜磊, 宁方飞. 低速叶型气动反问题设计方法[J]. 航空学报, 2011, 32(7): 1180-1188.
[12]	苗赢;孙兆伟. 星载GPS测量数据预处理方法研究[J]. 航空学报, 2010, 31(3): 602-607.
[13]	刘晨;王江峰;伍贻兆. 低马赫数流动中的预处理Euler方程的收敛特性[J]. 航空学报, 2009, 30(5): 842-848.
[14]	翟坤;李顺利;杨涤. 用户星天线空间扫描策略[J]. 航空学报, 2008, 29(5): 1308-1313.
[15]	肖中云;江雄;陈作斌;刘刚. 低速预处理方法在悬停旋翼模拟中的应用[J]. 航空学报, 2008, 29(2): 321-326.