简化模型(指数大气)
$\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}$
使用阿波罗、航天飞机、星舰预设对比再入加热效果。调节速度、角度和弹道系数,探索最大G值、热流密度和壁面温度。
$\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}$
大气密度随高度变化的简化模型(指数大气):
$$\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公里,速度11公里/秒)开始尝试。
最后,需注意此计算属于“局部”评估。飞行器整体的温度分布以及热量在结构内部如何传导(热传导分析),属于更复杂的CAE软件(热流耦合分析:CHT)的范畴。此工具可视为TPS(热防护系统)设计的第一步,用于“筛查”判断“飞行器哪个部位热环境最严苛”。
这种气动加热计算的思路,实际上不仅应用于航天器再入,还广泛适用于各类高速飞行器的设计。首先最直接关联的是高超音速飞行器(马赫5以上)的研发。例如,搭载超燃冲压发动机的实验机或未来的高超音速客机。即使在大气层内以高超音速巡航,机头或翼前缘也会产生堪比再入过程的气动加热。这里学到的热流密度公式 $\dot{q} \propto V^3$ 正是速度危险性的真实写照。
另一个领域是火箭上升段的气动加热。与再入过程相反,火箭发射穿越大气层时,在最大动压(Max Q)附近,箭体表面也会受热,尤其是顶端的整流罩和翼部。此外,研究陨石或太空碎片进入大气层的“天体物理学”和“太空碎片应对”领域,同样以相同的物理模型为基础。陨石能否未燃尽而坠落地面,恰恰取决于此工具涉及的“弹道系数”与“再入角”。
进一步拓宽视野,还会关联到等离子体物理学。速度极高时,被压缩的空气分子会分解、电离形成等离子体状态(电离气体)。这正是导致无线电中断的“通信黑障”现象的根源。此工具计算的高温环境,也可作为思考该等离子体层产生条件的切入点。
熟悉此工具后若产生“想了解更多!”的想法,可尝试进入下一阶段。首先,理解对“运动方程”进行时间推进求解至关重要。此工具的计算是瞬时性的,但实际现象是连续的。再入体的运动由重力、空气阻力(气动阻力)以及升力(如果存在)的平衡决定。通过数值积分(例如欧拉法)求解,即可计算出高度和速度随时间的变化,即“轨道”。使用Excel或Python编写一个利用重力加速度 $g$ 和阻力 $D = 0.5 \times \rho \times V^2 \times C_d \times A$ 的简单程序,是最好的练习方式。
其次,精细化大气模型。此工具使用的指数模型虽然简便,但实际大气随高度变化,其温度和成分都会改变。在程序中集成更贴近现实的“美国标准大气”等数据,可显著提升计算精度。同时,学习更复杂但精度更高的热流密度公式如“Fay-Riddell公式”等,能更接近实际工程分析方法。
最终,建议通过专业书籍或论文学习“热防护系统(TPS)的材料科学”。烧蚀材料(通过烧蚀带走热量的材料,阿波罗任务中使用)与可重复使用的辐射散热瓦(航天飞机使用)的根本区别是什么?星舰采用的不锈钢外壳如何通过辐射散热?当能够计算气动加热的“量”之后,另一个深邃的技术世界——如何与这些“热”作斗争——正等待着你的探索。