可压缩算法在燃油冷却数值研究中的应用

推进技术 ›› 2018, Vol. 39 ›› Issue (9): 2028-2034.

可压缩算法在燃油冷却数值研究中的应用

初敏,陈思员,史可天,艾邦成

中国航天空气动力技术研究院，北京 100074,中国航天空气动力技术研究院，北京 100074,中国航天空气动力技术研究院，北京 100074,中国航天空气动力技术研究院，北京 100074

发布日期:2021-08-15
作者简介:初敏，男，博士，研究领域为高超声速飞行器热防护设计。E-mail: chumin_caaa@126.com 通讯作者:艾邦成，男，博士，研究员，研究领域为高超声速飞行器热防护设计。

Application of Compressible Algorithm in Numerical Simulation of Fuel-Active-Cooling

China Academy of Aerospace Aerodynamics，Beijing 100074，China,China Academy of Aerospace Aerodynamics，Beijing 100074，China,China Academy of Aerospace Aerodynamics，Beijing 100074，China and China Academy of Aerospace Aerodynamics，Beijing 100074，China

Published:2021-08-15

摘要/Abstract

摘要： 为深入研究超燃冲压发动机的燃油主动冷却过程，在可压缩算法框架下，引入预处理技术、低耗散数值离散格式、一般流体状态方程、高压物性模型、燃油裂解模型等，发展了一套适用于求解一般流体低速可压缩问题的数值模拟方法。利用电加热管试验对数值方法进行了验证，并定量研究了PR-3燃油的冷却需求量及平板面板冷却槽道布置的影响。结果表明，数值结果与试验符合很好，证明当前数值方法在燃油冷却数值模拟中是准确有效的。针对350mm×90mm×6mm的平板面板冷却过程，需要约10g/s的燃油可确保1WM/m2热流下的防热安全。槽道布置对温度场影响明显，两侧进油方式比单侧进油的最大温差低约150K，均温性更好，更有利于降低结构热应力。

关键词: 燃油主动冷却；数值模拟；可压缩算法；预处理

Abstract: In order to study the fuel-active-cooling of scramjet， a numerical simulation method suitable to solve the low-speed compressible problems of general fluid was developed. Based on the framework of compressible algorithm， the method introduced precondition technology， low dissipation numerical scheme， general fluid state equation， physical property model of high pressure and fuel pyrolysis model. The results of electrically heated tube test were adopted to verify the physical models and algorithms. The mass flow rate of RP-3 fuel for cooling process of flat panel was analyzed quantitatively and the effects of the cooling channel arrangement on cooling process were studied. The results show that the physical models and algorithms are effective on fuel cooling. The mass flow rate of RP-3 should be more than 10g/s to ensure the thermal protection safety under 1MW/m2 heat flux when the size of flat panel is 350mm×90mm×6mm. The channel arrangement has obvious influence on the temperature field. The maximum temperature difference when there are two oil inlets is 150K lower than that when there is only one oil inlet. Better isothermal performance means it is beneficial to the decrease of structure thermal stress.

Key words: Fuel-active-cooling；Numerical simulation；Compressible algorithm；Precondition

初敏,陈思员,史可天,艾邦成. 可压缩算法在燃油冷却数值研究中的应用[J]. 推进技术, 2018, 39(9): 2028-2034.

CHU Min,CHEN Si-yuan,SHI Ke-tian,AI Bang-cheng. Application of Compressible Algorithm in Numerical Simulation of Fuel-Active-Cooling[J]. Journal of Propulsion Technology, 2018, 39(9): 2028-2034.

参考文献

[1] 范学军, 俞刚. 大庆RP-3航空煤油热物性分析[J]. 推进技术, 2006, 27(2): 187-192. (FAN Xue-jun, YU Gang. Analysis of Thermophysical Properties of Daqing RP-3 Aviation Kerosene[J]. Journal of Propulsion Technology, 2006, 27(2): 187-192.)
[2] 王夕. 超临界压力吸热型碳氢燃料热裂解及传热特性研究[D]. 北京:清华大学, 2013.
[3] 赵国柱, 宋文艳, 张若凌, 等. 超临界压力下正十烷流动传热的数值模拟[J]. 推进技术, 2014, 35(4): 537-543. (ZHAO Guo-zhu, SONG Wen-yan, ZHANG Ruo-ling, et al. Numerical Simulation on Flow and Heat Transfer of n-Decane under Supercritical Pressure[J]. Journal of Propulsion Technology, 2014, 35(4): 537-543.)
[4] 赵国柱, 宋文艳, 张若凌. 超临界压力下RP-3航空煤油吸热裂解反应的数值研究[J]. 航空学报, 2014, 35(6): 1513-1521.
[5] Jiang R P, Liu G Z, Wang X Q, et al. Thermal Cracking of Hydrocarbon Aviation Fuels in Regenerative Cooling Microchannels[J]. Energy and Fuels, 2013, 27(5): 2563-2577.
[6] Zhang R L, Zhao G Z, Le J L, et al. Numerical Study on Heat Transfer of Hydrocarbon Fuel with Thermal Cracking[R]. AIAA 2015-3621.
[7] 张强强. 微细通道内碳氢燃料传热与裂解过程的实验研究与CFD模拟[D]. 天津:天津大学, 2014.
[8] Bao W, Zhang S L, Qin J, et al. Numerical Analysis of Flowing Cracked Hydrocarbon Fuel inside Cooling Channels in View of Thermal Management[J]. Energy, 2014, 67(4): 149-161.
[9] Zhang S L, Qin J, Xie K L, et al. Thermal Behavior inside Scramjet Cooling Channels at Different Channel Aspect Ratios[J]. Journal of Propulsion and Power, 2016, 32(1): 57-70.
[10] Menter F R. Two Equation Eddy Viscosity Turbulence Models for Engineering Applications[J]. AIAA Journal, 2012, 32(8): 1598-1605.
[11] Yoon S, Jameson A. Lower-Upper Symmetric Gauss-Seidel Method for the Euler and Navier-Stokes Equations[J]. AIAA Journal, 1988, 26(9): 1025-1026.
[12] Edwards J R, Franklin R K, Liou M S. Low-Diffusion Flux-Splitting Methods for Real Fluid Flows at All Speeds[R]. AIAA 99-3327.
[13] Kim S K, Choi H S, Kim Y. Thermodynamic Modeling Based on a Generalized Cubic Equation of State for Kerosene/LOX Rocket Combustion[J]. Combustion and Flame, 2012, 159(3): 1351-1365.
[14] Cismondi M, Mollerup J. Development and Application of a Three-Parameter RK-PR Equation of State[J]. Fluid Phase Equilibria, 2005, 232(1): 74-89.
[15] Meng H, Yang V. A Unified Treatment of General Fluid Thermodynamics and its Application to a Preconditioning Scheme[J]. Journal of Computational Physics, 2003, 189(1): 277-304.
[16] Chorin A J. A Numerical Method for Solving Incompressible Viscous Flow Problem[J]. Journal of Computational Physics, 1967, 2(1): 12-26.
[17] Shuen J S, Chen K H, Choi Y. A Coupled Implicit Method for Chemical Non-Equilibrium Flows at All Speeds[J]. Journal of Computational Physics, 1993, 106(2): 306-318.
[18] 侯凌云, 董宁, 孙大鹏. 矩形冷却槽道内煤油热裂解过程数值研究[J]. 推进技术, 2014, 35(1): 128-132. (HOU Ling-yun, DONG Ning, SUN Da-peng. Numerical Study on Thermal Cracking Process of Kerosene in a Rectangular Cooling Channel[J]. Journal of Propulsion Technology, 2014, 35(1): 128-132.)
[19] Ward T A, Ervin J S, Striebich R C, et al. Simulations of Flowing Mildly-Cracked Normal Alkanes Incorporating Proportional Product Distributions[J]. Journal of Propulsion and Power, 2004, 20(3): 394-402.
[20] 童景山. 流体热物性学基本理论与计算[M]. 北京:中国石化出版社, 2008.　　　　　　　　　　　　　*　收稿日期:2017-09-08；修订日期:2017-11-02。作者简介:初敏,男,博士,研究领域为高超声速飞行器热防护设计。E-mail: chumin_caaa@126.com通讯作者:艾邦成,男,博士,研究员,研究领域为高超声速飞行器热防护设计。E-mail: aimen011@126.com(编辑:张荣莉)