动态负载平衡的二维非结构网格DSMC并行算法研究

基于PC-CLUSTER群机并行体系结构与消息传递库MPI并行环境,研究了二维非结构网格DSMC并行算法。提出一类非结构网格动态分区策略,保证各子区域的分子数量大致相等,实现计算进程间的动态负载平衡。利用MPI库函数构造了符合DSMC并行原理的单步通讯法。采用单控制多数据流(SPMD)以及Master/Slave并行模式,设计了二维非结构网格DSMC整体并行算法。在程序的编制过程中,充分展现了Fortran90高级语言的主要特性,引人动态数组、指针、链表及派生类型数据。最后对过渡流域高超声速绕流进行了并行计算,数值试验的结果表明本文设计的并行算法可以取得较为理想的加速比。     

低雷诺数下跨声速压气机转子失速工况时流动失稳触发过程研究

采用非定常数值方法模拟了低雷诺数条件下NASA Rotor 37跨声速压气机转子内部流动失稳机制.结果表明,随着压气机失速工况的推进,附面层径向涡不断扩大和增强,由附面层径向涡引发的叶顶附面层分离阻塞区不断向叶片前缘移动,直至与叶顶压力面前缘附近由激波和间隙泄漏流诱发的阻塞区相结合,使叶顶通道来流完全被阻塞,最终触发压气机流动失稳.      

MPI＋OpenMP算法在三维可压缩流场计算中的应用

在多核CPU集群并行体系结构下,采用MPI＋OpenMP的混合并行算法,对可压缩流场进行数值模拟,并在计算时间上与MPI算法进行比较。流场计算的控制方程为Euler方程,空间离散采用Jameson中心差分格式,时间离散采用R-K法,并行模式采用Master/Slave模式。通过对M6机翼和某弹丸采用多种并行方法进行流场的数值模拟,得出MPI＋OpenMP混合算法在一定条件下具有高效性的结论。     

基于Van Leer+AUSM混合格式超声速流场的并行数值算法研究

基于CLUSTER集群并行体系结构与消息传递库MPI并行环境,开展了超声速流场的并行数值算 法研究。计算网格采用多块对接粘性结构网格,控制方程为积分形式的Navier\|Stoke s(N-S)方程,空间离散采用基于Van Leer+AUSM混合格式的有限体积法,粘性项采用中心 格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,应用一方程S -A模型计算湍流的粘性。并行计算采用 分区并行的思想,设计了满足负载平衡要求的区域分解和通讯方法,编程上采用单控制流多 数据流(SPMD)模型以及Master／Slave并行模式。最后,通过对某弹丸和某制导火箭弹超声 速流动的数值模拟, 验证了并行算法的准确性、高效性和处理复杂几何外形的能力。
     

风扇前缘曲线相对前掠对风扇效率的影响

根据风扇前缘曲线相对前掠概念，将其应用于NASA67风扇叶片的改型中．并通过数值模拟研究了三种不同前缘曲线相对前掠程度对叶片通道内气流流动的影响。计算结果表明：前掠叶片改善了风扇的气动特性，降低了流动损失．进而提高了风扇效率。     

介质阻挡放电等离子体气动激励的动量特性

从诱导气流速度和体积力两个方面开展了介质阻挡放电等离子体气动激励的动量特性研究,对不同激励参数下的诱导气流速度进行了测量与结果分析;通过对体积力的实验测量与理论计算,对比研究了体积力的变化规律与主要影响因素.结果表明:等离子体气动激励可以增大气流速度,但随着气流速度的增大,等离子体气动激励的加速效果减弱;激励电压或激励频率增大,体积力均表现为线性增加,但激励电压增大时可以更好的增大体积力;与增大电极内间距相比,增大下层电极宽度可以产生更大的体积力.      

近空间高超声速流场通量分裂型DSMC-IP方法

针对近空间高超声速流场的特点,采用Van Leer格式对直接模拟蒙特卡洛-信息保存(DSMC-IP)方法质量守恒方程中的计算格式进行改进.将局部马赫数作为分裂通量的标准,并重构单元分界面两侧的左右输运通量,使得计算格式具有通量分裂的特点,解决了DSMC-IP方法在高超声速流场计算中的应用问题.将改进后的通量分裂型DSMC-IP方法引入非结构网格中,对二维近空间高超声速流场进行数值模拟.计算结果表明:通量分裂型DSMC-IP方法所得出的数值结果与实验值及参考值符合较好,明显降低了直接模拟蒙特卡洛(DSMC)方法所带来的统计耗散.当来流气体的稀薄程度增加时,其非平衡效应也更加明显,而通量分裂型DSMC-IP方法的计算结果与参考值相差均在10%以内,良好地反映了非平衡条件下的流场特征,验证了通量分裂型DSMC-IP方法的可行性和有效性.      

某中央翼燃油箱惰化流场的数值模拟及特性分析

在不同富氮气体入口参数下,对某大型运输机中央翼燃油箱的惰化过程进行了数值模拟,研究了在不同流量和体积分数富氮气体条件下燃油箱的惰化规律.计算结果显示:当富氮气体中氮气体积分数为95%、体积流量分别为0.021,0.028m3/s和0.042m3/s时,燃油箱内平均氧体积分数随体积置换次数的变化曲线几乎完全重合;当富氮气体体积流量为0.042m3/s、氮气体积百分比分别为92%,94%和98%时,燃油箱惰化率随时间的变化曲线重合.说明相同体积分数富氮气体条件下燃油箱惰化完成时所需富氮气体量与富氮气体流量无关,而富氮气体的流量直接决定了各舱室富氮气体的惰化率.通过对计算结果的分析,证明了利用缩比油箱去模拟真实油箱的惰化过程的合理性.为优化某运输机中央翼燃油箱惰化系统设计提供参考.      

高负荷压气机叶栅等离子体流动控制数值仿真与拓扑分析

为揭示叶栅等离子体流动控制的影响规律与作用机理,对等离子体气动激励前后高负荷压气机叶栅内部流动和拓扑结构进行了对比研究。结果表明:等离子体气动激励抑制叶栅流动分离的作用效果最明显的区域位于总压损失区域与主流区域的边界上;不同等离子体气动激励布局,对固壁面拓扑结构以及奇点总数的影响规律不同;吸力面流向激励通过增强附面层流体抵抗逆压梯度的能力,可以改善叶栅中间叶高流动特性;端壁横向激励通过抑制横向流动,抑制角区流动分离能力较强,并改变叶片展向的负荷分布;组合激励结合了吸力面流向激励和端壁横向激励的作用优势,因而提高叶栅气动性能、降低流动损失的效果最好。      

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

研究了二维热化学非平衡流动非结构网格DSMC方法实现的过程。提出了一种新型的高效搜索算法,该算法不仅可以跟踪模拟分子在网格之间的迁移,而且在搜索过程中可以准确判别分子与物面是否相互作用,避免了原有算法中分子表面反射非确定论判据。为加快流场的时间发展历程,设计了适合DSMC方法的动态局部时间步长技术,使其可以应用到定常和非定常流场的计算。利用Fortran90的动态分配内存技术编制了计算程序。最后对过渡流域高超声速圆柱绕流进行了数值试验,计算结果初步验证了该算法的可行性。