用于多级压气机非定常流动的谐波平衡方法研究

推进技术 ›› 2017, Vol. 38 ›› Issue (11): 2512-2521.

用于多级压气机非定常流动的谐波平衡方法研究

王子维¹,范召林²,江雄¹,陈逖¹

中国空气动力研究与发展中心计算空气动力研究所，四川绵阳 621000,中国空气动力研究与发展中心科技部，四川绵阳 621000,中国空气动力研究与发展中心计算空气动力研究所，四川绵阳 621000,中国空气动力研究与发展中心计算空气动力研究所，四川绵阳 621000

发布日期:2021-08-15
作者简介:王子维，男，硕士生，研究领域为叶轮机械CFD。

Investigation on Harmonic Balance Method for Simulation of Unsteady Flow in the Multi-Stage Compressor

Computational Aerodynamics Institute，China Aerodynamics Research and Development Center， Mianyang 621000，China,Science and Technology Department，China Aerodynamics Research and Development Center， Mianyang 621000，China,Computational Aerodynamics Institute，China Aerodynamics Research and Development Center， Mianyang 621000，China and Computational Aerodynamics Institute，China Aerodynamics Research and Development Center， Mianyang 621000，China

Published:2021-08-15

摘要/Abstract

摘要： 为了在保证计算精度的同时，进一步提高压气机流动非定常模拟的计算效率，在课题组自行开发的结构化有限体积解算器上实现了谐波平衡法。以1-1/2级轴流压气机为研究对象，采用谐波平衡法和相滞后法对其非定常流场进行了数值模拟，并采用谐波平衡法对该压气机的时序效应进行了初步研究。计算表明，随着谐波阶数的增加，谐波平衡法的计算结果逐渐接近于相滞后法的计算结果；在谐波阶数达到5阶以后，谐波平衡法的计算结果不再随谐波阶数的增加而变化，且与相滞后法的计算结果匹配；谐波平衡法在保证与相滞后法相同的计算精度的情况下，计算效率相比相滞后法高了近3倍；5阶谐波平衡法即可很好地模拟入口导叶-静叶时序效应。算例表明，谐波平衡法相比相滞后法具有更多的选择空间和更高的计算效率。

关键词: 轴流压气机；非定常流动；谐波平衡法；相滞后法

Abstract: In order to improve the simulational efficiency of the unsteady flow in the axial compressor without lossing the computational accuracy，the harmonic balance method were implemented into an in-house finite volume code. Taking a 1-1/2 stage axial compressor as an example，results of the harmonic balance method and the phase lag method were compared. The clocking effect of the compressor was also analysed using the harmonic balance method. The simulational results indicate that，the simulational result of the harmonic balance method grows closer to the phase lag method as the harmonic order of the harmonic balance method increases. When the harmonic order is set more or equal to 5，the computational result does not vary with the harmonic order and matches well with the computational result of the phase lag method. The harmonic balance method is as precise as the phase lag method with a factor of 3 speed-up. The harmonic balance method with the harmonic order of 5 can simulate the IGV-S1 clocking effect well. The results demonstrate that，compared with the phase lag method，the harmonic balance method has more computational options and higher computational efficiency.

Key words: Axial compressor；Unsteady flow；Harmonic balance method；Phase lag method

王子维,范召林,江雄,陈逖. 用于多级压气机非定常流动的谐波平衡方法研究[J]. 推进技术, 2017, 38(11): 2512-2521.

WANG Zi-wei¹,FAN Zhao-lin²,JIANG Xiong¹,CHEN Ti¹. Investigation on Harmonic Balance Method for Simulation of Unsteady Flow in the Multi-Stage Compressor[J]. Journal of Propulsion Technology, 2017, 38(11): 2512-2521.

参考文献

[1] Rai M M. Three-Dimensional Navier-Stokes Simulations of Turbine Rotor-Stator Interaction[J]. AIAA Journal, 1989, 5(3): 305-311.
[2] Arnone A, Pacciani R. Rotor-Stator Interaction Analysis Using the Navier-Stokes Equations and a Multigrid Method[J]. Journal of Turbomachinery, 1996, 118(4): 679-689.
[3] Clark J P, Stetson G M, Magge S S, et al. The Effect of Airfoil Scaling on the Predicted Unsteady Loading on the Blade of a 1 and 1/2 Stage Transonic Turbine and a Comparison with Experimental Results[R]. ASME 2000-GT-0446.
[4] Mayorca M A, De Andrade J A, Vogt D M, et al. Effect of Scaling of Blade Row Sectors on the Prediction of Aerodynamic Forcing in a Highly Loaded Transonic Compressor Stage[J]. Journal of Turbomachinery, 2011, 133(2).
[5] Trébinjac I, Benichou E, Buffaz N. Full-Annulus Simulation of the Surge Inception in a Transonic Centrifugal Compressor[J]. Journal of Thermal Science, 2015, 24(5): 442-451.
[6] 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, 4(5): 879-886.
[7] Ekici K, Hall K C, Kielb R E. Harmonic Balance Analysis of Forced Response in a Transonic Compressor[R]. AIAA 2009-1240.
[8] Sicot F, Dufour G, Gourdain N. A Time-Domain Harmonic Balance Method for Rotor/Stator Interactions[J]. Journal of Turbomachinery, 2012, 134(1).
[9] Ekici K, Hall K C. Nonlinear Analysis of Unsteady Flows in Multistage Turbomachines Using Harmonic Balance[J]. AIAA Journal, 2007, 45(5): 1047-1057.
[10] Guédeney T, Gomar A, Gallard F, et al. Non-Uniform Time Sampling for Multiple-Frequency Harmonic Balance Computations[J]. Journal of Computational Physics, 2013, 236: 317-345.
[11] 杜鹏程, 宁方飞. 跨声风扇周向畸变流动的谐波平衡法计算[J]. 推进技术, 2012, 33(3): 391-397. (DU Peng-cheng, NING Fang-fei. Numerical Simulation of Transonic Fan under Circumferential Inlet Distortion Using Harmonic Balance Method[J]. Journal of Propulsion Technology, 2012, 33(3): 391-397.)
[12] 施永强, 杨青真, 黄秀全, 等. 一种快速模拟振荡叶栅非定常流的数值方法[J]. 空气动力学学报, 2014, 32(4): 481-487.
[13] 张翔, 陈玉春, 黄秀全. 高低空环境下压气机叶片颤振特性研究[J]. 推进技术, 2015, 36(4): 572-578. (ZHANG Xiang, CHEN Yu-chun, HUANG Xiu-quan.Research on Flutter Characteristics of Compressor Blades in High and Low Altitude[J]. Journal of Propulsion Technology, 2015, 36(4): 572-578.)
[14] 马灿, 苏欣荣, 袁新. 用于非定常流动的谐波平衡方法研究[J]. 工程热物理学报, 2015, 36(4): 739-743.
[15] Gerolymos G A, Michon G J, Neubauer J. Analysis and Application of Chorochronic Periodicity in Turbomachinery Rotor/Stator Interaction Computations[J]. Journal of Propulsion and Power, 2002, 18(6): 1139-1152.
[16] Erdos J I, Alzner E, McNally W. Numerical Solution of Periodic Transonic Flow through a Fan Stage[J]. AIAA Journal, 1977, 15(11): 1559-1568.(编辑:史亚红)　　　　　　　　　　　　　*　收稿日期:2016-07-05；修订日期:2016-09-01。作者简介:王子维,男,硕士生,研究领域为叶轮机械CFD。E-mail: ziweigegeya@163.com