“仿生学”沟槽减阻仿真分析及机理研究

2021-05-18徐朋飞

航空发动机 2021年2期

陈璠，徐朋飞

（1.中国航空发动机研究院，北京101300；2.中国航发沈阳发动机研究所，沈阳110015）

0 引言

仿生学是指运用从生物界中所发现的机理和规律，来解决人类实际复杂问题的一门综合性交叉学科，“仿生学”在航空发动机领域最典型应用是沙丘驻涡火焰稳定器和耦合降噪仿生前缘非光滑形态翼型等。北京航空航天大学高歌教授[1]基于沙漠中新月形沙丘研制出的沙丘驻涡火焰稳定器，可提高燃烧效率和火焰稳定性，减小了流体阻力和振荡损失，大幅度提高航空发动机合格率；吉林大学孙少明博士[2]通过量化长耳鸮鸟翼形态与构型耦合消声降噪特征，建立的耦合仿生非光滑形态及特殊翼型构型，可有效延缓翼型绕流场边界层分离，减弱翼型表面流体压力脉动，进而减少声能产生；西安交通大学的刘小民等[3]通过提取苍鹰尾缘非光滑形态的降噪特征元素建立了仿生叶片结构模型，能够有效减小叶片表面的非定常压力脉动和尾迹涡引起的气动噪声。近年来的研究发现，部分生物如鲨鱼和海豚表皮具有沟槽结构可实现流动减阻。因此，在航空发动机叶片减阻和燃滑油系统流动减阻等方面，“仿生学”沟槽减阻方法具有潜在的应用价值。对此，国内外学者开展了一系列的研究工作。其中，在纵向沟槽研究方面，最早由NASA 开展了系统化试验研究[4]，确定了该结构能有效地减少表面摩擦阻力；Walsh[5-7]、Park 等[8]通过试验得出该结构可产生4%-8%不等的减阻率；Wang 等[9]发现V 型沟槽三角翼在6°攻角时减阻率高达40%；NASA 兰利中心在Learjet 型飞机上开展的飞行试验和李育斌等在运七飞机缩比模型上开展的试验均表明沟槽结构可使其减阻率最大达到8%[10-11]；Speedo公司[12]研发带有凹槽结构的泳衣可使总阻力减小4%，极大地提高选手的泳速。在横向沟槽研究方面，Kramer[13]制备的仿生海豚皮减阻率达50%；Scholle等[14]对低速流动下横向沟槽减阻性能进行了理论分析；De Angelis、Bechert 等[15-16]分析了横向沟槽减阻效果，研究了沟槽表面的速度场分布及表面边界层速度分布对表面减阻性能的影响；宋保维、刘占一等[17-18]通过数值模拟的方法研究了沟槽内产生的旋涡和沟槽间距对减阻效果的影响。由此可见，目前研究成果证实了“仿生学”减阻沟槽结构具有明显的减阻效果，但对于航空发动机高效率叶轮机械及管流部件的优化设计仍有差距，缺乏关键参数对于沟槽减阻效果的对应关系。

本文基于国内外对仿生学结构的减阻效果研究基础上，通过研究纵向沟槽及横向沟槽的减阻效果，利用数值模拟方法对流场结构、仿生学沟槽减阻问题进行机理分析。

1 模型与网格

1.1 模型简化及尺寸确定

便于对减阻效果进行量化，对减阻率进行定义

式中：R为减阻率；F为原构型所受阻力；F1为运用仿生结构的构型所受阻力。

针对单一沟槽结构，国内外学者进行了广泛地研究，根据丛茜等[19]的研究成果表明，刀刃型结构沟槽高应力区出现在刀刃尖端，减阻效果最理想，但该结构容易失稳。而本文研究的锯齿结构凭借其具有高稳定性而被广泛采用。刀刃型沟槽及锯齿型沟槽结构如图1 所示，其中s为沟槽宽度，h为沟槽深度。

纵向沟槽锯齿结构的无量纲沟槽深度h+和宽度s+定义为

图1 锯齿型沟槽（上）和刀刃型沟槽（下）

式中：v为流体黏度；uτ为壁面剪切速度

式中：τw为壁面剪切力；ρ为流体密度。

根据Walsh 等学者的试验研究结果，纵向沟槽的锯齿结构分别在（h+≤25，s+≤30）范围内具有减阻效果，且当h+=s+时减阻效果最佳。

本文探讨低速流对结构减阻的影响，同时考虑模型尺寸的可行性，单个沟槽的尺寸为

式中：U为来流速度。

选择来流速度为5 m/s 时，Re=1×106，计算流体为水，温度为293 K，根据计算结果设定单个沟槽s=h=0.1 mm，同时改变每个沟槽间的距离x，分别在x=0、0.05、0.10 mm时，计算探讨沟槽间距对减阻效果的影响。

本文建立了1段沿流向的3维空间计算域，如图2所示，设置流向尺寸为1750 mm，垂向尺寸为220 mm。为保证流体流动达到完全湍流状态，沟槽放置在沿流向尺寸为1050 mm的位置，并且垂向尺寸满足流场上下面互不影响。为了能在降低洞壁效应的同时控制计算量，将沟槽数量n设定为5，为此设定展向尺寸为1.2 mm；同样为使横向沟槽与纵向沟槽具有可比性，沟槽减阻区表面为1.2 mm×1.2 mm，沟槽形状如图3所示，具体尺寸见表1。

图2 3维计算域

图3 沟槽形状

表1 沟槽形状及尺寸

1.2 网格划分

采用GAMBIT6 软件对计算域生成的6 面体结构化计算网格如图4所示。网格第1 层厚度为10-5m，对沟槽附近局部加密，边界层网格按拉伸比例1∶1.1 拉伸 10 层，网格数量约120万。

图4 计算域网格

2 数值模拟

2.1 求解条件设置

流场介质为不可压缩黏性流体，流场为定常流动。流体温度为293 K，密度为998 kg/m3，黏性系数为1.003×10-3kg/（m·s）；采用速度入口，压力出口；下壁面设置为固壁条件，上壁面以及左右壁面均设置为对称边界条件。

2.2 数值方法验证

利用ANSYS CFX 软件进行流场仿真计算，并采用剪切应力输运k-ω模型进行数值模拟，通过在光滑平板计算得出的模拟结果与经验公式得出的理论摩擦阻力系数进行对比来验证数值方法，对比结果见表2。

表2 模拟结果与理论值对比

通过对比可知，在Re=1.18×107的流场中出现最大相对误差为0.88%，说明该模型能较好地模拟沟槽减阻性能。

2.3 数值模拟结果分析

选取雷诺数及沟槽构型作为变量进行CFD仿真计算，通过比较沟槽表面与平板的阻力求得减阻率。

不同尺寸沟槽构型如图5 所示。从图中可见横向及纵向减阻率随雷诺数的变化情况。MZ1、MZ2、MZ3 分别为构型 1、2、3 的纵向沟槽；MH1、MH2、MH3 分别为构型 1、2、3的横向沟槽。

图5 沟槽减阻率随雷诺数的变化

对图中曲线分析可得：

（1）所有构型的减阻率随雷诺数的增加而降低，且雷诺数在所有变量中对减阻率的影响占主导地位，当Re=5.08×106时达到最大，在Re=1.18×107时为最小。

（2）纵向沟槽在计算范围内的任一雷诺数下减阻效果均优于横向沟槽的，对于纵向沟槽，减阻率比横向沟槽高约1个百分点。

（3）3 种构型为纵向沟槽时，在不同雷诺数下减阻效果趋势一致，在同一雷诺数下，随着沟槽数量增加，沟槽减阻效果提高，但提升效果不大。

（4）从曲线斜率可见，对于纵向沟槽，随着雷诺数的增加，沟槽减阻率降低趋势逐渐明显。而对于横向分布的构型3，随着雷诺数的增加，其降低趋势减弱。

（5）3 种构型为横向沟槽时，曲线差异较大，沟槽构型的差异对横向分布影响较大，在横向沟槽减阻时通过优化沟槽结构得到较好的减阻效果。其中无沟槽间距的构型1 在Re=5.08×106时减阻率为4.14%，减阻效果优于其他2 种构型的，但在Re=1.18×107时减阻率接近0，减阻效果最差；3 种构型在横向分布的减阻率均值为构型2＞构型3＞构型1，沟槽宽度为0.1 mm，深度为0.1 mm，间距为0.1 mm的构型2减阻效果最好，说明适当增加沟槽间距可提升减阻效果。

3 机理分析

3.1 横向沟槽

M2 构型横向沟槽在Re=8.46×106下的速度矢量如图6 所示。从图中可见，在沟槽内顺时针旋转的涡形成涡垫，使上方流体直接从沟槽越过，不与壁面直接接触，可减小摩擦阻力，使顺时针涡顶部流向与主流的一致，对流动有促进作用，也对减阻效果有积极影响。M2 构型横向沟槽在Re=8.46×106下的速度云图如图7 所示。从图中可见，在沟槽内流体速度很慢，形成死水区，壁面边界层变厚，使壁面的速度梯度减小，流体耗散减弱，起到减阻效果。其他构型和雷诺数与之机理类似。

图6 M2构型横向沟槽在Re=8.46×106下的速度矢量

图7 M2构型横向沟槽在Re=8.46×106下的速度大小

分别从局部损失和沿程损失的角度，对几种构型在不同雷诺数下的流动减阻机理进行分析。结果表明：局部损失与壁面形状变化的剧烈程度有关；沿程损失与固体和流体接触的湿面积有关。沟槽间距段越长，沿程损失越大；沟槽数量越多，局部损失越大。通过对比不同构型，发现在低雷诺数工况下，构型1的减阻效果最优，其沟槽间距为0，但沟槽数量远大于其他2种构型的，说明在相对低雷诺数下以沿程损失为主。而在高雷诺数工况下，沟槽间距最大的构型2 的减阻效果最好，说明在相对高雷诺数下以局部损失为主。

3.2 纵向沟槽

Q 判据在流场中用于识别涡特征并可视化，代表流场中速度梯度张量的应变率、张量与涡张量的差值，用于判断流场中涡的情况。M2 构型纵向沟槽在Re=8.46×106下基于Q 判据的涡量如图8 所示。从图中可见，流体在向下游流动时产生方向相反的涡对，该涡对与纵向沟槽结构相互作用，在沟槽底部形成稳定的低速区，使流场稳定，同时产生离散涡，逐步削弱来流涡对强度，达到减阻效果。M2构型纵向沟槽在Re=8.46×106下的速度如图9 所示。从图中可见，纵向沟槽结构使底部形成速度死水区，从而使速度梯度大大减小，沟槽的存在使壁面的流场趋于稳定。其他构型和雷诺数与之机理类似。