预测悬停旋翼跨音速流动的一种组合自由尾流/CFD方法

一种新的组合自由尾流 / CFD方法用于悬停旋翼流动的 CFD解中以考虑实际尾流的作用 .用文中描述的尾流分析方法研究了螺旋尖涡的运动 .首先从广义尾流模型开始 ,用半经验公式模化了涡核对旋翼尾流的作用 ;然后在环量收敛和尾迹收敛的条件下完成了自由尾迹计算 ;最后应用 Jameson有限体积龙格 -库塔推进格式求解了欧拉方程 .所得结果与相关文献和实验数据进行了比较     

复杂旋翼流场的耦合欧拉-拉格朗日数值方法

针对影响旋翼流场求解精度的关键因素“桨叶复杂近体流动”和“尾迹涡畸变”,结合计算流体力学(CFD)方法和黏性涡方法,发展了一套适合于复杂旋翼涡流场分析的耦合欧拉-拉格朗日数值方法:为捕捉桨尖三维效应、激波等细节流场特征,在桨叶近体区域采用CFD方法对其进行求解;针对高雷诺数旋翼流场中桨尖涡的紧凑结构特点,引入黏性涡方法建立了高分辨率的尾迹求解模型;两计算域间的信息交换采用了集中涡源法和边界修正法.应用所建立的计算方法,以旋翼CFD标准验证试验(Caradonna-Tung旋翼)为算例,对尾迹影响明显的悬停状态进行了数值模拟,通过对比耦合边界处流场特征及桨叶表面压力系数分布,验证了方法的有效性.此外,还从旋翼尾迹捕捉精度、涡量耗散特征及计算时间等方面对不同计算方法进行了对比分析,结果表明耦合方法可充分发挥CFD和黏性涡方法各自的优点,在旋翼流场数值模拟方面具有独特的优势.     

一种新的自由涡尾迹计算方法

在使用自由涡理论进行旋翼气动计算时,尾迹几何结构的确定是计算旋翼气动载荷以及空间流场计算的关键。发展了一种新的自由涡尾迹结构计算方法,该方法能快速得到收敛的尾迹结构,能比较准确地计算桨叶载荷分布。算例与实验数据的比较验证了该方法的正确性。     

基于黏性涡模型的旋翼流场数值方法

 建立了一种适用于旋翼非定常流场特性分析的黏性涡数值方法。在该方法中:流场中的大尺度涡被离散为若干微小的涡元,通过求解涡量-速度形式的Navier-Stokes方程模拟涡元的输运等过程;黏性扩散效应采用高精度的粒子强度交换法进行计算,而桨叶附着涡以及新生涡环量采用了Weissinger-L升力面理论进行求解;为显著提高计算效率,在诱导速度及其梯度的计算中还引入了快速多极子算法(FMM)。应用上述方法,对悬停和前飞状态下的多个旋翼流场算例进行了计算,通过对比旋翼尾迹涡量特征和诱导速度分布等,验证了该方法的有效性。此外,还将本方法与旋翼计算流体力学(CFD)方法及传统的自由尾迹方法进行了比较,结果表明黏性涡方法在兼顾效率的同时,还能够更好地捕捉旋翼尾迹运动。     

前飞状态旋翼尾迹测量试验研究

利用PIV技术,测量前飞状态旋翼尾迹的横向速度分布及桨尖涡在横向剖面里的运动轨迹,得到了前飞状态旋翼两侧的尾迹边界及桨尖涡在运动过程中的耗散特性等,研究了风速及拉力系数对前飞状态旋翼尾迹和桨尖涡运动轨迹的影响。结果表明：前飞状态下,旋翼左右两侧尾迹的涡量值基本相当。旋翼尾迹沿径向急剧收缩,沿垂向逐渐下降,下降高度与流向距离几乎呈线性关系。风速、拉力系数的变化对桨尖涡及其运动轨迹有显著的影响。获得的测量结果为开展旋翼流动机理研究及提高CFD分析精度提供了试验依据。     

地面效应中悬停旋翼的自由尾迹计算

发展了一种具有良好收敛性的自由尾迹模型,用于对地面效应(IGE)状态下的旋翼悬停尾迹进行数值模拟。其中使用直涡元对桨尖涡丝进行离散,使用面元法模拟地面对流场的影响。针对以往类似研究中旋翼高度较低时自由尾迹迭代不易收敛的问题,引入了地下尾迹等比例修正并提出了环式地面面元分布等多项修改办法来改善计算模型的收敛性,同时这些改进也起到了削减计算量的作用。为了对新构建模型进行验证,以多个高度下的旋翼有地效状态为算例,对其尾迹几何结构进行了计算并与实验值进行了对比,结果显示二者符合得较好,证实了该模型具有较高的准确性。     

悬停和前飞状态倾转旋翼机的旋翼自由尾迹计算方法

将倾转旋翼机的旋翼尾迹问题区分为二类,即悬停和前飞等状态的定常尾迹以及过渡飞行状态的非定常尾迹.基于Weissinger-L升力面理论和卷起桨尖涡模型,建立了一个悬停和前飞状态的倾转旋翼自由尾迹分析方法.该方法不仅适合于单旋翼尾迹计算,也适用于双旋翼尾迹分析.分别以ATB旋翼和缩比的V-22旋翼为算例,计算了悬停状态的旋翼自由尾迹和诱导速度,并与可得到的试验数据进行了对比,验证了该方法的有效性.本文也对双旋翼的干扰尾迹进行了计算,表明了悬停和前飞状态不同的尾迹干扰影响.     

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

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

基于PC2B格式自由尾迹方法的旋翼响应特性分析

旋翼自由尾迹模型的建模和桨叶动态响应计算的问题是直升机空气动力学领域中富有关键意义的课题。本文结合桨叶气动模型、旋翼涡尾迹模型、桨叶挥舞动力学模型、"PC2B"尾迹求解算法和旋翼配平模型,建立了一个时间精确的旋翼自由尾迹和桨叶动态气动响应的分析方法,并针对旋翼总距突增时的旋翼载荷和桨叶动态响应进行了计算和分析,得出了一些有意义的结果。     

旋翼VPM/CFD耦合模型的建立及其在小前进比气动分析中的应用

为了提高旋翼计算流体力学（CFD）流场计算效率及克服其在尾迹捕捉上的不足,将旋翼黏性涡粒子方法（VPM）与CFD分析相结合,建立了一个新的旋翼VPM/CFD耦合气动分析模型。在该模型中,采用VPM分析以实现对旋翼尾迹中黏性涡的高效捕捉而不引入数值耗散,而采用CFD分析用于精确地模拟旋翼桨叶近体区域内复杂的流动现象,同时也为VPM分析提供一个较高精度的涡源模型。至于两者信息交换,则使用集中涡源法将CFD信息传递至VPM分析,而VPM计算得到的尾迹信息则通过边界修正施加至CFD域的远场边界上,从而可以鲁棒地实现CFD域与VPM域的耦合计算。在此基础上,对“Helishape 7A旋翼”小前进比前飞桨-涡干扰（BVI）状态进行了较为深入的气动分析,计算结果表明：与全CFD计算比较,建立的VPM/CFD耦合分析模型可以有效地避免旋翼尾迹区桨尖涡的数值耗散,从而更加可靠地捕捉桨-涡干扰状态下的桨叶非定常气动载荷脉动,同时对于本文算例,计算效率可以提高30%以上。     

基于非定常面元/黏性涡粒子混合法的旋翼/平尾非定常气动干扰

旋翼/平尾非定常气动干扰是导致直升机纵向“抬头（Pitch-up）”现象的主要原因。为在直升机设计阶段准确分析旋翼/平尾非定常气动干扰以及由此引起的低速纵向操纵特性变化,通过涡量等效原则和Neumann物面边界条件建立了适用于旋翼/平尾气动干扰分析的非定常面元/黏性涡粒子混合法。该方法耦合了考虑尾迹时变效应的非定常面元法、黏性涡粒子法及涡量镜面法,以准确模拟旋翼和平尾的非定常气动载荷、旋翼尾迹的非定常特性以及旋翼尾迹对平尾的气动干扰效应。首先通过计算NASA ROBIN（Rotor Body Interaction）旋翼尾迹几何和诱导速度分布,并与实验测量值、时间精确自由尾迹及CFD计算结果对比验证方法的准确性。相比于时间精确自由尾迹,本文方法计算精度更高。随后分析了旋翼/平尾非定常气动干扰对平尾向下气动载荷和气动导数的影响,并分析了平尾构型对旋翼/平尾非定常气动干扰的影响规律。分析表明：旋翼尾迹与平尾干扰导致低速状态的平尾载荷突增,气动导数反号;低平尾气动载荷突增较大,高平尾较小,但高速气动导数反号;前置平尾载荷突增量减小,但对应速度范围较宽;右旋直升机右平尾载荷突增量较小,但气动导数特性基本不变。     

旋翼桨尖涡模型及其在自由尾迹分析中的影响

 建立了一个旋翼桨尖涡初始涡核半径计算模型,用于消除旋翼自由尾迹分析中桨尖涡涡核半径依靠经验设定的现状。该方法利用涡量和涡矩守恒关系,考虑旋转效应的影响,并修正Vatistas涡量分布函数,得到了一个关于桨尖涡初始涡核半径与桨叶附着环量分布之间的非线性关系。将新建立的初始涡核模型同先前建立的涡核扩散与拉伸模型以微分方程初值问题的形式组成一个完整的桨尖涡模型,并将其应用于自由尾迹分析方法中,分别对两副模型旋翼的悬停尾迹进行分析。计算得到的桨尖涡初始涡核半径和展向位置与实验测量值的对比表明本文模型具有较高的准确性。此外,讨论了涡核湍流动量扩散效应在自由尾迹分析中的影响,分析表明该效应对自由尾迹分析的收敛性有重要影响。     

双旋翼自由尾迹计算分析中的加速算法

对快速多极展开算法(FMM)在共轴式和纵列式双旋翼自由尾迹互诱导速度计算中的应用进行了研究。介绍了FMM算法的基本原理,建立了FMM算法的数学计算模型,并对其加速特性和准确性进行了分析;将方法引入到基于直涡元的自由尾迹互诱导速度的计算分析中,推导了诱导速度的多极展开表达式。在此基础上,以某旋翼模型为例,通过对自由尾迹中涡元诱导速度的计算,验证了加速算法的加速度特性和准确性。使用方法对共轴式和纵列式双旋翼悬停和前飞状态下的自由尾迹几何形状进行了计算,分析了双旋翼之间的相互干扰。     

基于Navier-Stokes方程/自由尾迹/全位势方程的旋翼流场模拟混合方法

为充分考虑旋翼尾迹对流场的影响和减少尾迹的数值耗散，建立了一个基于Navier-Stokes方程／自由尾迹分析／全位势方程的旋翼流场求解的新的混合方法。该方法的求解域由三部分组成：一是围绕旋翼桨叶周围的粘性区域．采用可压Navier-Stokes方程来捕捉近场信息，包括激波及尾迹；二是离桨叶较远、粘性可以忽略的等熵流区域，以全位势方程来描述其流动；三是在无粘区域中嵌入自由尾迹模型，模拟桨尖涡从粘性区域进入势流范围的发展变化。为便于流场分区求解和信息传递，采用了嵌套网格方法，并给出了不同区域之间的信息传递方法。以两叶的Caradorma＆Tung模型旋翼和四叶的UH-60A直升机旋翼为算例，计算给出了旋翼桨叶表面的压强分布以及桨尖涡的位置，并与可得到的试验数据及无尾迹模型方法的计算结果进行了对比，表明本文的混合方法能够很大程度地减少旋翼尾迹的数值耗散。     

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

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

倾转状态的旋翼尾迹形状和挥舞响应分析

本文推导了旋翼在倾转时的桨叶挥舞运动方程,给出了桨叶升力面模型和尾迹时间推进预测—校正算法,并结合涡核模型和旋翼配平分析模型,建立了一个时间推进旋翼自由涡计算方法。通过计算的旋翼诱导速度与实验值的对比,验证了本文方法的有效性。然后利用该方法,计算了旋翼在倾转时的尾迹几何形状畸变和桨叶挥舞响应特性。     

悬停状态下旋翼尾迹测量试验研究

采用粒子图像测速(PIV)技术,对悬停状态下直升机旋翼尾迹进行了全面测量,获得了在不同试验条件时旋翼流场的涡量和速度分布。测量结果合理可信,为开展旋翼流动机理研究及提高CFD分析精度奠定了基础。     

悬停状态旋翼尾迹边界测量

介绍了用试验来定量确定悬停状态旋翼尾迹边界的方法，从物理本质分析了悬停状态旋翼尾迹的湍流特性，阐述了如何用热线风速仪来测量尾迹的湍流度并用湍流度来确定桨尖涡的位置和尾迹边界。揭示了悬停状态旋翼尾迹边界的一些特点，最后将试验结果与国外的相应研究结果作了比较，发现两者吻合得较好。     

旋翼地面效应的气动建模与特性

地面对旋翼气动特性影响明显,且导致旋翼流场更加复杂。为分析地面效应下的旋翼桨尖涡和流场变化特性,基于涡面和无滑移边界条件,求解第2类Fredholm方程获得地面涡面矢量分布,且将涡面矢量按涡扩散方程扩散到流体中,建立考虑黏性效应的地面气动模型,并耦合非定常面元/黏性涡粒子混合法以体现旋翼桨叶气动特性和旋翼尾迹的非定常效应,构建旋翼地面效应气动分析方法。通过计算Lynx尾桨地面效应下的性能和桨尖涡轨迹,并计算Maryland大学模型旋翼和NASA缩比旋翼地面效应下的垂向、径向速度分布,且与试验和CFD计算结果对比,验证了本文方法能较好捕捉地面效应下的旋翼尾迹变化特性和复杂速度场特性,且结果表明本文方法能较好模拟地面效应下旋翼桨尖涡的收缩、扩散、井喷、地面射流等物理现象。     

基于时间推进自由尾迹方法的俯仰状态旋翼尾迹和入流特性计算

推导了旋翼俯仰运动时的桨叶挥舞运动方程,给出了桨叶升力面模型和尾迹时间推进预测一校正算法,并结合涡核模型和旋翼配平分析模型,建立了一个时间推进旋翼自由尾迹计算方法.通过计算的旋翼入流值与实验值的对比,验证了计算方法的有效性,利用该方法,计算了旋翼在俯仰运动时尾迹几何形状的畸变特点,并分析了俯仰运动时旋翼的入流特性.