基于Gama-Theta转捩模型的气冷涡轮气热耦合研究

推进技术 ›› 2014, Vol. 35 ›› Issue (7): 932-940.

基于Gama-Theta转捩模型的气冷涡轮气热耦合研究

李　涛,边　昕,冯　宇,黄洪雁

哈尔滨工业大学能源科学与工程学院，黑龙江哈尔滨 150001;哈尔滨工业大学能源科学与工程学院，黑龙江哈尔滨 150001;哈尔滨工业大学能源科学与工程学院，黑龙江哈尔滨 150001;哈尔滨工业大学能源科学与工程学院，黑龙江哈尔滨 150001

发布日期:2021-08-15
作者简介:李　涛(1984—)，男，博士生，研究领域为气冷涡轮气热耦合计算。

Conjugate Heat Transfer Simulation of Air-Cooled Turbine

School of Energy Science and Engineering，Harbin Institute of Technology，Harbin 150001，China;School of Energy Science and Engineering，Harbin Institute of Technology，Harbin 150001，China;School of Energy Science and Engineering，Harbin Institute of Technology，Harbin 150001，China;School of Energy Science and Engineering，Harbin Institute of Technology，Harbin 150001，China

Published:2021-08-15

摘要/Abstract

摘要： 基于非结构化混合网格，采用有限体积法求解了N-S方程，对流项格式为AUSM+，梯度的求解方法为加权最小二乘法，求解器采用LUSGS隐式求解方法，并进行了预处理和网格重排序。与导热方程耦合，搭建了气热耦合平台，采用直接耦合方法，交界面互相传递温度，并且采用能保证通量守恒的面积加权类的插值方式实现数据传递，将[γ-Reθt]两方程转捩模型应用在流场程序中。通过MARKII叶片4311和5411工况的实验结果进行了验证，结果表明采用转捩模型对压力分布的影响不大，与实验值吻合较好，通过间歇因子的分布可以看出，[γ-Reθt]模型成功地预测了从层流到湍流的转捩过程，在层流区采用[γ-Reθt]模型计算的涡粘系数与真实流动情况更加吻合，由于涡粘系数对传热的影响很大，从而有效地提高了边界层的层流区和转捩区的传热计算的精度，得到的温度和换热系数与实验值吻合更好，但是该模型是和SST模型耦合在一起的，由于SST模型的局限性，所以对激波和边界层干扰区域的模拟产生了7%左右的温度误差。

关键词: 非结构化网格；直接耦合；气冷涡轮；气热耦合；转捩

Abstract: The finite volume method was adopted to solve the N-S equation based on unstructured meshes. The scheme of convective flux was AUSM+，and the method for calculating gradient was weighted least squares method. The LUSGS implicit method was applied with precondition and reordering of grids. The conjugate heat transfer(CHT)simulation platform，which coupled with the heat conduction equation，was established. The direct-coupled method was adopted，and the temperature was transmitted between interfaces by the interpolation method，which can maintain the conservation of flux. The[γ-Reθt]transition model was used in the calculation of the flow field in order to study the impact of transition on the CHT. The numerical results were compared with experimental data of the 5411 and the 4311 experimental condition of MARKII blade. It indicates that the transition model has little effect on pressure distribution and agrees well with the experimental ones. It can be seen from the distribution of intermittency that [γ-Reθt]transition model successfully predicts transition process from laminar to turbulent flow. The eddy viscosity calculated in laminar region agrees better with the real flow. Because the eddy viscosity has great influence on heat transfer simulation，the[γ-Reθt] transition model effectively improves the accuracy of heat transfer simulation in laminar and transition area，and the temperature and heat transfer coefficient agree better with the experimental results. However，because the model is coupled with the SST model，due to limitations of the model SST，a greater simulation error was produced in the shock and boundary layer interaction region.

Key words: Unstructured mesh；Direct-coupled method；Air-cooled turbine；Conjugate heat transfer；Transition

李　涛,边　昕,冯　宇,黄洪雁. 基于Gama-Theta转捩模型的气冷涡轮气热耦合研究[J]. 推进技术, 2014, 35(7): 932-940.

LI Tao,BIAN Xin,FENG Yu,HUANG Hong-yan. Conjugate Heat Transfer Simulation of Air-Cooled Turbine[J]. Journal of Propulsion Technology, 2014, 35(7): 932-940.

参考文献

[1] Bohn D, Tummers C. Numerical 3-D Conjugate Flow and Heat Transfer Investigation of a Transonic Convection-cooled Thermal Barrier Coated Turbine Guide Vane with Reduced Cooling Fluid Mass Flow[R].ASME 2003-GT-38431.
[2] Kusterer K, Bohn D, Sugimoto T, et al. Conjugate Calculations for a Film-Cooled Blade Under Different Operating Conditions[R].ASME 2004-GT-53719.
[3] Montomoli F, Adami P, Della Gatta S, et al. Conjugate Heat Transfer Modelling in Film Cooled Blades[R].ASME 2004-GT-53177.
[4] Rigby D L, Lepicovsky J. Conjugate Heat Transfer Analysis of Internally Cooled Configurations[R].ASME 2001-GT-0405.
[5] Dunn M G, Convective Heat Transfer and Aerodynamics in Axial Flow Turbines [J].Turbomach, 2001, 123(4):637–686.
[6] Facchini B, Magi A, Greco A D. Conjugate Heat Transfer Simulation of a Radially Cooled Gas Turbine Vane[R]. ASME 2004-GT-54213.
[7] 郭兆元, 王强, 黄家骅,等. 涡轮叶片气热耦合数值模拟方法 [J].推进技术, 2009, 30(2):159-164. (GUO Zhao-yuan, WANG Qiang, HUANG Jia-hua,et al. Conjugate Heat Transfer Technique Applied to a Turbine Blade[J].Journal of Propulsion Technology, 2009, 30(2):159-164.)
[8] 周驰, 王强, 郭兆元, 等. 气冷涡轮气热耦合数值模拟 [J].推进技术, 2009, 30(5):566-570. (ZHOU Chi, WANG Qiang, GUO Zhao-yuan. Coupled Heat Transfer Analysis of Air Cooling Gas Turbine [J]. Journal of Propulsion Technology, 2009, 30(5):566-570.)
[9] 王强, 郭兆元, 周驰, 等. 考虑转捩的跨声速气冷涡轮叶片气热耦合计算 [J].航空动力学报, 2009, 24(13):2703-2710.
[10] 孙润鹏, 朱卫兵, 徐凌志, 等. 应用Transition k-kl-w转捩模型对内冷叶片气热耦合的数值模拟 [J].推进技术, 2012,33(2): 274-282.(SUN Run-peng, ZHU Wei-bing, XU Ling-zhi, Numerical Simulation of Coupled Heat Transfer for an Innercooled Blade by Applying Transition k-kl-ω Transition Model[J].Journal of Propulsion Technology, 2012, 33(2):274-282.)
[11] Menter F R, Langtry R B, Likki S R, et al. A Correlation-Based Transition Model Using Local Variables-Part Ⅰ : Model Formulation [R]. ASME 2004 -GT-53434.
[12] 曾军, 卿雄杰. 涡轮叶栅外换热系数计算 [J]. 航空动力学报, 2008, 23(7):1198-1204.
[13] Misaka T, Obayashi S. Application of Local Correlation-Based Transition Model to Flows Around Wings [R]. AIAA 2006-918.
[14] Blazek J. Computational Fluid Dynamics Principles and Applications [M]. UK: Kidington, Oxford OX5 IGB, 2001.
[15] 商立英. 非结构化网格上非稳态导热问题的求解 [D].南京:南京理工大学, 2006.
[16] 宁方飞, 徐力平. GMRES算法在二维定常无粘流计算中的应用 [J].计算物理, 2000, 17(5).
[17] 刘鑫, 陆林生. 拼接网格通量守恒插值算法研究 [J].计算机应用与软件, 2012, 29(2).
[18] 王强. 提高气冷涡轮气热耦合计算精度的计算方法研究[D].哈尔滨:哈尔滨工业大学, 2009.
[19] Hylton L D, Milhec M S, Turner E R, et al. Analytical and Experimental Evaluation of the Heat Transfer Distribution Over the Surface of Turbine Vanes[R].NASA-CR-168015, 1983.(编辑:张荣莉)收稿日期:2013-09-30；修订日期:2013-11-28。作者简介:李涛(1984—),男,博士生,研究领域为气冷涡轮气热耦合计算。E-mail:hitaolittle@126.com