高超声速弹头气动热工程算法与数值传热

论文摘要

采用有限体积法,在非结构化网格上离散稳态导热方程,耦合迭代计算法向扩散项和切向扩散项,以满足高精度离散格式的要求。运用可视化语言Fortran95模块化完成了数值传热计算程序的编制,采用C++编写网格文件的接口程序,并对来自与Tecplot和GAMBIT的网格文件重新进行读取与编号,实现了网格前处理的兼容性。在标准大气环境下,对球锥型弹头在5～8马赫高超声速来流时,采用工程算法计算弹头表面热通量。基本公式为费-里德尔（Fay-Riddell）平衡边界层驻点热流密度公式,李斯（Lees）层流热流分布公式。在二维零攻角来流情况下,考虑了两种流动状况,一种是完全层流无转捩；另一种为采用Batt边界层转捩准则来确定边界层转捩点。完全湍流区采用湍流热流计算公式确定其热流值,转捩过渡区热流由层流传热与湍流传热的加权平均值计算得到。在三维有攻角来流情形下,弹体周向热流不再对称,结合上面的基本公式与经过适当修正的等价锥法,分析攻角对层流传热的影响。然后,采用一种简化的轴对称比拟法计算相同攻角来流条件下的弹体热流分布,并与等价锥法以及实验值相互比较,并得到了令人满意的一致性。最后,在验证了自编气动热程序的正确性后,将上述工程算法得到的弹体表面热流值,线性化后作为第三类边界条件耦合到弹体热分析程序中,对弹体的温度场进行了数值模拟。计算结果表明,湍流边界层传热会导致弹头高温区扩大,表面温度上升；随着攻角的增大,迎风面的热流值会增大,而且温度也会上升；而背风面的热流值会随着攻角的增大而减小,温度也会相应的有所降低。上述所得结论将为弹箭气动热设计和热防护提供了有价值的参考。

论文目录

• 摘要
• ABSTRACT
• 目录
• 图表清单
• 1 绪论
• 1.1 科学意义与背景
• 1.2 国内外研究现状
• 1.3 高超声速流动基本特征
• 1.4 本文主要研究工作
• 2 高超声速气动热工程算法
• 2.1 弹头表面压力分布工程计算
• 2.2 气动热工程算法基本公式
• 2.2.1 驻点对流热通量
• 2.2.2 圆球头部热通量分布
• 2.2.3 锥体部热通量分布
• 2.3 高超声速激波后空气输运参数计算
• 2.4 本章结束语
• 3 非结构化网格数值传热计算方法
• 3.1 网格描述
• 3.2 几何计算
• 3.3 方程离散
• 3.3.1 内部控制体离散
• 3.3.2 边界控制体离散
• 3.4 三维稳态传热算例
• 3.5 本章结束语
• 4 边界层转捩对弹头温度场的影响
• 4.1 湍流热流计算
• 4.2 边界层转捩准则
• 4.3 转捩过渡区热流计算
• 4.4 数值计算算例和分析
• 4.5 本章结束语
• 5 攻角对弹头温度场的影响
• 5.1 等价锥体法求解有攻角的弹体热流分布
• 5.1.1 球形头部等效热流计算
• 5.1.2. 锥体部分等效热流的计算
• 5.2 算例的计算与分析
• 5.3 热流的轴对称比拟方法
• 5.4 工程算法程序
• 5.5 本章结束语
• 6 总结与展望
• 6.1 全文总结
• 6.2 工作展望
• 致谢
• 参考文献
• 攻读硕士学位期间发表的论文