基于k-kl-w湍流模型的无人机翼型流动特性研究.docxVIP

基于k-kl-w湍流模型的无人机翼型流动特性研究 近年来，太阳能无人机越来越受到重视。1-2。由于临近空间空气稀薄, 飞行器雷诺数较低, 同时飞行器的大展弦比特征使机翼三维效应相对较弱, 翼型的二维气动力特性成为影响飞行器性能的最主要因素。因此低雷诺数翼型的设计与分析显得愈发重要。 不同于高雷诺数状态, 低雷诺数条件下, 流体流动常常处于层流状态, 其抗拒逆压梯度的能力弱, 容易产生分离、转捩等流动现象, 从而对翼型气动性能产生非常大的影响。大量实验表明, 随着雷诺数的降低, 翼型的气动性能急剧下降。Lissaman［3］对低雷诺数翼型的概念、力学基础、实验研究和理论设计等方面进行了综述。Drela Dark［4］描述了层流边界层分离的相关现象, 并分析了不同设计参数对翼型气动特性的影响。白鹏等［5］用数值方法对低雷诺数对称翼型绕流特征进行了研究。李建华等［6］指出了低雷诺数前缘层流分离和翼尖涡对小展弦比机翼的气动特性起的决定性作用。尽管很多学者对低雷诺数流动进行了一系列的理论、实验和数值模拟研究, 但这些研究大都是基于小尺寸飞行器低空状态采用传统湍流模型进行的数值模拟, 通常只针对某一典型雷诺数的流动进行细致分析。因此, 对于高空低速条件下大展弦比无人机翼型较大雷诺数范围的流动特性, 还需要采用更精细化的转捩模型来进一步开展详细的研究。 本文采用求解k-kl-w湍流模型的雷诺平均N-S方程, 数值研究了某太阳能无人机前翼设计翼型在一定雷诺数范围内的气动力参数的变化趋势和流场结构的变化情况。模拟了该翼型在大迎角条件下升力的非线性特征, 并分析研究了其流场结构和产生失速现象的机理, 以及提供高升力的原因。 1 数值研究方法 1.1 从流场转浚特性出发 分离气泡是低雷诺数条件下翼型表面流动的重要特征, 其生成机理与流场的转捩特性有很大关系, 采用转捩模型可以较好地捕捉这一流动特征［7］。本文采用Transition k-kl-w湍流转捩模型, 此模型包含3 个输运方程［8］: 1.2 数值求解及分析 本文采用Gridgen软件对设计翼型进行几何建模, 计算网格为O型结构化网格, 图1 是最终计算网格的局部放大图。翼型的前、后缘处网格进行了局部加密, 计算域的外场为90 倍弦长。 数值求解采用二维雷诺平均N-S方程的有限体积方法; 控制方程中对流通量项采用二阶精度Roe迎风格式离散; 粘性通量项采用中心差分格式离散;时间离散采用隐式LU-SGS方法。 1.3 模型求解结果 本文以SD7037 原始翼型作为验证算例, 参考美国伊利诺伊大学 ( UIUC) 低湍流亚声速风洞实验, 选取雷诺数Re= 3. 0×105状态, 攻角a = -6° ~16°。 为排除因网格条件造成的计算误差, 首先进行网格效应分析。分别构造第一层网格y+为5. 0、1. 0、0. 5 和0. 1 的4 种网格结构进行验证计算, 网格量呈倍数增加, 计算得到的翼型升阻力系数对比如图2 所示。 可以看出y+为1. 0、0. 5 或0. 1 时的计算结果及趋势基本相同, 网格疏密对于气动力计算影响很小。考虑到计算条件及耗费机时, 选用第一层网格y+=0. 5 的网格结构进行接下来的分析计算。 计算得到的SD7037 翼型升阻力系数随攻角变化曲线与实验结果［9］对比如图3 所示, k-kl-w湍流模型及SA湍流模型数值模拟的6°迎角压力分布对比及表面摩擦力分布对比如图4、图5 所示。 由图3 可以看出, k-kl-w湍流模型的计算结果与实验值吻合较好, 而SA湍流模型数值计算得到的阻力系数与实验值相差较大, 由图4、图5 可以看出, SA湍流模型并没有模拟出流动分离与再附等现象。 这是因为SA湍流模型是从经验和量纲分析出发, 其中大部分参数值都是通过一些特殊条件下的实验结果获得, 并不一定适用于所有状态下的流动问题, 对于本文雷诺数Re= 3. 0×105条件下的流动, 该模型没有进行转捩模化或低雷诺数修正, 而且其经验参数涡粘系数产生项的增长系数Cb1取值过大, 导致计算边界层厚度过大, 摩擦阻力增加［10］。因此SA湍流模型并不适用于数值模拟低雷诺数条件下近壁面处的复杂流动。 由图3a) 可以看出当攻角增大时, k-kl-w湍流模型计算误差有所增大, 升力计算值相对较大。这是因为该模型在较大迎角下捕捉翼型上表面分离泡时, 对前缘转捩区域的判断提前, 对后缘分离区［11］的预估过小, 致使数值模拟精度降低。而在压力分布中则表现为: 靠近前缘存在压力平台区 ( 如图4所示) , 后缘分离区比SA湍流模型小, 压力分布包线包围的面积大, 升力系数计算值相对较大。 综合来看, 本文通过与实验公布的气动力参数对比验证了该数值方法的准确性, 另外借鉴相关的直接数值模拟的压力