考虑地面效应的翼型动态特性数值模拟

为了研究考虑地面干扰的飞行器非定常气动特性,基于滑移网格技术,通过求解二维非定常Euler方程,对不同近地高度下的NACA0012翼型的动态振动特性进行了数值模拟。分析了存在地面干扰时的升力系数以及俯仰力矩系数迟滞环变化,并结合强迫振动法,计算了翼型各个近地高度下的纵向组合动导数。结果表明,地面效应不仅对定常流场产生影响,更显著地影响了非定常气动力及力矩,近地高度越小,这种效应越明显,使得升力系数的迟滞环面积变小,幅值减小,而俯仰力矩系数的迟滞环变化不规律,力矩系数的动导数随着近地高度减小而增大(绝对值减小),纵向的动态稳定性产生损失。地面效应干扰在飞行器非定常气动研究中应该引起重视。     

基于谐振舵面的跨声速抖振抑制探究

跨声速抖振引起的非定常脉动载荷会造成飞行器结构疲劳甚至引发飞行事故,所以跨声速抖振的控制研究逐渐成为航空领域的热点。采用基于Spalart-Allmaras(S-A)湍流模型的非定常雷诺平均方程开展了基于谐振舵面的跨声速抖振抑制研究。首先验证静止NACA0012翼型的抖振边界和频率特性,然后分别从舵偏平衡位置、舵偏幅值、频率以及相角等角度研究了谐振舵面的控制效果。舵偏平衡位置等效于减小了翼型的有效迎角;幅值和频率对抖振抑制效果影响较大,当舵面振荡频率与抖振频率接近时发生共振现象;相角对控制效果有一定影响,在270°相角附近,升力系数幅值减小了60%。在合适的舵偏幅值、频率以及相角组合下,谐振舵面有可能成为跨声速抖振的有效开环控制策略。     

旋翼翼型非定常动态失速响应的计算

基于旋翼非定常翼型气动模型,给出了计算分离流和深度失速状态下的翼型非定常升力、俯仰力矩的数值计算方法。该方法采用半经验指数响应公式,利用数值离散方法来求解翼型的非定常法向力和俯仰力矩。分别计算了NACA0012和SC-1095翼型上的非定常气动载荷,并与可得到的试验结果进行了对比,验证了方法的有效性。文中还讨论了缩减频率和马赫数对动态失速响应的影响;然后,这个模型被改进以适用于后掠流下的翼型动态失速响应计算,分析了后掠角对翼型动态失速响应的影响。最后,得出了一些结论。     

零质量射流对翼型增升影响的数值模拟研究（上）

针时安装了激振器的NACA0012二维翼型，采用数值模拟的方式时零质量射流技术进行研究。流场解算器采用了我院计算所自主研发的ENSMBV1．2计算平台。添加了零质量射流模块。计算结果给出了激振器的存在时翼型升力、阻力、俯仰力矩系数的影响。并且刘划了流场图谱。分析了其流场结构和流动机理。     

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

翼型尾缘低频大功率合成射流的动态载荷特性研究——B：动态气动载荷特性分析

针对NACA0015翼型,设计了适用于风洞研究的尾缘低频大功率合成射流致动器,对翼型尾缘合成射流作用下的非定常气动特性进行了风洞实验研究。研究表明：尾缘合成射流与横流的相互作用能够有效改变作用在翼型上的气动载荷;单侧喷口喷/吸时,喷冲程气动力系数响应幅值约为吸冲程幅值的3倍;双侧喷口同时工作时,升力和力矩系数的幅值并不是单侧喷口单独工作时喷气幅值的简单叠加,而是处于单侧喷气幅值和喷吸幅值和之间;升力和力矩响应的幅值与喷流动量系数的平方根之间存在近似的线性关系;在动量系数不变时,升力和力矩系数响应的幅值会随减缩频率的增加而减小;给定合成射流器行程,升力和力矩响应幅值与合成射流频率之间近似呈线性关系。     

超临界翼型Gurney襟翼增升技术实验研究

通过在二元翼型风洞中进行测力实验,研究了不同高度Gurney襟翼对超临界翼型气动力和力矩的影响规律。实验结果表明：在亚声速条件下,Gurney襟翼同样可以明显增加翼型的升力系数,使整个升力曲线向上平移,并使翼型低头力矩增加。高度为翼型弦长0.5%的Gurney襟翼可以带来超临界翼型的最大升阻比。同Gur-ney襟翼对NACA 0012翼型气动特性改变的对比表明,其在超临界翼型上带来的升力系数增量要大于在NACA0012翼型上的效果,但是带来的低头力矩增量较小。     

基于CFD方法的机翼阵风响应研究

通过引入＂网格速度＂方法模拟阵风条件,求解非定常欧拉方程实现了不同展弦比0012平直机翼阵风响应的数值模拟。首先采用该方法对NACA0006翼型迎角阶跃型阵风的气动力响应进行了计算,计算结果与文献结果、理论值吻合良好。进一步对展弦比分别为5、10的0012平直机翼在迎角阶跃型、One-Minus-Cosine型阵风作用下的气动力响应过程进行了模拟分析。研究结果表明,当展弦比增大时,阵风作用下机翼的升力系数响应会增加,机翼翼根部位气动响应幅值大于翼尖部位响应特性。     

Gurney襟翼改善翼型动态失速特性研究

采用CFD数值方法,研究了NACA0012翼型、加装传统Gurney襟翼及改进不对称Gurney襟翼后翼型的动态失速特性。给出了传统襟翼对翼型动态失速特性的影响,并带来较大动态低头力矩的不足,基于传统襟翼的不足,提出了改进的不对称Gurney襟翼方案。研究表明,不对称Gurney襟翼可较好改善翼型的动态失速特性,在增加动态升力的同时,俯仰低头力矩明显减小,可能是直升机旋翼的较理想翼型。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

Gurney Flap绕流气动特性研究

采用C-H型多块结构网格和SST k-ω湍流模型求解二维定常雷诺平均NS方程,对NACA4412单段翼型后缘附近和NACA63(2)-215B多段翼型襟翼后缘附近的Gurney Flap绕流进行了数值模拟。讨论了局部网格对计算结果的影响,重点研究了Gurney Flap位置和高度对气动特性的影响。计算结果表明:对于NACA4412单段翼型,Gurney Flap位置对气动性能的影响较小,但随着高度的增加,零升攻角和失速攻角减小,最大升力系数增大,同时阻力系数和力矩系数增加;对于NACA63(2)-215B多段翼型,Gurney Flap对气动性能的影响类似于NACA4412单段翼型情形。     

NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

大振幅下薄翼和流线型箱梁的气动迟滞研究

回顾了机翼和箱梁的气动迟滞效应研究现状,介绍了大振幅下获取薄翼和流线型箱梁自激气动力的风洞试验。研究结果表明：薄翼在超过失速角的振幅条件下,升力并没有出现明显的失速现象,但气动力的高阶谐波分量显著增加;大振幅条件下,流线型箱梁的气动力高阶谐波分量也比较显著,并以第2和第3阶谐波分量为主。此外,在大振幅条件下,流线型箱梁的气动力矩迟滞曲线可出现“8字环”,即气动力在一个振动周期内既做了正功也做了负功;其中,迟滞曲线“8字环”的气动正功部分随着振幅和折算风速的增大而增加。该现象可导致桥梁也出现类似于失速机翼的极限环震荡。最后,基于不同振幅下流线型箱梁的力矩迟滞曲线,简要讨论了大跨度桥梁在颤振后状态可能出现的振动形式和气动稳定性。     

基于遗传算法的旋翼翼型综合气动优化设计

结合直升机飞行特性及旋翼工作气动环境,提出了适用于中型运输直升机旋翼翼型优化设计的目标函数及约束条件.在此基础上采用遗传算法,以"黑鹰"(UH-60A)直升机旋翼翼型SC1095为基础翼型,进行多目标、多状态以及多约束条件下的旋翼翼型优化设计.优化翼型同SC1095翼型相比弯度稍大,最大厚度略有增加.结果表明:优化翼型气动力特性在满足约束条件的情况下,零升力矩系数减小了57.2%,最大升力系数增加了3.7%,阻力系数减小了7.5%.同时,非定常动态失速状态下的气动特性也有一定的改善,阻力系数和力矩系数的发散范围均有明显的减小.此外,悬停状态计算结果表明:采用优化翼型构成的旋翼与SC1095旋翼相比,具有更高的悬停效率及更低的力矩系数.      

无网格法耦合RNGk-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

 基于移动最小二乘无网格方法,耦合RNG (Re-Normalisation Group) k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM (Advection Upstream Splitting Method)⁺-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP (Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNG k-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。