多段翼型绕流NS方程数值模拟.pdfVIP

第十届全国计算流体力学会议论文集 多段翼型绕流NS方程数值模拟+ 李风蔚鄂秦陈海昕桑为民肖志祥 西北工业大学空气动力学研究室西安710072 摘要 采用结构与非结构网格对多段翼型绕流进行雷诺平均NS方程数值模拟，获得了计算的压力分布与 实验吻台茛好的结果，说明本文的网格生成技术与流场求解方法是实用、有效的。计算结果的流动显示对进一 步认识多段翼型复杂流动的机理也有一定的帮助。 引言 多段翼型的应用是飞机获得高升力以改善起飞、着陆性能的主要途径。由于其构形与流动 的复杂性，在进行流场理论计算时会遇到一些困难．如缝道、空腔、拐角处流动的模拟等。本 文采用多块网格生成技术．生成贴体及与边界正交的高质量的结构化网格。采用推进面及推进 层方法，生成质量较高的非结构化网格。为了较好的模拟缝道、空腔的流动，应该考虑粘性的 影响。本文采用BL湍流模型、流场分区求解和有限体积法，求解雷诺平均NS方程，研究具有 多缝道的多段翼型绕流。分别对GAW-1两段、某三段及四段翼型进行了计算。结果表明，计算 与实验结果吻合良好，说明本方法可以较好地模拟多段翼型绕流，为多段翼型的设计提供一个 有效的分析手段。 1网格生成 计算两格是空间离散的基础．是流场数值求解的前提和决定流场求解精度的一个关键因 素。本文采用目前比较流行的结构化和非结构网格，对多段翼型绕流进行数值模拟。 1．1多段翼型结构网格生成 对多段翼型这样复杂的外形，应采用多块网格技术分区生成网格。对每个分区分别采 用求解椭圆型方程与代数方法相结合生成贴体、与边界正交的高质量的网格”。…。对网格疏密 及光滑性也进行了成功的控制。网格的分区示意图见图2，生成的GAW-I两段、某三段、四段 翼型网格分别见图1、3、4。 1．2多段翼型非结构网格生成 应用推进层方法“1，生成物面附近大长宽比的粘性计算的网格．应用推进面方法”’．生成 推进层以外的空间的计算网格。最后对网格进行光顺优化以进一步提高网格的质量。本文采用 Laplacian光顺和对角线交换方法，使网格结构达到比较满意的状态，图lO至图12分别给出 GAW—I两段、某三段、四段翼型网格． ‘航空基金资助 2NS方程求解 积分形式非定常NS方程为： 盖5u岫+』F媚2意p8S 采用J∞eson中心有限体积法进行空间离散．应用四步Runge-Kutta方法进行时间推进并获得 定常解。物面边界采用无滑移条件，即u--v--w=O：物面法向压力梯度为零，即OP!知=0；绝热 壁条件，即／，T／an=0。远场边界应用无反射越界条件。为抑制解的波动和顺利地捕捉激波， 对显示的J鲫eson有限体积法必须加入人工粘性。本文对常用的人工粘性差分格式及边界条件 进行了改进”1，实践证明这大大改善了计算的稳定性并在物面减少了人工粘性对真实粘性的干 扰，提高了粘流计算的精度。此外，针对太长宽比粕性网格的特点，将网格长宽比对人工粘性 的影响也加以考虑”3。这将进一步减少人工粘性对物理粘性的干扰，并使方法更具有适用性。 进一步改善了计算的收敛性。另外，还采用了变系数残值光顺技术及传统的当地时间步长以加 快迭代的收敛。 湍流模型采用Baldwin—Lomax两层代数模型。在物面拐角及缝道处，网格中心点有不同的 物面距离，对每一个网格我们取其较近的物面距离计算湍流粘性系数m。 3算例计算及结果分析 圈5给出GAW—l翼型带30％襟翼的计算结果。襟翼偏角为40。，Re=2．2×10。。计算的压力 分布与实验结果吻台的还是比较好的，只是在襟翼上表面接近后缘处与实验结果有一些偏差。 从实验结果看，压力分布出现了拐折，说明流动出现了分离。计算结果虽也算出了分离，但稍 偏后(表明分离的流线图本文没有给出)。这可能是由于湍流模型的限制。 xi06。 图6至图8给出某三段翼型的计算结果。前缘缝翼和襟翼的偏角均为30。，Re=9．0 计算的压力分布与实验吻合良好。从图8的流线圈可清楚地看到缝翼背风面、主翼后缘空腔中 的驻涡以及流动在这-此IIZ涡外缘与物面形成的缝道问加速的情况。从图中可看出襟翼背风面的 旋涡并非由谍翼物面流动分离而产生，困为从流线圈上找不到分离点。但从速度矢量圈(本文 未给出)可发现，从主