三维非结构网格DSMC并行算法及应用研究

基于PC-CLUSTER机群并行体系结构与消息传递库MPI并行环境,研究了三维非结构网格DSMC并行算法。提出一类基于结构背景网格上的非结构网格动态分区策略,保证各子区域的分子数量大致相等,实现计算进程间的动态负载平衡。利用MPI库函数构造了两类符合DSMC并行原理的通讯法:单步通讯法与多步通讯法。采用单控制多数据流(SPMD)以及Master/Slave并行模式,设计了三维非结构网格DSMC整体并行算法。给出了跟踪模拟分子在四面体网格间迁移运动的详细计算过程。最后对全尺寸航天飞机高超声速绕流进行了并行模拟,验证并行算法的有效性。     

一种改进的区域覆盖星座优化设计方法

为缩短区域覆盖星座优化设计的计算时间，对网格点仿真法进行了改进。在统计时间内，以解析法计算卫星覆盖时刻集，并采用遗传算法优化星座参数。仿真计算结果表明，本文改进方法的计算时间仅为优化前的百分之一。     

HTPBⅠ型裂纹粘弹性有限元分析

在简化模型的基础上，推导了复合固体推进剂中Ⅰ型拉伸裂纹的裂尖区域应力、应变与位移的关系式，并进行了粘弹性有限元方法分析。结果表明，有限元分析能获得具一定精度的结果，并与理论分析基本吻合。     

变形铝合金中微量合金元素研究

分析了微量合金元素在变形铝合金中的作用，论述了微量合金元素与杂质的关系。指出微量合金元素与杂质是不能混淆的两个概念。前者的作用不仅与其种类和含量有关，还与加工方法有关。     

二维非结构网格Euler方程的LU-SGS算法

研究了将二维非结构三角形网格重新排序后，用LU-SGS隐式算法计算Euler方程的方法。分析了网格不平衡情况下该方法的可行性和应用效果。计算了不同翼型在不同流动条件下的二维Euler方程，并与四步Runge-Kutta显式方法进行了比较，结果验证了此方法具有良好的计算和收敛效果，以及在粘性计算方面的潜力。运用此方法对某预研型号的头部外形进行了设计和选择。该方法克服了以往隐式方法大量耗费内存的弱点，达到了计算耗时短和占用内存少的统一。     

复合材料的基体混杂——互穿聚合物网络

本文介绍了高聚物基复合材料基体研究的一种最新进展——互穿聚合物网络(IPN),IPN可望满足高级结构复合材料对于基体性能的苛刻要求,是与常用的共混、共聚、后反应、填充与强化等完全不同的新的基体改性方法。本文综述了IPN目前的研究水平,并初步介绍了这种改性的主要机制。     

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

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

Rankine涡作用下激波变形和分叉的数值研究

数值模拟了激波和Rankine涡相互作用，重点研究了波涡相互作用引起的激波结构的变化过程。应用非结构化适应网格下的二阶精度Godunov型的PLM格式来求解Elder方程，对一平面激波与不同涡强Rankine涡相互作用下瞬态激波的运动、发展进行了研究。漩涡与激波间相对强度的大小对激波结构的影响明显，弱相互作用诱导激波变形，强相互作用诱导激波变形并产生分叉。计算结果表明，应用非结构化适应网格的PLM格式所得到的数值解较好地反映了瞬态激波结构的有关特征和信息。     

结构/直角切割网格下时域有限体积法在计算电磁中的应用研究

 采用有限体积法直接求解时域麦克斯韦方程组，研究了二维圆柱、翼型和三维球体、立方体、双球体等典型完全导电目标的电磁散射问题。空间离散使用近似黎曼解构建网格通量，并针对二维结构网格、三维自适应直角切割网格分别给出重构方法，时间方向采用二步龙格库塔法。计算与理论、实验或文献结果符合良好。     

三维贴体坐标系下燃烧室热态流场的大涡模拟

大涡模拟三维贴体坐标系下环形燃烧室火焰简热态紊流瞬态流场。利用椭圆方程方法生成三维贴体网格，计算中采用k方程亚网格尺度模型估算亚网格紊流粘性；亚网格EBU燃烧模型估算化学反应速率；热通量辐射模型估算辐射通量。并在非交错网格系下采用sIMPLE算法和混合差分格式求解离散方程，利用壁面函数处理固壁边界条件。计算结果与实验结果的比较表明，采用大涡模拟方法能更真实反映环形燃烧室火焰简内紊流化学反应流气流结构和燃烧过程。