一类二维非结构网格DSMC方法的实现策略及其应用

研究了一类二维非结构网格DSMC方法的实现策略。在数据结构方面，设计了局部化的数据组织方式，节约了内存与计算时间。发展了一种跟踪模拟分子迁移的算法，该算法仅需少量的逻辑运算与代数运算，不仅可以快速跟踪模拟分子在网格之间的迁移，而且可以准确判别分子与物面是否相互作用，搜索过程中的附带信息给出了分子与物面碰撞的精确时间与位置，避免了重新计算。引入碰撞距离的思想，既减少存储又保证正确的模拟结果。在程序编制过程中，我们充分展现了Fortran90高级语言的主要特性，引入动态数组、指针、链表以及派生类型数据，编制了计算程序。最后对过渡流域高超声绕流进行了数值试验。     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法，前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率，而后者保证了网格间边界附近流场求解精度；本文非定常Euler方程的求解采用有限体积法，通量计算采用Osher格式，时间积分采用双时间步长推进。为验证本文提出的方法，分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率，保证了网格边界处的精度，具有较高的工程应用价值。     

非结构聚合多重网格法流场数值模拟研究

 在非结构网格上应用聚合多重网格技术来加速 Euler方程的收敛过程。设计了一种高效率的网格聚合方法。采用四重网格,在每一层网格上采用有限体积法进行计算。通过对多段翼型绕流的数值求解验证了多重网格加速收敛的高效性。     

基于非定常低速预处理和DES的三角翼数值模拟

发展出一套基于低速预处理技术采用SA-DES和SST-DES模型,对三维低速非定常流动进行数值模拟的方法.采用非结构混合网格有限体积方法求解非定常流场,时间离散采用基于LU-SGS隐式格式的双时间步长方法.对非定常低速预处理方法进行了推导和分析,并应用该方法对大攻角三角翼的非定常流场进行了数值模拟,对不同站位压强分布与试验值和参考文献值进行了对比,吻合较好.     

跨声速欧拉方程的一种隐式有限元解法

本文用有限体积Ｇａｌｅｒｋｉｎ法求解了跨声速二维欧拉方程。求通量运用了Ｏｓｈｅｒ格式，通过外插把原来只有一阶精度的逆风格式改进成二阶精度。并采用基于一性化流量的隐式格式和非结构化网格，使程序具有较高的收敛效率。     

折叠式群伞充气接触的膜索非线性有限元算法

降落伞由折叠到打开的充气过程和多束降落伞捆绑在一起充气的过程都伴随着一个高度非线性的柔性结构接触问题,因此模拟降落伞群伞流固耦合问题必须首先解决降落伞群伞非线性结构接触的数值模拟问题。对于降落伞这类柔性瞬变非线性织物结构,数值模拟织物系统的接触现象以及预测接触结果对降落伞工作状况的影响对于降落伞群伞设计具有很重要的指导意义。基于三维降落伞膜索非线性有限元编程模拟技术设计了一种针对单伞和群伞非线性动力系统的接触搜索算法,推导设计了针对接触随机性的接触切线刚度矩阵计算方法。针对降落伞群伞非线性有限元计算的庞大数值计算量,设计编写了一种基于消息传递接口(MPI)通信协议的膜索结构并行接触模拟FORTRAN程序,测试了该数值模拟程序的并行计算效率。针对高度折叠的C-9降落伞群伞充气问题,使用PC-Cluster计算机群进行了数值仿真模拟,验证了该接触非线性有限元程序的计算效果,预测了C-9降落伞群伞的接触过程,并分析接触现象对降落伞群伞充气工作的影响。     

旋翼干扰流场数值模拟中的新型Blend低速预处理方法

针对定常低速预处理在双时间步长内收敛性差、旋翼干扰流场数值模拟中非定常低速预处理自动退化为无预处理的问题,在非定常低速预处理方法的基础上,采用量级分析方法研究了对流项和黏性项量级,通过组合定常和非定常预处理特征值矩阵,建立了一种新型Blend低速预处理方法。通过圆柱非定常扰流、前飞直升机旋翼/机身干扰非定常流场和直升机/舰船干扰非定常流场的数值模拟,对各低速预处理方法的收敛性和预测精度进行了细致研究。结果表明,本文方法解决了在旋翼干扰流场数值模拟时非定常低速预处理自动退化为无预处理的问题,相比定常低速预处理大幅提高了收敛效率。     

三维机翼气动结构多学科优化方法

采用Euler方程、连续伴随优化方法、自由形面变形技术（free-form deformation,FFD）以及一种四边形线性壳单元模型对三维机翼气动结构多学科优化方法开展了研究。目标函数为阻力系数和机翼翼盒结构质量的加权和,设计变量为FFD控制点以及10段蒙皮的厚度。气动方面升阻力系数对设计变量的梯度由伴随方法求得,结构方面翼盒质量对于蒙皮厚度的梯度可直接计算,应力的梯度采用有限差分方法获得。对跨声速大展弦比机翼进行气动结构多学科优化,并将该计算结果与机翼气动单学科优化设计后的静气弹计算结果进行对比,结果表明:该气动结构多学科优化方法能够在满足约束条件下,实现阻力系数和翼盒质量同时降低,保证优化结果的合理性。      

非结构多重网格流场数值模拟研究

本文在非结构网格上应用多重网格技术来加速Euler方程的收敛过程。此方法在一系列粗网格和细网格之间进行值的传递，并且粗网格和细网格之间结构上无关。在每一重网格上用有限体积法进行计算。采用本文的方法数值求解了NACA0012翼型的绕流。