三维烧蚀内部热响应数值计算研究

对热解型碳化复合材料三维烧蚀内部热响应数值计算关键技术进行了研究。采用碳化层—热解面—原始材料层模型,将热解气体与碳化层之间的对流换热处理为源项,通过有限元法建立移动边界条件下温度场求解方程组,采用Gauss-Seidel迭代法计算热解气体质量流量和温度场。同时,研究和分析了三维烧蚀移动边界处理方法以及动网格生成方法。由于每个时间步都需要网格重划,烧蚀热防护数值计算对存储效率和计算效率要求较高,本文研究了内部热响应计算中影响存储效率和计算效率的主要因素,并提出了相应的压缩存储方案和求解方案。计算结果表明,移动边界处理方法准确合理;存储方案的存储效率较高;保持刚度矩阵和形函数矩阵正定对称性可以加快温度场计算的收敛速度。

再入钝锥体烧蚀热防护内部热响应的数值仿真

研究了烧蚀热防护系统内部热响应的计算模型和计算方法.采用碳化层-热解面-原始材料层模型,建立碳基材料内部热响应物理模型和数学模型,利用有限元法分析和计算再入目标热防护系统轴对称内部热响应.着重研究和分析了轴对称烧蚀过程中热解气体质量流率计算方法和传热机制.将热解气体与碳化层之间的对流换热处理为源项,通过保证刚度矩阵和形函数矩阵的正定对称性可以加速温度场计算收敛.计算表明:热解气体的质量流量主要由厚度方向构成,占80%以上;头部驻点附近最大烧蚀厚度接近10mm,需要采用抗烧蚀能力强的碳-碳材料,身部烧蚀量小于2mm,可以采用密度较小的碳-酚醛材料.

热解气体流动的二维烧蚀热防护数值仿真研究

研究了热解型碳化复合材料二维烧蚀的计算模型和计算方法。采用碳化层—热解面—原始材料层模型;将热解气体与碳化层之间的换热作为热载荷处理;计算了热解气体在碳化层中的流动和传热过程,获得了热解气体质量流量分布从而为烧蚀材料温度场计算提供边界条件。计算结果表明热解气体质量流量在碳化层明显成二维分布,不同位置的流量相差很大,并且对碳化层温度和烧蚀率有很大影响,因此,内部热响应计算需要考虑热解气体在碳化层中的流动和换热过程。