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

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

汽车超车过程的空气动力特性研究

采用稳态和瞬态两种方法,对两辆车超车过程中外部流场的空气动力特性进行了数值模拟研究。采用网格变形和局部网格重构相结合的动网格技术,并耦合刚体动力学方程和PISO算法实现了汽车超车过程的三维瞬态数值模拟。通过稳态和瞬态计算,探讨了超车过程中两车纵向相对位置对气动特性的影响以及气动力对汽车行驶稳定性的影响。研究发现超车过程中被超车的阻力、侧向力和横摆力矩都有最大值,并找出了其变化规律。通过对稳态与瞬态计算结果的比较,可以看出:本文采用的瞬态模拟方法能够充分捕捉超车过程中的动态效应,优于传统的稳态模拟。     

混合式CRP吊舱推进器水动力性能数值模拟及试验

建立混合式CRP吊舱推进器数值模型并置于数值空泡水洞内,结合RANS方程和SST k-ω湍流模型,运用滑移网格方法对混合式CRP吊舱推进器在均匀流场中水动力性能进行了数值模拟,得到了敞水性能预报结果,并分析了混合式CRP推进器的非定常力;同时,利用空泡水洞、吊舱动力仪及长轴动力仪对混合式CRP吊舱推进器敞水性能进行了试验研究,通过3组重复性试验表明了试验方法的稳定性;并将3组试验平均值与将数值模拟所得的敞水性能结果进行了比较,彼此吻合较好。结果表明,本文所用试验方法和数值计算方法对吊舱推进器模型敞水特征的预报具有较高的可信度,能达到工程应用要求的要求;和国外学者在空泡水洞内得到的不附带吊舱相同对转桨敞水试验数据进行比较,分析了吊舱对前后桨敞水性能的影响。结果表明,吊舱对前桨的敞水性能基本没有影响,但对后桨的敞水性能影响显著。     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

基于结构网格点的Euler方程无网格自适应求解算法

提出了一种自适应无网格Euler方程求解算法。构造了基于无网格点云结构的布点加密技术,实现了借助压力梯度变化准则的流场局部加密;借助流场结构化剖分技术,给出初始无网格节点分布,并采用四步显式Runge—Kutta时间格式推进计算,求解了绕翼型的典型流动问题。数值算例表明,本文提出的方法能有效地提高流场的分辨率,如捕捉激波等流动特征清晰明了。     

三维非结构网格的生成及其应用

采用阵面推进法生成了高质量的三维四面体非结构网格,该方法的特点是在生成单元的同时引入结点,与其他预布结点的方法相比,其好处是能更好地控制网格单元尺度在整个计算域的合理２并提高网格质量。用Ｄｅｌａｕｎａｙ方法连接存储有控制参数的结点生成了背景网格。网格生成过程中采用了合理的数据结构,有效地提高了程序运行效果,最后给出了Ｍ６机翼的数值算例。     

基于网格重生成的涡轮叶片变形传递方法

在涡轮叶片多学科优化中需要在气动模型和结构模型之间进行耦合信息(气压、温度和变形等数据)传递,来考虑学科耦合的性质.采用修正几何外形和网格重生成的方法实现了叶片结构变形向气动模型传递.传递时在叶片几何模型上增加若干外形控制线,叶片的外形由这些控制线的形状来决定.根据叶片结构变形,调整控制线节点位置坐标,以修改叶片外形控制线的形状,可以实现变形向叶片几何模型的传递.进一步,重新生成气动网格,可以得到结构变形后叶片的气动网格模型.      

滚转状态下卷弧翼火箭弹气动特性的数值模拟

采用滑移网格技术对卷弧翼火箭弹的滚转特性进行数值模拟,并对模拟所得的马格努斯力与风洞实验数据进行分析比较。误差分析证明了滑移网格技术模拟火箭弹滚转特性的准确性;模拟结果与风洞实验数据对比表明,采用细长体外形效果显著,升阻力系数受滚转特性影响较小,因而在本次模拟的火箭弹滚速下,可不考虑滚转来模拟火箭弹的升力系数和阻力系数,得出火箭弹滚转运动对气动参数的影响。     

基于RBF插值边界条件的笛卡尔网格方法

介绍了一种基于RBF(Radial Basis Function,径向基函数)插值处理物面边界条件的笛卡尔网格方法.整个流场采用有限体积方法计算,但是对于物面附近切割的网格计算域不完整,因此通过RBF插值来得到物体内部实单元上的值,从而使得流场计算顺利的进行推进.相比于其他插值方法来说,RBF插值特别适合于多维离散点的插值,不会因为插值点的位置关系而产生病态的系数矩阵,而且使用比较少的网格点就可以得到同最小二乘方法相当的精度.将这种方法应用于圆柱和翼型绕流中,得到了满意的结果,证明了方法的有效性,还对边界条件的各种不同形式对结果的影响作了对比.