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

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

非结构同位网格SIMPLE类算法收敛性能比较

以非结构网格的SIMPLE算法为基础,将算法扩展为SIMPLEC算法.利用30°角的斜方腔流动计算成果,分析了非结构同位网格的SIMPLE/SIMPLEC算法的收敛性能;比较了因网格的非正交而引入的非正交项的取舍对该算法收敛性能的影响;并采用显式校正步法对SIMPLEC算法进行了显式校正.比较表明,在非结构同位网格SIMPLEC算法中可忽略非正交项,但有必要对压力作亚松弛.显式校正步法可显著地加速在非结构网格上求解N-S方程的收敛性能,而且在不同的松弛因子组合下,均有较好的收敛速率.     

子母弹分离过程的数值模拟方法

以分区拼接网格数值模拟方法为基础，通过气动方程和弹道方程联立求解数值模拟了子弹从母弹中分离并穿越母弹激波的过程，研究了子母弹干扰流场的气动特性，分析研究了子弹各初始分离状态（如初始抛撒速度、初始姿态角、初始姿态角速度等）对分离过程的影响。     

用Newton子迭代方法计算前飞旋翼粘性绕流

通过求解Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场，为了同时保证计算的时间精确性和计算效率，时间推进格式采用了双时间推进方法，在该方法中，子迭代过程由十分高效的LU-SSOR方法完成，且使子迭代过程成为Newton子迭代，空间上应用中心平均的有限体积法进行离散，为了模拟前飞桨时间的相对运动，网格布局采用了运动嵌套网格方法，应用本文方法对-悬停流场进行了数值计算，计算结果与实验吻合较好，尽管缺乏实验数据的验证，对-有升力前飞状态的数值模拟结果是可信的。     

关于非结构网格生成中的背景函数

背景信息对于用阵面推进法生成非结构网格是一个关键的环节.文章基于众多文献中提到的背景信息提供方法,提出一种更为实用的背景信息提供技术.此方法不仅避免了不必要的机时,内存的消耗,又使网格具有光滑性,网格的生成具有灵活性,方便性的特点.基于生成的网格采用有限体积法研究无粘高超音速化学非平衡流的流场分布,生成的网格能够成功的运用于高超音速化学非平衡流的数值计算。     

柔性扑翼非定常流场的数值计算方法

提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据扑翼的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。非定常可压缩Navier-Stokes方程在非结构的动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代求解。几个扑翼算例的结果表明,该方法充分利用了Delaunay图映射网格变形方法的高效率,同时也发挥了嵌套网格处理大幅运动的优势;用于既有整体大幅扑动又有局部小变形的柔性扑翼流场计算,可取得令人满意的精度和效率。     

带稳定板装置弹射座椅偏航稳定性能研究

针对国内弹射救生领域现状,设计了一种适合于国内某型座椅的弹射救生系统稳定板。运用计算流体力学(CFD)数值计算方法,对改装稳定板前后的人椅系统进行了高速流场下的静态和动态特性分析,并与风洞试验结果进行对比,两者吻合较好。改装稳定板前后的数值计算对比表明:在不对俯仰和滚转性能产生不利影响的情况下,稳定板的改装对人椅系统的偏航性能有比较明显的改善。六自由度动态计算比较真实地模拟了弹射出舱的过程,也进一步验证了稳定板对改善偏航稳定性能所起到的作用。数值计算完全可以应用于弹射座椅的气动研究,同时稳定板的设计对工程上的应用具有一定的参考价值。     

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。      

进口速度分布对短突扩压器性能的影响

在任意曲线坐标系下,对环型燃烧室的三维流场进行了数值计算,研究了四种不同的进口速度分布对扩压器性能的影响。计算结果表明,进口速度大小和方向都对扩压器性能有重要影响,进口速度分布均匀以及方向贴近扩压器壁面,有利于提高扩压器的性能。计算中采用标准ｋ－ε双方程紊流模型,采用控制容积法进行离散,在非交错网格体系下用ＳＩＭＰＬＥ法求解     

面向深空探测的30cm离子推力器寿命预估

为评价30cm离子推力器的寿命,提出了有限寿命考核结合栅极仿真模型的推力器寿命预估方法,利用粒子-蒙特卡洛（PIC-MCC）方法建立了栅孔溅射腐蚀模型,开展了30cm离子推力器寿命预估研究,分析了栅孔刻蚀速率、单孔电场及离子引出特性,给出了30cm推力器寿命预估值。结果表明,栅孔直径仿真与实测值一致性较好,误差在20%以内,基于每个寿命小节栅孔实测值对模型的修正是有效的;根据仿真结果,10000h寿命考核后,减速栅最先失效,30cm离子推力器在最大功率（3kW）工况下预估寿命为37540h,能够满足小天体探测任务对电推进系统的长寿命需求。