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

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

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

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

直升机桨—涡干扰预测方法研究

本文简要介绍了桨—涡干扰(BVI)预测的一些情况(预定尾迹法、自由尾迹法、等环量曲涡元法和CFD/CSD耦合法等)。总的来说,在桨—涡干扰预测方面已取得了很大进步,但要达到准确预测桨—涡干扰尚需作进一步努力。     

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

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

机翼尾迹脱落涡噪声的理论研究

Ｍｉｃｈｅｌ等人１９９８年应用平面传声器阵列对飞机过顶噪声进行的测量研究首次发现，机翼尾迹脱落涡噪声是某些类型飞机重要的噪声源。为发展一种预测这种噪声源的理论预测模型，应用ｖｏｎ Ｋａｒｍａｎ涡街模型模拟二维机翼下游尾迹脱落涡，尾迹涡的强度和脱落频率应用这个模型进行计算。基于Ｈｏｗｅ后缘噪声理论，并结合尾迹模型，本文发展了一种预测脱落涡噪声声压级和指向特征的气动声学模型。对６架现代商用飞机的机翼尾迹脱落涡噪声的计算表明，本文理论模型预测的涡脱落频率、声压级以及噪声的指向性等与实验测量结果有较好的一致性。     

高负荷对转压气机尾迹涡级间迁移影响因素研究

为探究高负荷对转压气机上游转子尾迹脱落涡级间迁移的影响因素,截取高负荷对转压气机60%叶高叶型,保持叶型几何、总压比不变的前提下,分别取12%,20%,25%,30%,35%,45%下游转子弦长的轴向间隙进行了非定常数值模拟。研究发现,在每个轴向间隙下,尾迹脱落涡迁移至下游转子前缘平面的周向位置均不随时间发生变化;当轴向间隙改变时,尾迹脱落涡迁移路径随之发生变化。基于理论分析和数学推导,建立了高负荷对转压气机尾迹脱落涡级间迁移轨迹的无量纲准则关系式,并结合通流计算与激波模型,发展了预估上游转子尾迹脱落涡迁移轨迹的方法,为设计阶段考虑非定常效应提供指导。     

翼型近尾迹流中大尺度涡结构的理论模型

根据流动稳定性理论相干结构的相关理论 ,提出了一种二维波模型来描述翼型近尾迹流动中的大尺度涡结构 ,通过采用对近尾迹流动区域分区处理的方法计算了近尾迹区域的大尺度涡结构。计算所得的流线图 ,涡量分布以及雷诺应力分布具有与实验结果相似的分布规律。说明该二维模型能够较好的描述近尾迹流动中这种大尺度涡结构的运动学和动力学特性      

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

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

尾流撞击效应对轴流压气机下游叶排涡面的影响

针对某高负荷高压压气机末级,通过研究上游转子尾迹涡面周期性扫掠对下游静子涡面的影响,从涡面非定常作用的角度分析流场时空结构.结果表明,流场性能的改善与涡面时空结构相对应,尾迹涡面同分离涡面的相互作用可显著降低流动损失;与上游尾迹涡面发生耦合的是叶片尾缘旋涡交替脱落频率而非吸力面前缘分离涡不稳定频率.      

有/无尾迹作用下低压涡轮叶栅分离边界层转捩的大涡模拟

采用经过大量算例验证的可压缩大涡模拟求解器对雷诺数为60154、马赫数为0.402的低压涡轮叶栅T106D-EIZ进行了细致模拟,计算了定常来流和周期性尾迹来流两种工况.对计算结果的分析表明:定常来流工况下,叶片吸力面后部出现大尺寸的层流分离泡,分离剪切层的转捩过程受Kelvin-Helmholtz （K-H）不稳定性控制;尾迹来流工况下,由于来流尾迹的周期性扫掠,时均分离泡尺寸变小,叶栅总压损失降低.对相位平均和瞬态流场的分析表明,尾迹引起的逆射流使分离点后移,形成卷起涡结构,逆射流掠过卷起涡的过程中与其发生强烈的相互作用,产生大量气动损失,而后卷起涡破碎,流动转捩为湍流.      

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

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

确定性涡方法的精确度研究

本文以无界流场中圆形涡扩散问题为模型问题，对几种为典型的确定性涡方法进行了讨论，分析了方法的参数条件，数值结果表明，确定性涡方法比随机涡方法有更高的精度。     

动态尾涡间隔计算工具的开发与实现

基于建立的尾涡消散模型、运动模型、遭遇模型,以Delphi 7.0为平台,采用面向对象的设计思想,开发了动态尾涡间隔计算工具,包括前机诱导尾涡流场计算、后机遭遇响应过程分析、安全间隔计算等功能模块。根据输入的飞机参数、气象数据、机场数据,实现对尾涡间隔的实时、高效计算。基于开发的计算工具,以B777-200和B737-800为前后机进行了验证计算。     

Numerical investigation of full scale coaxial main rotor aerodynamics in hover and vertical descent

The original free vortex wake model was used for numerical investigation. Calculation of the aerodynamic characteristics in hover and vertical descent modes in the range of vertical descent speed of 0–30 m/s including the Vortex Ring State (VRS) area was performed. The calculations were carried out under the condition of variable blade pitch angle values providing a fixed time-average thrust value. Visualization data of free vortex wake shapes, flow structures, and velocity fields were obtained and analyzed. The time-dependences of the rotor’s thrust and torque coefficients were obtained and analyzed. The obtained data allows determining the boundaries of the VRS area by various criteria such as rotor thrust and torque pulsations, growth of rotor power consumption relative to the hover, growth of rotor induced velocities relative to the hover, and growth of the required rotor blade pitch angles values. The results of the study are compared with experimental and calculated data of other authors and can significantly supplement the available results of experimental and computational studies in this area.     

虚拟体法数值模拟双圆柱强迫振荡系统

用虚拟体方法研究了均匀流及尾流中圆柱横向振荡的流动特征.单圆柱结果与相关试验吻合,得到了"锁定"区域,发现涡量跳跃现象.对小间距双圆柱的数值模拟发现,在"锁定"区域附近随激励频率的增加流场特征会出现若干显著的变化,包括涡量的跳跃,涡模态的改变,正负涡量带的形成和消失.在大间距情况下发现下游圆柱的振荡对上游尾涡脱落没有明显的影响,下游圆柱尾流场没有发现"锁定"现象.     

基于气弹耦合的旋翼桨涡干扰噪声预测方法

建立了基于气动/弹性耦合的旋翼桨涡干扰（BVI）气动和噪声分析方法。气动模型包括修正Beddoes尾迹模型和CFD模型，噪声计算采用基于声学类比法推导出的FW H（Ffowcs Williams Hawkings）方程，弹性桨叶动力学建模采用有限元方法。应用所建立的方法，对刚性的OLS（operational load survey）旋翼桨涡干扰状态的气动和噪声特性进行了计算，对比了两种气动模型在研究桨涡干扰问题的有效性；以弹性的HART Ⅱ旋翼为研究对象，分析了桨叶弹性、时间步长对桨涡干扰气动载荷和噪声的影响。结果表明:进行桨涡干扰计算时所采用的时间步长不宜超过2°。CFD方法由于固有的数值耗散，计算出的OLS旋翼噪声声压峰值仅为试验值的60%，而修正Beddoes尾迹模型能够避免数值耗散，且具有高效率的优势。考虑桨叶气动弹性能够提高旋翼桨涡干扰噪声的预测精度。      

反弹自消散系统尾涡控制机制研究

为了研究加快尾涡消散的机理,以获得更合理的尾流间隔,进行了壁面反弹实验的数值模拟计算。实验采取在模拟机翼后缘安置反弹面的方法,诱发反弹二次涡与主涡相交出现不稳定性,以加快尾流衰减。数值计算运用N-S方程,对雷诺应力项采用Realizableκ-ε涡粘模型进行计算,分析经二次涡干涉后尾涡轴向涡量衰减、涡核下沉运动等参数的变化,以探究反弹二次涡对飞机主涡消散的演变影响。实验结果表明,反弹涡对主涡衰减有明显的促进作用。该实验为人工干预缩减机场尾流间隔的研究奠定了理论基础。     

Planar laser scattering visualization of streamwise vortex pairs in a Mach 6 flow

A series of cross-sectional flow fields of Counterrotating Vortex Pairs(CVPs) generated by a large-scale ramp vortex generator is observed using an ice-cluster-based Planar Laser Scattering(PLS) method in a shock tunnel with a nominal flow Mach number of 6. Combined with a numerical simulation, two streamwise CVPs with opposite rotating directions are identified in the wake flow of the vortex generator with an absence of a boundary layer, namely, a Primary CVP(PCVP) and a Secondary CVP(SCVP). Th...