波浪条件下地效翼型气动力的环量控制研究
Circulation Control of Airfoil Aerodynamic Force Under Ground Effect of Wavy Wall
通讯作者: 孙建红,男,教授,博士生导师,电话(Tel.):025-84896099;E-mail:jhsun@nuaa.edu.cn.
责任编辑: 王一凡
收稿日期: 2021-09-30
基金资助: |
|
Received: 2021-09-30
作者简介 About authors
刘浩(1995-),男,河北省保定市人,博士生,从事地效飞行器和流动控制研究.
波浪海面气流环境相对复杂,严重影响了地效飞行器巡航过程中的稳定性与安全性.针对波浪条件下的地效翼型,应用数值方法分析了不同海况波浪形状和不同飞行攻角对翼型气动特性的影响,进一步研究了定常吹气和周期吹气方法对地效下的环量控制翼型气动力影响规律.计算结果表明:波浪海面地效下,翼型升力系数随波浪周期性变化,波浪的波高较高、波长较短,翼型攻角较大时升力系数波动的幅度较大;通过与波浪相对运动相同周期的吹气环量控制方法,可以有效地控制波浪海况下翼型升力系数的波动,增强地效飞行器的稳定性与安全性.
关键词:
The interaction of airflow and sea waves seriously affects the flight stability and cruising safety of ground effect vehicles. The influence of different sea states and different angles of attack were analyzed numerically on the aerodynamic characteristics of the airfoil under ground effect of wavy wall. The influence of the constant blowing and periodic blowing methods was further studied on the aerodynamic force of the airfoil under ground effect. The simulation results show that the lift coefficient of the airfoil changes periodically with the wave under the wavy ground wall. The amplitude of the lift coefficient fluctuation is larger with the increasing of wave height and angle of attack, or the decreasing of wavelength. Applying the circulation control method for periodic blowing in the same period as the relative motion of the waves can effectively weaken the fluctuation of the airfoil lift coefficient under wavy sea conditions and improve the flight stability and safety of ground effect vehicles.
Keywords:
本文引用格式
刘浩, 孙建红, 孙智, 陶洋, 王德臣, 刘光远.
LIU Hao, SUN Jianhong, SUN Zhi, TAO Yang, WANG Dechen, LIU Guangyuan.
地面效应就是指飞行器近地飞行时,受地面边界的影响,气动特性发生改变的现象.充分利用地面效应,贴近地面或水面巡航飞行的高速运载工具一般称为地效飞行器.由于陆地建筑物、人类活动以及地形地貌的不利影响,地效飞行器研究和应用主要集中在水面.20世纪60年代开始,苏联先后研制并服役了多款大型地效飞行器,包括起飞质量达540 t的“里海怪物”.然而受当时技术条件的限制,复杂环境下的安全性与稳定性成为地效飞行器发展的瓶颈.特别在海面,宽广的环境便于高速飞行,但是由于海面气象环境复杂,如气流与波浪环境,地效飞行器与波浪的相对运动会造成飞行器升力剧烈变化,造成颠簸和飞行不稳定等问题[1-2],为地效飞行器操纵性能与安全性能带来挑战.
对于光滑平面地效工况下翼型的气动特性研究相对广泛.Hsiun等[3]通过数值模拟的方法研究了雷诺数、攻角和相对飞行高度对翼型地面效应下气动力的影响.Ahmed等[4]应用移动带进行了地面效应的风洞试验,认为下翼面与地面之间气流通道的形状是影响翼型气动性能的主要因素.Qu等[5]进一步根据地效对翼型升力的影响,将攻角与飞行高度分为3个区间,即一个正地效区和两个负地效区.刘浩等[6]通过数值模拟的方法研究了不同的特征参数对地面效应下翼型气动特性的影响.姜裕标等[7-8]通过实验和数值计算的方法研究了地面效应对射流增升翼型性能的影响.秦绪国等[9-10]分别用流体体积(Volume of Fluid,VOF)方法和固定壁面模拟波浪,研究了波长和波高对气动力平均值和波动幅度的影响规律.米百刚等[11]通过动网格的方法对二维翼型在自由来流与地面效应的非定常运动流场进行了数值计算分析.Hu等[12]研究了正弦波浪条件下攻角对翼型气动力的影响.当前波浪海面下对于高速地效飞行器的气动性能的研究仍然不够全面,且无有效方法解决地效飞行器稳定性与安全性问题.环量控制是一种有效的气动力控制方法,可以用于地效飞行器在波浪海面的增稳方法研究.
本文应用嵌套网格的方法,对波浪海况下地效翼型非定常流场进行了数值模拟,分别研究了不同波浪形状、攻角条件下翼型的气动特性,以及定常吹气和周期吹气方法对地效下的环量控制翼型气动力的影响.通过周期吹气的环量控制方法有效地减轻了波浪条件下地效翼型升力的波动[13],为高海况下地效飞行器设计与抗浪性能的评估提供了参考.
1 数值方法与验证
不可压非定常流动的雷诺平均Navier-Stokes(URANS)方程:
式中:ρ为密度;ui、uj为雷诺平均速度分量;t为时间;p为压力;xi、xj为空间方向;u'i、u'j为脉动速度分量;σij为应力张量;-ρ
图1
图1
数值计算的嵌套网格与背景网格
Fig.1
Overlapped grids and background grids for numerical calculation
图2
图2
NACA4412翼型压力系数分布
Fig.2
Pressure coefficient distribution of NACA4412 airfoil
2 计算结果与分析
表1 典型地效飞行器机翼参数
Tab.1
飞行器名称 | 巡航速度/(km·h-1) | 平均气动弦长/m |
---|---|---|
里海怪物 | 555 | 19.12 |
鹈鹕 | 490 | 30.50 |
小鹰 | 231 | 9.67 |
表2 不同算例边界条件
Tab.2
算例 | 海况 等级 | 浪高/ m | 风级 | 波速/ (m·s-1) | 周期/ s | 波长/ m | 风速/ (m·s-1) |
---|---|---|---|---|---|---|---|
Case1 | 3 | 0.7 | 4 | 7.2 | 3.8 | 18 | 4.7 |
Case2 | 4 | 1.4 | 5 | 10.0 | 5.4 | 34 | 6.3 |
Case3 | 1.4 | 4 | 7.2 | 3.8 | 18 | 4.7 | |
Case4 | 0.7 | 5 | 10.0 | 5.4 | 34 | 6.3 |
图3
图3
固壁波浪与水面波浪下翼型升力系数曲线
Fig.3
Lift coefficient curves under waves of solid wall and water surface
2.1 波浪形状对翼型气动特性的影响
图4为不同波浪形状海面下NACA4412翼型在两个周期内的升力系数曲线.图中最下方的曲线表示对应时间下翼型尾缘与波浪的相对位置.在翼型与波浪的相对运动中,升力系数随波浪周期波动,会导致飞行器产生周期性的纵向浮沉响应.不同波浪形状海面下,翼型升力系数波动的幅度和趋势不同.波浪形状的波长较短时(Case1、Case3),翼型的升力系数变化趋近于正弦曲线;而波浪形状的波长较长时(Case2、Case4),翼型的升力系数上升的过程比较平缓,下降比较迅速.
图4
图4
NACA4412翼型在不同波浪形状海面下的升力系数曲线
Fig.4
Lift coefficient curves of NACA4412 airfoil under different waves
4种波浪形状海面下(Case1~4)翼型升力系数波动的幅度分别为7.04%、6.37%、15.25%、2.15%,说明波浪海面地效下波浪形状的波长和波高都对翼型升力系数波动的幅度有显著影响.Case3属于波高较高、波长较短的算例,波高高放大了有效迎角的变化,波长短意味着相对运动的波动频率更高.因此,Case3的算例出现翼型升力系数剧烈变化,进而影响地效飞行器的安全性.
图5(a)和5(b)分别为不同波浪形状海面条件翼型在升力系数最大和最小时周围流场的压力系数分布.从图中可以看出,波浪海面与翼型的相对运动主要对翼型下翼面附近流场的压力产生影响,对上翼面的影响较小.
图5
图5
不同波浪形状海面下翼型周围流场压力系数分布
Fig.5
Pressure coefficient distribution of flow fields around the airfoil under different waves
图6
图6
平面地效和波浪地效的气流通道
Fig.6
Air flow channel under ground effect of plane and wave
Case2和Case4翼型在位于迎波面位置时升力系数最大,在位于背波面位置时升力系数最小.主要是由于翼型驶向迎波面的过程中,波浪与下翼面的相对运动对空气产生冲压作用,使气流的动能转换为压力势能,翼型下翼面压力升高;翼型在背离波浪运动时,气流的压力势能转换为动能,使下翼面的压力降低.
图7显示了翼型前缘和尾缘到波面的相对高度hl、ht和气流通道面积s.基于前述分析,为了验证波浪位置和气流通道形状与翼型升力系数的关系,取迎波系数α=ht/hl表示翼型处于迎波或背波位置,气流通道形状参数β=s/[0.5(ht+hl)c]表示气流通道的“扩张—收缩”或“收缩—扩张”形状.以迎波系数α和气流通道形状参数β为自变量,K0、K1和K2为待定系数,采用一阶多项式拟合翼型在3、4级海况正弦波浪海面运动过程中升力系数的变化,即
图7
表3 拟合曲线的待定系数和确定系数
Tab.3
算例 | K0 | K1 | K2 | R2 |
---|---|---|---|---|
Case1 | 0.726 8 | 0.010 9 | 0.519 0 | 0.997 8 |
Case2 | 0.879 3 | 0.093 4 | 0.254 7 | 0.993 9 |
Case3 | 0.675 2 | 0.080 7 | 0.510 9 | 0.994 5 |
Case4 | 0.991 7 | 0.123 1 | 0.096 9 | 0.997 5 |
图8
图9为不同波浪形状海面下翼型升力系数最大(CL,max)和最小(CL,min)时表面压力系数分布.从图中可以看出,波高较高的波浪形状使翼型上翼面的前缘吸力峰峰值变化更加剧烈;波浪形状的波长主要影响翼型下翼面的压力变化,对上翼面几乎没有影响.对于相同的波高,翼型在短波长的波浪形状表面运动时下翼面压力变化更加剧烈.波长较长、波高较矮的波浪形状对地效翼型的气动特性影响最小.
图9
图9
不同波浪形状海面下,升力系数最大和最小时翼型表面压力系数分布
Fig.9
Distributions of pressure coefficient on airfoil surface with maximum and minimum lift coefficients under different waves
2.2 攻角对波浪海况下翼型气动特性的影响
为了进一步研究攻角对地效翼型在波浪海况下的气动特性的影响,采用4级海况波浪形状边界条件,仅改变翼型飞行攻角,进行对比研究.
图10为不同攻角下 NACA4412翼型升力系数随时间变化曲线,翼型以不同攻角(AOA)在相同波浪形状海面上运动时,升力系数波动的相位和趋势相同.图11为不同攻角下翼型升力系数最大值、平均值(CL, ave)、最小值曲线,可以看出在4级海况的地效条件下,NACA4412翼型升力系数随攻角增大而增大,同时翼型升力系数波动的幅度明显增加,当攻角大于16° 时翼型开始失速,升力系数开始下降,波动幅度减小.图12为不同攻角下NACA4412翼型俯仰力矩系数最大值(CM, max)、平均值(CM, ave)、最小值(CM, min)曲线(取距翼型前缘0.25c为力矩中心,其负值表示低头力矩).随着攻角增大,翼型的低头力矩增加,俯仰力矩系数波动的幅度相应增大.较小的机翼攻角有利于减小地效飞行器在波浪海面下巡航过程中的气动力波动.
图10
图10
不同攻角下NACA4412翼型升力系数随时间变化曲线
Fig.10
Lift coefficient curves of NACA4412 airfoil over time at different angles of attack
图11
图11
不同攻角下NACA4412翼型升力系数最大值、平均值、最小值曲线
Fig.11
Curves of maximum, average and minimum lift coefficients of NACA4412 airfoil at different angles of attack
图12
图12
不同攻角下NACA4412翼型俯仰力矩系数最大值、平均值、最小值曲线
Fig.12
Curves of maximum, average and minimum pitching moment coefficients of NACA4412 airfoil at different angles of attack
2.3 环量控制对翼型气动特性的影响
图13
定常吹气和周期吹气是两种常用的环量控制方法.从吹气动量系数方程可知[20],通过控制吹气出口速度可以调节吹气动量系数:
式中:ρj为气流密度;vj为吹气速度;ρ∞为自由来流密度;Sj为吹气口面积;v∞为自由来流速度;Sw为机翼参考面积.二维翼型的吹气口面积和参考面积分别为吹气口高度hj.周期吹气采用正弦周期吹气的方式,最大吹气速度为100 m/s,其吹气速度随时间变化如图14所示,吹气周期与4级海况波浪运动周期相同.
图14
图15为定常吹气不同吹气动量系数下的翼型升力系数,可以看出在定常吹气的地面效应条件下,翼型升力系数与吹气动量系数成正比例关系,吹气动量系数越大,升力系数相应增大.
图15
图15
不同吹气动量系数的定常吹气下环量控制翼型的升力系数
Fig.15
Lift coefficients of circualtion control airfoil under steady blowing with different blowing momentum coefficients
图16
图16
Cμ=0, 0.002时翼型尾缘流线与流场速度分布
Fig.16
Streamlines and velocity distribution of airfoil trailing edge at Cμ=0, 0.002
图17
图17
Cμ=0, 0.002时翼型前缘滞止流线
Fig.17
Stagnant streamlines of airfoil leading edge at Cμ=0, 0.002
升力的瞬时响应特性是非稳态流动控制的关键.Friedman等[22]和姜裕标等[23]分别研究了升力对于脉冲吹气环量控制的非稳态延迟响应,对于正弦周期吹气的升力响应问题的研究较少.图18为周期吹气环量控制翼型与对应时刻相同速度的定常吹气升力系数对比,可以看出周期吹气的翼型最大升力系数显著小于定常吹气,最小升力系数略大于定常吹气.应用周期吹气的方法产生升力增量的响应过程存在延迟,约为0.1倍的吹气周期(即0.1Tw).拟合曲线为定常吹气升力系数曲线在0.2Tw时间内的动态平均曲线,拟合曲线与周期吹气升力系数曲线基本一致.因此可以认为对于当前环量控制翼型,当吹气气流离开吹气口作用于周围流场之后,在0.1Tw时刻对翼型气动力影响最大,在0.2Tw之后对翼型气动力的影响基本消失.
图18
图18
周期吹气环量控制与对应时刻相同速度的定常吹气环量控制的升力系数
Fig.18
Lift coefficients of periodic blowing and steady blowing with the same velocity at corresponding time
式中:tL,max为翼型最大升力的时刻;0.75为相位差;0.1为周期吹气产生升力增量的响应延迟.
图19
图19
最大和最小吹气速度时周期吹气环量控制翼型的压力系数分布
Fig.19
Pressure coefficient distributions of circulation control airfoil at maximum and minimum blowing velocity
图20
图20
最大和最小吹气速度时周期吹气环量控制翼型的流场速度分布
Fig.20
Velocity distributions of circulation control airfoil flow field at maximum and minimum blowing velocity
图21为环量控制翼型有、无施加周期吹气条件下的升力系数曲线,可以看出通过周期吹气的环量控制方法可以明显减小升力系数波动的幅度,使曲线趋于平稳,但难以完全消除升力系数的波动.在4级海况条件下无吹气翼型的升力系数波动幅度为7.45%, 施加周期吹气的翼型升力系数波动幅度为2.49%,通过与波浪相对运动相同周期的吹气环量控制方法可以显著改善地效飞行器的稳定性与安全性.
图21
图21
环量控制翼型有、无施加吹气条件下的升力系数曲线
Fig.21
Lift coefficient curves of circulation control airfoil with and without blowing
3 结论
本文主要研究了不同波浪形状海面、攻角以及环量控制方法对地效飞行器的二维翼型在波浪海面下气动特性的影响,得到结论如下:
(1) 波浪海面条件下,地效翼型升力系数随波浪周期性波动变化,升力系数的变化与迎波系数和气流通道形状参数成正相关关系.波浪形状的波高较高、波长较短时,翼型升力系数波动的幅度较大,对地效飞行器的安全性造成威胁.
(2) 在4级海况的波浪地效下,攻角小于失速攻角时,翼型升力系数和俯仰力矩系数波动的幅度随攻角增大而增大.
(3) 在平面地效条件下,定常吹气的环量控制翼型升力系数与吹气动量系数成正比例关系;周期吹气的方法产生升力增量的响应过程存在延迟;通过与波浪相对运动相同周期的吹气环量控制方法,可以有效地控制波浪海况下翼型升力系数的波动,为适用于复杂海况的地效飞行器设计提供了参考.
参考文献
大型水陆两栖飞机增升装置特殊设计综述
[J]. ,
Review of high-lift devices design for amphibious aircraft
[J]. ,
水陆两栖飞机的关键技术和产业应用前景
[J]. ,
Key technologies and industrial application prospects of amphibian aircraft
[J]. ,
Aerodynamic characteristics of a two-dimensional airfoil with ground effect
[J]. ,DOI:10.2514/3.46949 URL [本文引用: 1]
Aerodynamics of a NACA4412 airfoil in ground effect
[J]. ,DOI:10.2514/1.23872 URL [本文引用: 2]
Airfoil aerodynamics in ground effect for wide range of angles of attack
[J]. ,DOI:10.2514/1.J053366 URL [本文引用: 1]
地面效应下的不同翼型亚声速气动特性分析
[J]. ,
Aerodynamic characteristics of different airfoil in subsonic flow with ground effect
[J]. ,
地面效应对射流增升翼型性能影响实验研究
[J]. ,
Experimental investigation on blowing control airfoil influenced by ground effect
[J]. ,
定常吹气对无缝襟翼翼型地面效应影响的数值模拟
[J]. ,
Numerical simulation of effect of steady blowing slot-less flap airfoil in ground effect
[J]. ,
翼型波浪水面巡航地面效应数值模拟
[J]. ,
Numerical simulation on aerodynamics of airfoil flying over wavy water surface
[J]. ,
地效飞行器近波浪地面大迎角飞行分离流动数值研究
[J]. ,
Numerical research on separated flow around a WIG craft flying over wavy ground at high incidence angle
[J]. ,
近地、水面时的飞行器动态稳定特性数值模拟
[J]. ,
Numerical simulation of aircraft dynamic stability characteristics flying over ground and water surface
[J]. ,
Airfoil aerodynamics in proximity to wavy ground for a wide range of angles of attack
[J]. ,DOI:10.3390/app10196773 URL [本文引用: 1]
Zonal two equation k-ω turbulence models for aerodynamic flows, AIAA-93-2906
[R]. :
Flying-hot-wire study of flow past an NACA 4412 airfoil at maximum lift
[J]. ,DOI:10.2514/3.61127 URL [本文引用: 1]
NACA4412翼型低速绕流数值计算中湍流模型对比
[J]. ,
Comparison of turbulence models for numerical simulation of low-speed flow around NACA4412 airfoil
[J]. ,
环量控制技术研究
[J]. ,DOI:10.7527/S1000-6893.2015.0282 [本文引用: 1]
未来军/民运输机的高性能要求促使近年来环量控制技术正成为研究的新热点。本文简单介绍了环量控制研究的进展;深入讨论了包括二维环量控制翼型标模和CCA/OTW(Circulation Control Airfoil/Over the Wing)实验、半模型子系统实验和三维翼身融合体全机实验等可供CFD验证用的NASA实验研究。在2个尺寸相近的风洞中对同一二维标模的实验结果表明,源于切向吹气的最大升力系数C<sub>L</sub><sub>max</sub>在中等缝道出口高度时可达8~9。数据对比表明此实验结果可供计算流体力学(CFD)验证用。二维CCA/OTW实验表明,发动机位置前移可大幅增大失速迎角和C<sub>Lmax</sub>;CCA后缘吹气噪声的低频部分强度与速度的8次方成正比,高频部分与速度的6次方成正比。半模型子系统的FACT-MAC跨声速实验不仅可研究高雷诺数效应,且可提供2种飞行状态的数据。初步结果表明,与无射流的低速数据相比,在α=25°时C<sub>L</sub>增大约33%,跨声速时在非设计状态下射流可有效地使激波诱导的分离再附,在保持原有强度下激波位置可后推5%的弦长。三维全机CCW/OTW的实验数据尚在整理分析中,但初步结果已表明,应用前缘吹气可将失速迎角增大至25°,C<sub>Lmax</sub>增大至6,正确安排OTW位置可增大升力线斜率等。
Study of the circulation control technology
[J]. ,DOI:10.7527/S1000-6893.2015.0282 [本文引用: 1]
<p>The circulation control technology is being a new reseach highlight recently, due to the high performance needs of the fufure military/commercial transports. The research progress of the circulation control is introduced briefly and the NASA experimental research plan, including 2D CC airfoil benchmark and CCA/OTW(circulation control airfoil/over the wing)experiments, half-model subsystem experiments, and 3D hybrid wing body aircraft experiments, intended for CFD code validation is discussed in detail in the present paper. The experimental data of the same model in two similar sized wind tunnels shows that when the slot exit height is moderate <em>C</em><em><sub>L</sub><sub>max</sub> can be reach 8-9 due to the tangential blowing. Comparison of the 2 tunnels' data indicates that the data can be used for CFD validation. 2D CCA/OTW experiments show that both stall angle of attack and <em>C</em><em><sub>L</sub><sub>max</sub> can be increased largely by the position forward movement of the engine. The lower frequency intensity of the trailing blowing noise is proportional to the velocity's 8th power and the higher frequency's is 6th. Results of the half model subsystem FACT-MAC transonic experiments can be used not only to study the effect of high Reynolds number, but also to provide the results of two flight regimes. The preliminary results indicate that comparing with no blowing the increase of <em>C</em><em><sub>L</sub></em> is about 33% at <em>α</em>=25° at low speed, and the separation induced by shock wave can be effectively reattached and the shock wave's position is moved aft 5% chord while keeping its strength in off-design condition at transonic speed. Being analysed and compared now, the results of 3D large scale full span CCW/OTW aircraft tests will be presented in the upcoming NASA TM. The published data illustrates that with leading edge blowing the wing stall angle of attack can be greater than 25ånd <em>C</em><em><sub>L</sub><sub>max</sub> can be increased to over 6. And the lift curve slope is increased with the application of OTW.</em></em></em></p>
Numerical investigation of circulation control applied to flapless aircraft
[J]. ,DOI:10.1108/AEAT-10-2019-0208 URL [本文引用: 1]
Active circulation control on the blunt trailing edge wind turbine airfoil
[J]. ,DOI:10.2514/1.J056223 URL [本文引用: 2]
超临界翼型的双射流环量控制研究
[J]. ,
Circulation control of double jet flow on supercritical airfoil
[J]. ,
Lift build-up on circulation control airfoils
[J]. ,DOI:10.2514/1.C033841 URL [本文引用: 1]
/
〈 | 〉 |