基于LBM-LES方法的多段翼型流场数值模拟的研究.pdfVIP

第十六届全国计算流体力学会议  基于LBM-LES 方法的多段翼型流场数值 模拟研究 石岩∗，司海青 ，王兵   （南京航空航天大学 民航/飞行学院 飞行技术系，江苏 南京，210016）   摘 要 本文采用LBM与LES 结合的方法，数值模拟研究了二维翼型绕流流场，其中，粒子速度模型采 用D2Q9模型，松弛时间中对应的湍流粘性系数由Smagorinsky模型计算。首先，运用本文研究的LBM-LES 方法对NACA0012翼型进行数值模拟，验证其有效性；然后，为进一步研究其模拟复杂物体绕流的能力，本 文以三段翼型30P30N为研究对象进行数值研究，计算出了给定攻角条件下的气动特性，并与实验值进行比   较。计算结果表明，LBM-LES 方法的结果能够和实验值一致，它是一种有效、可行且较新的数值模拟方法。 关键词 LBM-LES；D2Q9模型；Smagorinsky模型；30P30N  0 引言 增升装置在飞机起飞与着陆过程中发挥着至关重要的作用。虽然相对于飞行时间，使用 它们的时间很短暂，它们却对于有效荷载、燃油消耗、机体噪声(尤其是缝翼的非定常流动 更是产生噪声)具有强烈的影响。  上世纪 70 年代，A.M.O.Smith[1]深入研究了高升力系统流动的复杂现象，为高升力空气 动力学的研究奠定了基础。然后，越来越来的学者涉足到高升力的研究领域，例如，国际上 定期召开的 AIAA CFD High-Lift Prediction Workshop[2]，学者之间互相交流关于高升力研 究工作的进展。  流动的数值模拟方法大致可以分为宏观方法、微观方法。传统的计算流体力学 （Computational Fluid Dynamics，CFD）是将流体视为不间断的整体，尽管采用不同的研 究体系，但它们都是基于连续模型之上的。随着研究的发展，经典力学在某些方面的不足逐 渐显示出来。比如：微米尺度和跨尺度流动，传统方法难以描述或根本不适用。而 LBM 从不 同的角度建立起宏观与微观、连续与离散之间的联系，通过对大量离散粒子的统计分析就可 [3,4] 得出流体运动的宏观特征 。  LBM 常用的方法为 LBGK 逼近法，又称为单松弛因子法。虽然具有编程简单、实施容易 的优点，但由于其对边界条件依赖性太强，对处理高雷诺数时计算稳定性不是很好，难以收 敛。其根本原因是由于 Re增大，根据公式Re UL / υ，可知随着雷诺数的提高，粘度减小，                                                                基金项目：国家自然科学基金面上资助 作者简介：石岩（1990-01），男，江苏徐州，研究生，研究方向为飞行安全与模拟，安全科学与工程，email: symon90@126.com 1    第十六届全国计算流体力学会议  流场中的最大速度梯度引起计算的不稳定，虽然通过加密网格能够提高雷诺数，但这额外的 增加了计算能源与效率，而且可能带来压缩性误差。所以，需要在标准 LBM 中引入大涡模拟 的思想，即本文的 LBM-LES 方法，由函数分布的非平衡态项求得变形率张量，由变形率张量、 Smagorinsky 常数、滤波宽度来确定湍流粘性系数，即松弛