基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

倾转四旋翼UAV过渡状态下旋翼/机身气动干扰数值分析

通过数值模拟的方法对十字构型倾转四旋翼UAV在倾转过渡状态下前旋翼/机身/后旋翼间的气动干扰进行了计算与分析。分别建立了孤立前旋翼、前后旋翼、前旋翼-机身、前后旋翼-机身组合时的计算模型。在不同旋翼倾转角下,研究了以上组合模型对前旋翼、后旋翼以及机身的气动力和力矩的影响并分析了相互间的气动干扰情况,得出了前旋翼是整机气动干扰影响的主导因素,揭示了旋翼/机身间气动干扰形成的原因,为该构型旋翼布置与控制系统的设计提供参考。     

纵列式直升机双旋翼流场及性能试验

通过加装一台套单旋翼试验台,与原有的单旋翼试验台组成了双旋翼试验台.基于测力天平和PIV技术,通过改变前后两旋翼的水平和轴向间距,测量了纵列式双旋翼不同气动布局干扰状态下的流场特性和旋翼性能,并与单旋翼情况进行了比较.分别给出了悬停和前飞状态下的旋翼速度场和涡量场分布,对比了不同纵向和轴向间距对双旋翼性能的影响.结果表明:悬停时,后旋翼性能和单旋翼基本一致;前飞时,后旋翼性能比单旋翼差,且随前进比的增加,后旋翼性能与单旋翼性能的差距会更大.     

基于目标压力分布的旋翼先进气动外形反设计分析方法

基于雷诺平均Navier-Stokes（RANS）方程、运动嵌套网格、目标压力分布及余量修正方法,构建了一套直升机旋翼桨叶先进气动外形反设计方法。为避免由桨叶气动外形变化导致的网格畸变,发展了一套基于Poisson方程求解的旋翼桨叶结构化贴体正交网格的快速、自动化生成方法,为提高运动嵌套网格的生成质量和通用性,采用剖面间网格插值与桨尖翻折相结合的方法,同时建立了基于“Top Map”和“Inverse Map”相结合的洞边界划定和贡献单元搜寻的新方法。基于Navier-Stokes方程和双时间法建立了旋翼非定常流场模拟方法,通量求解采用Roe-MUSCL格式,并使用低速预处理法来克服前飞旋翼流场收敛中遇到的刚性问题。在计算流体力学（CFD）方法基础上,基于旋翼翼型压力系数余量联立各方位角处的反设计MGM（Modified Garabedia-McFadden）超定方程组,并依据激波分离、失速等约束设置了各方位角处的反设计权重系数,创建了基于MGM超定方程组最小二乘解的旋翼气动外形（翼型）设计方法。应用所建立的方法,分别针对多目标、多状态和前飞时的旋翼（翼型）气动外形进行反设计分析,验证了本文方法的有效性。最后,将该方法拓展应用到旋翼桨尖气动外形设计中,设计得到与UH-60A直升机旋翼气动特性相似的矩形桨叶外形。     

共轴旋翼高速直升机雷达散射截面计算及雷达吸波材料影响分析

为准确高效地预估共轴旋翼高速直升机的雷达散射特性,结合雷达吸波材料(Radar absorbing material,RAM)在隐身设计中的应用,开展了共轴旋翼高速直升机雷达散射截面(Radar cross section,RCS)特性及涂覆型RAM对其影响的研究。首先,基于计算涂覆目标表面散射的物理光学法(Physical optics,PO)和计算涂覆边缘绕射的等效电磁流法(Method of equivalent current,MEC),建立了计算RCS的高频方法,并通过涂覆了RAM的金属球和直升机矩形桨叶算例验证了其有效性。在此基础上,研究双旋翼、尾部螺旋桨、平垂尾和机身在鼻锥、侧向和尾追3个典型方位的雷达散射特性和强散射源分布,并采用局部涂覆RAM的方法进行隐身设计。研究表明:尾部螺旋桨、共轴旋翼桨毂及其整流罩部位、机身上曲率较大的鼻锥和尾部以及曲率较小的侧面护板是机身的重要强散射部位。在强散射部位涂覆RAM能有效降低高速直升机各方位双站RCS的均峰值,显著提升高速直升机隐身性能的效果。     

变来流速度下旋翼翼型非定常气动特性分析

为分析变来流速度状态下的旋翼翼型气动特性,提出了利用翼型平移来模拟来流速度变化的数值方法.在此方法基础上,采用基于隐式LU-SGS(lower upper symmetric Gauss-Seidal)方法的非定常雷诺平均N-S(Navier-Stokes)(RANS)方程,模拟了SC1095旋翼翼型在定迎角 变来流速度及变迎角 变来流速度状态下的非定常气动特性.通过对比分析发现:翼型在变速度-定迎角状态下会表现出明显的非定常现象,产生了前缘分离涡,气动特性会出现明显的迟滞效应及波动现象,脉动速度越大,非定常效果越明显.并且基准速度越大,翼型气动特性的峰值越大;翼型迎角越大,非定常涡出现的也越早.考虑直升机旋翼翼型实际工作环境,在变速度-动态失速状态下,翼型最大迎角处的气动力会得到一定程度的削弱,在小迎角下的气动力得到一定程度的增强,且脉动速度越大,翼型的非定常特性也越强.      

倾转旋翼机多部件对机翼气动干扰的分析及优化

为研究倾转旋翼机各部件对机翼的气动干扰效应,建立了一套适用于倾转旋翼机流场CFD求解方法,提出并生成了一套由三棱柱/四面体组成的适用于倾转旋翼机多种飞行状态的非结构混合网格系统.以三维Navier-Stokes方程为主控方程,采用动量源方法进行倾转旋翼的模拟,并引入了高效的OpenMP并行加速技术等,提出了一套新型的动量源项添加及搜索方法.计算分析了倾转旋翼机旋翼/机身/短舱对于机翼气动特性的干扰影响,得出了一些对设计有指导意义的结论.采用基于径向基函数（RBF）的代理模型方法,根据机翼不同展向位置的干扰程度对各段翼型配置进行气动优化.优化结果表明:考虑气动干扰作用,在巡航速度下优化后的全机升阻比增长达到了36.78%.      

使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟

为减少Jameson二阶中心差分有限体积法导致的旋翼尾迹的数值耗散,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,建立了一个三维雷诺平均N-S方程数值模拟悬停状态旋翼流场的方法。为了充分考虑旋翼尾迹对流场的影响,采用周期性边界条件和由动量理论导出的远场边界条件。为进一步减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。然后,进行了算例计算,给出了桨叶表面的压强分布,与可得到的试验数据及二阶中心差分方法的计算结果进行了对比,并针对几种不同桨尖形状的旋翼悬停流场进行计算,数值结果显示:后掠桨尖可减弱超临界流动。此外,还计算和分析了旋翼下方不同位置上的涡量分布。计算结果表明,本文方法能够有效地减少旋翼尾迹的数值耗散。      

纵列式双旋翼悬停状态气动干扰特性参数影响分析

建立了一个纵列式直升机双旋翼气动干扰特性分析的自由尾迹迭代方法。在该方法中,考虑了双旋翼以及旋翼与尾迹之间的相互干扰影响,将旋翼尾迹模型、桨叶气动力模型以及旋翼配平模型进行耦合求解。同时给出了一个适合于双旋翼干扰计算特点的配平方法。计算了悬停时干扰状态下的双旋翼诱导速度分布以及旋翼性能,并与试验数据进行对比,验证了方法的有效性。应用上述方法,对比分析了纵列式双旋翼与单旋翼的性能,结果表明:悬停时纵列式前、后旋翼的性能都比单旋翼时的要差。文中进一步系统地分析了悬停状态双旋翼纵向间距和轴向间距对气动干扰特性的影响,结果显示:当纵向间距为1.3R时,纵列式双旋翼相比两单独旋翼需要额外的附加功率为8.5%,且随纵向间距的增大,附加功率减小并出现负值,当纵向间距为1.85R时,附加功率最小。     

基于混合CFD方法的旋翼桨-涡干扰噪声计算

建立了一个基于Euler方程和自由尾迹分析的新的、高效的混合计算方法,用于旋翼桨-涡干扰(BVI)气动和噪声特性的研究.在该方法中,仅针对单片桨叶进行计算,通过求解Euler方程来捕捉流场中旋翼近尾迹的影响以及桨叶表面的非定常载荷,流场空间中旋翼远尾迹的位置由自由尾迹模型给出;提出了广义网格速度方法以将旋翼远尾迹的影响...