平流层飞艇具有驻留时间长、效费比高、生存能力强和高度适中等特点,作为搭载侦察、预警、通信等电子设备的平台可以发挥其特有优势,在民用和军用领域均有很大现实意义[1]。但平流层飞艇螺旋桨在低雷诺数飞行环境下,叶素边界层转捩、分离等问题造成螺旋桨性能和效率降低[2],提高螺旋桨低雷诺数下飞行效率和宽工况飞行能力已成为重要研究课题。当前国内外对平流层螺旋桨增效的研究主要集中在优化翼型气动外形、采用桨梢小翼等领域,获得了有价值的研究成果,但是螺旋桨增效能力有限[3, 4, 5, 6, 7]。
表面介质阻挡放电(Surface Dielectric Barrier Discharge-SDBD)是等离子体流动控制技术的重要研究方向,能够有效抑制翼型流动分离,发挥减阻增升作用,提高螺旋桨低雷诺数下飞行效率和飞艇宽工况飞行能力;同时具有结构简单、质量轻、功耗低和响应快等特点,已成为流动控制技术的热点研究领域[8, 9]。
目前,在翼型等离子体流动控制方面,主要通过优化激励器几何参数和激励电源参数,探讨翼型在不同的来流速度和迎角以及低雷诺数条件下抑制流动分离的控制效果[10, 11, 12],其中有美国圣母大学和我国南京航空航天大学等科研单位开发了等离子体流动控制的实验无人机。
针对平流层高效螺旋桨,基于螺旋桨叶素理论,国防科技大学雷光新、刘巍等人[13]给出了平流层螺旋桨设计流程;中国电子科技集团二十七研究所姜 波、翟彬等人[14]建立了计算NACA-0015翼型螺旋桨气动性能的模型,并结合实验数据验证了模型的正确性。平流层螺旋桨等离子体流动控制技术方面,程钰锋等人[15, 16]通过数值仿真开展了螺旋桨等离子体流动控制增效的研究,但是当前还没有见到通过实验方法开展平流层螺旋桨等离子体流动控制研究的公开文献,难以推广螺旋桨在平流层的应用。实验开展高空平流层等离子体螺旋桨的研究需要低密度风洞,虽然可以实现,但是存在技术难度高、效率低以及成本高等问题。因此,迫切需要综合采用现有地面实验手段模拟平流层螺旋桨等离子体流动控制,定量评估效果。
作者根据刘沛清等[17, 18]开展的地面常规风洞进行平流层螺旋桨实验的可行性研究,在车学科等人[19]发展的地面条件下模拟高空等离子体流动控制效果的方法基础上,基于螺旋桨叶素理论,采用实验手段开展平流层螺旋桨等离子体流动控制增效的研究。
1 平流层螺旋桨地面实验原理 1.1 平流层螺旋桨叶素运动参数确定叶素理论又称孤立桨叶理论[20],如图 1所示,沿切线方向将桨叶切成若干个称之为叶素的微小段,认为绕过每个叶素的运动气流是二维的,所以叶素之间、桨叶之间互不影响。
图 1中,r为叶素展向中点距桨心距离;R为螺旋桨半径;ξ为叶素相对半径,即ξ=r/R;br为叶素弦长;θr为叶素安装角,dr为叶素微小段宽度。螺旋桨在平流层以转速ns、相对来流前进速度V0运动时,根据运动相对性,其r处叶素运动分析如图 2所示。
在不考虑螺旋桨轴向干涉速度和环向干涉速度的情况下,忽略其干涉合成速度V1影响,几何合成速度W0可以近似代替实际入流速度W1,则几何入流角度Φ0可以近似等于实际入流角度Φ1。根据叶素运动分析,其周向速度U0、几何合成速度W0、几何入流角度Φ0和叶素迎角α可以表示为:
不同半径r处叶素的几何合成速度W0和迎角α分别由Wr、αr表示,根据式(1)~(4)可知,桨叶叶素的运动状态可以由来流速度Wr、迎角αr 2个参数确定。因此,在平流层以转速ns、相对来流速度V0运动的螺旋桨可以转化成不同半径r处迎角为αr的叶素在速度为Wr的来流中的运动。
1.2 地面实验参数确定螺旋桨相同转速在不同高度工作时,叶素绕流的雷诺数也会随之变化。在平流层20km低气压、低密度条件下,螺旋桨雷诺数相当于低空飞行时的1/10,叶素近壁面层流难以克服逆压梯度而容易发生流动分离造成螺旋桨性能下降。文献[18]认为,密度基本不会对螺旋桨的相似参数产生影响。根据螺旋桨雷诺相似准则,地面二元翼型风洞实验模拟平流层各叶素微段流动时对应雷诺数相等,即:
式中:ρh、uh为平流层大气密度、动力粘度;ρg、ug为实验气体密度、动力粘度;b为实验翼型弦长;Wrg为地面实验风速。根据式(5)可得到地面实验风速Wrg,实验翼型迎角与平流层叶素迎角αr相同。
文献[19]认为,根据等离子体射流雷诺相似准则,通过调整激励器参数和激励电源参数,可以由地面SDBD等离子体放电模拟不同高度下放电的流动控制效果。因此,实验基于射流雷诺相似准则,根据高空平流层等离子体射流雷诺数确定地面激励电源频率和电压。
1.3 螺旋桨推力评估通过二元翼型风洞实验系统,将采集的实验数据经过编程计算得到分别垂直和平行于来流方向的升力dL系数CL、阻力dD系数CD,如图 3所示。
根据翼型受力分析可知,叶素沿桨盘竖直向上的分量dT即螺旋桨产生的推力,则dL、dD和dT可表示为:
式中:pr为实验测得来流空气动压。对桨叶叶素推力dT积分可得螺旋桨单桨叶推力T:
式中:r0为螺旋桨桨盘半径。本实验主要探讨施加等离子体流动控制对螺旋桨推力的变化效果,采用计算得出的推力指标定量评估平流层螺旋桨性能的改变情况。
2 实验方案实验模拟的平流层飞艇螺旋桨基于文献[18]的设计,叶素剖面为S1223翼型,螺旋桨半径R= 3250mm。实验在装备学院的低湍流度风洞中进行,湍流度0.02%,满足平流层低湍流度要求,风洞试验段尺寸为0.8m×1.0m×4.0m,实验S1223翼型弦长200mm,展长790mm,采用环氧树脂加工。实验舱段和翼型如图 4所示,翼型上下表面分别布置18和17个测压孔,测压设备采用西北工业大学研制的压力采集系统,由DS-32电子压力扫描器、DS-104微压校准单元和气压 传感器3部分组成,测量误差0.1%。
分别在翼型上表面前缘、下表面后缘处布置1个SDBD激励器,激励器暴露电极和植入电极宽度分别为5mm、10mm,电极长度500mm,重合区长度230mm,电极间隙0mm,两电极间通过粘贴3层Kapton胶带作为介质阻挡层,其等效模拟的平流层20km螺旋桨激励器的暴露电极、植入电极宽度为10mm和60mm,电极间隙0mm[19]。翼型上表面激励器暴露电极靠近翼型前缘,距离前缘5.0mm;下表面激励器暴露电极靠近翼型后缘,距离后缘5.0mm。根据上述激励器铺设方式可知上表面等离子体射流方向与来流方向相同,下表面等离子体射流与来流方向相反[20],翼型上表面1~5号测压孔由于激励器的遮挡无法正常测压。
激励电源采用HFHV30-1高频高压交流电源,输出电压±15kV,输出频率1~50kHz,上、下表面激励器采用并联方式,使得激励器工作状态时刻相同。电压采用安捷伦N2771B高压探头测量,电流采用皮尔森电流线圈6595测量,电压、电流测量结果使用安捷伦DSO3024A示波器进行显示和记录。
实验方案设定平流层20km处飞艇螺旋桨以300和600r/min 2种转速工作,由于平流层平均风速为10m/s[21],同时又要满足飞艇驻空悬停和机动能力,螺旋桨的前进速度V0分别取5、10和20m/s,根据前文所述实验原理可得地面二元翼型风洞试验参数,如表 1所示。
文献[19]通过低压舱平板射流实验,得出平流层20km、激励电压峰-峰值13.6kV时,等离子体沿壁面切向的射流雷诺数为36~40。在本实验中,进一步考虑了等离子体对低气压稀薄气体的压缩作用导致密度的变化,沿壁面射流的雷诺数约为79~112,由于激励电源频率对雷诺数影响很小,实验频率采用10kHz,主要是为了更好模拟平流层激励器放电诱导的流场,实验电压峰-峰值根据射流雷诺数分别为6和7kV。
(a)螺旋桨转速300r/min | ||||||
0.30 | 5.48 | 38.948 | 5.69 | 38.794 | 6.46 | 38.531 |
0.35 | 7.52 | 36.701 | 7.73 | 36.567 | 8.53 | 36.329 |
0.40 | 9.31 | 34.448 | 9.51 | 34.329 | 10.29 | 34.114 |
0.45 | 10.76 | 32.191 | 10.95 | 32.085 | 11.67 | 31.889 |
0.50 | 11.82 | 29.832 | 11.99 | 29.736 | 12.64 | 29.556 |
0.55 | 12.47 | 27.481 | 12.62 | 27.393 | 13.18 | 27.227 |
0.60 | 12.70 | 25.118 | 12.83 | 25.038 | 13.32 | 24.884 |
0.65 | 12.55 | 22.504 | 12.65 | 22.430 | 13.07 | 22.287 |
0.70 | 12.04 | 20.430 | 12.13 | 20.361 | 12.47 | 20.227 |
0.75 | 11.22 | 18.234 | 11.29 | 18.170 | 11.57 | 18.044 |
0.80 | 10.13 | 16.048 | 10.18 | 15.988 | 10.41 | 15.869 |
0.85 | 8.80 | 14.242 | 8.84 | 14.185 | 9.01 | 14.073 |
0.90 | 7.25 | 12.445 | 7.28 | 12.391 | 7.41 | 12.285 |
0.95 | 5.48 | 10.948 | 5.50 | 10.897 | 5.58 | 10.796 |
0.975 | 4.50 | 10.449 | 4.52 | 10.399 | 4.58 | 10.301 |
(b)螺旋桨转速600r/min | ||||||
ξ /(r·R-1) | V0=5m/s | V0=10m/s | V0=20m/s | |||
Inflow velocity /(m·s-1) | Attack angle /(°) | Inflow velocity /(m·s-1) | Attack angle /(°) | Inflow velocity /(m·s-1) | Attack angle /(°) | |
0.30 | 10.85 | 39.028 | 10.96 | 38.948 | 11.38 | 38.794 |
0.35 | 14.93 | 36.770 | 15.04 | 36.701 | 15.47 | 36.567 |
0.40 | 18.52 | 34.508 | 18.63 | 34.448 | 19.03 | 34.329 |
0.45 | 21.44 | 32.245 | 21.53 | 32.191 | 21.91 | 32.085 |
0.50 | 23.57 | 29.881 | 23.65 | 29.832 | 23.99 | 29.736 |
0.55 | 24.87 | 27.525 | 24.94 | 27.481 | 25.24 | 27.393 |
0.60 | 25.35 | 25.159 | 25.41 | 25.118 | 25.66 | 25.038 |
0.65 | 25.05 | 22.542 | 25.10 | 22.504 | 25.31 | 22.430 |
0.70 | 24.04 | 20.465 | 24.09 | 20.430 | 24.26 | 20.361 |
0.75 | 22.41 | 18.267 | 22.44 | 18.234 | 22.59 | 18.170 |
0.80 | 20.23 | 16.079 | 20.26 | 16.048 | 20.37 | 15.988 |
0.85 | 17.58 | 14.271 | 17.60 | 14.242 | 17.69 | 14.185 |
0.90 | 14.48 | 12.472 | 14.50 | 12.445 | 14.56 | 12.396 |
0.95 | 10.95 | 10.974 | 10.96 | 10.948 | 11.00 | 10.897 |
0.975 | 9.00 | 10.474 | 9.01 | 10.449 | 9.04 | 10.399 |
300r/min(ns=5r/s)时不同相对半径ξ处叶素小微段推力dT如图 5所示。
图 5中无控制指激励器关闭状态,在不同前进速度V0下,峰-峰值6和7kV 2种模式与无控制时的叶素推力进行比较可以看到,SDBD诱导的等离子体射流均有明显增效作用。
相对半径ξ介于0.30~0.45区间时,施加激励与无控制的推力曲线重合度较高,这主要是此区间叶素迎角相对较大,出现了严重的流动分离,超出了等离子体抑制流动分离的能力范畴,因此该区间等离子体对叶素推力作用效果不明显;ξ介于0.45~0.90区间段,该区间推力dT明显提高,是等离子体增效的主要区间,具体在ξ=0.45处开始出现增效,其中ξ~(0.50,0.80)推力dT增加幅度最为显著,ξ~(0.80,0.90)处的增加幅度相对减小;ξ介于0.90~0.975区间段,此范围叶素已接近螺旋桨桨尖部位,迎角和弦长都比较小,几乎没有发生流动分离,等离子体不发挥增效作用,因此推力曲线重合。
随着激励电压增大,等离子体射流对螺旋桨增效幅值的分布规律大致相同。7kV时增效作用略好,主要体现在ξ介于0.50~0.65区间段,说明该区间等离子体体积力和射流速度的增大可以更有效地抑制流动分离。
图 6给出了螺旋桨在实验电压6和7kV时等离子体射流对螺旋桨推力的增效情况。可以看出,300r/min转速时在3种前进速度V0下等离子体射流均可显著增大螺旋桨推力,不同V0时其增效情况略有差异,其中7kV时的流动控制效果要好于6kV,螺旋桨推力分别增加了6.6%~10.8%、9.2%~10.9%。
3.2 600r/min时螺旋桨推力分析螺旋桨转速600r/min(ns=10r/s)时,在2种激励电压下的各叶素微段推力dT如图 7所示。
可以看出,600r/min时等离子体作用范围主要介于ξ~(0.35~0.80)区间。其中,仅在V0=10、20m/s 2种速度时,ξ~(0.45,0.55)区间的等离子体增加了叶素微小段推力,而其余工况和叶素区间均出现了使推力减小的负作用。其中,ξ介于0.65~0.8叶素区间螺旋桨推力减小幅度最大,此区间同时也是叶素推力下降最为剧烈的一段,说明该处流场处于临界状态,等离子体射流对流场的改变存在不确定性。峰-峰值6和7kV 2种激励电压下桨叶推力分布曲线重合程度很高,表明螺旋桨转速600r/min时,流场气体动量很大,等离子体射流对流场的影响能力减小。
与螺旋桨300r/min时施加等离子体射流的效果比较,两种转速下等离子体对螺旋桨推力产生影响的叶素分布区间长度大致相同,相对半径比区间长度均为0.45左右。但600r/min时等离子体在ξ=0.35处开始产生气动影响,相对于300r/min时ξ=0.45,其气动作用区间整体向桨根移动了大约0.1的相对半径长度。
图 8为螺旋桨600r/min时,施加等离子体射流后对螺旋桨推力的作用情况。结合叶素推力分布图 7可知:激励器在6和7kV诱导等离子体射流时,均只在不同V0下的一小段叶素区间略微增加了叶素推力,但总作用效果减小了螺旋桨推力,减小幅值0.52%~2.1%。其中,6kV比7kV放电时使推力减小的幅度更大,对应分别减小了1.4%~2.1%、0.52%~1.7%。
综上可知,螺旋桨转速不同时等离子体射流产生的作用效果也会不同,由于不同转速螺旋桨流场的不同,等离子体射流对推力存在增效和减效2种情况;在增效和减效2种情况时的等离子体射流,当提高激励电压,相比之前较低电压均有利于提高螺旋桨性能;螺旋桨前进速度变化会使等离子体的控制作用发生轻微变化,但其影响程度与转速相比差距较大,不是主要影响因素,说明飞艇在以较低的相对来流速度前进时,应主要针对螺旋桨转速这一参量寻求最佳的等离子体增效方案。
4 结 论(1) 提出了一种模拟平流层螺旋桨等离子体流动控制效果的地面实验方法,该方法基于等离子体射流雷诺相似准则和螺旋桨雷诺相似准则,根据螺旋桨叶素理论,通过地面实验分析平流层螺旋桨在低雷诺数条件下的气动特性。
(2) 平流层螺旋桨SDBD激励器采用微秒正弦连续式激励模式时,等离子体射流对平流层螺旋桨推力的作用效果主要与转速有关;在同一转速下,螺旋桨不同前进速度对等离子体控制效果的影响很小;在螺旋桨300r/min可增效转速,提高激励器电压有助于增大螺旋桨推力。
(3) 本文实验条件下,模拟了平流层20km SDBD激励器在峰-峰值13.6kV,频率10kHz时的等离子体射流对螺旋桨的气动效果,具体当螺旋桨以5,10和20m/s的速度前进时,在转速300和600r/min 2种情况下,等离子体放电时使螺旋桨推力分别提高了6.6%~10.9%和减小了0.52%~1.7%。
[1] | 郭建国, 周军. 临近空间低动态飞行器控制研究综述[J]. 航空学报, 2013, 35(2): 320-331. Guo Jianguo, Zhou Jun. Review of the control of low dynamic vehicles in near space[J]. Acta Aeronautica et Astronautica Sinica, 2013, 35(2): 320-331. |
[2] | Gambler D E, Arena A. Automated dynamic propeller testing at low Reynolds numbers[R]. AIAA-2010-853, 2010. |
[3] | Li Y H, Wu Y, Zhou M, et al. Control of the corner separation in a compressor cascade by steady and unsteady plasma aerodynamic actuation[J]. Experiment in Fluids, 2010, 48: 1015-1023. |
[4] | Ma R, Liu P Q. Optimization design study of low-Reynolds-number high-lift airfoils for the high-efficiency propeller of low-dynamic vehicles in stratosphere[J]. Science China Technological Sciences, 2010, 53(10): 2792-2807. |
[5] | Selig M S. Wind tunnel aerodynamic tests of six airfoils for use on small wind turbines[R]. AIAA-2004-1188, 2004. |
[6] | Benard N, Moreau E. Effects of altitude on the electromechanical characteristics of dielectric barrier discharge plasma actuators[R], AIAA-2010-4633, 2010. |
[7] | 许成杰, 杨旭东, 朱敏临. 近空间桨梢小翼螺旋桨布局气动增效研究[J]. 航空计算技术, 2011, 41(5): 61-64. Xu Chengjie, Yang Xudong, Zhu Minlin. Aerodynamic synergistic mechanism of high altitude proplet propeller configurations[J]. Aeronautical Computing Technique, 2011, 41(5): 61-64. |
[8] | Wang J J, Choi K, Feng L H, et al. Recent development in DBD plasma flow control[J]. Progress in Aerospace Sciences, 2013, 62: 52-79. |
[9] | 王勋年, 王万波, 黄勇, 等. 介质阻挡放电等离子体对翼型流动分离控制的实验研究[J]. 实验流体力学, 2011, 25(4): 9-14. Wang Xunnian, Wang Wanbo, Huang Yong, et al. Investigation of flow separation control on an airfoil using DBD plasma actuators [J]. Journal of Experiments in Fluid Mechanics, 2011, 25(4): 9-14. |
[10] | Rizzetta D, Visbal M. Plasma flow control simulations of a low-Reynolds number low-aspect-ratio wing[J]. Computers and Fluids, 2012, 70: 95-114. |
[11] | 张攀峰, 王晋军, 施威毅, 等. 等离子体激励低速分离流动控制实验研究[J]. 实验流体力学, 2007, 21(2): 35-39. Zhang Panfeng, Wang Jinjun, Shi Weiyi, et al. Experimental study on the separation control by plasma actuator in subsonic flow[J]. Journal of Experiments in Fluid Mechanics, 2007, 21(2): 35-39. |
[12] | Benard N, Jolibois J, Moreau E. Lift and drag performances of an axisymmetric airfoil controlled by plasma actuator[J]. Journal of Electrostatics, 2009, 67: 133-139. |
[13] | 雷光新, 刘巍, 杨涛. 平流层飞艇螺旋桨初步设计方法[J]. 导弹与航天运载技术, 2014, (2): 1-4. Lei Guangxin, Liu Wei, Yang Tao. Primary design method of stratosphere airship propeller[J]. Missiles and Space Vehicles, 2014, (2): 1-4. |
[14] | 姜波, 翟彬, 靳磊. 基于叶素和动量组合理论的直升机性能计算[J]. 电光系统, 2011, (1): 45-48. Jiang Bo, Zhai Bin, Jin Lei. Performance calculation of helicopter based on blade element and momentum theory[J]. Electronic and Electro-optical Systems, 2011, (1): 45-48. |
[15] | 程钰锋, 聂万胜, 车学科. 等离子体提高螺旋桨桨根翼型气动性能的仿真研究[J]. 核聚变与等离子体物理, 2012, 32(3): 265-270. Cheng Yufeng, Nie Wansheng, Che Xueke. Numerical study on plasma flow control to enhance the aerodynamic characteristic of the aerofoil on propeller root region[J]. Nuclear Fusion and Plasma Physics, 2012, 32(3): 265-270. |
[16] | 程钰锋, 聂万胜. 离子体提高平流层螺旋桨气动性能的数值分析[J]. 核聚变与等离子体物理, 2012, 32(4): 372-378. Cheng Yufeng, Nie Wansheng. Numerical analysis of the effect of plasma flow control on enhancing the aerodynamic characteristics of stratospheric screw propeller[J]. Nuclear Fusion and Plasma Physics, 2012, 32(4): 372-378. |
[17] | 刘沛清, 马蓉, 段中喆, 等. 平流层飞艇螺旋桨地面风洞试验[J]. 航空动力学报, 2011, 26(8): 1775-1781. Liu Peiqing, Ma Rong, Duan Zhongzhe, et al. Ground wind tunnel test study of the propeller of stratospheric airships[J]. Journal of Aerospace Power, 2011, 26(8): 1775-1781. |
[18] | 刘沛清, 马利川, 段中喆, 等. 平流层飞艇螺旋桨相似准则分析与验证[J]. 北京航空航天大学学报, 2012, 38(7): 957-961. Liu Peiqing, Ma Lichuan, Duan Zhongzhe, et al. Study and verification on similarity theory for propellers of stratospheric airships[J]. Journal of Beijing University of Aeronautics and Astronautic, 2012, 38(7): 957-961. |
[19] | 车学科, 聂万胜, 侯志勇, 等. 地面实验模拟高空等离子体流动控制效果研究[J]. 航空学报, 2015, 36(2): 441-448. Che Xueke, Nie Wansheng, Hou Zhiyong, et al. Study on high altitude plasma flow control simulation through ground experiment[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(2): 441-448. |
[20] | 冯立好, 王晋军, Choi K S, 等. 离子体环量控制翼型增升的实验研究[J]. 力学学报, 2013, 45(6): 815-821. Feng Lihao, Wang Jinjun, Choi K S, et al. Experimental investigation on lift of a plasma circulation control airfoil[J]. Chinese Journal of Theoretical and Applied Mechanics, 2013, 45(6): 815-821. |
[21] | 刘沛清. 空气螺旋桨理论及其应用[M]. 北京: 北京航空航天大学出版社, 2006: 59-62. |
[22] | Liu P Q, Ma L C, Duan Z Z. Aerodynamics properties and design method of high efficiency-light propeller of stratospheric airships[R]. IEEE-2011-8041, 2011. |