大气层再入计算器 返回
航空航天工程模拟器

大气层再入/气动加热计算器

使用阿波罗、航天飞机、星舰预设对比再入加热效果。调节速度、角度和弹道系数,探索最大G值、热流密度和壁面温度。

预设方案
再入参数
再入速度 V₀
km/s
再入角度 γ
°
弹道系数 β
kg/m²
头部半径 R_n
m
辐射率 ε
计算结果
最大G值 (g)
最大热流密度 (MW/m²)
壁面温度 (K)
累积热量 (MJ/m²)
Traj
Vel
理论与主要公式

$\rho(h) = \rho_0 e^{-h/H}$(H=7000 m)

$\dot{q}\approx C_h \rho^{0.5}V^3 R_n^{-0.5}$

平衡壁温: $\varepsilon\sigma T_w^4 = \dot{q}$

什么是大气层再入与气动加热

🙋
“气动加热”是什么?听起来好吓人,飞船真的会被烧化吗?
🎓
简单来说,就是飞船以超高速冲进大气层时,前方的空气被剧烈压缩和摩擦,产生几千度高温,像被巨大的喷灯持续烧烤一样。比如阿波罗飞船返回地球时,隔热罩表面温度超过了2700°C,比火山岩浆还热!在实际工程中,我们最关心的就是这股热流有多强,以及飞船表面温度有多高。你可以试着拖动上面模拟器的“再入速度”滑块,从7km/s调到11km/s,看看热流密度会飙升多少倍,就能直观感受到速度的威力了。
🙋
诶,真的吗?那为什么飞船要做成圆圆的钝头,而不是像子弹一样尖尖的?尖的不是阻力更小吗?
🎓
这是个绝佳的工程悖论!简单来说,钝头设计(头部半径大)恰恰是为了“隔热”。尖头会让激波紧贴表面,所有热量都传给飞船;而钝头能让脱体激波远离,把大部分热量挡在外面,让热空气从旁边流走。你可以马上在模拟器里试试:把“头部半径”从0.1米调到10米,你会发现最大热流密度反而会下降!工程现场常见的就是航天飞机和阿波罗那种大钝头,这就是“以退为进”的智慧。
🙋
原来如此!那“弹道系数”又是什么?听起来好专业。
🎓
别被名字吓到,简单理解就是飞行器“有多重、有多钝”。弹道系数 $\beta = m/(C_d A)$,质量大、阻力小(横截面积小)的物体,$\beta$值就高,像一块石头,空气很难让它减速。改变参数后你会看到,高$\beta$值的飞行器(比如某些返回舱)会像石头一样快速砸进大气层深处,减速和加热都集中在很短时间,峰值很高;而低$\beta$值的飞行器(比如航天飞机)则像滑翔机一样,在高空慢慢飘着减速,加热过程更平缓。试试切换“阿波罗”和“航天飞机”预设,对比一下它们的弹道系数和热流曲线吧!

物理模型与关键公式

大气密度随高度变化的简化模型(指数大气):

$$\rho(h) = \rho_0 e^{-h/H}$$

其中,$\rho_0$是海平面参考密度(约1.225 kg/m³),$h$是高度,$H$是大气标高(约7000米)。这个公式意味着,高度每增加7公里,大气密度就减少到原来的约37%。

驻点热流密度的工程估算公式(Chapman公式简化形式):

$$\dot{q}\approx C_h \rho^{0.5}V^3 R_n^{-0.5}$$

$\dot{q}$是热流密度(W/m²),$C_h$是热传导系数,$\rho$是当地大气密度,$V$是飞行速度,$R_n$是头部曲率半径。这个公式清晰地告诉我们:热流与速度的三次方成正比(速度翻倍,热流变8倍!),与头部半径的平方根成反比。

辐射平衡壁面温度:

$$\varepsilon\sigma T_w^4 = \dot{q}$$

当表面温度足够高时,它以辐射形式散发的热量($\varepsilon\sigma T_w^4$)等于气动加热带来的热量($\dot{q}$),温度就稳定了。$\varepsilon$是表面辐射率(“黑不黑”),$\sigma$是斯特藩-玻尔兹曼常数。这就是为什么飞船隔热罩要做成高辐射率的“黑色”。

现实世界中的应用

载人飞船返回:阿波罗指令舱、神舟飞船、龙飞船的返回舱都采用大钝头+烧蚀隔热罩的设计。工程师利用这些公式精确计算热防护系统(TPS)的厚度,确保宇航员在经历几分钟的“火球”考验后安全回家。

航天飞机轨道器:航天飞机采用低弹道系数的升力体设计,实现更平缓的再入和高达2000°C的机翼前缘温度。其表面覆盖了数万块可重复使用的硅基防热瓦,对热流和壁温的预测精度要求极高。

星际探测器进入:如“好奇号”火星车进入火星大气时,速度高达5.9 km/s。由于火星大气稀薄,减速和加热过程更为剧烈,需要结合本文的模型进行特殊修正,以设计能承受极端条件的隔热罩。

高超声速飞行器:正在研发中的高超声速滑翔飞行器(如HTV-2),其再入走廊极窄,对气动加热和热管理是核心挑战。精确的热流预测直接关系到飞行轨迹规划和热防护材料的选型。

常见误解与注意事项

开始使用此工具时,有几个需要特别注意的要点。首先,“热流密度大的位置 ≠ 飞行器温度最高的位置”。虽然鼻锥尖端确实承受着极高的热流密度,但热量输入的“强度”由热流密度决定,而部件实际升温到多高温度则取决于材料的热容和散热能力(导热系数)。例如,即使热流密度很高,如果结构能快速将热量传导至内部或向后耗散,表面温度反而可能保持在较低水平。相反,如果使用隔热性好的材料,热量可能积聚导致温度缓慢上升,从而损坏内部设备。模拟器中的“平衡壁面温度”终究只是“经过足够长时间后,输入热量与辐射散热达到平衡状态”的理论值。

其次,参数设置的陷阱。“速度”与“高度”并非独立变量。在实际再入过程中,随着高度下降,空气阻力会使速度急剧降低。此工具为简化起见,允许输入某一瞬间的“条件”进行计算,但在实际工程中,必须通过“轨道计算”追踪整个“轨道”的时间变化。例如,高度70公里处秒速7公里与高度40公里处秒速7公里的状态,大气密度相差百倍以上,热流密度值也截然不同。使用工具时,请养成结合实际考虑高度与速度组合的习惯。建议从“阿波罗”预设的初始值(高度120公里,Velocity11公里/秒)开始尝试。

最后,需注意此计算属于“局部”评估。飞行器整体的温度分布以及热量在结构内部如何传导(热传导分析),属于更复杂的CAE软件(热流耦合分析:CHT)的范畴。此工具可视为TPS(热防护系统)设计的第一步,用于“筛查”判断“飞行器哪个部位热环境最严苛”。

使用指南

  1. 选择再入速度(单位:km/s,范围7.8-12.0),对应低轨卫星、阿波罗(11.2)或星舰(9.5)预设
  2. 输入比热比γ(空气取1.4,高温电离气体取1.25-1.3)和动压系数β(0.85-0.95)
  3. 点击计算获得最大G值、热流密度、壁面温度、累积热量四项关键指标

具体计算示例

阿波罗13号再入工况:V0=11.2km/s,γ=1.3(高温空气),β=0.92,RN=1.5m。计算结果:最大G值8.6g,最大热流密度3.8MW/m²,壁面温度(TPS烧蚀层)2100K,累积热量52MJ/m²。相比之下,国际空间站返回舱(V0=7.8km/s)热流仅0.8MW/m²,壁面温度1200K。

实务注意事项

  1. 钝体再入时驻点热流与(ρV³)^0.5成正比,速度增加10%热流增加约50%;星舰采用铜合金TPS,允许壁面温度达2700K
  2. 低轨再入(V0≤8.0km/s)以动压为主(Q~100kPa);深空再入(V0≥11km/s)热流主导,需考虑辐射冷却和催化复合
  3. 实际工程需耦合轨迹计算(高度-速度曲线)和形状系数,钝球体β=0.85-0.90,椎体β=0.70-0.80