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

针对安装了激振器的NACA0012二维翼型,采用数值模拟的方式对零质量射流技术进行研究.流场解算器采用了中国航空工业空气动力研究院计算所自主研发的ENSMB V1.2计算平台.计算结果给出了湍流模型对零质量射流数值模拟技术的影响,并且分析了激振器参数对翼型气动力特性的影响规律.     

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

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

基于零质量射流技术的旋翼翼型流场主动控制研究

旋翼翼型经常工作在大迎角来流条件下,流动分离普遍存在。通过在具有分离流动的翼型表面施加零质量射流控制,可以实现对流动分离的抑制,达到增升减阻、推迟失速的目的。本文通过对某旋翼翼型的流动控制的计算,采用数值模拟方法研究了零质量射流主动流动控制机理,并对影响零质量射流控制效果的射流喷角、动量系数和无量纲频率三个控制参数进行了分析比较,得出了一些有意义的结论。     

零净质量射流的数值模拟

对零净质量射流致动器的流场和有零净质量射流时翼型的绕流流场进行了数值模拟,分析了动网格与几何守恒率、边界条件、湍流模型对零净质量射流致动器流场计算的影响,分析了零净质量射流的速度幅值和驱动频率对翼型增升效果的影响.研究结果表明,在复杂计算中可以对能动部件的边界作赋值处理;随着射流速度幅值的增加,翼型的平均升力系数和阻力系数都要增加;射流频率对升力的影响呈非线性.     

零质量射流推迟翼型失速的数值模拟

 通过在具有分离流动的翼型表面局部地引入零质量射流干扰，可以较小的能量代价实现对翼型绕流的全局性主动控制，从而达到抑制分离、推迟翼型失速的目的。本文发展了零质量射流与翼型绕流干扰的数值分析方法。在格心有限体积法的基础上，空间离散采用AUSM+up格式，时间推进运用含牛顿型LU-SGS子迭代的全隐式双时间法，且引入了预处理方法和多重网格法等技术。以TAU0015翼型失速主动控制为算例，通过在翼型上表面前缘附近引入零质量射流控制，使翼型失速特性得到改善。将计算结果与实验结果做了比较，验证了本文发展的数值模拟方法及程序的正确性。还以NACA0015翼型为例研究了射流动量系数、射流频率、射流偏角等关键参数对失速控制效果的影响，得出一些关于如何选择控制参数的结论。     

翼型风洞侧壁干扰的数值模拟研究

运用Navier-Stokes数值模拟对翼型模拟试验对风洞侧壁干扰进行模拟，将简单代数湍流模型扩展用于机翼／风洞侧壁结合区流动，分析了风洞实验侧壁干扰问题的形成机理和影响翼型实验侧壁干扰的各种因素，如翼型展长、风洞侧壁边界层厚度及侧壁边界层抽吸等，对实验结果的影响，得出了一些有用的结论。计算格式空间采用中心有限体积离散，时间采用多步Runge-Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

薄翼型阵风响应的数值模拟

 采用有限体积法来数值求解非定常Euler和N-S方程，并在非定常计算中通过引入“网格速度”来计入阵风的影响，实现了对刚性薄翼型阵风响应的有效数值模拟。采用该方法来对薄翼型迎角突然增大的阵风、陡边型阵风这两种典型阵风情况下的翼型气动力响应进行了计算，结果表明：计算结果与现有的一些理论结果都符合较好，与文献中计算结果也吻合。同时通过计算也表明厚翼型的阵风响应需根据工况开展进一步的计算分析。     

翼型表面粗糙度对结冰冰型的影响研究

本文采用数值模拟方法对NACA0012翼型前缘结冰进行模拟，在同一大气环境中，对不同表面粗糙度情况下的翼型进行结冰模拟计算，并将计算所得的各条件下的冰型进行对比分析，从而研究机翼表面粗糙度对结冰的影响。     

风力机叶片翼型的结冰数值模拟研究

研究了风力机叶片翼型的结冰数值模拟方法,及翼型结冰后对其气动性能的影响。求解雷诺平均N-S方程,引入k-ωSST湍流模型封闭方程,获得风力机叶片流场;采用拉格朗日法计算风力机叶片翼型周围的水滴运动轨迹,同时考虑了多尺寸水滴分布的影响,获得翼型表面的局部水收集系数分布;根据质量守恒和能量守恒原理,计算翼型表面的各项热流,获得翼型表面的结冰速率和结冰冰形;考虑到翼型结冰是时间的动态函数,采用多时间步长法完成结冰数值模拟。计算了风力机叶片翼型在不同环境条件和气象条件下的结冰冰形,同时模拟了风力机翼型结冰后周围流场的变化,并与干净翼型的气动特性进行了对比。结果表明,环境温度较高时形成的明冰对翼型气动性能的影响较大,结冰导致翼型升力下降,升阻比减小,最大减小幅度达到61%,同时结冰后的翼型会提前进入失速区,导致桨叶气动性能恶化。     

Gurney Flap绕流气动特性研究

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

OA212翼型主动流动控制的数值模拟研究

采用数值模拟的方法,探讨了基于零质量射流的主动流动控制技术对OA212旋翼翼型动态失速的控制效果和控制特性.以积分形式雷诺平均Navier-Stokes(N-S)方程为控制方程,采用格心有限体积法进行求解.空间离散采用AUSM~+-up格式,时间推进采用含牛顿型LU-SGS子迭代的全隐式双时间法,且引入了预处理方法和多重网格方法加速收敛.通过在喷口上施加非定常边界条件来模拟射流对翼型绕流的影响.研究了不同类型射流、不同位置射流以及不同控制参数(频率、相位、偏角、动量系数等)对动态失速控制效果的影响.研究表明:零质量射流和传统的定常射流均可减小动态失速迟滞环的回线面积,但在提高最大升力方面零质量射流明显优于定常射流;在12%c和62%c处施加组合零质量射流的控制效果最为明显.     

关于跨音速零升力翼型绕流的松弛和两种自修正风洞的收敛性

本文采用横向小扰动而纵向大扰动速势方程,计算了跨音速零升力翼型的绕流。在线松弛的数值实验中,φ_γ的差分式用简单迭代和φ_(xx)的差分式用改进迭代时,稳定性较好。此结论与文献的线化理论分析相符。 本文用混合差分法数值模拟,证明了基于两个控制面上的静压和基于一个控制面及翼面上的静压的跨音速零升力翼型自修正风洞的收敛性。对前一种方案,NACA0012翼型,M∞=0.9,RAE104翼型,M∞=0.8,对后一种方案,NACA0012翼型,M∞=0.72,0.8,在迎角为零和风洞高度与翼弦之比为3时,均能收敛到无洞壁干扰的自由流。     

翼型和机翼阵风响应的数值模拟

采用有限体积法来数值求解非定常Euler和NS方程,并在非定常计算中通过引入"网格速度"来计入阵风的影响,实现了对刚性翼型和机翼阵风响应的有效数值模拟。首先采用该方法来对翼型攻角突然增大的阵风、陡边型阵风这两种典型阵风情况下的翼型气动力响应进行了计算,结果表明:计算结果与现有的一些理论结果都符合较好,与文献中计算结果也吻合,NS方程的计算结果与Euler方程计算结果在变化趋势上是一致的,但具体的数值有差别。最后针对攻角突然增大的阵风,对不同马赫数下不同展弦比平直机翼的阵风响应进行了数值模拟,得出的计算结果与文献中的计算结果一致。     

二元翼型结冰数值模拟研究

对二元翼型前缘结冰数值模拟进行了研究。采用中心有限体积法求解N—S方程；四阶龙格一库塔法求解水滴运动轨迹方程；逐次二分法寻找水滴撞击极限，然后计算局部水滴收集率；分析翼型表面的传质、传热项，建立相应质量、能量方程并迭代求解，得到翼型表面结冰质量；根据角度二分法生成结冰冰形。预测的NACA0012结冰冰形与实验及Lewice模拟结果吻合很好，结冰过程中的特征参数与冰型及相应冰形特征吻合。     

电磁力控制翼型体绕流的数值研究

为考察包覆电磁激活板的翼型体绕流特性,采用数值模拟的方法研究了低雷诺数下电磁力对翼型体流场结构和绕流过程的影响,探讨了不同迎角和不同电磁力大小对翼型体尾流涡街形态和流体边界层结构的影响。计算结果表明,作用于翼型体表面的电磁力可以明显改变流体边界层的结构,抑制边界层分离,并有效提高翼型体的升力。     

多段翼型增升小翼纳维-斯托克斯方程分析

装有扁平板增升小翼的多段翼型流场数值分析已经完成。研究的小翼装在接近主翼段后缘，高度范围是0．25％～1．25％参考翼型弦长。二维数值模拟在嵌入式结构拓扑网格上采用不可压缩纳维-斯托克斯（N-S）方程求解分析。新的网格细化提高了网格边界重叠部位的计算精度。在固定雷诺数条件下，对装有窄缝襟翼的二段翼型研究了各种小翼的影响。计算和实验结果都表明，相对未装小翼的基本翼型，主翼段凹陷处装有小翼提高了最大升力和升阻比。计算得出的流线分布揭示，正是小翼的存在引起了襟翼流场分离的减弱。也研究了三段翼型在一个雷诺数范围内的流场情况。对优化的襟翼来说，计算和实验所得的雷诺数影响结果一致。当襟翼从最有利的位置移动时，数值分析表明，小翼有助于翼型重新优化，大概能起到所优化襟翼效益的1％量级的作用。     

转捩对风力机翼型和叶片失速特性影响的数值模拟

采用基于k-ωSST湍流模型的Gamma-Theta转捩模型对S809翼型和NREL Phase VI叶片进行了考虑转捩的气动力数值模拟,研究了转捩对其气动特性特别是失速特性的影响。首先对S809翼型在迎角0°～30°范围内开展了数值模拟,比较了转捩模拟和全湍流模拟获得的翼型升力系数和流场特征,发现前缘层流分离泡的存在显著影响了翼型的失速特性;然后对NREL Phase VI叶片开展了类似的数值模拟,结果显示转捩对叶片失速特性和翼型失速特性产生影响的作用方式是相似的。经过对以上数值模拟结果加以分析认为,转捩对翼型和叶片失速特性的影响主要通过前缘层流分离泡的作用体现出来,前缘层流分离泡的存在使翼型和叶片更早地进入深失速。     

后缘隔板对充气式机翼翼型气动力影响分析

以充气式机翼为背景,研究了安装后缘隔板充气式翼型的气动特性。以翼型E387为基础生成的钝后缘翼型为研究对象,采用大涡模拟方法对安装中部隔板的充气式翼型进行了流场数值计算,并与未安装隔板的情况进行了对比。总结了后缘中部隔板对充气式翼型气动特性影响的变化规律,为充气式机翼设计提供参考。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构，层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟，研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型，基于结构化拼接网格，对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性，雷诺数对翼型气动特性具有显著影响。随雷诺数降低，时均分离泡外形增大、位置后移，平均阻力系数增大，特别是在较低雷诺数下，翼型升阻力系数随时间出现振荡现象。进一步研究表明，造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低，流动黏性增大，导致分离剪切层速度梯度减小，流动发生转捩及再附位置后移，直至翼型表面不再发生转捩和再附。     

基于Euler方程的二维高速海面效应数值分析

以近海面飞行器的海面效应为研究背景,通过数值求解可压缩Euler方程模拟了波形壁对翼型高亚声速气动力特性的影响.计算网格布局采用了重叠网格方法,避免了翼型与波形壁之间网格生成的困难.详细分析了波形壁边界条件,给出了正确的通量计算公式.时间推进方法采用了双时间方法,子迭代过程由LU-SGS方法完成.通过数值分析,给出了飞行高度、波长和波幅等因素对翼型气动力特性的影响.