低动能来流下背负式进气道非定常流动特性分析

2019-05-05刘志敏闫盼盼张群峰黎星佐孙超

北京航空航天大学学报 2019年4期

刘志敏，闫盼盼，张群峰,*，黎星佐，孙超

(1. 北京航空航天大学航空科学与工程学院，北京 100083； 2. 北京交通大学土木建筑工程学院，北京 100044； 3. 中国航空工业集团有限公司沈阳飞机设计研究所，沈阳 110035)

良好的隐身性能可以提升战斗机的生存能力和战斗力，已成为评价新一代飞行器的重要指标。背负式S弯进气道可以大幅降低进气系统的雷达反射面积，提升战斗机隐身性能，在无人机设计中得到了广泛的应用。针对背负式S弯进气道的设计及气动特性研究，国内外已开展了一些研究工作。李大伟等[1-2]设计了一种适用于无人机的背负式S形进气道，随后探究了采用流场控制技术改善背负式S弯进气道出口流场分布，降低畸变指数。郁新华等[3]利用风洞试验探究了背负式进气道在不同迎角和侧滑角时的气动特性。谭慧俊和郭荣伟[4-5]设计了一种背负式无隔道进气道，其选取的进口鼓包能够有效地隔除机身上表面的附面层。Shi和Guo[6]设计了一种锯齿形入口背负式S弯进气道，并对不同迎角和偏航角下进气道流动特性进行了分析。Zhang等[7]采用数值模拟方法探究了无人机背负式进气道总压恢复系数和流量系数随迎角的变化规律。Murakami[8]给出了一种适用于超声速无人机的背负式进气道设计方案，得到了较好的总压恢复和流量特性。张乐等[9]对比研究了三角形、矩形和梯形进口形状对进气道气动及隐身性能的影响。孔德英等[10]提出了一种背负式S弯进气道辅助进气门设计方案。徐诸霖等[11]基于五孔探针的测量方法对S弯进气道旋流畸变进行了评估。张航等[12]对背负式狭缝进口进气道进行了设计和仿真研究，得到了进气道的流动特性和工作特征。但上述研究大多侧重于高动能来流时进气道工作特性，针对低动能来流条件下进气道气动特性的研究相对较少，某些无人机在低动能来流条件下进行地面测试时，进气道内出现强烈的压力脉动，导致进气道结构疲劳破坏，严重影响进气道的性能。

低速时飞机冲压能量降低，基本上依赖发动机本身抽吸能力使气流加速到压气机入口所需速度。由于机身的影响，气流主要由两侧和上部吸入进气道，而且上下气流流动极不均匀，对于追求高隐身特性的无人机来说，采用尖唇缘大后掠的进口形式，在唇口绕流处容易出现分离，导致进气道性能下降，因此针对唇口后掠角度对进气道分离流动的影响需要深入探究。

本文采用改进的延迟分离涡模拟(Improved Delayed Detached-Eddy Simulation,IDDES)方法、高质量的计算网格和高精度计算格式对飞翼布局无人机背负式进气道低动能来流时流场特性进行数值模拟，研究唇口分离流动对进气道性能及其内部压力脉动特性的影响，并通过改型设计对唇口分离流动进行控制改进，降低分离带来的不利影响，提升低动能来流时进气道性能。

1 数值计算方法

1.1 改进的延迟分离涡模拟方法

本文研究问题涉及显著的分离流动，IDDES方法求解该类问题较为适合，因此选用其进行湍流模拟[13-15]，公式如下：

(1)

(2)

(3)

(4)

lLES=CDESΔ

(5)

(6)

(7)

Δ=min{max[cwd,cwhmax,hwn],hmax}

(8)

(9)

fdt=1-tanh[(8rdt)3]

(10)

(11)

fB=min{2exp(-9α2),1.0}

(12)

α=0.25-d/hmax

(13)

式中：ρ为密度；k为湍动能；t为时间；U为速度矢量；μ和μt分别为分子黏性及湍流黏性；Pk为湍动能生成项；ω为湍流耗散比；β=0.075；hwn为垂直壁面方向的网格步长；cw为经验常数，取0.15；hmax为hwn的最大值；F1为经验混合函数；CDES为比例系数；σk、σω、σω2为待定常数；fe为尺度混合函数；Cμ为常数，取0.07；d为网格间距；νt为湍流动力黏度;κ为卡门常数，取0.4。

1.2 离散格式

采用基于格心的有限体积法进行计算，黏性项选用中心差分格式进行离散；对流项的空间离散采用具有二阶精度的Roe格式。选取修正的Venkatakrishnan[16]限制器保证二阶精度插值且具有TVD性质，同时又具有较小的数值耗散。非定常计算采用Jameson[17]提出的双重时间步法进行时间推进，即在控制方程中引入虚拟时间项，利用物理时间步求解真实解，而每一物理时间步通过虚拟时间迭代达到收敛。内迭代收敛判据为残差下降2个量级或内迭代步达到20步。

1.3 本征正交分解方法

背负式S弯进气道在低动能来流时，唇口附近和S弯进气道内部流场存在很强的非定常特性，本征正交分解(Proper Orthogonal Decomposition, POD)方法可以用来从大量的非定常流场数据中提取相干结构，其基本思想是将原始数据分解为一系列POD模态的线性叠加。常用的POD处理方法有经典方法和快照法。当数据结构的空间点数大于采样时间点数时，快照法更加适用，本文选用快照法。Sirovich[18]于1987年提出的快照法处理步骤如下：首先以相同的时间间隔得到截面A的N份脉动速度场数据，假设截面A上存在M个节点，将得到的非定常数据按照下列方式排列。

(14)

由此计算得到自协方差矩阵C=XTX，并求解其特征值及特征向量。

CAi=λiAi

(15)

将特征值按照由大到小的顺序排序，使得λ1>λ2>…>λN=0，特征值的大小代表该模态所含能量的多少。POD模态φi按照下列公式计算：

(16)

继续求解得到各阶POD模态系数ai=φiun，某时刻速度场即可以按下式展开为N个POD模态的叠加：

(17)

2 模型网格及计算条件

2.1 计算模型

计算模型为飞翼布局无人机前机身加背负式S弯进气道，模型尺寸为7 m×4 m×1.5 m，计算域尺寸选取为80 m×40 m×40 m。原型进气道几何外形及安装位置如图1 (a)所示，进气道进口形状为三角形，唇口后掠角为40°，进气道唇口宽1 m，最大高度为0.47 m。改进型进气道几何外形及安装位置如图1 (b)所示，进气道唇口后掠角由40°增加为60°。

图1 背负式S弯进气道计算模型Fig.1 Calculation model of dorsal S-shaped inlet

2.2 网格划分

在计算域内生成非结构的Trim网格，在近壁区采用棱柱体网格，为保证壁面y+值接近于1，壁面垂向第一层网格尺寸设置为3×10-6m，进气道内网格尺度为8 mm，唇缘附近进行局部加密，计算模型总网格数为1 600万。图2为对称面网格分布。

图2 对称面网格分布Fig.2 Grid distribution on symmetry plane

2.3 计算条件

本文主要研究进气道在低动能来流时气动特性，来流速度设为1 m/s，静压为87 900 Pa，静温设为290 K。进气道出口背压设为0.84倍环境静压，壁面均采用无滑移壁面条件，计算时间步长选取为5×10-5s。

3 计算结果与分析

3.1 流场特性分析

为了分析进气道入口流场特性，选用Q判据和进气道唇口附近流线图相结合来识别流场中的分离涡结构。图3给出了背负式进气道内Q值等值面分布图，图4为背负式进气道唇口附近流线图，流线的颜色表示马赫数的大小。可以看出，背负式S弯进气道受机身的影响，进气道下部的气流沿着前机身上表面平顺吸入，进气道两侧和上部气流以较大的转折角吸入，两股气流形成较大的气流夹角，在相互诱导下，在唇口底部沿着唇口边缘卷绕，在左、右两侧形成分离涡。

两侧的分离涡(见图3中红色引线)由进气道唇口底部一直延伸到顶部位置，且受进气道抽吸作用的影响，分离涡并不紧贴唇口而是向进气道内部倾斜。分离涡向进气道顶部发展过程中不断加速，涡核区为高速区而具有很高的能量。分离涡发展到进气道顶部附近时发生破裂，下游流场表现为脉动很强的湍流涡团，使得进气道内呈现强烈的非定常特性。

图3 背负式进气道内Q值等值面分布 (Q=1×107 s-2)Fig.3 Q iso-surface distribution inside dorsal inlet (Q=1×107 s-2)

采用1.3节介绍的快照法对进气道入口截面流态进行POD分析，截面位置如图5所示，共选取2 000个快照进行计算，图5给出了前四阶POD模态对应的速度矢量图及涡量云图，可以看出速度矢量图与涡量云图分布基本重合。从图中还可以发现，进气道上部存在较强的涡结构，该旋涡结构即为进气道入口分离涡。其中POD第1阶模态(模态1)和第2阶模态(模态2)分别对应进气道左侧和右侧的分离涡结构。

图6给出了改进后背负式进气道内Q值等值面分布图。图7为改进后背负式进气道唇口附近流线图。可以看出，唇口后掠角增加使得进气道唇口气流与主流间夹角减小，因此唇口分离涡仅分布于唇缘附近，两侧的分离涡(见图6中红色引线)向进气道内部发展很短的距离后强度便大幅衰弱，并未形成向进气道内部延伸的高强度分离涡，同时进气道顶部区域的涡结构也大幅减小。

图4 背负式进气道唇口附近流线图Fig.4 Streamline near dorsal inlet lip

图5 背负式进气道入口前四阶POD模态Fig.5 The first four order POD modes of dorsal inlet

采用同样的方法对改进型进气道入口截面流态进行POD分析，图8给出了前四阶POD模态对应的速度矢量图，云图表示涡量分布。从前四阶POD模态可以看到紧贴进气道两侧存在条形结构，其代表两侧唇缘的分离涡引起的低速流动区域，从涡量云图可以看出，其强度明显低于原型进气道。改进后进气道入口截面各阶POD模态均未出现高强度的涡结构，表明进气道唇口分离流动得到了有效改善，其强度大幅降低。

图6 改进型背负式进气道内Q值等值面分布 (Q=1×107 s-2)Fig.6 Q iso-surface distribution inside improved dorsal inlet (Q=1×107 s-2)

图7 改进型背负式进气道唇口附近流线图Fig.7 Streamline near improved dorsal inlet lip

图8 改进型背负式进气道入口前四阶POD模态Fig.8 The first four order POD modes of improved dorsal inlet

上述分析表明，进气道唇口后掠角是影响分离涡产生、发展和破裂的重要因素。唇口后掠角由40°增加为60°能够显著改善唇口分离流动。

3.2 流量特性及出口畸变特性分析

图9为原型进气道内横截面流向速度分布云图。可以看出，进气道入口速度具有很强的不均匀性，唇口两侧分离涡的涡核区域流动速度达到超声速，这是由于气流被吸入后，进入涡核内区域的气流能够在压力驱动下持续加速进而形成了局部超声速区。分离涡破碎后该高速区消失，形成大范围的低速尾迹区，加剧了顶角分离区的强度和规模，使得进气道入口有效流通截面积减小，导致进气道质量流量下降。

图10为改进后进气道内横截面流向速度分布云图。对比图9可以看出，入口截面速度分布较为均匀，分离涡内未出现高速区，大面积流场区域流速有所增加。分离涡破碎过程中没有产生大规模的涡结构，进气道顶角分离区的范围减小，进气道有效流通截面积增大，质量流量增加。

图11为0.62～0.70 s进气道质量流量随时间变化曲线。可以看出，改进后进气道质量流量由48.69 kg/s提高到51.41 kg/s。

采用面平均紊流度和稳态周向总压畸变指数来表征进气道出口畸变特性。面平均紊流度εav计算公式如下[19]：

图9 背负式进气道内部横截面流向速度分布云图Fig.9 Streamwise velocity distribution contour of interior sections of dorsal inlet

式中：n为测点数；εi表征每一个监测点上的总压脉动特性，即

其中：P(t)为测量点得到的随时间变化的总压值；Tu为脉动气流取样时间；Pav为时间Tu内该测点的总压平均值。

本文选取的总压监测点分布如图12所示。

图10 改进型背负式进气道内部横截面流向速度分布云图Fig.10 Streamwise velocity distribution contour of interior sections of improved dorsal inlet

图11 背负式进气道质量流量随时间变化曲线Fig.11 Variation curves of mass flux with time of dorsal inlet

图12 背负式进气道出口截面总压监测点分布Fig.12 Distribution of total pressure monitor points on outlet section of dorsal inlet

低压区范围：

θ-=θ2-θ1

综合畸变指数W按照下式进行计算：

从总压脉动特性来看，原型进气道面平均紊流度εav为0.034，改进后进气道唇口分离流动得到改善，流场内非定常特性减弱，进气道出口面平均紊流度εav降低为0.031。改进后进气道出口的综合畸变指数W由0.074降低为0.067，下降了9.5%，出口流场畸变减小，进气道性能得到提升。

图13 σr沿周向分布Fig.13 Circumferential distribution of σr

工况εavΔσ0θ-/(°)W原型0.0340.041350.074改进型0.0310.0361280.067

3.3 压力脉动特性分析

原型进气道唇口两侧高强度分离涡具有很高的能量，在顶部位置破碎后形成大规模的涡结构并随气流向S弯进气道下游运动。这些旋涡加剧了S弯进气道内部流场非定常特性，使得进气道内部出现剧烈的压力脉动。图14为快照法获得的原型进气道前四阶POD模态对应的速度矢量图及涡量分布云图。可以看出，第1阶模态流动结构主要集中在唇口顶部附近，这对应于3.2节提到的进气道入口高强度分离涡破碎引起的进气道顶部区域强涡流场。第2阶模态和第3阶模态分别给出了进气道内部的主要含能结构，这些含能结构直接导致了进气道内的强压力脉动。

图15为改进型进气道前四阶POD模态对应的速度矢量图及涡量分布云图。可以看出，改进后进气道入口顶部的强涡流动结构消失，第1阶模态和第2阶模态分别对应于进气道内部的含能结构。从涡量云图可以看出，其幅值显著小于原型进气道，表明改进后进气道内部含能结构强度下降，因此改进后进气道内压力脉动幅值显著降低。

图16和图17分别给出了数值模拟计算得到的原型和改进后进气道下表面最低点附近脉动压力功率谱密度分布及声压级频谱特性，图16还给出了原型进气道高空台试验的实测值。由图16可以看出，针对原型进气道，数值模拟得到的3个频率值100 Hz、266 Hz和296 Hz与试验测量得到的结果符合较好，这也说明本文选取的数值模拟方法可以较为精确地模拟进气道内部非定常压力脉动。3个频率中，296 Hz为主频，对应最大声压级幅值约为145 dB。由图17可以看出，改进后进气道内部脉动压力的主频率变为200 Hz，其对应的声压级最大值下降了8 dB，为137 dB，其余各阶频率的声压级幅值同样出现一定幅度降低。