风力机标模非定常数值模拟中的影响因素研究

基于自主研发的“亚跨超 CFD 软件平台”(TRIP3.0)，采用刚性运动网格技术和动态拼接网格技术，开发了针对旋转类机械的非定常求解模块。本文开展了 NREL Phase VI 风力机标模非定常数值模拟中的影响因素研究，影响因素主要包括不同时间步长的影响、不同湍流模型的影响、刚性动网格技术和动态拼接网格技术的影响三个方面。本文数值模拟采用的控制方程为雷诺平均 N-S 方程，采用有限体积法离散控制方程，离散方程的时间方向采用“双时间步”方法求解，空间方向无粘项离散采用 MUSCL-Roe 格式，湍流模型采用 SA 和 SST 湍流模型，并引入多重网格和并行计算技术加速求解。数值模拟结果表明：时间步长、湍流模型、网格方案等因素主要影响风力机叶片吸力面的流动结构，进而影响吸力面的压力分布，而对压力面的流动结构和压力分布基本没有影响；采用刚性运动网格结合 SA 湍流模型所得到的压力分布更接近实验值。     

CFD技术在带动力飞机气动设计中的应用

应用数值模拟的方法,研究了飞机带动力发动机的流场模拟,基于此计算分析了其性能特性.采用基于雷诺时均的N-S方程,动力方程中的雷诺应力通过SST二方程湍流模型和自动壁面处理计算;同时,采用守恒隐式有限体积法.离散方程使用有限高精度对流程序类似于Barth和Jesperson进行计算.质量流动通过Rhie和Chow算法压力-速度耦合进行计算.     

超音速气流绕局部透气凹角流动的数值模拟

以具有压力分裂形式的简化Ｎ Ｓ方程为控制方程,数值模拟了超音速来流条件下的激波 边界层干扰被动控制（ｐａｓｓｉｖｅｃｏｎｔｒｏｌｏｆｓｈｏｃｋ ｂｏｕｎｄａｒｙｌａｙｅｒｉｎｔｅｒａｃｔｉｏｎ）。模拟是以预先给定激波前吹气和激波后吸气的流量来实现的。为了定性地确定吹气或吸气对激波 边界层干扰的影响,首先计算了单独吹气和单独吸气两种情况。数值计算时采用了多重扫描法对控制方程差分离散,以反映亚音速区压力对流场的椭圆性影响。     

用多重网格法模拟三维流动

采用多重网格法模拟三维流动。控制方程在三维直角或圆柱坐标系中通过有限体积方法进行离散 ,SIMPLE算法对速度与压力进行耦合 ,亚松弛策略提高离散方程的收敛速度。在压力方程的求解过程中 ,界面速度应用由Rhie、Chow和Majumdar等人提出的动量方程插值法来剔除不合理的压力波动。三维流动算例表明采用多重网格法模拟三维流动可以得到良好的计算精度和收敛性。     

Sajben跨声速扩压器分离流动中湍流模式数值研究

为研究湍流模式对激波/湍流边界层干扰内流流动的影响,提高数值计算准确度,使用SA,SST k-ω,非线性EASM k-ω,Gao-Yong四个湍流模式对Sajben扩压器内激波/湍流边界层干扰流动进行了数值计算。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。文中展示了四个湍流模式计算得到的壁面压力、速度剖面、摩阻系数等分布。计算值与实验值符合很好,四个湍流模式总体上能够较好地模拟扩压器内激波/湍流边界层干扰复杂分离流动。Gao-Yong湍流模式对分离区内的压力、速度型的模拟更加准确,而非线性EASM k-ω模式对分离再附点位置计算最理想。     

基于粘性非结构网格的三维增升装置N-S方程数值模拟

应用基于粘性非结构网格求解N—S方程方法数值模拟了某型飞机三维增升装置的绕流。在物面附近粘性主导作用区域内，采用三棱柱型非结构网格，在其他流场区域，采用传统四面体非结构网格的混合非结构网格生成技术。湍流模型采用Spalart—Allmaras一方程模型，运用格心格式有限体积法，对控制方程和湍流模型方程进行数值离散。用四步龙格-库塔方法作显式时间推进，采用当地时间步长、隐式残值光顺等加速收敛措施．求解某型飞机全机带增升装置构型的粘性绕流流场，其数值计算与风洞试验结果显示了良好的吻合度。     

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

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

静止叶栅三维粘性数值模拟

应用Ｎ－ Ｓ方程压力修正算法和ｋ－ ε两方程湍流模型,对端弯叶片和基准（无端弯）叶片的平面叶栅进行了三维粘性数值模拟。在采用简单交错网格的情况下,通过对计算控制体通量的速度进行修正的方法解决了Ｎ－ Ｓ方程压力修法算中遇到的压力场和速度场呈锯齿状分布的问题。用实验测量值对本数值模拟方法的精度和可靠性进行了检验     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

导弹垂直发射系统三维超声速射流的数值模拟

对导弹垂直发射系统(VLS,Vertical launching system)内的三维湍流流场进行了数值模拟.采用了单元中心有限体积法和三阶迎风TVD差分格式离散Navier-Stokes方程,采用双时间步和预处理过程提高计算收敛速度,引入对角化的双时间ADI算法求解控制方程.仿真结果准确地模拟了发射系统内的激波形成过程,真实地反映出VLS中的复杂流动.研究表明,压力室底部射流中心的波后压强以大约1 200Hz的频率振荡.