喷流流动控制增升效应数值模拟

2020-03-11朱晓军欧东斌陆志良

空气动力学学报 2020年1期

朱晓军, 李锋, 欧东斌, 周凯, 陆志良

(1. 中国航天空气动力技术研究院, 北京 100074; 2. 南京航空航天大学航空学院, 江苏南京 210016)

0 引言

飞机设计初级阶段必须考虑飞机的起降性能,而增升系统的设计是影响起降性能的最主要因素，它不但决定着飞机结构设计、操纵系统设计、载荷的大小和飞机重量，并且对增升系统的气动效率有着直接的影响[1]。对于大型运输机，无论是高亚声速飞机，还是声速或者是超声速飞机，增升系统性能都是影响设计成败的一个重要因素[2]。因此可以采用一些流动控制技术来改善飞机增升系统的气动性能，从而可以缩短起飞-着陆距离，提高起降性能[3]。

流动控制技术通常可分为被动流动控制和主动流动控制两种方式。主动流动控制[4-7]在增升方面的应用吸引了大量学者的关注，通过对机翼下表面进行喷气控制就是一个很好的应用，可以改变机翼上下表面的压力分布，从而达到增加升力的目的。目前对喷气控制效果研究主要实验和数值模拟两种手段[8-9]。为了能够准确模拟喷气的增升效果，首先要能够准确、细致地模拟带有喷气控制的流场[10-11]。由于带有喷气控制的流场非常复杂且带有强烈的干扰，数值模拟面临很多困难，因此目前大部分研究是以风洞实验为主。随着计算机技术的快速发展，使得数值模拟翼型大迎角绕流成为了可能[12-13]。在此基础上可以对翼型喷气流场进行数值模拟，研究喷气增升的机理。

国外Krathapalli等[14]最先是在一个基于NACA0018翼型的有限翼展模型的下表面0.5倍弦长处开一口，通过风洞实验的方法测出，机翼升力系数和弦向截面表面压力分布随喷气速度的变化。国内北京空气动力研究所李锋等[15]在NACA0018翼型的基础上通过求解 N-S方程，格式为Beam-worming格式的改进型，研究了喷气口位置、压强和喷气速度对翼型升力系数的影响。但鲜有人对三维机翼和飞机的喷气控制进行数值模拟研究。

本文基于NACA0012翼型和机翼对飞机下表面喷流增升效应做了全方面的研究。目前，对于喷气控制广泛应用的限制主要集中在两个方面。第一是要能获得稳定的气源，第二是在机翼全展向开口会对机翼结构安全性产生重大影响[16]。通过从发动机引气的方法可以解决气源问题，然而第二个问题至今都没有得到很好地解决。本文提出只在机翼展向的某一段开口进行局部喷气控制来避免由于全展向开口带来的结构安全问题，并对喷气控制进行了数值研究，找出最合理的局部喷气展向分布，研究了不同喷气速度对机翼气动特性的影响。

1 喷气控制数值模拟方法

1.1 控制方程和湍流模型

本文采用不可压缩的雷诺时均N-S方程组[17]，可表示为：

(1)

动量方程：

(2)

采用有限体积法进行空间离散，对流项采用二阶迎风格式，扩散项采用中心差分格式，压力速度耦合采用基于压力的Simple算法处理，离散代数方程组采用Gauss-Seidel代数法求解[18]。湍流模型采用Spalart-Allmaras模型。SA模型是一个相对简单的一方程模型，求解了一个有关涡黏性的运输方程。该模型比较适合具有壁面限制的流动问题，对于有逆压梯度的边界层问题能够给出一个很好地计算结果,常常用于空气动力学问题中，例如飞行器、翼型等绕流流场分析[19]。

1.2 边界条件

1.3 物理模型及网格划分

二维模型采用NACA0012翼型，翼型弦长为1 m。计算网格采用C-H型网格，如图1所示。计算域四周边界条件距翼型表面距离均为20倍弦长，喷气口附近加密后的网格总数是43 586。最内层网格高度为10-5m，y+≈1。

三维模型采用基于NACA0012翼型的有限翼展模型，如图2所示。翼根弦向截面(z=0)是对称面，半翼展的参考面积是3 m，计算域半展长是12 m，所以翼稍(z=3 m)处于流场之中，这样就考虑了翼尖涡的影响。

图1 二维模型局部网格Fig.1 Local 2D mesh

2 二维模型数值模拟结果和分析

2.1 算法验证及增升机理研究

NACA0012是最简单的对称低速模型，几十年来人们对它进行了大量的风洞实验[20-21]和数值模拟。本算例采用Gregory和O’Reilly的实验设置和实验数据，此实验是1970年在美国 NPL(National Physical Laboratory)实验室4.0 m×2.7 m的低速风洞中进行的。实验来流速度为55 m/s，模型弦长为0.76 m，雷诺数为Re=2.88×106，实验设备和细节可参考文献[20]。Gregory公布的实验数据仅有三个不同迎角下吸力面的压力系数，升力和阻力系数采用文献[21]提供的数据。

图2 三维模型网格Fig.2 3D mesh model

表1 不同迎角下的升力系数Table 1 Lift coefficient at different angles of attack

表2 不同迎角下的阻力系数Table 2 Drag coefficient at different angles of attack

翼型上下表面压力分布也是验证数值模拟准确性的一个重要参数。图3所示为上述条件下NACA0012翼型数值模拟和实验得到的上下表面压力分布。由图可以看出数值模拟得到的结果和实验结果相差不大，在可接受的范围之内。

接下来对喷气增升机理展开讨论。从流场分布上来看，图4是喷口和翼型附近流线图，从图中可以看出在喷气口上游产生一个分离泡，喷气口后游产生一对反向旋转的旋涡。喷气口的存在导致翼型后缘驻点移动到翼型的压力面上。Kutta条件的改变必然引起翼型后缘上下表面出现压差。对涡的存在使得上表面流速增加，压力变小，相对应的下表面流速减小，压力变大。

图3 升力系数比较 Fig.3 The comparasion of lift coefficients

图4 喷口附近流线图Fig.4 Streamline near the jet

图5是有无喷气控制时的翼型周围压力云图。从图中可以看出，没有喷气控制时翼型上下表面压力基本一样，采用喷气控制后翼型下表面的压力要大于翼型上表面的压力，对图3中有喷气控制翼型的压强系数积分后有正升力(CL=0.675 32)。

从外形上来看在翼型下表面进行喷气控制效果和在下表面加装Gurney襟翼类似，相当于增加了翼型的有效弯度，从而增加了绕翼型的典型环量，进而提高了翼型的升力系数。

(a) 无喷气控制

(b) 有喷气控制

2.2 不同喷气位置的影响

基于上述研究，为了研究不同喷口弦向位置对增升效果的影响，只改变吹气口的弦向位置，其余计算条件和2.1节一样。得到的结果如图6所示，图7为不同喷气口位置时的Cp曲线图。

图6 升力系数随喷口位置变化Fig.6 CL Change with the position of the jet

图7 压强系数随喷口位置变化Fig.7 Cp vs the position of the jet

从图中可以看出在喷气控制作用下，下翼面喷气口前产生较大正压，喷气口后形成负压。并且随着喷气口后移，可以增大下翼面的正压区，较少负压区，这样有利于升力系数的增加。

2.3 不同喷气速度的影响

为了研究不同喷气速度对增升效果的影响，只改变喷气速度，其他计算条件均不变。图8和图9是升力系数和压力系数随喷气速度变化曲线。

图8 升力系数随喷气速度变化Fig.8 CL vs the velocity of the jet

图9 压力系数随喷气速度变化Fig.9 Cp vs the velocity of the jet

3 有限翼展模型数值模拟结果与分析

基于上述研究，本节将喷气口的弦向位置固定于0.7c～0.74c之间，研究不同展向位置和喷气速度的喷气控制对机翼气动特性的影响。当不采用喷气控制时，经数值模拟得到：CL=-0.000 016，CD=0.009 59。采用三维模型进行计算时得到的阻力系数稍微大于二维模型的计算结果，随着迎角增加，三维模型的升力系数会逐渐小于二维模型的计算结果。这是由于对于三维模型在下表面进行喷气控制时可以增加升力，但是湍流摩擦阻力和压差阻力也是增加的，而且三维模型会受到翼尖涡的影响，诱导阻力使得升力减小，升力越大诱导阻力也就越大。

将翼根到翼稍平均分为3段，每段的区域分别编号为1到3，每段展向长度均为1 m。本节按喷气口位置和喷气速度分为以下6种情况进行讨论：