应用不同湍流模型计算叶栅流场的结果与比较

采用有限体积TVD-LW格式求解N-S方程,使用了B-L模型、高雷诺数k-ε模型和低雷诺数k-ε模型等多种湍流模型。并且用具体算例说明本文所用程序的可靠性,比较了3种的湍流模型的特点。B-L模型最大的优点是简单、使用方便。k-ε方程适应性强,分为低雷诺数k-ε方程和高雷诺数k-ε方程两种。高雷诺数k-ε方程比低雷诺数k-ε方程精度略低,但所需计算时间和存储量都比后者小,网格也可以粗些。k-ε模型对粘性层、分离区有比较高的分辨能力,高雷诺数k-ε模型比较适于一般数值计算。      

旋翼翼型非定常动态失速特性的CFD模拟及参数分析

构建了一套基于运动嵌套网格技术和可压缩RANS方程的旋翼翼型非定常流动特性模拟的高效、高精度的CFD方法。首先,发展了基于Poisson方程求解的围绕翼型的粘性贴体正交网格生成方法,并提出了基于最小距离法(MDM)改进策略的运动嵌套网格生成方法,克服了弹簧法可能导致网格畸变的不足;其次,为准确模拟由湍流分离和气流再附引起的气动力的迟滞效应,基于RANS方程、双时间方法和高阶插值格式,建立了旋翼翼型非定常气动特性分析的高精度数值方法,并采用能够较好捕捉气流分离现象的S-A湍流模型;再次,针对旋翼后行桨叶动态失速时桨叶剖面来流速度较低、迎角较大的特点,为解决低来流速度时L-B半经验模型在旋翼翼型非定常动态失速计算中的局限性,并克服可压缩方程对低速流场计算收敛困难和精度低的问题,建立了基于Pletcher-Chen低速预处理方法、FAS多重网格法和隐式LU-SGS方法相结合的高效数值方法。应用发展的方法,分别针对NACA0012、SC1095旋翼翼型静态和轻度、深度动态失速进行计算,精确捕捉了气动力迟滞效应以及翼型前缘脱体涡的产生、对流和脱落过程,验证了本文方法的有效性;最后,着重针对NACA0012动态失速状态,开展了振荡参数对旋翼翼型非定常动态失速特性影响的分析,研究结果表明翼型迎角平均值、振幅及减缩频率的变化均能引起迟滞效应的改变并使得气动力峰值发生有规律的前、后移现象等。     

纵列式直升机双旋翼气动特性分析

由于纵列式直升机前后旋翼之间会产生复杂的干扰现象,本文针对两旋翼之间的气动干扰问题进行了流场分析。首先,建立了基于动量源模型和N-S（Navier-Stokes）方程的数值研究方法,选用k?ω SST两方程湍流模型,采用基于压力隐式求解器,并对所建立求解方法进行了算例验证。随后,对悬停状态,前向来流状态,后向来流状态的纵列式双旋翼直升机前后旋翼干扰流场分别进行了数值模拟,并结合计算结果与流场特性,对干扰流场所产生的气动现象进行了分析。结果表明,针对本文算例,悬停状态前后旋翼性能均降低,后旋翼对前旋翼的性能影响较大。前向来流状态后旋翼升力损失大于悬停状态,后向来流状态前旋翼升力损失更为明显,最低处旋翼效率损失近半。     

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.      

滑移流区内不同Kn数对流动及传热影响的数值分析

针对微尺度二维通道气体滑移流区的流动阻力和传热特性,采用滑移模型和SIMPLER数值进行了模拟分析.研究了定热流情况下不同Kn数时通道内压力、速度和温度场的分布以及流动阻力和换热的特性.数值分析结果表明随着Kn数的增大,速度滑移和温度跳跃的程度都增加,并且影响着流动阻力的下降,但两者对于换热系数的影响趋势是相反的:速度滑移有利于提高换热系数,而温度跳跃则降低换热系数,但总体上增加了微通道中的热阻.     

直升机典型城区环境流场计算与试验

首先,针对典型城区模型,在黏性非结构直角网格上求解了三维Navier Stokes（N-S）方程,并通过数值模拟计算了典型城区环境流场。其次,采用粒子图像测速（Particle image velocimetry,PIV）技术开展了城区流场测量试验研究。最后,通过对比试验与计算结果分析了典型城区流场特性。     

民机翼梢小翼气动特性数值模拟分析研究

针对某民机翼身组合体加装不同翼梢小翼——融合式翼梢小翼和翼尖涡扩散器进行对比分析。用软件ICEM生成原始翼身组合体及加装不同小翼的翼身组合体的点对点对接多块网格技术生成高质量的数值计算网格,运用Roe三阶迎风偏置通量差分裂方法和隐式近似因子分解方法求解雷诺平均N-S方程。对两种不同翼梢小翼在巡航点进行数值模拟,得到合理的翼稍小翼几何参数;对比升阻特性数据,两种翼梢小翼都可以提高升力系数,减少阻力系数。相对原始翼身组合体,融合式翼梢小翼升阻比增加9%,翼尖涡扩散器升阻比增加6%;但翼尖涡扩散器在力矩特性上有较大优势。所得结论对民用飞机翼稍小翼设计工作具有较大的工程应用价值。     

基于 N-S 方程的飞机风挡雨刷气动力计算

风挡雨刷在飞行过程中易出现被气流吹离复位位置、飘起现象,阻碍飞行员视野,严重影响飞行安全。选取某型机风挡雨刷作为优化对象,采用了结构、非结构混合网格的办法,在机身壁面附近生成结构化网格,在雨刷机构表面生成非结构化网格,用 Interface 面将两部分结合生成雨刷最终计算模型。利用 Fluent 软件,基于 N-S 方程和 S-A 湍流模型,采用时间平均法对不可压缩流的参数进行时均化,计算了不同飞行状态下,风挡雨刷在复位位置以不同角度停放及增加扰流板时所受气动力,对风挡雨刷的外形及停放角度进行优化,比较分析了优化前后雨刷所受气动力。仿真计算结果表明,雨刷合理的停放角度、在刷臂上加装扰流板能有效引导风挡雨刷迎风面气流的方向,从而减小雨刷迎风面所受气动力。通过试飞验证了数值计算的可靠性及优化方案的可行性,为风挡雨刷预紧力的设置及复位位置的确定提供依据。     

多面体网格在CFD中的应用

介绍了多面体网格技术在计算流体力学(CFD)中的应用.考虑到粘性计算的需要,将混合网格技术引入多面体网格中,发展了一套多面体网格生成及求解技术,并用于M6机翼、DLR-F6复杂外形的气动特性模拟中.计算结果表明,多面体网格技术具有准确度较高、效率较高及适用性广等优点.     

非定常N－S方程的线化及其差分求解

在微幅简谐振荡的假设下，对计算平面内非定常Ｎ－Ｓ方程进行了线比，并在ＭａｃＣｏｒｍａｒｋ于１９８０年提出的半隐格式的基础上，发展了一种适用于求解计算平面内线化Ｎ－Ｓ方程的差分格式．用该差分格式求解计算平面内的线化Ｎ－Ｓ方程时，可以避免求解五对角方程；且时间步长可以取得较大，计算效率高．