基于边界层理论的球头驻点热流计算方法

doi:10.7527/S1000-6893.2024.30448

流体力学与飞行力学

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

前一篇 |

基于边界层理论的球头驻点热流计算方法

田润雨¹^,², 龚红明¹^,²(), 常雨¹^,², 孔小平¹^,²

^1.飞行器流体物理全国重点实验室，绵阳 621000
^2.中国空气动力研究与发展中心超高速空气动力研究所，绵阳 621000

收稿日期:2024-03-26 修回日期:2024-06-04 接受日期:2024-06-21 出版日期:2024-06-26 发布日期:2024-06-25
通讯作者: 龚红明 E-mail:gh_ming@163.com

Calculation method for heat flow at stagnation point of spherical head based on boundary layer theory

Runyu TIAN¹^,², Hongming GONG¹^,²(), Yu CHANG¹^,², Xiaoping KONG¹^,²

^1.National Key Laboratory of Aerospace Physics in Fluids，Mianyang 621000，China
^2.Hypervelocity Aerodynamics Institute，China Aerodynamics Research and Development Center，Mianyang 621000，China

Received:2024-03-26 Revised:2024-06-04 Accepted:2024-06-21 Online:2024-06-26 Published:2024-06-25
Contact: Hongming GONG E-mail:gh_ming@163.com

摘要/Abstract

摘要：

准确获知球头驻点热流，对于飞行器防热设计至关重要。本文基于平衡空气球头驻点层流自相似边界层假设，开展了边界层方程相关推导，采用四阶龙格-库塔法数值求解坐标变换后的常微分边界层方程组，建立正向逼近打靶法，以寻找方程最优解，从而建立了通过求解边界层方程组获得球头驻点热流的方法，将该热流简称为BLES热流，在与试验值的对比中获得了较好的对比效果；以此方法为基础，针对高度10~60 km、不同速度、不同壁面温度的182个工况，开展了Fay-Riddell平衡流条件球头驻点热流公式计算偏差分析，发现在大部分工况下，壁面温度接近边界层外缘温度时，其热流计算值偏差明显增加；变形驻点热流表达式，将其中的无量纲参数拟合为多个无量纲参数组合的形式，得到了多参数拟合的新的球头驻点热流计算公式，并在多个工况下开展了热流计算效果对比研究，获得了优于Fay-Riddell平衡流条件球头驻点热流公式的应用效果。

关键词: 边界层理论, 层流自相似, 球头, 驻点热流, 平衡流, 气动热, BLES热流

Abstract:

In aircraft thermal protection design， it is very important to accurately know the heat flux at the stagnation point of spherical head. Based on the assumption of self-similar boundary layer at the stagnation point of equilibrium air， the boundary layer equations are derived. The fourth-order Runge-Kutta method is used to numerically solve the ordinary differential boundary layer equations after coordinate transformation， and the forward approximation shooting method is established to find the optimal solution for the equations. Thus， by solving the boundary layer equations， a method for obtaining the heat flux at the stagnation point of the spherical head is established， which is referred to as Boundary Layer Equations Stagnation （abbreviated as BLES） heat flux in this paper. The results obtained are consistent with the experimental values. Using this method， the calculation deviation of heat flux formula of spherical stagnation point under Fay-Riddell equilibrium flow condition at 182 working conditions of 10–60 km height and different velocity and wall temperature is analyzed. It is found that in most working conditions， the calculation deviation of heat flux increases obviously when the wall temperature approaches the outer edge temperature of boundary layer. New heat flux formulas for spherical stagnation point are obtained by fitting the dimensionless parameters into a combination of several dimensionless parameters. The heat flux calculation results under several working conditions are compared， obtaining better results than those of the spherical stagnation point heat flux formula under the Fay-Riddell equilibrium flow condition.

Key words: boundary layer theory, laminar flow self-similarity, spherical head, stagnation heat flux, equilibrium flow, aerodynamic heat, BLES heat flux

中图分类号:

V211.19

田润雨, 龚红明, 常雨, 孔小平. 基于边界层理论的球头驻点热流计算方法[J]. 航空学报, 2025, 46(2): 130448.

Runyu TIAN, Hongming GONG, Yu CHANG, Xiaoping KONG. Calculation method for heat flow at stagnation point of spherical head based on boundary layer theory[J]. Acta Aeronautica et Astronautica Sinica, 2025, 46(2): 130448.

图/表 13

图 1

图 2

表 1

表 2

图 3

图 4

图 5

图 6

表 3

表 4

图 7

表 5

图 8

参考文献 28

1	FAY J A， RIDDELL F R. Theory of stagnation point heat transfer in dissociated air［J］. Journal of the Aerospace Sciences， 1958， 25（2）： 73-85.
2	万云博，马戎，王年华，等. 基于混合网格多维梯度重构的热流预测方法研究［J］. 力学学报， 2018， 50（5）： 1003-1012.
	WAN Y B， MA R， WANG N H， et al. Accurate aero-heating predictions based on mul-ti-dimensional gradient reconstruction on hybrid unstructured grids［J］. Chinese Journal of Theoretical and Applied Mechanics， 2018， 50（5）： 1003-1012 （in Chinese）.
3	胡守超，庄宇，李贤，等. 高超声速气动热标模 HyHERM-I试验［J］. 航空学报， 2022， 43（S2）： 727804.
	HU S C， ZHUANG Y， LI X， et al. Hypersonic aero-heating environment research model HyHERM-I： Experiment［J］. Acta Aeronautica et Astronautica Sinica，2022， 43（S2）： 727804 （in Chinese）.
4	PRANDTL L. Applications of modern hydrodynamics to aeronautics： Report No.116［R］. Washington， D.C.： NACA， 1923.
5	孔祥志，韩宇峰. 高超声速平板边界层中波纹粗糙壁峰值热流的变化规律及预测［J］. 航空学报， 2023， 44（12）： 127769.
	KONG X Z， HAN Y F. Patterns and prediction of surface peak heat flux in hypersonic flat plate boundarylayer over wave wall［J］. Acta Aeronautica et Astronautica Sinica， 2023， 44（12）： 127769 （in Chinese）.
6	RESHOTKO E， COHEN C B.Similar solutions for the compressible laminar boundary layer with heat transfer and pressure gradient［J］.Technical Report Archive & Image Library， 1954， 22（3325）.
7	LEES L. Laminar heat transfer over blunt-nosed bodies at hypersonic flight speeds［J］. Journal of Jet Propulsion， 1956， 26（4）： 259-269.
8	TAUBER M E. A review of high-speed， convective， heat-transfer computation methods［J］. NASA STI/Recon Technical Report A， 1989， 90： 8_1-8_63.
9	DORODNITSYN A A. Laminar boundary layer in compressible fluid［J］. Doklady Akademii Nauk， 1942， 34： 213-219.
10	PARK S H， NEEB D， PLYUSHCHEV G， et al. A study on heat flux predictions for re-entry flight analysis［J］. Acta Astronautica， 2021， 187： 271-280.
11	田润雨，龚红明，常雨，等. 高焓流场球头外形气动热试验研究［J］. 空气动力学学报， 2024， 42（1）： 1-12.
	TIAN R Y， GONG H M， CHANG Y， et al. Experimental study on aerodynamic heat of sphere heads in high enthalpy flow［J］. Acta Aerodynamica Sinica， 2024， 42（1）： 1-12 （in Chinese）.
12	DE FILIPPIS F， SERPICO M. Air high-enthalpy stagnation point heat flux calculation［J］. Journal of Thermophysics and Heat Transfer， 1998， 12（4）： 608-610.
13	ZUPPARDI G， VERDE G. Improved fay-riddell procedure to compute the stagnation point heat flux［J］. Journal of Spacecraft and Rockets， 1998， 35（3）： 403-405.
14	LEE S， YANG Y， KIM J G. Evaluation of Fay and Riddell formula under hypersonic flight conditions［J］. International Journal of Numerical Methods for Heat & Fluid Flow， 2023， 33（1）： 14-41.
15	OLIVIER H. An improved method to determine free stream conditions in hypersonic facilities［J］. Shock Waves， 1993， 3（2）： 129-139.
16	曲齐齐. 高超声速钝头体表面热流数值模拟研究［D］. 天津：天津大学， 2017： 24-31.
	QU Q Q. Numerical simulation of surface heat flux of hypersonic blunt body［D］.Tianjin： Tianjin University， 2017： 24-31 （in Chinese）.
17	GUPTA R N， LEE K P， THOMPSON R A， et al. Calculations and curve fits of thermodynamic and transport properties for equilibrium air to 30000 K： NASA Reference Publication 1260［R］. Washington， D.C.： NASA， 1991.
18	ANDERSON J D Jr. Hypersonic and high-temperature gas dynamics［M］. 2nd ed. Reston： AIAA， 2006.
19	SRINIVASAN S， TANNEHILL J C， WEILMUENSTER K J. Simplified curve fits for the thermodynamic properties of equilibrium air： NASA Reference Publication 1181［R］. Washington， D.C.： NASA， 1987.
20	柳军. 热化学非平衡流及其辐射现象的实验和数值计算研究［D］. 长沙：国防科学技术大学， 2004： 53-101.
	LIU J. Experimental and numerical study on thermochemical non-equilibrium flow and its radiation phenomenon［D］.Changsha： National University of Defense Technology， 2004： 53-101 （in Chinese）.
21	李海燕. 高超声速高温气体流场的数值模拟［D］. 绵阳：中国空气动力研究与发展中心， 2007： 38-90.
	LI H Y. Numerical simulation of hypersonic high temperature gas flow field［D］. Mianyang： China Aerodynamics Research and Develop-ment Center， 2007： 38-90 （in Chinese）.
22	刘国庆. 应用数值分析［M］. 北京：化学工业出版社， 2020： 201-220.
	LIU G Q. Applied numerical analysis［M］. Beijing： Chemical Industry Press， 2020： 201-220 （in Chinese）.
23	GUPTA R N， YOS J M， THOMPSON R A， et al. A review of reaction rates and thermodynamic and transport properties for an 11-species air model for chemical and thermal nonequilibrium calculations to 30000 K： NASA Reference Publication 1232［R］. Washington， D. C.： NASA， 1990.
24	CHASE M W Jr， CURNUTT J L， DOWNEY J R Jr， et al. JANAF thermochemical tables， 1982 supplement［J］. Journal of Physical and Chemical Reference Data， 1982， 11（3）： 695-940.
25	MCBRIDE B J， ZEHE M J， GORDON S. NASA glenn coefficients for calculating thermodynamic properties of individual species： NASA/TP-2002-211556［R］. Washington， D. C.： NASA， 2002.
26	IRIMPAN K J， MENEZES V. Stagnation heat flux estimation in spherically blunt axisymmetric hypersonic models［J］. Journal of Aerospace Engineering， 2023， 237（6）：1369-1375.
27	CHADWICK K， CHADWICK K. Stagnation heat transfer measurement techniques in hypersonic shock tunnel flows over spherical segments： AIAA-1997-2493［R］. Reston： AIAA， 1997.
28	MACLEAN M， MARINEAU E， PARKER R， et al. Effect of surface catalysis on measured heat transfer in an expansion tunnel facility： AIAA-2012-0651［R］. Reston： AIAA， 2012.

E-mail：hkxb@buaa.edu.cn

关于我们

期刊社服务

专业学科

封面文章

友情链接

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

步长	壁面温度760 K工况			11个工况
步长	迭代启动次数	平均迭代步数	耗时/s	总耗时/s
0.005	123	1 060.21	92	906
0.01	123	535.20	46	476
0.02	122	264.72	23	249
0.04	151	130.65	15	123

风洞	速度u_∞/（m·s^-1）	静压p_∞/Pa	密度ρ_∞/ （10^-3kg·m³）	壁温T_w/K	球头直径D/cm	试验值/（W·cm^-2）	BLES热流/（W·cm^-2）	相对试验值偏差/%	FR热流/（W·cm^-2）	相对试验值偏差/%
IITB-ST激波风洞^［26］	1 828.00	79.0	2.572	300.00	5	31.8	26.92	-15.35	26.68	-16.10
Calspan 96英寸激波风洞^［27］	1 922.43	47.66	4.619	296.22	1.27	85.99	81.40	-5.34	80.67	-6.19
					1.905	70.01	66.46	-5.07	65.87	-5.91
					3.175	59.58	51.48	-13.60	51.02	-14.37
					7.62	41.09	33.23	-19.13	33.93	-17.43
	2 150.21	369.92	15.055	296.11	1.27	232.24	222.50	-4.19	222.78	-4.07
	2 306.82	1 036.76	28.29	295.83	1.27	400.09	394.30	-1.45	393.98	-1.53
	3 565.86	40.80	0.840	294.78	1.27	291.41	288.48	-1.01	282.64	-3.01
CARDC2 m激波风洞B	2 623.14	959.0	15.062	286.88	3	290.66	295.96	1.82	295.82	1.78
	2 980.00	538.0	7.013	298.00	3	304.25	309.63	1.77	302.64	-0.53
	3 118.67	703.0	8.743	287.69	3	409.40	403.16	-1.52	390.96	-4.50
LENS XX^［28］	5 732.00	180.155	1.16	300.00	7.62	651.00	648.96	-0.31	639.21	-1.81
CARDC高焓膨胀风洞^［11］	6 260.00	18.4	0.44	282.1	5	659.70	638.83	-3.16	613.98	-6.93
	7 750.00	14.3	0.221	281.0	2	1 525.10	1 368.06	-10.30	1 347.06	-11.67
	7 990.00	17.5	0.177	290.0	2	1 574.60	1 340.51	-14.87	1 337.27	-15.07

拟合式	3	4	5	6
P₁	0~2	0~2	0~1	0.7~0.9
P₂	无	无	无	-0.5~0
P₃	-2~2	-2~2	-1~1	0.4~1
P₄	-2~2	-2~2	-1~1	0~0.5
P₅	无	无	-1~1	-0.5~0
P₆	无	-2~2	-1~1	0~0.45

拟合式	FR公式	3参数	4参数	5参数	6参数
总偏差	7.173 16	2.989 87	2.984 06	2.672 90	2.365 93
P₁	0.76	0.87	0.89	0.91	0.77
P₂	无	无	无	无	-0.26
P₃	0.4	0.68	0.72	0.77	0.68
P₄	0.4	0.14	0.17	0.14	0.15
P₅	无	无	无	-0.23	-0.13
P₆	0.4	无	0.04	0.31	0.19

拟合式	FR公式	3参数	4参数	5参数	6参数
总偏差	10.940 1	5.536 83	5.663 09	5.207 14	4.432 80
P₁	0.76	0.87	0.89	0.91	0.77
P₂	无	无	无	无	-0.26
P₃	0.4	0.68	0.72	0.77	0.68
P₄	0.4	0.14	0.17	0.14	0.15
P₅	无	无	无	-0.23	-0.13
P₆	0.4	无	0.04	0.31	0.19

基于边界层理论的球头驻点热流计算方法

Calculation method for heat flow at stagnation point of spherical head based on boundary layer theory

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 13

参考文献 28

相关文章 15

编辑推荐

Metrics

本文评价

[1]	刘宇, 赵淼, 何青松. 充气式减速结构壁面变形对热化学非平衡流场的影响[J]. 航空学报, 2025, 46(1): 630727-630727.
[2]	刘宏康, 陈坚强, 向星皓, 赵雅甜. 改进k-ω-γ转捩模式对不同雷诺数下HIAD的转捩预测[J]. 航空学报, 2023, 44(6): 126868-126868.
[3]	杨国涛, 岳振江, 刘莉. 基于自适应采样的高超声速飞行器气动热全局快速预示[J]. 航空学报, 2023, 44(6): 127391-127391.
[4]	马崇立, 刘景源. 网格对高超声速钝头体表面热流数值模拟结果的影响[J]. 航空学报, 2023, 44(5): 126710-126710.
[5]	曾品棚, 陈树生, 李金平, 贾苜梁, 高正红. 减阻杆与环形喷流组合构型钝头降热数值模拟[J]. 航空学报, 2023, 44(22): 128407-128407.
[6]	刘重晓, 王江峰. 滑移流区化学非平衡流中气动热环境的数值模拟[J]. 航空学报, 2023, 44(16): 127980-127980.
[7]	张振康, 徐万武, 李智严, 叶伟. 钝锥体辐射平衡温度不确定度量化分析[J]. 航空学报, 2023, 44(15): 528673-528673.
[8]	孔祥志, 韩宇峰. 高超声速平板边界层中波纹粗糙壁峰值热流的变化规律及预测[J]. 航空学报, 2023, 44(12): 127769-127769.
[9]	廖东骏, 石安华, 黄洁, 简和祥, 谢爱民, 王宗浩. 来流速度5~7 km/s时CO₂流场中CO离解现象的验证[J]. 航空学报, 2022, 43(S2): 106-114.
[10]	胡守超, 庄宇, 李贤, 江涛. 高超声速气动热标模HyHERM-Ⅰ试验[J]. 航空学报, 2022, 43(S2): 233-248.
[11]	李昊歌, 杨华, 杨雨欣, 陈伟芳. 高超声速升力体迎风面精细化降热优化设计[J]. 航空学报, 2022, 43(S2): 124-137.
[12]	聂春生, 袁野, 马伟, 曹占伟, 于明星. 主动引射气体参数对平板空气舵气动热影响[J]. 航空学报, 2022, 43(S2): 170-179.
[13]	尧少波, 蒋励剑, 赵文文, 卢铮, 吴昌聚, 陈伟芳. 耦合聚类的数据驱动稀薄流非线性本构计算方法[J]. 航空学报, 2022, 43(S2): 40-53.
[14]	傅杨奥骁, 丁明松, 刘庆宗, 江涛, 石润, 董维中, 高铁锁. 轨控系统热喷干扰效应数值模拟[J]. 航空学报, 2022, 43(9): 125941-125941.
[15]	吕俊明, 李飞, 李齐, 程晓丽. 火星大气高温光谱建模与非平衡辐射热流预测[J]. 航空学报, 2022, 43(3): 626551-626551.