三维非结构网格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格式所得到的数值解较好地反映了瞬态激波结构的有关特征和信息。     

"神舟二号"大气密度探测器的探测结果Ⅱ.在太阳和地磁扰动期间高层大气密度的变化

选用了Sz—2大气密度探测器在2001年2—4月间的探测数据,进行日照和阴影区域热层大气密度对太阳和地磁活动程度的响应变化的探讨．结果表明,日照区大气密度峰值主要随F10.7值而变,在地磁扰动期间,阴影区大气密度对扰动的响应更明显,通常响应变化开始于高纬度地区,然后向低纬度地区推移．     

武汉上空Na层季节和夜间变化的观测研究

根据武汉大学激光雷达2001年3月到2002年10月的观测数据,对武汉上空(30.5°N,114.4°E)的背景Na层结构进行了研究,并与其他观测站尤其是Urbana(40°N,88°W)的观测结果进行了比较。发现武汉上空11月的Na层柱密度最大,是4月的2.5倍,质心高度的极小值出现在5月和11月,极大值在4月、9月和12月;全年平均Na层均方根宽度为4.47km.Na层各参数呈现出周期性变化规律,其中Na层均方很宽度呈现出明显的半年周期性变化规律。夜间Na层柱密度呈现总体增加的趋势,质心高度在入夜后立即迅速上升,在2000LT以后上升速度趋缓,凌晨0500LT以后出现下降;Na层均方根宽度在夜间约有1.3km的起伏变化,它入夜后迅速减小,在午夜0035LT时出现一低谷,然后随时间逐渐变宽。     

单流体二维太阳风结构的磁流体力学模拟

在添加动量项的条件下，对单流体二维磁流体力学方程组进行模拟，得到了子午面上的太阳风结构，结果表明，添加动量项的形式及其被加入的位置对远区太阳风速度和质子密度有重要的影响。本文在适当的区域加入合适的动量项得出了远区太阳风速和质子密度与Ulysses观测基本符合的结果。文中给出了较合适的动量添加区域为3.5-10Rs(Rs为太阳半径）。