直升机旋翼的瞬态飞行地面效应流场模拟

辛 冀, 马成江, 李 攀, 陈仁良

(1. 中国直升机设计研究所, 江西 景德镇 330001; 2. 南京航空航天大学 航空宇航学院, 江苏 南京 210016)

直升机旋翼的瞬态飞行地面效应流场模拟

辛 冀1, 马成江1, 李 攀2,*, 陈仁良2

(1. 中国直升机设计研究所, 江西 景德镇 330001; 2. 南京航空航天大学 航空宇航学院, 江苏 南京 210016)

发展了一种用于分析地面效应中瞬态飞行状态旋翼流场的理论模型，采用CB3D时间步进自由尾迹格式求解旋翼的尾迹结构，并采用面元法模拟地面对旋翼流场的影响，针对离散数值计算中出现的部分尾迹涡线落入地面下方的非物理现象，引入了“等体积”修正方法对地面下方尾迹节点的位置进行修正。新模型计算结果显示，计算所得到的地面效应中“环流”和“地面涡”等特殊流动现象的位置和涡强与试验图像符合良好，可以反映瞬态飞行旋翼地面效应中的旋翼气动力变化机理。

旋翼；地面效应；瞬态飞行；时间步进自由尾迹；面元法

0 引 言

近地飞行是直升机特有的飞行模式。当直升机离地较近时，其旋翼会处于地面效应(In Ground Effect，IGE)状态中，此时由于地面的干扰，流场会表现出与无地效(Out of Ground Effect，OGE)状态不同的某些特殊现象，需要专门研究。

在IGE悬停状态下，地面的出现可以降低旋翼桨盘平面内的诱导速度[1]，因而有利于提高直升机的悬停性能，这一方面的研究较多[2-3]。而当旋翼在IGE状态下由悬停转入小速度前飞时，文献[4-6]中的试验显示该过程中需用功率会呈现出增大的趋势，旋翼拉力和桨毂力矩随前进比的变化方式相比于OGE状态也呈现出不规则的特点。

针对IGE状态下的旋翼流场模拟，学者们开发过多种数值仿真模型，文献[7-8]等采用自由尾迹和镜像法对旋翼地面效应进行模拟；文献[9-11]则采用面元法模拟地面对旋翼流场的影响。近年来关于旋翼IGE状态的相关文献[12-14]对于IGE瞬态飞行状态关注较少。

对于瞬态运动状态下的孤立旋翼，已发表的时间步进自由尾迹分析格式有PC2B[15-16]、CB2D[17]等。但它们在用于模拟更为复杂的旋翼地面效应状态时，暴露出数值稳定性偏低的问题；自由尾迹模型在离散数值计算中，还容易出现部分涡线节点落入地面下方的“非物理”现象[11，18]。

针对这些问题，本文将引入稳定性更高的时间步进自由尾迹格式，以及对“落入”地面下方的尾迹涡线进行位置修正的方法，建立起一个新的旋翼地面效应流场分析模型，对IGE加速前飞时的旋翼气动力变化特点进行分析。

1 计算模型

1.1旋翼时间步进自由尾迹模型

如图1所示，旋翼桨叶使用升力面模型进行建模，桨叶后缘拖出的近尾迹采用Weissinger-L模型建模[19]。因桨尖涡对旋翼流场的影响起主要作用，因此忽略桨叶内侧涡片的影响[20]。

本文引入一种CB3D格式[21]分析旋翼尾迹结构，如式(1)所示：

根据文献[21]的分析，式(1)是一种3阶显式格式，并且消去了二阶误差项中的反阻尼项，数值稳定性更强。

1.2地面干扰模型

本文在地面上布置均布源面元网格以模拟地面效应中地面对旋翼流场的影响，单个均布源面元对空间某点诱导速度计算和控制点坐标求解可参考文献[19]。各源面元的强度依据面元控制点处的流动不可穿透边界条件求解得到。为减小地面离散引起的“噪声”误差[11]，本文根据瞬态飞行旋翼的尾迹特点，以图2所示的环式网格加矩形正交网格相结合的方式布置地面网格。这种网格除了与前飞旋翼的尾迹形状更加契合之外，还能比较方便地调整网格的密度分布，可在不降低计算精度的前提下减小模型的计算量。

自由尾迹方法在离散的时间和空间上求解旋翼尾迹结构时，很容易发生涡元运动至地面下方的非物理现象。针对瞬态飞行的特点，本文引入一种“等体积修正”法[18]，该方法在低速流场不可压假设下，修正原理如图3所示，假设尾迹节点位于一个壁厚为小量δr的薄壁中空“圆管气团”上，“圆管”的下底与地面重合。假定进入地面下方尾迹节点的径向位置满足精度要求。据此根据体积守恒关系:

δr·2πrlast·hlast=δr·2πrorigin·hnew

求解圆管的新高度hnew。

式中，hnew是修正后节点所在“气团”的高度，hlast和rlast是上一步迭代结束时的“气团”厚度和半径，rorigin是这一步迭代结束、但未修正时的“气团”半径，zground是地面处的轴向坐标，znew是修正后的节点轴向坐标。

2 计算结果及分析

算例计算将针对旋翼稳态和瞬态飞行地面效应两种状态开展，以对新建立的模型进行验证。

2.1IGE匀速飞行状态下的旋翼流场计算

以2010年Nathan所进行的试验[5]为参考算例，计算得到四种前进比下的自由尾迹结构如图4所示。

各前进比下流场涡量图的计算与试验结果对比如图5所示。

由图5可见，计算可以捕捉到地面流场中旋翼特殊的卷起涡流结构，并且其形态和强度与试验值相符良好，验证了本文所建立模型的准确性。

2.2IGE瞬态飞行状态下的旋翼气动力分析

由图6可见，本文模型对IGE匀速和匀加速前飞两种情况下旋翼气动力随前进比变化趋势的计算较为准确，在前进比μ小于0.04时，由于尾迹涡线波动剧烈[4-6]，导致计算值与试验值的偏差略大。

由图7可以看出，匀加速前飞时，在时变的自由来流速度下，没有足够的时间形成匀速前飞对应前进比下那种充分发展的旋翼流场，导致旋翼前部(面对自由来流一侧)的特殊涡结构，都会在时间上延迟形成、延迟消失，位置也更加靠前。由于这一原因，造成了图6中匀速前飞时的旋翼气动力变化趋势相比于瞬态飞行更加靠前。

3 结 论

本文引入了稳定性更高的时间步进自由尾迹格式，以及对“落入”地面下方的尾迹涡线进行位置修正的方法，建立起一个新的旋翼地面效应流场分析模型，得到如下结果：

1) 新建立的旋翼地面效应自由尾迹分析模型可以捕捉到有地效前飞中旋翼前部的特殊涡流结构,且其位置和形态与试验结果符合良好。

2) 新模型计算得到的IGE加速飞行状态下旋翼力随前进比的变化趋势同试验值符合良好，验证了新模型对IGE瞬态飞行状态模拟的有效性。

3) IGE匀加速前飞状态下，由于地面附近特殊涡流动的延迟形成，使得旋翼力随前进比的变化趋势相比于IGE匀速前飞状态出现了一个延迟。

[1]Hayden J S. The effect of the ground on helicopter hovering power required[C]//American Helicopter Society 32nd Annual Forum Proceedings, Washington DC, 10-12 May, 1976.

[2]Light J S. Tip vortex geometry of a hovering helicopter rotor in ground effect[J]. Journal of the American Helicopter Society, 1993, 38(2): 34-42.

[3]Lee T E, Leishman J G, Ramasamy M. Fluid dynamics of interacting blade tip vortices with a ground plane[J]. Journal of the American Helicopter Society, 2010, 55(2): 022005.

[4]Curtiss H C Jr, Erdman W, Sun M. Ground effect aerodynamics[J]. Vertica, 1987, 11(1): 29-42.

[5]Nathan N D, Green R B. Wind tunnel investigation of flow around a rotor in ground effect[C]//American Helicopter Society Aeromechanics Specialists’ Conference, San Francisco, CA, 20-22 January, 2010: 362-380.

[6]Ganesh S, Komerath N. Study of ground vortex structure of rotorcraft in ground effect at low advance ratios[R]. AIAA 2006-3475, 2006.

[7]Pulla D, Conlisk A T. The long time structure of the rotor wake in ground effect[R]. AIAA 2005-1408, 2005.

[8]Lee C S, He C J. A free wake/ground vortex model for rotors at low speed in ground-effect flight[C]//Proceedings of the 51st Annual American Helicopter Society Forum, Fort Worth, TX, USA, 9-11 May, 1995: 841-850.

[9]Itoga N, Nagashima T, Iboshi N. A new numerical method for predicting aero-mechanical behavior of a rotor hovering at close proximity to inclined flat surface[C]//Proceedings of the 55th Annual American Helicopter Society Forum, Montreal, Can, 25-27 May, 1999: 511-520.

[10]Griffiths D A, Ananthan S, Leishman J G. Predictions of rotor performance in ground effect using a free-vortex wake model[J]. Journal of the American Helicopter Society, 2005, 50(4): 302-314.

[11]Griffiths D A. A study of multi-rotor interference and ground effect using a free-vortex wake model[D]. America, Maryland: University of Maryland, 2001.

[12]Verna A, D’Andrea A. Unsteady simulations of rotorcrafts in ground effect using a fast panel-fast vortex formulation[C]//35th European Rotorcraft Forum Proceedings, Hamburg, Germany, 22 September, 2009: 188-202.

[13]Phillips C, Brown R E. Eulerian simulation of the fluid dynamics of helicopter brownout[J]. Journal of Aircraft, 2009, 46(4): 1416-1429.

[14]Syal M, Govindarajan B, Leishman J G. Mesoscale sediment tracking methodology to analyze brownout cloud developments[C]//American Helicopter Society 66th Annual Forum Proceedings, Phoenix, AZ, 11-13 May, 2010: 645-674.

[15]Bhagwat M J, Leishman J G. Rotor aerodynamics during maneuvering flight using a time-accurate free-vortex wake[J]. Journal of the American Helicopter Society, 2003, 48(3): 143-158.

[16]邓景辉, 谭剑锋, 张呈林. 基于时间步进自由尾迹的旋翼总距突增气动响应分析[J]. 空气动力学学报, 2011, 29(2): 217-221.

[17]Li P, Chen R L. Rotor unsteady aerodynamics model using an efficient free-vortex method[J]. Aircraft Engineering and Aerospace Technology, 2012, 84(5): 311-320.

[18]陈仁良, 辛冀, 李攀. 一种新的尾迹-地面干扰修正方法[J]. 南京航空航天大学学报, 2012, 44(5): 694-699.

[19]Katz J, Plotkin A. Low-speed aerodynamics(2nd edition)[M]. Cambridge, UK: Cambridge University Press, 2001.

[20]Bagai A, Leishman J G. Rotor free-wake modeling using a relaxation technique-including comparisons with experimental data[J]. Journal of the American Helicopter Society, 1995, 40(3): 29-41.

[21]辛冀, 李攀, 陈仁良. 基于三阶显式格式的旋翼时间步进自由尾迹计算与验证[J]. 航空学报, 2013, 34(11): 2452-2463.

Flowfieldsimulationforahelicopterrotorintransientflyingwithgroundeffect

XIN Ji1, MA Chengjiang1, LI Pan2,*, CHEN Renliang2

(1.ChinaHelicopterResearchandDevelopingInstitute,Jingdezhen330001,China; 2.CollegeofAerospaceEngineering,NanjingUniversityofAeronauticsandAstronautics,Nanjing210016,China)

A new aerodynamic model was developed for the flow field simulation of a rotor in transient flight with ground effect. In the new model, a stable and efficient time-stepping free-wake algorithm CB3D is used to determine the wake geometry, and the ground is modeled with a panel method. In consideration of unphysical phenomenon that a part of the wake vortices may be generated under the ground in the numerical prediction, a constant volume rectifying method is incorporated to rectify the position of the unphysical wake vortices. The analysis results of the new model show that, for a rotor at various advance ratios with ground effect, the predicted recirculation and ground vortex flow regimes agree well with the experimental images with respect to the location and vorticity, which could provide an insight for the rotor transient flying in ground effect.

rotor; ground effect; transient flight; time-stepping free-wake; panel method

V212.4

A

10.7638/kqdlxxb-2015.0066

0258-1825(2017)05-0650-05

2016-05-22;

2016-09-28

中央高校基本科研业务费专项资金资助(NS2014011)

辛冀(1988-)，男，博士，高级工程师，研究方向：直升机空气动力学. E-mail：xinji@nuaa.edu.cn

李攀*, E-mail：lipan@nuaa.edu.cn

辛冀, 马成江, 李攀, 等. 直升机旋翼的瞬态飞行地面效应流场模拟[J]. 空气动力学学报, 2017, 35(5): 650-654.

10.7638/kqdlxxb-2015.0066 XIN J, MA C J, LI P, et al. Flow field simulation for a helicopter rotor in transient flying with ground effect[J]. Acta Aerodynamica Sinica, 2017, 35(5): 650-654.