在非结构化网格中,采用有限体积方法离散非稳态扩散方程。对于非稳态项采用全隐格式,对于扩散项,采用法向扩散项(一次扩散项)和切向扩散(二次扩散项)进行相互交叉迭代的算法。在此基础上运用可视化Fortran77/90/95语言模块化完成了二维和三维程序的编制,并且与通用商业软件(Tecplot和GAMBIT)的网格生成进行了对接。在弹头热源区内,采用单组法求解核反应材料与反射层中的热中子通量密度,然后运用类似方法确定核反应材料两区中热中子产生的热功率分布。该种分布分为两种情况,一种是热功率是空间位置函数(径向为贝塞尔函数分布,轴向为余弦分布)；另一种是温度的一阶线性分布函数(发生相应的正负反应)；并对其进行数值传热计算。以上计算模拟结果显示,加设反射层后不仅提高了比威力,而且有利于堆芯小型化的设计。最后,对圆球锥形弹头头部外流场的气动热计算,选择工程算法计算热弹头部驻点热通量、半球和锥体部分外表面热流通量,并将其线性化作为第三类边界条件附加在弹体表面。在弹头内部设有不同半径的球形热源区,它具有径向正弦函数分布的释热率。将气动热的影响和内部余热的影响相耦合,数值模拟弹头内部温度场的分布,得到高温区和动态温度场分布规律。计算模拟结果表明,随着内部热源项半径的增大,弹头内部的高温区有所增大。热源区产生的余热和外部气动加热是弹头热环境分布的重要影响因素。这些结果可为含高热源弹头的优化设计和小型化以及热防护设计提供有益的参考。