热推力器层板换热芯流动与传热研究-research on flow and heat transfer of heat exchanger core of hot thruster laminate.docx

国防科学技术大学研究生院硕士学位论文 国防科学技术大学研究生院硕士学位论文 万方数据 万方数据 目 录 摘 要 i ABSTRACT ii 第一章 绪论 1 1.1 研究工作的背景及意义 1 1.2 发展现状 3 1.2.1 微推进与太阳能热推进 3 1.2.2 层板与层板冷却 5 1.3 本文的主要工作 8 第二章 层板换热芯连续流动的数值模拟 10 2.1 层板换热芯物理模型 10 2.1.1 层板换热芯外形 10 2.1.2 层板换热芯工作状态 11 2.2 流动分区概述 11 2.3 连续流基本方程 12 2.3.1 直角坐标系下的控制方程 12 2.3.2 无量纲形式的控制方程 13 2.3.3 NS 方程的坐标变换 14 2.3.4 NS 方程的离散 15 2.4 层板换热芯内部流场的二维数值模拟 17 2.4.1 二维计算模型 17 2.4.2 网格划分 17 2.4.3 结果分析 17 2.5 层板换热芯内部流场的三维数值模拟 19 2.5.1 计算模型与定解条件 19 2.5.2 结果分析 19 2.6 小结 21 第三章 并行化计算方法在层板换热芯流场计算中的应用 31 3.1 并行的定义 31 3.2 并行算法的性能指标 32 3.3 并行思想及方法验证 34 3.4 层板换热芯流场的三维并行计算 35 第 I 页 第 I 第 II 页 万方数据 3.5 小结 36 第四章 层板换热芯内微小尺度流动仿真 43 4.1 微小尺度传热及流动特点 43 4.2 DSMC 方法的原理 45 4.2.1 基本思想 45 4.2.2 一般步骤 45 4.3 非结构网格在 DSMC 中的使用方法 47 4.3.1 基本思想 48 4.3.2 程序流程 48 4.4 DSMC 中碰撞的计算 49 4.4.1 碰撞对抽样方法 49 4.4.2 碰撞计算 51 4.4.3 碰撞中的能量交换 52 4.5 DSMC 中的并行计算方法 53 4.5.1 并行计算的网格读入 53 4.5.2 并行计算的信息交界面处理 53 4.5.3 并行计算结果的后处理 53 4.6 边界条件 54 4.6.1 来流边界条件 54 4.6.2 对称边界条件 55 4.6.3 真空边界条件 55 4.6.4 固壁边界条件 55 4.6.5 信息交界面边界条件 56 4.7 DSMC 方法的验证 57 4.8 DSMC 方法并行算法数值验证算例 57 4.8.1 计算条件描述 57 4.8.2 计算结果分析 57 4.9 层板换热芯流场的 DSMC 数值模拟 58 4.9.1 计算模型及条件 58 4.9.2 网格生成与应用 58 4.9.3 结果分析 59 4.10 小结 60 结 束 语 69 第 PAGE 第 PAGE IV 页 万方数据 致 谢 70 参考文献 71 作者在学期间取得的学术成果 76 表 目 录 表 1.1 层板应用及时间 6 表 3.1 单双进程计算π的结果对比 33 国防科学技术大学研究生院硕士学位论文 国防科学技术大学研究生院硕士学位论文 万方数据 万方数据 图 目 录 图 1.1 Lamilloy 层板结构示意图 5 图 1.2 Transply 多层壁层板结构示意图 5 图 1.3 层板制作过程 6 图 2.1 单个层板 10 图 2.2 层板内流动方向 10 图 2.3 计算区域外形 22 图 2.4 计算区域网格 22 图 2.5 流场马赫数（截面 0.1mm） 22 图 2.6 流场马赫数（截面 0.15mm） 22 图 2.7 流场压力（截面 0.1mm） 22 图 2.8 流场压力（截面 0.15mm） 22 图 2.9 流场温度（截面 0.1mm） 23 图 2.10 流场温度（截面 0.15mm） 23 图 2.11 流场密度（截面 0.1mm） 23 图 2.12 流场密度（截面 0.15mm） 23 图 2.13 流场流线图（截面 0.1mm） 23 图 2.14 流场流线图（截面 0.15mm） 23 图 2.15 两种截面马赫数对比 24 图 2.16 两种截面压力对比 24 图 2.17 两种截面温度对比 24 图 2.18 两种截面密度对比 24 图 2.19 计算外形 24 图 2.20 计算网格 25 图 2.21 流场马赫数（z=0.01mm,截面 0.1mm×0.1mm） 25 图 2.22 流场马赫数（z=0.01mm,截面 0.1mm×0.15mm） 25 图 2.23 流场马赫数（z=0.07mm,截