三维热化学非平衡流动非结构网格DSMC方法及其应用

研究了过渡流域三维热化学非平衡流动DSMC方法实现的过程.以四面体非结构网格为基本单元,提出一种新型的高效搜索算法,该算法不仅可以快速跟踪模拟分子在网格之间的迁移,而且可以准确判别分子与物面是否相互作用,避免了原有算法中分子表面反射的非确定论判据.设计了适合三维DSMC方法的动态局部时间步长技术,节约了计算时间.将碰撞距离的思想引入到非结构网格上来,有效地消除了网格尺度小于三分之一平均自由程的限制.利用Foaran90的动态分配内存技术编制了适用于任意外形的通用计算程序.最后对全尺寸航天飞机高超声速绕流进行了数值试验,计算结果验证了该算法的可行性及高效性.     

轴对称降落伞稳定下降阶段的流场特性

文章根据降落伞的特点,通过3点假设(伞衣的薄膜假设、降落伞的轴对称假设和流场的定常 假设),将三维复杂流动问题,转化成二维轴对称问题,以节约计算时间。然后,建立轴对称降落伞的流体力 学计算模型,利用数值模拟手段,求解RNG(Renormalization Group)k-epsilon湍流模型下的N—S方程组,获 得与有关单位试验相吻合的计算结果。分析发现带顶孔的轴对称降落伞绕流和亚临界状态下带中心孔的 圆球绕流,在尾流区拓扑结构上几乎完全一致。结果表明建立的轴对称降落伞模型,能够揭示降落伞流场 的本质特性,为进一步研究降落伞流场的流动机理和流固耦合问题打下基础。     

激波在狭缝中绕射过程研究

针对固体火箭发动机药柱存在裂纹的情况，对运动激波在狭缝中绕射传播的非定常流场进行了计算，得到了清晰的流场变化图谱，从流场图中可分辨出激波在狭缝中的传播过程。对不同深宽比的狭缝进行研究，发现流场结构存在明显差别，在深宽比较小时狭缝内压力波振荡较小，当深宽比超过一定值时狭缝内的压力波产生剧烈振荡。研究结果表明，该方法可很好地描述激波在狭缝中的传播过程。     

2D碳/碳复合材料CVI过程的数值模拟研究

本文根据３Ｄ碳／碳复合材料的结构特征及ＣＶＩ工艺的特点建立了孔隙模型和动力学模型，并利用该模型对其ＣＶＩ过程进行了模拟与分析。模拟结果与实验结果的对比表明该数学模型是合理的。     

喷管内横流流场的数值模拟

采用显式的原始ＭａｃＣｏｒｍａｃｋ时间分裂格式和改进的ＴＶＤ ＭａｃＣｏｒｍａｃｋ格式分别求解了二维轴对称喷管内横流流场的Ｎ－Ｓ方，得到了与理论分析结果相近似的流场结果图谱，并对以上两方法得到的结果进行了比较，结果表明：ＴＶＤ的方法具有更高的精度。文中还介绍了一种基于ＴＴＭ方法的改进的网格生成方法，用于生成贴壁性和正交性都比较好的网格，实现了完全由边界网点分布来确定内部网格点的分布。     

斜切反喷管流场的数值模拟

为了在级间分离期间提供反推力，许多固体火箭发动机前端都装有一组斜切反喷管。由于反喷管的气动型面具有许多尖点，并且在超声速区有一个台阶，喷管内存在一系列激波，而且亚声速区和超声速区互相混杂。本文用时间相关法数值模拟了反喷管流场，控制议程用Ｍａｃ Ｃｏｒｍａｃｋ显格式求数值解，边界参数采用物理边界条件和从双特征线方法推导来的有效的特征方法为计算，为于固壁边界点，在不同的区域采用不同的方法对计算方程组求     

再入体激波层和近尾流热辐射实验研究

本文介绍了在弹道靶设备上对再入速度下模型的热辐射进行的实验研究工作。铝、高硅氧和聚碳酸酯三种材料的模型在速度为4.5～6.0(km/s),压力为600—8000Pa条件下做了实验,在258nm—1100nm光谱范围内完成了九个波段的辐射测量,给出了三种材料模型的辐射光谱分布曲线。结果表明辐射量除了受环境条件影响外,还强烈地依赖于模型材料;由烧蚀产物引起的辐射通量主要在可见光和近红外区(465nm以上)。估算了烧蚀产物引起的辐射加热率,对头部和近尾流辐射通量做了比较。     

石英晶体微量天平温度热效应分析

通过对在 KM3模拟器中进行的天平热效应专项试验测量数据分析,利用热力学定律,计算了天平热效应参数,并给出了在实际测量中进行实时温度预估的公式,为天平温度补偿中的温度控制设计提供了基础。     

数值模拟高超音速化学非平衡流的隐式格式

针对高超音速化学非平衡流提出了两种隐式格式.一种是近似因子分解算法对角化形式在隐式TVD(全变差递减)格式上的推广,另一种可视为隐式TVD格式的LU分解形式.数值试验表明,本文给出的算法保持了TVD格式的高分辨率特性,同时计算量大大减少,是高超音速化学非平衡流数值模拟的有效方法.     

部分开放多极磁场中的日冕磁绳灾变

采用球坐标下二维三分量理想MHD模型，研究部分开放多极背景磁场中日冕磁绳的灾变现象．背景磁场由含3个闭合双极场的冕流和带赤道中性电流片的开放场构成，磁绳位于中心双极场的下方，其特性由环向磁通和轴向磁通表征．对给定的环向磁通，存在轴向磁通的一个临界值；对给定的轴向磁通，也存在环向磁通的一个临界值．在该临界值以下，磁绳附着于太阳表面，系统处于平衡状态；该临界值一旦被超越，磁绳将脱离太阳表面向上喷发，说明部分开放多极背景磁场中的日冕磁绳系统存在灾变现象．本文算例表明，灾变点对应的磁能阈值超过对应部分开放场(中心双极场开放，两侧的双极场仍维持闭合)能量约15％，其超过部分可为日冕物质抛射一类太阳爆发提供能源．