振荡翼型非定常流动数值模拟研究

进一步完善了非定常雷诺平均N-S方程求解程序.并借助运动网格生成方法,对周期性俯仰运动下NACA0012翼型的非失速绕流和动态失速问题进行了数值模拟研究,分析了振荡翼型的迟滞效应以及俯仰运动各个参数对迟滞效应的影响.数值模拟结果与实验结果对比表明,本文发展的数值模拟方法在模拟非失速绕流时是可靠的;结果还表明,随着振幅、折合频率的增加迟滞效应增强,平均攻角对迟滞效应无影响;保持无量纲参数--振幅和折合频率的乘积--恒定时,迟滞效应相同.在模拟动态失速问题时.本文方法在深度失速下结果与实验符合较好,而在轻度失速下符合较差.经分析表明这是不同条件下非定常流动对湍流模型和转捩的敏感程度不同所致.     

双三角翼俯仰振荡运动的流场特性数值模拟

 建立了适用于双三角翼大迎角非定常分离流场模拟的数值方法,研究双三角翼俯仰振荡时的动态流场特性,给出动态流场结构和气动力性能随迎角的变化规律,重点考察了减缩频率、转轴位置、平均迎角和振幅等参数对动态流场迟滞效应和气动力曲线迟滞环的影响。研究结果表明:俯仰振荡到相同大迎角时上仰和下俯的流场存在明显差异;减缩频率对气动力迟滞效应的影响相对大于转轴位置;平均迎角的变化导致双三角翼背风区流场结构呈现不同流态,而振幅的大小决定这些流态的数目,事实上俯仰运动时如果跨越的流态数目越多则流场结构的动态响应滞后现象就越显著。通过数值分析,有利于提高对双三角翼在俯仰振荡运动条件下的非定常特性和流场滞后效应等非线性现象的认识。     

基于气/固耦合非定常流动的叶栅颤振分析

采用二维N-S方程／结构振动方程组耦合数值方法，分析计算了二维非定常气动力强迫振动条件下振动能量及气动力对叶栅所做的功，分别根据振动能量和气动力做功分析了叶栅的颤振特性。NACA0012和PROF叶栅在不同折合频率下气动参数随时间的变化表明，折合频率是影响叶栅颤振的重要因素。气固耦合方法得到的NA-CA0012叶栅颤振折合频率比非耦合的结果高，而对PROF叶栅则相反。这表明叶栅结构动力参数对其发生颤振时的折合频率影响很大且很复杂，要准确预测颤振应该考虑叶栅结构动力与气动力的耦合因素。     

环量控制翼型气动力迟滞效应研究

为研究环量控制翼型在非定常运动中的气动力迟滞效应,通过求解环量控制翼型的强迫俯仰振动数值模拟数据,研究了不同减缩频率、振幅和射流动量系数对环量控制翼型气动特性的影响规律,分析了不同动量系数对非定常流场变化的作用机理。研究结果表明:减缩频率和迎角振幅越大,气动力迟滞效应越明显;小迎角下动量系数越大,迟滞效应越微弱,而大迎角下动量系数越大,迟滞效应越明显,并表现出提前动态失速的特征;气动力迟滞效应主要受翼型上表面压力分布影响,与后缘涡和上表面分离涡关系密切。所得结论对于环量控制技术在无人飞行器控制和操纵方面的应用具有一定的实际意义。     

80°后掠三角翼强迫俯仰、自由滚转双自由度耦合运动特性数值研究

耦合求解NS方程和刚体动力学方程数值研究了80°后掠三角翼强迫俯仰、自由滚转双自由度运动特性。内容包括:对比双自由度俯仰/滚转运动和单自由度滚转及俯仰运动特性的差异,研究动态俯仰角对三角翼滚转运动振幅和频率的影响以及三角翼双自由度运动条件下的气动力特性。结果表明:在动态俯仰角的影响下,细长三角翼不能形成极限环形式的振荡;俯仰运动的迟滞效应造成滚转振幅的变化落后于俯仰角的变化;俯仰耦合滚转运动将会导致飞行器升力降低,横向稳定性变差。     

错频对叶片的气动弹性稳定性影响

采用求解动网格下的非定常可压缩Navier-Stokes方程模拟了振荡叶片下的气动弹性问题,研究了叶片的频率错频对叶片气动弹性稳定性的影响.通过数值模拟平面叶栅的气动弹性第10标准算例,验证了气动弹性的数值模拟方法,计算了不同叶片间相位角和折合频率下的气动阻尼系数,研究了叶片振动频率改变对叶片气动弹性稳定性的影响.计算结果说明:频率错频是提高气动弹性稳定性的有效方法,其主要作用是减小叶片间振动的耦合效应和叶片间相位角的影响,并且随着错频量增大叶片稳定性增强;通过模拟三维环形叶栅的气动弹性第4标准算例,计算了气动阻尼系数随错频量变化的规律,验证了错频量和气动弹性稳定性增强的规律.      

非定常运动下的激波/边界层干扰分离特性研究

采用高精度数值方法求解非定常雷诺平均N-S方程和单自由度强迫俯仰运动方程,研究二维压缩拐角模型在等速抬头/低头和周期性俯仰振动等运动方式下的非定常激波/边界层干扰问题,考察角速度、振幅和频率等参数对分离区非定常变化规律的影响。结果表明,压缩拐角等速抬头运动时分离区减小、等速低头运动时分离区增大,周期性俯仰振动时分离点、再附点以及分离区大小的变化都与模型振动频率相同。与相同来流条件下定姿态计算对比分析表明,压缩拐角壁面运动产生的非定常迟滞作用使分离区大小变化减慢,俯仰振动振幅对分离区大小变化影响较大、而对变化的相位影响较小,而振动频率的改变对分离区大小变化的相位影响很大、对分离区大小影响较小。     

飞行器自激振荡的流动物理与动力学机制

全局亚迭代耦合求解流体动力学方程和刚体动力学方程,研究轴对称飞行器单自由度俯仰失稳运动的非定常特征.高超音速锥-柱-裙飞行器的单自由度俯仰失稳运动发展为极限环形式,周期性运动伴随着波系结构的非定常变化;自激振荡源于静稳定构型的振动回复机制和尾裙激波流动迟滞效应所构成的阻尼机制.基于第二拉格朗日方程和虚功原理,导出能够描述迟滞现象的参数化非线性动力学模型.多尺度近似分析获得参数化运动特征:自激振动过程是拟简谐运动;静稳定构型才能出现自激振荡周期运动;平衡点阻尼是决定运动稳定特性的分叉参数;振幅特性与阻尼非线性相关,频率特性与刚度非线性相关.基于数值结果的参数辨识,非线性模型的理论分析、重构都与数值结果高度一致,从而有效地佐证了自激振荡建模研究的合理性.     

二维叶栅气固耦合颤振分析

引入叶片结构动力学方程,通过数值求解气动/结构耦合方程体系,研究了气固耦合条件下叶栅内的二维非定常粘性流动,并初步进行了叶栅的颤振分析。针对两种叶型分别计算了不同折合频率下的流场,通过比较振动叶栅相关参数随时间变化的曲线,表明折合频率是决定叶片颤振与否的重要因素。      

压气机线性振荡叶栅气弹试验研究（一）：非定常气动响应

 为了加强对叶轮机三维气动弹性机理的理解,为数值模拟方法提供试验验证数据,开展了压气机线性叶栅气动弹性试验研究。本研究分为两部分,本文为第1部分,研究压气机线性振荡叶栅的气弹稳定性,第2部分重点辨别叶尖间隙对此振荡叶栅气弹稳定性的影响。组建了线性叶栅试验台,其试验段由7个真实可控扩散叶型组成,中间叶片由凸轮机构驱动模拟三维一阶弯曲振动。使用外置压力传感器在6个叶高处测量叶片表面的定常压力分布和3个折合频率下叶片表面非定常压力分布。由管传递函数修正连接管造成的非定常信号失真。运用影响系数法构造协调叶栅任意叶片间相位角下的非定常气动响应结果。试验结果表明：非定常气动响应具有完全的三维特征;折合频率和叶片间相位角对叶栅的气动弹性稳定性具有决定性影响。     

边界层抽吸和脉动反压作用下隔离段内流动特性研究

为了解上游边界层抽吸控制和下游周期脉动反压作用下隔离段内流动特性,采用非定常数值模拟和理论分析相结合的方法,对来流Ma=2情况下的隔离段内激波串动态演化特性、激波串形态结构变化以及激波串演化迟滞现象进行了研究。结果表明,在脉动反压和边界层抽吸作用下,激波串在上游抽吸狭缝与下游隔离段出口之间周期振荡,振荡频率与脉动反压一致。在振荡过程中,首道激波串形态在规则反射与马赫反射以及马赫反射与弧形激波(包含正激波)之间相互转换。边界层抽吸将激波串固定在抽吸狭缝位置,有效提高了隔离段抗反压性能,脉动频率越大,可承受的瞬态反压峰值越大。在一个振荡周期内,激波串向上移动速度较向下移动更快,且在上下移动过程中形态变化存在迟滞现象。     

三角翼俯仰振荡中若干参数影响的实验研究

本文研究一个 6 0°三角翼以不同振荡频率 ,不同振幅从不同起始迎角进行俯仰振荡的纵向气动力特性 ,研究发现跨越具有不同时间尺度的流态的振荡会产生明显的气动力迟滞现象 ,这一迟滞现象随减缩频率的增大有一个极限值 ;在同一种流态下振荡 ,气动力迟滞现象不明显 ,这一结果对充分利用非定常气动力有指导意义。     

Nonlinear aeroelastic analysis of airfoils: bifurcation and chaos

Different types of structural and aerodynamic nonlinearities commonly encountered in aeronautical engineering are discussed. The equations of motion of a two-dimensional airfoil oscillating in pitch and plunge are derived for a structural nonlinearity using subsonic aerodynamics theory. Three classical nonlinearities, namely, cubic, freeplay and hysteresis are investigated in some detail. The governing equations are reduced to a set of ordinary differential equations suitable for numerical simulations and analytical investigation of the system stability. The onset of Hopf-bifurcation, and amplitudes and frequencies of limit cycle oscillations are investigated, with examples given for a cubic hardening spring. For various geometries of the freeplay, bifurcations and chaos are discussed via the phase plane, Poincaré maps, and Lyapunov spectrum. The route to chaos is investigated from bifurcation diagrams, and for the freeplay nonlinearity it is shown that frequency doubling is the most commonly observed route. Examples of aerodynamic nonlinearities arising from transonic flow and dynamic stall are discussed, and special attention is paid to numerical simulation results for dynamic stall using a time-synthesized method for the unsteady aerodynamics. The assumption of uniform flow is usually not met in practice since perturbations in velocities are encountered in flight. Longitudinal atmospheric turbulence is introduced to show its effect on both the flutter boundary and the onset of Hopf-bifurcation for a cubic restoring force.     

垂直交叉双圆柱绕流强迫振动数值模拟

采用虚拟体法数值模拟了平面单圆柱绕流强迫振动,以及空间垂直交叉双圆柱绕流下游圆柱强迫振动流场。验证了单圆柱强迫振动中的锁定状态以及相位突变现象,从而证实了该数值方法模拟振动流场的可靠性。研究了在雷诺数Re=150、间距为5倍圆柱直径、下游圆柱按正弦曲线振动时,对尾流场的影响。下游圆柱两端尾流场在振动的作用下,涡的横向间距增大,而中心尾流场由于受上游圆柱尾流的影响而保持原先的状态。下游圆柱在锁定区的振动使尾流场变得稳定,表现在流向二次涡结构的减少,以及尾流场的速度场频谱趋向单一化。     

振动鸭翼复杂流场测量

研究鸭式布局飞机模型振动鸭翼对翼面涡流场的非定常干扰影响，进行了有无鸭翼、鸭翼不同偏角和不同振动频率，不同振动平均偏角及不同模型攻角下的主翼面涡流场静动态流动显示和翼面及立尾上压力分布测量。分析上述参数对主翼涡大小和强度、主翼涡位置和破散特性，压力分布特性的影响及其造成该现象的上下洗效应，涡系干扰和动态迟滞特性等复杂流动机理。     

三角翼过失速非定常洞壁干扰修正

分别以固壁条件和洞壁附近实测的压力分布,模拟各类实壁和透气壁试验段的洞壁边界条件,利用N-S方程数值求解模型在风洞中的绕流场,得出洞壁干扰对跨音速模型绕流的影响和气动力修正结果。初步研究结果表明,该方法能较有效地模拟模型在跨音速风洞中的绕流场,GBM04A模型在0.6m风洞中的试验结果经洞壁干扰修正后与无干扰参考结果吻合较好