一种适合于旋翼前飞非定常流场计算的新型运动嵌套网格方法

针对旋翼流场计算中关键的运动嵌套网格技术,提出了一个新的网格单元关系快速判断方法。该方法通过计算特殊网格点在背景网格上的位置来建立桨叶网格在背景网格中所占据范围的关系,并用于嵌套网格洞单元的识别,又进一步将该方法推广用于Inverse Map的建立,在此基础上,建立了一套用于前飞流场计算的运动嵌套方法。通过多种前飞旋翼测试算例,分别评估了背景网格、桨叶网格尺寸改变对计算时间的影响,同时验证了该方法的可靠性,结果表明构建网格嵌套关系所用时间随桨叶网格和背景网格尺寸的变化幅度较小,且具有良好的鲁棒性,能够满足旋翼前飞状态计算的需要。最后,在基于前飞非定常N-S/Euler流场求解器中,采用该方法模拟了UH-60A和7A旋翼悬停/前飞状态的气动特性,通过计算时间的对比体现出该方法的高效性,并且流场计算结果与试验数据对比表明该方法可以有效地用于旋翼前飞非定常流场的计算。     

非结构嵌套网格的直升机旋翼/机身前飞流场数值模拟

建立了一个适用于前飞状态直升机单独旋翼、机身以及旋翼/机身组合体流场的数值计算方法及模型.在该方法中,使用非结构运动嵌套网格描述桨叶之间以及桨叶与机身之间的空间位置变化关系,控制方程采用惯性坐标系下的非定常Navier-Stokes(N-S)主控方程,空间方向上使用二阶迎风格式,时间方向上使用Lower-upper symmetric Gauss-Seidel(LU-SGS)格式.对于求解中遇到的旋翼/机身间距离过近所导致的贡献单元生成困难问题提出了一种新的解决方法.应用所建立的方法,分别对单独旋翼、机身及旋翼机身组合体的前飞绕流场进行了数值模拟与分析,计算结果与试验结果的对比表明本文建立的方法对直升机的前飞绕流场的计算和分析是很有效的.      

基于非结构嵌套网格和逆风格式的旋翼悬停流场数值模拟

为更好地计及桨叶间相互影响,同时捕捉桨尖附近的流场细节,发展了一套非结构运动嵌套网格方法,并针对关键的贡献单元搜寻问题提出了一种高效的搜索方法.为了可方便地扩展到前飞状态的计算,同时为提高旋翼流场涡量的捕捉精度和效率,采用双时间法求解惯性坐标系下的非定常N-S主控方程,空间方向上使用二阶迎风格式,在此基础上,将不同精度格式的计算结果与试验值进行对比,并讨论了桨尖形状变化对旋翼气动特性的影响,结果表明本文建立的CFD方法对旋翼悬停流场的模拟是很有效的.     

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。     

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器,用来精确模拟复杂的旋翼近地流场,为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中,控制方程采用惯性坐标系下的非定常N-S方程,空间方向上采用二阶迎风格式,并用背景网格的一个面模拟地面作用,以方便地面边界条件的处理。应用所建立的模型,首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟,以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益,并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上,得出了一些有意义的结论。     

基于非结构网格CFD技术的旋翼气动噪声计算方法研究

将基于非结构网格技术的旋翼流场CFD计算方法与基于FW-H和Kirchhoff方程的声学方法相结合,建立了一套既适合于直升机旋翼厚度、载荷和桨-涡干扰噪声,又适合于跨声速高速脉冲噪声的综合计算模型。为提高旋翼流场及桨叶表面气动载荷计算的精度,主控方程的求解采用了三维可压非定常的N-S方程,网格划分则使用非结构运动嵌套网格方法。在噪声计算中,通过FW-H方法计算旋翼的厚度噪声和载荷噪声,并选取能够包含流场非线性区域的旋转面作为Kirchhoff积分面,由Kirchhoff方法计算包含四极子项的高速脉冲噪声。应用该模型,以AH-1旋翼为算例,计算了不同飞行状态下的旋翼气动噪声,并与可得到的试验结果进行比较,验证了方法的有效性。然后,着重对两种声学方法对计算结果的影响进行了对比研究,并分析了旋翼厚度噪声、载荷噪声和四极子噪声的特性。     

共轴式双旋翼悬停流场和气动力的CFD计算

考虑到共轴式双旋翼流场特征高度复杂、桨叶承受非定常气动载荷的特点,为更好地预测共轴式双旋翼的气动特性,把非结构嵌套网格方法和网格的自适应技术相结合,发展了一套适合于共轴式双旋翼流场数值计算的求解器.在该求解器中,采用非结构嵌套网格方案来模拟桨叶之间存在的相对运动,自适应网格技术用来捕捉尾迹对流场和气动特性的影响,求解惯性坐标系下的非定常N-S主控方程来模拟流场的非定常特性.应用该求解器,首先计算了有试验结果可供对比的一副试验旋翼的诱导速度场分布,在此基础上,计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性,并与单旋翼的计算结果进行了对比和分析,得到了一些有意义的结论.     

基于非结构嵌套网格的涵道螺旋桨数值模拟

基于非结构嵌套网格方法,通过求解可压缩Euler方程对涵道螺旋桨系统的非定常流场进行了数值模拟,并研究了桨尖间隙、螺旋桨位置、涵道尾扩张角等参数的影响规律.将计算区域分为相互重叠的旋转区域和静止区域,旋转区域包含螺旋桨和旋转的桨毂,静止区域包含涵道和部分静止的桨毂,采用嵌套网格算法真实地模拟了螺旋桨和涵道之间的相对旋转运动,得到了细致的螺旋桨流场.计算值和实验测量值的对比结果相一致.参数影响研究表明,减小桨尖间隙、在一定范围内后移螺旋桨或增加尾扩张角,均可以提高涵道螺旋桨的气动性能.     

Euler方程的自适应多重非结构网格计算

介绍一种基于四面体网格的多重网格技术,用于求解三维Euler方程定常解,具有十分明显的加速收敛效果。求解Euler方程的基本算法采用有限体积中心差分格点格式。多重网格技术采用一组几何上相互独立的粗细不同的网格。流动变量、余量和修正量在相邻粗细网格间的相互传递通过线性插值实现。采用自适应技术对网格进行加密,可以获得较优的网格,并作为一套细网格,进行多重网格计算。数值实验表明本文自适应多重网格法十分有效。     

三维非结构自适应多重网格技术

发展了一种基于几何外形的非结构自适应多重网格技术。根据流场参数的变化梯度确定加密边,由加密准则对粗网格进行自适应剖分得到分布合理的较密网格。通过预先生成的初始极密表面网格来将边界网格的加密点投影到边界上,使得边界保持初始外形。为了提高流场求解速度,将自适应加密后的一系列网格作为多重网格,采用了多重网格算法进行计算。M6机翼的流场计算表明该方法可以准确捕捉激波位置,并且表面的压力系数分布与实验值相吻合。     

高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

共轴刚性旋翼直升机着舰流场计算分析

以CFD软件为基础,建立了适用于共轴刚性(双)旋翼直升机着舰飞行时的气动计算模型.应用这一模型,首先计算了Robin旋翼的悬停状态算例,计算结果与试验数据进行了对比.然后着重分析了共轴刚性旋翼直升机在着舰与着陆时的流场异同,并进一步分析了着舰飞行时不同风向角对直升机俯仰、滚转和偏航力矩的影响,获得了一些对直升机着舰飞行有指导意义的结果.     

Experiments on unsteady vortex flowfield of typical rotor airfoils under dynamic stall conditions

A new experiment for airfoil dynamic stall is conducted by employing the advanced particle image velocimetry(PIV) technology in an open-return wind tunnel. The aim of this experimental investigation is to demonstrate the influences of different motion parameters on the convection velocity, position and strength of leading edge vortex(LEV) of airfoil under different dynamic stall conditions. Two different typical rotor airfoils, OA209 and SC1095, are measured at different free stream velocities, oscillation frequencies, and angles of attack. It is demonstrated by the measured data that the airfoil with larger leading edge radius could notably decrease the strength of LEV. The angle of attack(Ao A) of airfoil can obviously influence the dynamic stall characteristics of airfoil,and the LEV would be effectively inhibited by decreasing the mean pitch angle. In addition, the convection velocity of LEV is estimated in this measurement, and the results demonstrate that the influence of airfoil shape on convection velocity of LEV is limited, but the convection velocity of LEV would be increased by enlarging the oscillation frequency. Meanwhile, the convection velocity of LEV is a time variant value, and this value would increase as the LEV convects to the trailing edge of airfoil.     

基于图像处理思想的自适应网格技术

为了提高激波和旋涡等复杂流动现象的捕捉精度和计算效率,基于图像处理的思想,研究了非结构网格的自适应加密和稀疏算法.将流场计算网格与图像像素类比,借鉴图像处理的方法,提出了一种流场自适应区域的分割算法.该方法能自动对自适应区进行分层,中间区域网格可以达到原来多次自适应的效果,而且提高了生成网格的质量.将算法嵌入到流场求解器中,根据流场求解结果对网格进行自适应,并使用算例进行了验证.     

旋翼干扰流场数值模拟中的新型Blend低速预处理方法

针对定常低速预处理在双时间步长内收敛性差、旋翼干扰流场数值模拟中非定常低速预处理自动退化为无预处理的问题,在非定常低速预处理方法的基础上,采用量级分析方法研究了对流项和黏性项量级,通过组合定常和非定常预处理特征值矩阵,建立了一种新型Blend低速预处理方法。通过圆柱非定常扰流、前飞直升机旋翼/机身干扰非定常流场和直升机/舰船干扰非定常流场的数值模拟,对各低速预处理方法的收敛性和预测精度进行了细致研究。结果表明,本文方法解决了在旋翼干扰流场数值模拟时非定常低速预处理自动退化为无预处理的问题,相比定常低速预处理大幅提高了收敛效率。     

用于旋翼气动特性计算的一种时间推进自由涡新方法

建立了一个新的用于旋翼尾迹和诱导速度分析的时间推进自由涡计算方法.在该方法中,推导了一个新的二步二阶预测-校正差分算法--"D2PC"算法,并对采用该算法离散表示的涡线主控方程进行修正,以提高尾迹的求解精度;基于Vatistas公式给出涡核模型,并考虑了实际的耗散影响;采用更适合于旋翼气动特性计算的简化升力面桨叶气动模型;同时在自由涡方法中结合了挥舞动力学模型和旋翼配平模型.通过计算的旋翼尾迹几何形状和下洗速度与实验值的对比,验证了本文建立的自由涡方法的有效性.然后利用该方法,计算和分析了在悬停和前飞状态旋翼的尾迹几何形状和下洗诱导速度.     

适于旋翼CFD模拟的高效预定边界运动嵌套网格方法

通过引入新的预定边界嵌套策略,并结合提出的"逆向边界"和"Local Direct-Map"(LDP)技术,建立了一套预定边界嵌套网格方法。"逆向边界"通过预估计算对边界进行动态调整克服了透视图方法很难明确避开流动非线性区的缺点;LDP方法解决了传统的基于Inverse-map透视图嵌套方法中分辨率与计算效率矛盾的问题。对不同操纵特性下相同嵌套网格的洞边界分布特性进行检验,结果表明预定边界嵌套方法保持了高鲁棒性并且嵌套效率提高了16.7倍。为进一步验证对非线性流场模拟的有效性以及对旋翼复杂运动的适应性,建立了一套适用于旋翼非定常流场计算的可压雷诺平均N-S方程数值模拟方法,分别对悬停和前飞状态下的C-T旋翼和UH-60A旋翼进行了数值分析。模拟结果与试验值吻合较好,表明了预定边界嵌套网格方法能够有效地用于旋翼非定常流场和气动特性的数值模拟分析。     

基于伴随方程方法的非结构网格自适应技术及应用

对基于伴随方程的网格自适应技术的基本原理进行了研究,构造并实现了基于流场伴随方程的非结构网格自适应探测器,利用局部加密网格的办法,建立了基于Euler方程的无粘伴随自适应能力。对于粘性流动问题,提出首先应用无粘计算和伴随自适应技术从基准网格自动生成具有较好网格分辨率的非结构自适应网格,然后采用层推进法自动生成粘性计算网格的求解策略。采用RAE2822二维翼型和ONERA-M6三维机翼对所建立的无粘伴随自适应方法进行了验证,并应用所提出的方法和策略,对第二届国际涡流试验项目(VFE-2)中的65°尖前缘三角翼大攻角涡流场进行了数值模拟。对比计算结果表明,所建立的非结构网格伴随自适应方法是可行的和有效的。