共轴双桨涵道无人飞行器非定常气动特性

基于滑移网格技术和SST(shear stress transport) k -ω湍流模型,建立了模拟共轴双桨涵道无人飞行器(UAV)流场的CFD(计算流体动力学)数值方法,并通过计算风洞实验算例验证其有效性。数值模拟了共轴双桨涵道UAV在飞行过程中的动态流场,分析了涵道、飞行速度、螺旋桨转速、攻角等因素对其非定常气动特性的影响规律。计算结果表明:所建立的CFD数值计算方法适于模拟共轴双桨涵道UAV动态流场;涵道的存在显著削弱涵道螺旋桨的桨尖涡、后缘脱体涡和尾流收缩,具有较弱的桨-涡干扰和涡-涡干扰现象,明显减小共轴双桨涵道UAV的需用功率;随前飞速度增大,共轴双桨涵道UAV的升力和阻力同时增大;随螺旋桨转速增大,共轴双桨涵道UAV升力增大,而阻力减小;随攻角增大,共轴双桨涵道UAV的阻力增大,而升力先增大后减小。      

涵道螺旋桨与孤立螺旋桨气动特性的 数值模拟对比

基于非结构动态嵌套网格方法,对涵道螺旋桨与孤立螺旋桨的气动特性进行了非定常Euler方程数值模拟对比研究和分析.将计算区域分成旋转区与静止区两个子域,生成相互重叠的嵌套网格,有效地解决了螺旋桨与涵道以及静止部分桨毂之间的相对旋转问题.计算结果与实验结果相一致,验证了计算方法的正确性.涵道入口前缘形成的负压区是产生涵道附加拉力的原因所在.对比计算结果表明:与孤立螺旋桨相比,涵道的存在改善了螺旋桨桨尖区域的绕流特性,减小了桨尖损失;相同转速情况下,涵道螺旋桨产生更大的拉力,而所需功率略小,涵道螺旋桨系统具有更高的气动效率.      

一种涵道螺旋桨性能预测的CFD方法

为适应高性能涵道螺旋浆精细化设计需求,基于雷诺平均Navier—Stokes控制方程和非结构网格旋转／静止滑移面技术,开展了涵道与螺旋桨动／静部件干扰的非定常流及性能预测数值模拟方法研究。提出了．静止与旋转域的多块网格布局策略及转／静滑移面技术,研究了涵道壁而与螺旋桨桨梢间隙的网格布局形式,有效解决了螺旋桨旋转流及螺旋...     

悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。     

倾转涵道倾转过渡阶段的非定常气动力

为研究倾转涵道动力装置在倾转过渡阶段的非定常气动力,使用基于滑移网格技术的非定常计算方法,利用基于内场信息和叶素理论的压力盘模型模拟风扇螺旋桨,通过求解Navier-Stokes(N-S)方程,对涵道风扇俯仰拉起过程进行数值模拟.结果表明:在倾转过程中涵道风扇非定常气动力的迟滞特性明显,俯仰角速度延迟了流场分离,增加了涵道风扇的升力和阻力;低速飞行时,涵道风扇在整个倾转过渡阶段气动性能优良;在高速大迎角飞行时,涵道风扇气动性能恶化,俯仰力矩曲线紊乱,不利于进行倾转过渡飞行.      

涵道尾桨气动特性及翼型CFD分析

以Jameson/TVD混合格式对低速时的OAF和NACA63a312翼型的升阻特性进行了计算。该方法的特点是采用有限体积法进行离散,并用龙格-库塔时间推进法进行求解。通过将计算的结果与NASA公布的试验曲线进行对比,验证了该混合格式用于旋翼气动力计算的可行性。为了进一步研究旋翼的拉力和功率,本文建立了以微圆环动量理论和叶素理论为基础的计算模型,并用该模型计算了涵道尾桨在不同条件下的拉力和功率情况,计算的结果具有一定的代表性。      

直升机涵道尾桨与孤立尾桨的气动特性对比研究

本文推导了涵道尾桨的叶素理论公式，将文献［１］的滑流理论与叶素理论公式组合，建立了计算直升机悬停和侧风状态下涵道尾桨气动特性的方法。     

不同形状涵道风扇推进特性数值分析

以分布式涵道风扇推进为背景,使用多重参考系（MRF）和给定力分布的动量源方法（MSM）求解雷诺平均N-S（RANS）方程,对不同形状涵道风扇推进特性进行了数值模拟分析。对桨盘与不同形状涵道壁面之间的相互作用原理进行了研究,进而对不同外形分布式涵道风扇进行了分析。结果表明:纯圆形机匣推进特性最佳,纯方形最差,由方转圆的机匣推进特性居中;非圆机匣圆角的存在会诱使机匣内壁出现分离,产生干扰阻力,且圆角半径越小,影响越显著;非圆机匣影响风扇进口面积和桨尖涡的大小,从而从桨盘效率和唇口吸力两方面影响涵道推进效率。      

旋翼桨-涡干扰气动特性计算及参数影响研究

建立了一个适用于旋翼桨-涡干扰(BVI)气动特性分析的数值方法。在该方法中,控制方程采用惯性坐标系下的非定常Euler方程,以适合于分析BVI流场的特点;为便于前飞流场分区求解和信息传递,使用了嵌套网格方法;对于流场中涡线的模拟,建立了一种适用于有限体积格式的涡引入方法——广义网格速度法,以简化计算。应用上述方法对旋翼BVI流场进行了计算,并与可得到的试验数据进行对比,验证了方法的有效性。通过对比桨叶弦向位置压强的变化得出,在旋翼BVI过程中,气动载荷主要来自桨叶10%弦长内的前缘部分的压强突变。文中还进一步分析了涡强、干扰距离和干扰角度对BVI气动特性的影响。本文的BVI数值计算表明:当涡接近桨叶前缘时,升力达到最大;而涡位于后缘位置时,诱导速度改变桨叶环量分布,破坏了后缘的库塔条件,但随着涡的远离,桨叶表面环量开始重构。     

涵道风扇布局参数对其气动特性的影响研究

基于Navier-Stokes方程开展了不同气动布局参数下的涵道风扇非定常气动特性数值模拟研究。首先,结合滑移网格方法建立了适用于涵道风扇流场求解的数值模拟方法,并基于NASA涵道风扇试验模型开展了气动性能计算验证。随后,基于所建立的CFD方法开展了涵道风扇气动布局参数对其气动特性的影响研究,揭示了悬停状态下桨尖间隙和桨盘/涵道轴向相对位置对涵道风扇流场与气动性能的影响规律。结果表明：较大的桨尖间隙对桨尖涡的抑制作用减弱,导致涵道风扇的气动性能大幅降低,同时,桨尖涡强度的增加极大地改变了桨叶尖部的压强分布;当桨盘由基准位置向涵道出口方向移动时,涵道内的诱导速度分布发生改变,涵道唇口处气流速度降低,唇口的负压峰值下降,导致涵道拉力减小,当桨盘位置由基准位置向涵道入口移动时,桨尖涡的卷起对唇口处流动的影响加剧,导致涵道拉力明显降低。     

应用多块对接结构网格方法的直升机涵道尾桨气动特性分析

基于湍流Navier-Stokes方程,建立了一套直升机涵道尾桨流场及气动特性的CFD分析方法.在该方法中,针对直升机涵道尾桨流场的特点以及动量源方法对网格系统的要求,提出了一种多块对接网格生成方法,分块贴体网格采用求解Poisson方程的方法生成.通过结合涵道尾桨的运动方式、几何特征及气动特征,建立一个包含动量源项的N-S方程的涵道尾桨流场计算方法和迭代流程;为较好地捕捉涵道壁附近存在的旋涡及分离现象,采用了一方程S-A湍流模型.通过对ONERA M6机翼的绕流特性的计算分析,以及对涵道尾桨在悬停状态下的拉力、诱导速度、载荷分布等的计算分析,验证了该CFD方法的有效性.在此基础上,采用该方法进行直升机悬停、侧飞和前飞状态下涵道尾桨流场与气动特性的数值分析,得到了关于涵道尾桨流场和气动特性的一些有意义的结论.     

共轴刚性旋翼气动干扰及操纵特性分析

建立了一套基于高效配平策略和嵌套网格技术的共轴刚性旋翼气动数值分析方法。流场求解采用雷诺平均Navier Stokes方程,空间离散采用Roe monotone upstream centered schemes for conservation laws（MUSCL）格式,湍流模型为Spalart Allmaras模型;配平计算中引入“差量修正”策略。分析了共轴刚性旋翼模型在不同飞行状态下的气动干扰现象和操纵特性。结果表明:所建立的配平方法在同等精度下,比传统方法可显著节省计算时间;共轴刚性旋翼间距较小,其气动干扰明显强于普通共轴旋翼;悬停状态下,由于受上旋翼下洗流影响,下旋翼总距较大但拉力较小;小速度前飞时气动干扰引起下旋翼气流不对称性增大,双旋翼周期变距差也变大;采用适当的升力偏置量,可以提高共轴刚性旋翼气动效率。      

共轴双旋翼非定常流场干扰特性

直升机共轴双旋翼相遇过程的非定常气动干扰在其周期性干扰中最为强烈。为了研究该过程中的干扰特征,建立了一个基于非定常雷诺平均Navier-Stokes方程的气动干扰数值方法,双旋翼的反转运动采用运动嵌套网格进行模拟。为了探究双旋翼相遇时的气动干扰机理,以两个相对运动的双翼型系统来模拟双旋翼特定展向截面的非定常相遇过程。分析了双旋翼和双翼型相遇时的气动特性和流场特征,并对双翼型系统进行了参数影响研究,结果表明:双旋翼上下桨叶相遇时,上下旋翼拉力均会出现先增后减的变化趋势且波动幅值分别为其对应总拉力的30%和22%,双翼型系统的升力波动趋势与之相似;间距增大使上下翼型间的气动干扰减小,且上下翼型升力波动对间距变化的敏感时间不同;翼型相对厚度增大使双翼型升力波动幅值增大,并出现二次波动;有附加来流时,相对来流速度较小的翼型升力波动的幅值与范围更大。      

典型直升机旋翼翼型气动特性试验研究

采用高升阻比特性的翼型是提高直升机旋翼气动性能的关键。对在典型的OA309旋翼翼型基础上开发的CH309翼型,进行了低、高速气动特性风洞试验研究。试验分别在FL14风洞和NF-6风洞中进行,采用表面测压和尾耙型阻测量技术。试验结果表明：CH309翼型的总体性能优于OA309翼型。试验结果为直升机旋翼设计时翼型的选取提供了参考。     

基于Navier-Stokes方程的旋翼前飞状态气动特性数值模拟研究

建立了一个适用于旋翼前飞状态非定常流场和气动载荷计算的数值方法。在该方法中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上将三阶逆风格式(MUSCL)与通量差分分裂方法相结合。为较好地模拟流动分离,采用了一方程的Spalart-Allmaras湍流模型,时间方向则采用双时间法推进求解。为计入桨叶的旋转、周期性挥舞和变距运动,采用了运动嵌套网格方法。此外,为了更真实地反映旋翼桨叶的实际运动,发展了一个旋翼配平分析模型及求解方法。应用所建立的方法,首先计算了具有先进气动外形的四片桨叶的UH-60A直升机旋翼悬停状态的流场及气动性能,验证了该方法对悬停状态旋翼气动性能数值模拟的有效性。在此基础上,着重对Caradonna模型旋翼无升力前飞状态和AH-1G直升机旋翼桨-涡干扰前飞状态的气动特性(表面压强分布、扭矩、气动载荷)进行了计算,数值结果和试验值吻合良好,表明本文建立的CFD方法对旋翼前飞状态流场的模拟和气动载荷的计算是有效的。     

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

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

OA212翼型主动流动控制的数值模拟研究

采用数值模拟的方法,探讨了基于零质量射流的主动流动控制技术对OA212旋翼翼型动态失速的控制效果和控制特性.以积分形式雷诺平均Navier-Stokes(N-S)方程为控制方程,采用格心有限体积法进行求解.空间离散采用AUSM~+-up格式,时间推进采用含牛顿型LU-SGS子迭代的全隐式双时间法,且引入了预处理方法和多重网格方法加速收敛.通过在喷口上施加非定常边界条件来模拟射流对翼型绕流的影响.研究了不同类型射流、不同位置射流以及不同控制参数(频率、相位、偏角、动量系数等)对动态失速控制效果的影响.研究表明:零质量射流和传统的定常射流均可减小动态失速迟滞环的回线面积,但在提高最大升力方面零质量射流明显优于定常射流;在12%c和62%c处施加组合零质量射流的控制效果最为明显.     

共轴刚性双旋翼/机身干扰流场数值模拟

采用滑移网格技术求解Navier-Stokes (RANS)方程的方法,研究了共轴刚性双旋翼/机身的干扰问题。通过Caradonna-Tung旋翼、Robin直升机、Maryland直升机旋翼/机身干扰和Harrington 2共轴双旋翼等算例,验证了所提出的旋翼流场数值模拟方法的正确性。在此基础上,以Maryland机身为原型,分析了不同桨距的共轴刚性双旋翼与机身之间的干扰特性。结果表明:所提出的数值模拟方法能够很好地模拟共轴刚性双旋翼/机身的气动干扰特性;由于机身对于共轴刚性双旋翼下洗流场的阻滞作用,旋翼的悬停效率增加5%左右,并且随着拉力系数的增大使得悬停效率的增量更加明显;旋翼的悬停效率增加主要来源于下旋翼0°方位角附近的桨叶升力系数的增大,并且拉力系数的增量由桨根向桨尖方向逐渐减小。