高升力装置数值模拟研究.pdfVIP

高升力装置的数值模拟研究  吴宗成 朱自强 陈泽民 （北京航空航天大学航空科学与工程学院，北京 100191 ） 摘 要： 本文采用数值模拟方法对高升力装置进行了计算研究，数值方法为求解非定常雷诺平均 N-S 方程，空间离 散采用有限体积方法，无粘通量矢量的离散选用 AUSM+迎风格式，并用 MUSCL 技术获得二阶精度和保单调性，四 步龙格-库塔时间推进，湍流模型采用 B-L 模型和 S-A 模型。数值模拟了典型的二维三维高升力装置及高升力翼身组 合体构型，结果表明了本文数值方法对三维复杂流场的较强模拟能力。流场分析发现，三维空穴分离存在向“剪刀 差”方向的明显展向运动，形成特殊的三维空穴螺旋运动型分离模式。 关键词：高升力装置，N-S 方程 Abstract: The flows around high-lift configuration are calculated and investigated by using numerical simulation method. The unsteady Reynolds-averaged Navier-Stokes equations are solved by using finite volume method, in which the inviscid flux vector is discretized by using AUSM+ scheme and MUSCL technique for second-order accuracy and monotonicity, four-step Range-Kutta method is used for time marching, Baldwin-Lomax and Spalart-Allmaras turbulence models are used for turbulent viscous flow. Several typical multi-element airfoils, wing and high-lift wing-body are simulated. The numerical results indicate that the present method for three dimensional complicated flow is effective. The evident spanwise flow is observed in the three dimensional cove separation region and the especial volute kinetic separation pattern is formed. Key words: high-lift device, N-S equations, 采用计算流体力学的方法准确、可靠、高效地预测三维高升力装置的气动特性，是一个极具挑 战性的课题，也是目前民机研制中必须依靠风洞试验的一个问题。多段翼型绕流流场几乎完全由粘 性绕流所决定，包含了众多的复杂流动现象，层流分离、湍流分离、层流向湍流转捩、激波和边界 层干扰、强逆压梯度、汇流边界层等；三维高升力构型，由于襟翼端面与常规部分机翼端面的结合 处是几何不连续的，形成“剪刀差”构型，此处流场十分复杂，另外，还有接触线转捩和再层流化 等，这些流动现象都是流体力学中的艰难课题，同时，又十分重要，国际国内均投入大量人力物力 开展研究，目前仍是计算流体力学应用的国际前沿问题。 设计部门关于高升力构型的迫切需求，其一，升力和阻力由于缝襟翼缝隙和重叠量以及偏角的 改变而产生的增量的预测能力，这将对高升力装置的试验成本和设计周期产生巨大影响[1] ，目前前 缘和后缘装置的配置优化，包括偏角、缝隙和重叠量，需要长时间大量的风洞试验。其二，雷诺数 的反向效应[2][3]，预测最大升力系数对雷诺数的反向效应的能力将是CFD计算发展的里程碑，目前， 预测附着线转捩、再层流化和再转捩的能力还远远不够。不断提高和改进数值模拟的精确程度使其 成为增升装置的有力设计工具对今后民机设计是十分重要的，可大大地节约成本，降低风险和提高 周期。在设计过程中，先用CFD技术进行大量的数值选型试验，最后再经过相对较少量的风洞