基于流动显示的翼尖涡不稳定频率测量-薛栋.pdf

《基于流动显示的翼尖涡不稳定频率测量-薛栋.pdf》由会员分享，可在线阅读，更多相关《基于流动显示的翼尖涡不稳定频率测量-薛栋.pdf（7页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、2016年 4月第42卷第4期北京航空航天大学学报Journal of Beijing University of Aeronautics and AstronauticsApril 2016V0142 No4http：bhxbbuaaeducn jbuaabuaaeducnDOI：1013700jbh1001-596520150594基于流动显示的翼尖涡不稳定频率测量薛栋，潘耕，李广超+(北京航空航天大学航空科学与工程学院流体力学教育部重点实验室，北京100083)摘 要：翼尖涡涡核振荡频率的准确测量是翼尖涡控制技术得以有效实施的重要前提。采用流动显示方法，研究了椭圆机翼翼尖涡在低雷诺数条件

2、下的不稳定特性。分别采用单点谱分析和动力学模态分解技术，从流动显示图像序列中提取了涡核振荡的短波不稳定模态的频率，2种方法得到的频率相对误差最大不超过5。研究结果表明：涡对的空间运动通常展现出长波与短波模态的耦合，涡核的高频短波振荡耦合在低频长波摆动中，以前者为主要含能模态；短波不稳定性的无量纲振荡频率随雷诺数的增大而增大、随机翼攻角的增大而减小。关 键 词：翼尖涡；流动显示；短波不稳定性；频率；动力学模态分解中图分类号：0357；V2117文献标识码：A 文章编号：10015965(2016)04_0837-07固定翼飞机在飞行过程中，其上下翼面存在压力差，使得气流在机翼两侧翼尖处翻卷形成一

3、对反向旋转的翼尖涡，翼尖涡是飞机尾流中主要的相干结构。在无外加干扰的情况下，翼尖涡的强度在100倍机翼展长范围内不会发生明显的衰减，其携带的旋转能量会引起作用范围内后续飞机的飞行速度、高度、航向、倾斜角及其他飞行特性发生显著变化，影响其飞行安全，由此制定的飞机尾流安全标准决定了飞机起降频率、影响机场运营效率。对翼尖涡流场进行深入研究，进而发展促使尾涡失稳耗散的流动控制技术，具有重要的应用价值。近年来，翼尖涡不稳定特性及控制是流体力学研究的热点问题。Crow等旧31最早对反向旋转的飞机尾涡系统进行了系统的研究，提出了长波不稳定性的观点。他们认为若仅依靠黏性耗散，旋涡强度不会发生显著衰减；而在外界

4、扰动的作用下，涡管会发生失稳变形，随着涡管向下游发展，其变形幅度逐渐变大，最终2根相邻涡管对接、破裂并形成涡环，涡强迅速减小，即：涡对经历了“稳定阶段一不稳定线性放大阶段一非线性发展阶段”，最终形成小尺度湍流并被黏性耗散掉。基于Crow等的研究，Widnall等H剖提出了另一种不稳定模态，即短波不稳定性，其波长远小于长波不稳定性，通常由一个涡对另一个涡的应力作用触发而产生，表现为涡核的微幅振荡。Leweke和WilliamsonMl的研究表明，工程实际中，反向旋转涡对的空间形态通常表现为长波和短波不稳定2种模态的相互耦合。不稳定模态通常表现出特定的时间尺度特征。传统的展向脉冲射流1、活动元件振

5、荡。81和合成射流等翼尖涡控制技术，主要是通过在机翼翼尖附近施加周期性扰动，在特定频率和强度下，翼尖涡涡核的空间位置、涡核半径和涡量峰值等会发生显著变化。为保证施加的外部扰动与翼尖涡的不稳定模态发生作用，需要明确知道原始流收稿日期：2015-09-10；录用日期：20151030；网络出版时间：2015-1119 10：36网络出版地址：WWWcnkinetkcmsdetaiL112625V201511191036010html基金项目：航空科学基金(2013ZC51030)通讯作者：Tel：01082318069 Email：liguangchaobuaaeducn引用格式：薛栋，潘种，李广

6、超基于流动显示的翼尖涡不稳定频率测量J北京航空航天大学学报，2016，42“)：837843XUE DPAN CLI G CFrequency measurement of wingtip vortex instability by flow visualization ej 3Journal of Beijing Universit),of Aeronautics and Astronautics，2016，42(4)：837-843(in Chinese)万方数据838 北京航空航天大学学报 2016年场中涡对的不稳定时间尺度，即不稳定模态频率。基于上述原因，本文主要通过定量化流动显示技术

7、研究翼尖涡在低雷诺数下的不稳定特性。基于翼尖涡流场灰度信息，首先采用模式识别方法追踪涡核中心位置的时间演化；其次使用动力学模态分解(Dynamic Mode Decomposition，DMD)方法直接从流动显示灰度场中提取涡核的振荡频率信息，并与模式识别方法进行对比验证。本文工作将为后续翼尖涡不稳定演化的研究提供必备的研究手段。1 实验技术和数据后处理方法11实验技术实验在北京航空航天大学低速回流式水洞中进行，该水洞实验段长12m，横截面为12ml m的长方形。水洞自由来流速度可在050 cms的范围内无级调节，实验段湍流度不超过08。该水洞进行过大量条带稳定性方面的实验研究，流场品质满足旋

8、涡稳定性实验要求。9。实验模型为三维椭圆机翼，如图1所示，机翼截面为低雷诺数翼型E387，展长b=6 cm，翼根弦长c=35 cm。机翼模型由3D打印制作，模型表面进行了抛光处理。实验装置示意图如图2所示，其中石、Y和z分别代表流向、法向和展向方向。以机翼模型根弦长后缘作为坐标原点。模型通过连接杆固连在攻角装置上，机翼攻角可在一20。一20。范围内调整。图l 椭圆机翼模趔Fig1 Elliptical wing n1图2实验装置示意图Fig2 Schematic of experimental setup传统的流动显示方法，例如染色液方法，需要在模型表面开孔，或者引入导管，不可避免地会对流场引

9、入额外的扰动。为避免外界扰动对旋涡稳定性的影响，本实验使用了染色剂表面涂层方法来定性显示翼尖涡涡核的形态，具体方法是：将少量罗丹明染料和芦荟胶混合均匀，然后均匀涂抹于机翼模型翼尖处；模型置于水洞中后，胶状混合物逐渐溶解于水，染料被释放到周围流体中。由于旋涡的卷携作用，染料被逐渐卷入涡核中心，从而显示出涡核附近的结构。需要说明的是，表面涂层非常薄，厚度小于01 am，在低雷诺数条件下涂层并不影响机翼表面的流动。涂层可稳定持续释放荧光染料超过20 min，足以保证实验图像的采集。流动显示的照明光源为半导体激光器，提供波长为532 nm的连续激光，输出能量为8 W，通过透镜组产生厚度约为1 mm的平

10、面片光，垂直来流方向照亮实验段横截面。片光下游60 cm处竖直安装一个45。倾斜的光学平面反射镜。分辨率为2456像素2 058像素的CCD相机(相机型号IMPERX B2520)通过平面反射镜反射记录激光片光照亮的流动图像，相机配备45 mm尼康定焦镜头。图像采样频率为14 Hz。根据不同来流速度，曝光时间间隔在210 ms范围内调整。相机拍摄的物理平面的视野范围为9 cm12 cm(y-z平面)，图像空间分辨率约为21像素mm。图3为采用上述流动显示方法得到的垂直于旋涡轴线的翼尖涡涡对图像。从图中可以看出，荧光染料逐渐卷入旋涡中心，轨迹清晰，涡核中心为高浓度的染料，涡对中问区域为机翼后缘脱

11、落的剪切层，翼尖涡保持层流状态。幽3反向旋转翼尖涡对流动显不Fig3 Flow visualization of counterrotatingwingtip voaex pair研究类似翼尖涡这种拉格朗日流动结构，通常有2种观测方法：一种是拖曳实验；另一种是在固定截面定点观测。相比拖曳实验，定点观测有以下3点优势：可针对涡对运动进行固定站位的长时间统计观测；可避免拖曳实验中拖车启动停止对流场的扰动；由于涡对的诱导作用，翼万方数据第4期 薛栋，等：基于流动显示的翼尖涡不稳定频率测量839尖涡卷起后会以速度w=F2wb。(J1为翼尖涡环量，b。为涡核中心间距)下沉，进行拖曳实验需要保证足够的视野

12、范围，这对实验观测设备提出很高的要求。当前实验采用定点观测方法对垂直于流向的翼尖涡切面进行长时间观测统计。在固定流向站位观测时，涡对的平均法向高度并不会改变，由于涡对自身的不稳定性，涡对仅会在平均的涡核中心空间位置附近波动。观测平面布置在距离机翼后缘50 CB处，总共进行了4个攻角(2。，4。，6。，8。)、2个流速(5 cms，10 cms)的实验观测，基于翼根弦长和自由来流速度的特征雷诺数Re一1 750，3 500。12模式识别技术模式识别技术是辨识并追踪流动相干结构的一种常用数据后处理方法，主要的应用包括边界层中展向涡识别。和圆柱尾迹中卡门涡街的追踪121等。模式识别技术通常针对某一类

13、物理场数据，例如速度场、涡量场及灰度场数据。首先需要定义或者构建数学模型对需要识别的流动结构做初始估计，然后基于互相关计算在相邻2帧数据之间进行匹配分析，由此可以得到流动结构的空间位置。本文基于流动显示的翼尖涡灰度场数据，利用模式识别技术连续追踪翼尖涡涡核中心在不同时刻的空间位置，具体操作如下：首先，在第1帧图像上，选定大小为100像素100像素并包含翼尖涡涡核中心的矩形区域作为初始模型，如图4中小正方形框包含的区域；然后，在下一帧图像中进行互相关匹配，找到初始模型在下一帧图像上的空间位置，以此类推，通过在连续2帧图像间进行匹配，可以得到翼尖涡涡核区域在不同时刻的空间位置；最后，通过计算涡核区

14、域灰度值权重中心，可以得到涡核中心的空间位置坐标，如图4中“+”所示。需要说明的是，为了节省计算图4模式识别技术追踪翼尖涡涡核中心Fig4 Illustration of pattern recognition technique fortracking wingtip vortex core center时间，根据翼尖涡振荡的幅值，人工设定一更大的区域作为互相关计算区域，如图中点划线正方形框所示，而不是在整个图像上进行互相关匹配。由于涡对的形态随时间变化相对比较稳定，相邻2帧之间无量纲化的互相关系数均大于09，基于灰度场的涡核中心识别方法被证明是可行的。13 DMD技术DMD技术是Schmi

15、d等3“。最近提出的一种针对流体动力学系统的模态分解技术。DMD技术基于瞬时流场数据，能够分离、提取复杂流场中暗含的动力学信息而无需知道流动控制方程“，这对于理解复杂流动机理提供了极大的便利。以速度场为例简单介绍DMD原理。更多有关DMD分解原理可参见文献14一15。DMD技术通过一系列矩阵运算将时空变化的速度项H(戈，t)分解成2个独立项，其中一项与时间变化有关，即时间系数项(a，(t)，另外一项与空间分布(咖。(z)有关，即H(戈，t)=：ai(t)咖。(戈) (1)置DMD方法通过时空解耦，将流场分解为不同的空间模态以及时间系数，每一模态对应的时间项暗含模态频率、增长(衰减)率以及能量信

16、息，因此，DMD可以定量地刻画流场中存在的不同尺度流动结构的增长衰减及频率特性。DMD方法已被大量应用到复杂流场的定量分析中，如Gurney襟翼尾迹涡脱落频率测量5|、高速列车尾流频谱分析钊和边界层旁路转捩扰动分析3等。本文中，将DMD技术应用到翼尖涡图像的灰度场数据中，提取翼尖涡涡核振荡的频率信息。2实验结果与讨论21 基于模式识别的不稳定频率辨识采用11节介绍的流动显示方法，对垂直于流向的翼尖涡截面图像进行连续纪录，纪录时长为72 S，采集的总的图像序列为1 000帧，图3为一个典型的翼尖涡对图像。基于翼尖涡灰度图像序列，采用模式识别技术提取涡核中心在y-z平面的空间位置。图5(a)、图5

17、(b)给出了攻角d=8。、自由来流速度V。=5 cms工况下，左右涡核中心在展向和法向方向的空间位置随时间的变化过程(纵坐标的长度单位采用机翼展长无量纲化)。从图中可以看出，有2种明显的不稳定波动，即长波不稳定和短波不稳定，同时存在于翼尖万方数据840 北京航空航天大学学报涡的运动模态中。Breitsamter。按照不稳定波长来划分不稳定模态，长波：Ab。2,rr，短波：叮T2Ab。=4，A为波长。高频的涡核振荡属于短波不稳定性，低频长周期的涡对摆动属于长波不稳定性。与Leweke和Williamson1类似，当前实验中同样观测到长波与短波不稳定模态的耦合(如图5所示)，高频的涡核振荡耦合在低

18、频涡对摆动中。相比于短波，长波的波动振幅明显更大，振幅大约为01b。在展向方向(z方向)，左右涡对的长波摆动近似呈反对称；而在法向方向(Y方向)则呈近似对称。翼尖涡涡核中心平均间距为072b。理论计算给出的椭圆机翼翼尖涡涡心间距为db=叮T41(d为涡对间距)，与实验测量结果存在一定差异，实验测量结果略小于理论值。ts(b)涡_=：中心法向分布图5翼尖涡涡核中心空间位置分布Fig5 Spatial position distribution of wingtipVOrtex core center当前实验重点关注短波不稳定特性。通过计算每一瞬时涡核中心展向间距d：和法向间距d，在原有数据点上减

19、去展向和法向半间距，剔除掉长波的影响，得到短波不稳定模态下涡核中心的空问位置分布，如图6所示，短波的振幅最大不超过机翼展长的4。图7给出了涡核中心短波振荡的放大图。在展向方向，如图7(a)所示，涡核中心存在约百r2的相位差；而在法向方向，如图7(b)所示，涡核中心几乎同相位。008景0-008O 10 20 30 40 50 60 70tfs(a)涡核中心展向分布“S(b)涡核中心法向分布图6翼尖涡涡核中心短波不稳定性空间位置分布Fig6 Spatial position distribution of wing-tipvortex core center for shortwave inst

20、ability图7翼尖涡涡核中心短波不稳定性放大图Fig7 Zoomin view for shortwave instability ofwingtip vortex core center从图6、图7可以看出，左右翼尖涡涡核中心振幅并不完全对称，且振幅随时间发生变化。这样的非对称性在实验观测中普遍存在。通常来讲，对于这种类似机翼翼尖涡的稳定性实验，测量的主要误差来源有以下3个方面：自由来流中的扰动、机翼支撑装置的扰动及机翼后缘分离剪切层的扰动，这样的误差在当前实验中不可避免。尽管如此，实验并未观测到转捩等现象，翼尖涡保持良好的层流状态，而且上述的非对称性并不影响涡核振荡频率和相位的测量。对

21、短波模式下翼尖涡涡核中心空间位置进行谱分析，可得涡核中心短波振荡的不稳定频率。图8(a)、图8(b)分别给出了左右翼尖涡展向和法向短波不稳定性的能谱特性。如图8所示，展向方向上，左右翼尖涡振荡的特征频率均为0952 Hz；法向方向为0966 Hz。此外，均出现能量集中的高次谐波。22基于DMD的不稳定频率辨识21节通过模式识别技术追踪每一时刻翼尖涡涡核中心的空间位置，然后分离出短波不稳定模态，进而进行快速傅里叶变换(FFT)，得到的翼尖涡特征频率精确，是一种单点的局部分析方法，但是实现起来比较复杂。相比于单点分析方法，DMD技术是一种全局特征辨识与分离方法。本文将DMD技术直接应用到灰度图像的

22、分析上。DMD技术将瞬时场数据分解为不同的模态，并且OOO万方数据第4期 薛栋，等：基于流动显示的翼尖涡不稳定频率测量 84l可以得到不同模态的频谱特性。图9给出了相同工况下翼尖涡涡核中心振荡的DMD能谱特性。横坐标，代表某一模态的特征频率，纵坐标l|a。lI是该模态时间系数的范数，代表这一模态所具有的能量。如图9所示，能量峰值所对应的主频为正=0924 Hz，即涡核振荡的短波频率。其他的高频能量峰值对应2倍频和3倍频等高阶谐波一丑删趟I 左翼涡弋孓N 。愀频率Hz(a)展向能谱图8涡核中心空间位置坐标能谱Fig8 Energy spectrum of spatial displacement

23、 ofvortex eore Center，=0924HzC C00吣嗵 0汹图9 DMD能谱Fig9 DMD energy spectrum对比图8和图9，DMD得到的主频与单点谱分析得到的主频存在一定误差，误差的主要来源为本实验的图像采集频率过低。在对实验数据进行DMD分解之前，Pan等副对DMD算法进行了系统的误差分析，研究结果表明，DMD对数据分解的精度依赖于输入的数据质量，在一定的范围内，数据的采集频率越高，误差越小。对于当前实验，由于CCD相机采集频率所限，本实验最大采样频率为14 Hz，但是DMD分解得到的频率已经非常接近单点分析的结果，相对误差最大不超过5。23翼尖涡不稳定频率

24、的影响规律采用上述单点谱分析方法和DMD技术，分别计算了4个攻角(2。，4。，6。，80)、2个自由来流速度(V。=5，10 cms)下翼尖涡涡核振荡的无量纲频率，表1和图10分别给出了2种方法的计算结果(单点谱分析方法给出的是左右翼尖涡法向无量纲频率)，结果表明无量纲频率随雷诺数增大而增大；随机翼攻角增大而减小，且近似呈线性减小。表1 翼尖涡涡核振荡无量纲频率Table 1 Non-dimensional frequency of wingtipvortex core oscillation图10翼尖涡涡核中心振荡无量纲频率对比Fig1 0 Comparison of nondimensio

25、nal frequency ofwingtip vortex core oscillation3 结 论本文通过水洞实验，采用流动显示技术对垂直于流向的翼尖涡截面流场进行了观测，研究了翼尖涡涡对的不稳定频率特性；分别采用模式识别方法和动力学模态分解技术提取了涡核中心振荡的特征频率。研究结果表明：1)在低雷诺数下，翼尖涡涡对的不稳定运动通常表现出短波与长波2种模态，高频的涡核振荡耦合在低频涡对摆动中。2)在当前的实验条件下，产生了层流翼尖涡对，通过模式识别方法识别出涡核中心进而得到涡核不稳定振荡的频率。将动力学模态分解技术应用于流动显示图像，也可提取出不稳定频率信息，与单点分析方法相比，DMD方

26、法具有更强的万方数据842 北京航空航天大学学报 2016拄适用性，且精度有保证。3)从现有的8个实验工况来看，涡核振荡的无量纲频率随机翼攻角增大而减小，随雷诺数增大而增大。参考文献(References)1BREITSAMTER CWake vortex characteristics of transport air-craftJProgress in Aerospace Sciences，2011，47(2)：89-1342CROW S CStability theory for 8 pair of trailing vorticesJAIAA Journal，1970，8(12)：21

27、72-21793CROW S C，BATE E RLifespan of trailing vortices in a turbu-lent atmosphereJJournal of Aircraft，1976，13(7)：476-4824WIDNALL S E，BLISS D，TSAI C YInstability of short-waveson a vortex ringJJournal of Fluid Mechanics，1974，66(Part1)：35-475WIDNALL S E，BLISS D，ZALAY ATheoretical and experi-mental stu

28、dy of the stability of a vortex pairMAircraftwake turbulence and its detectionBoston：Springer，1971：3053386LEWEKE T，WILLIAMSON C H KExperiments on longwavelength instability and reconnection of a vortex pairJPhysicsof Fluids，2011，23(2)：465-4747HEYES A L，SMITH D A RSpatial perturbation of a wing-tipvo

29、rtex using pulsed span-wise jetsJExperiments in Fluids，2004，37(1)：1201278BREITSAMTER C，ALLEN ATransport aircraft wake influenced by oscillating winglet flapsJJournal of Aircraft，2009，46(1)：1751889GUO H，BORODULIN V I，KACHANOV Y S，et a1Nature ofsweep and ejection events in transitional and turbulent b

30、oundarylayersJJournal of Turbulence，2010(11)：N3410HU Y，WANG J JDual leading-edge vortex structure for flowover a simplified butterfly modelJExperiments in Fluids，2011，50(5)：1285-12921 1SCARANO F，BENOCCI C，RIETHMULLER M LPattern rec-ognition analysis of the turbulent flow past a backward facing1213st

31、epJPhysics of Fluids，1999，11(12)：3808-3818GIRALT F，FERRI；J AStructure and flow patterns in turbulent wakesJPhysics of Fluids A：Fluid Dynamics，1993，5(7)：1783-1789SCHMID P J，SESTERHENN JDynamic mode decompositionof numerical and experimental dataJJournal of Fluid Mechanics，2010，656(10)：5-28SCHMID P JA

32、pplication of the dynamic mode decompositionto experimental dataJExperiments in Fluids，2011，50(4)：11231130PAN C，YU DWANG JDynamical mode decomposition ofGurney flap wake flowJTheoretical and Applied MechanicsLetters，2011，1(1)：012002MULD T W，EFRAIMSSON G，HENNINGSON D SFlowstructures around a high-spe

33、ed train extracted using proper or。thogonal decomposition and dynamic mode decompositionJComputersFluids，2012，57(4)：87-97HE G S，PAN C，WANG J JDynamics of vortical structures incylinderwall interaction with moderate gap ratioJJournal ofFluids and Structures，2013，43(7)：100-109PAN C，XUE D，WANG J JOn th

34、e accuracy of dynamic modedecomposition in estimating in stability of wave packetJExperiment in Fluids，2015，56(8)：115作者简介：薛栋男，硕士研究生。主要研究方向：实验流体力学。Tel：01082318069Email：xuedongl202buaaeducn潘种男，博士，副教授。主要研究方向：湍流、转捩，流动不稳定性。Tel：01082318069E-mail：panchongbuaaedu cn李广超男，博士，讲师。主要研究方向：流动控制，飞机防、除冰。Tel：01082318

35、069E-mail：liguangehaobuaaeducn】HM墙万方数据第4期 薛栋，等：基于流动显示的翼尖涡不稳定频率测量843Frequency measurement of wing-tip vortex instability by flow visualizationXUE Dong，PAN Chong，LI Guangchao+(Fluid Mechanics Key Laboratory of Education Ministry，School of Aeronautic Science and EngineeringBeijing University of Aeronau

36、tics and Astronautics，Beijing 100083，China)Abstract：The accurate measurement of the dominant frequency of wingtip vortex core oscillation isessential for controlling the wingtip vortexIn the present study，the instability of wing-tip vortex，which isgenerated from an elliptical wing，was investigated u

37、nder low Reynolds number via flow visualization techniquePointwise spectrum method and dynamic mode decomposition were applied to extract the dominantfrequency of the shortwave instability of the wing-tip vortex core from the flow visualization image sequence，and the maximum relative error of these

38、two methods is less than 5The results show that shortand longwave instability modes are developed simultaneously among the vortex pair；the highfrequency vortex coreoscillation is coupled with the slow sidetoside movement of vortex tube，and the former is the main energymode；the nondimensional frequen

39、cy of the shortwave instability increases with the increase of Reynoldsnumber，and decreases with the angle of attack increasingKey words：wingtip vortex；flow visualization；short-wave instability；frequency；dynamic modedecompositionReceived：2015-0910；Accepted：2015-10-30；Published online：2015-1119 10：36URL：WWWcnkinetkcmsdetail112625V201511191036010htmlFoundation item：Aeronautical Science Foundation of China(2013ZC51030)十Corresponding authorTel：0 1 0823 1 8069 E-mail：liguangchaobuaaeducn万方数据