网格变形在俯仰振荡翼型中的应用

为解决航空器结构变形问题,在研究相关非定常流场基础上,建立了包含结构变形的非定常流场数值求解方法,运用无限插值法（TFI）生成任意时刻的非定常动态网格,通过与刚性动网格比较,证明了方法的可行性,并以NACA0012翼型为算例,对翼型俯仰振荡运动中表面柔性变形问题进行了数值模拟,发现俯仰振荡中叠加变形形式为一阶振型时,可以明显地提高升力、增加动稳定、提高升阻比。     

细长旋成体大迎角俯仰振荡的数值研究

为研究俯仰振荡对细长旋成体流场与气动特性的影响,利用雷诺平均方法与大涡模拟的混合方法相结合的脱体涡模拟方法对大迎角下细长旋成体俯仰运动进行了数值模拟研究。通过对背风面分离涡的强度和位置、细长体绕流形态、截面压强系数和侧向力系数的观察和分析,发现施加小振幅、特定频率的俯仰振荡对大迎角细长旋成体背风面流场有明显的控制能力,可以改变流场结构,使背风面非对称流场趋于对称。同时,研究了振荡参数对耦合作用的影响,发现振荡中振幅和频率的增加,加强了俯仰振荡对非对称涡的抑制作用。研究成果对充分利用非定常气动力具有指导意义。     

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

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

旋翼振荡翼型与旋涡干扰问题数值模拟研究

以直升机旋翼桨/涡干扰问题为研究背景,构建了二维振荡翼型-动态旋涡相互作用模型,基于混合笛卡尔网格生成与动态自适应技术,建立了能够准确捕捉流场特征的非定常自适应混合笛卡尔网格（UAHCG）流场数值模拟方法。通过数值模拟振荡翼型绕流问题和旋涡-静态翼型相互作用问题,验证了方法的有效性。最后,针对振荡翼型-动态旋涡相互作用问题,从翼型振荡减缩频率和俯仰运动方向的变化,揭示了旋涡与振荡翼型相互作用过程中产生的复杂流动现象和非定常气动力变化。     

三角翼受迫俯仰滚转运动气动特性的数值模拟研究

基于动网格的技术,利用有限体积法对N—S方程进行数值离散,选用湍流模型,对65°大后掠三角翼的受迫俯仰滚转运动进行了数值模拟,在此基础上,对俯仰滚转耦合运动的气动力特性和流场结构进行了分析。计算结果表明：俯仰滚转耦合运动时,三角翼上表面涡分布的非对称性将产生横侧方向的偏航力矩和滚转力矩,滚转力矩和偏航力矩随着滚转振幅角和滚转缩减频率的增大而增大,但对法向力影响不大。     

基于时间谱方法的振荡翼型和机翼非定常黏性绕流数值模拟

 传统的双时间方法在非定常计算中长时间的过渡迭代推进求解导致其计算效率相对较低,针对周期性非定常流动问题的流动特征,发展了一种基于离散傅里叶变换的高效时间谱方法,用于求解振荡翼型和机翼的非定常黏性绕流。在时空耦合的雷诺平均Navier-Stokes(RANS)方程的求解中,对流项的离散应用了Roe的通量差分格式,物理时间项的离散方法为时间谱方法,伪时间推进采用了隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式。考虑到湍流的时空耦合效应,时空耦合的Spalart-Allmaras一方程湍流模型的物理时间项同样采用时间谱方法进行离散。为了进一步提高计算效率,当地时间步长和多重网格技术等加速收敛的措施均被采用。算例对俯仰振荡NACA0012翼型和Lann机翼的周期性非定常流场进行了数值计算。结果表明:对于周期性非定常流场的数值模拟,相比于传统的双时间方法,用时间谱方法近似物理时间项,不仅能够提高流场的计算精度,而且更能够大幅度提高计算效率。     

利用DES方法对大迎角细长旋成体气动特性受俯仰振荡影响的数值研究

采用DES(Detached-Eddy Simulation)对大攻角下细长旋成体俯仰振动的气动特性进行了数值模拟,分析俯仰振动对细长旋成体流场与气动力特性的影响.通过对背风面分离涡的强度和位置,流场形态,截面压强系数和侧向力系数的观察与分析,发现施加小振幅、特定频率的俯仰振动对大迎角细长旋成体背风面流场有明显的控制能力,将改变流场结构,使背风面非对称流场趋于对称.     

非失速二维振荡叶栅非定常流动数值模拟研究

进一步完善了非定常雷诺平均N-S方程求解程序,借助快速网格生成方法,对周期性俯仰运动下NACA0012叶栅的非失速绕流问题进行了数值模拟研究,并分析了振荡叶栅的迟滞效应以及俯仰运动各个参数对迟滞效应的影响.数值模拟结果与实验结果对比表明,本文发展的数值模拟方法是可靠的;结果还表明,随着振幅,折合频率的增加迟滞效应增强,平均攻角对迟滞效应无影响;保持无量纲参数--振幅和折合频率的乘积--恒定时,迟滞效应相同.     

于拼接网格的栅格舵动态特性数值模拟研究

采用基于拼接混合网格的数值计算方法,模拟了栅格舵在不同马赫数下的静态、动态气动特性,计算结果与试验的一致性说明了所采用的数值方法的有效性。通过模拟发现栅格舵的气动焦点随马赫数先向上游移动然后再向下游移动;栅格舵在不同速域下的非定常俯仰迟滞现象均不明显。     

细长旋成体大迎角非定常运动数值研究

利用DES方法对处于非对称流态下的细长旋成体进行了俯仰和偏航振荡的数值模拟,观察细长体气动特性尤其是背风面非对称分离涡的变化.计算结果表明在大迎角非对称多涡系情况下,固定频率和振幅的非定常运动可以改变流场结构和气动力,对细长体背风面流态有巨大的影响.俯仰振荡对非对称分离涡有明显的控制作用,抑制流场的非对称性,使分离涡趋于对称:而固定频率的偏航振荡则破坏背风面分离涡的稳定性,激励非对称背涡产生随时间变化的周期性脱落.所进行的非定常运动与细长体流场耦合作用研究在国内外研究尚少.     

稳定性导数计算方法研究

稳定性导数的准确计算对于飞行器的操稳特性具有重要意义.应用计算流体力学(CFD)方法中的非结构化动网格技术建立能够模拟飞行器做周期性俯仰运动的强迫振荡法,以国际动导数标模Finner导弹为验证算例,获得不同马赫数下俯仰力矩系数的迟滞环曲线,进而计算Finner导弹在不同马赫数下的静稳定性导数和动稳定性导数.结果表明:本文计算得到的俯仰静稳定性导数与试验数据非常接近;在亚音速和超音速范围内,动稳定性导数计算结果与试验值很接近,但在跨音速范围内,本文计算结果与试验曲线的规律性一致,但误差较大.     

舵面矢量喷流对现代战斗机气动特性的影响

在双三角翼、椭圆机身的先进战斗机气动布局中配置舵面矢量推力装置,系统研究了该舵面矢量推力装置引起的纵向矢量喷流和横向矢量喷流对机翼绕流及其气动特性的影响,同时还研究了舵面矢量喷流的落压比NPR、舵面的纵向和横向偏角D_p、D_y以及纵向舵面的宽度W_ｐ和舵面的配置等因素对喷流干扰效应的影响。     

垂直下降状态下螺旋桨气动特性及涡环特征

螺旋桨在垂直下降状态下极易进入涡环状态,从而引起复杂的空气动力学特性。采用黏性涡粒子方法（viscous vortex particle method, VVPM）,对APC thin electric 10X7螺旋桨在不同垂直下降率下的气动特性和涡环特征开展研究。结果表明：螺旋桨气动特性容易受到下降率的影响,当下降率增加时,拉力呈现先略增大再减小后增大的趋势,同时尾迹压缩效应逐渐增强,涡环中心逐渐上移。在此过程中,涡环先后经历了悬停状态、涡环前期、涡环初期、深度涡环和类风车流动状态,涡环特征从管状螺旋下洗涡逐渐演变成螺旋涡与环形涡并存,再到完全涡环结构,最后涡系结构移动至桨盘上方。涡场拓扑结构的演变直接影响螺旋桨的气动特性,当桨盘附近尾迹为管状螺旋涡结构时,螺旋桨气动特性变化不明显;随着下降率的增加,逐渐演变出的环形涡会严重影响螺旋桨的气动特性,表现为拉力骤降,尾迹压缩效应增强,回流区包裹桨盘;随着下降率进一步增加,涡系结构到达桨盘上方,螺旋桨进入类风车状态,气动力逐渐增大。     

三角翼受迫俯仰滚转耦合运动的气动特性研究

基于刚性动网格的技术,选用B-L湍流模型,利用有限控制体积法对N-S方程进行数值离散,对76°大后掠三角翼的受迫俯仰滚转耦合运动进行了数值模拟,在此基础上,对俯仰滚转耦合运动的气动力特性和流场结构进行了分析.计算结果表明:俯仰滚转耦合运动时,三角翼上表面的涡分布的非对称性将产生横侧方向的偏航力矩和滚转力矩,滚转力矩和偏航力矩随着滚转振幅角和滚转缩减频率的增大而增大,但对法向力影响不大.     

涡流冷却推力室流场特征与性能仿真

针对2000N气氢/气氧涡流冷却推力室,采用三维全尺寸计算模型开展了仿真研究,得到了流场速度分布特点,验证了涡流冷却推力室内具有双向涡旋结构,内外涡流分界面约占涡流冷却推力室圆柱段半径的86%,燃烧区域约占涡流冷却推力室圆柱段半径的70%.分析表明:外层涡流主要受来流速度与涡流冷却推力室几何参数影响,内层涡流在黏性、燃烧等作用下室压、密度稳定.侧壁温度平均为388K,比冲效率达92%以上,仿真结果与试验对比一致.      

卡门涡街与翼型的相互作用

本文利用保角变换,用离散涡方法求解了卡门涡街与翼型的相互作用。较详细地分析了它作为偶极子声源的特性,并模拟了在涡翼相互作用过程中尾迹的发生和演化。计算结果表明:在相互作用时,尾迹所吸收的声能可以是正值也可以是负值;升力和声偶极子强度出现了周期性脉动,这种脉动与翼型的振动和噪声有直接的关系。     

Kàrmàn涡街流场的定性分析研究

本文在采用Ｒａｎｋｉｎｅ涡模型和冻结Ｏｓｅｅｎ涡模型消除Ｋａｒｍａｎ涡街流场的速度奇性后,对Ｋａｒｍａｎ涡街流场进行了拓扑结构分析。证明Ｋａｒｍａｎ涡街流场的自由驻点表现为中心和鞍点。随着空间多数的变化,流动结构可具有三种形式。