基于非结构重叠网格的二维ns方程求解与应用研究-空气动力学学报.pdfVIP

第26卷第3期 空气动力学学报 V01．26．No．3 2008年09月 ACTAAERODYNAMICASINICA Sep．，2008 文章编号：0258．1825【2008}03．0405．07 基于非结构重叠网格的二维N．S方程求解与应用研究 田书玲，伍贻兆，夏健 (南京航空航天大学空气动力学系，江苏南京210016) 摘要：本文发展了非结构重叠网格技术。通过构造任意单元类型网格重叠的网格同边界确定算法、宿主单元搜 索算法以及网格问插值算法，把非结构重叠网格推广应用于求解Ns方程，并结合湍流模型模拟二维湍流流场。N- s方程的求解采用有限体积法，通量计算采用Osher格式，湍流模型采用sA模型。为验证本文方法的正确性和有 效性，模拟了RAE2822翼型和某三段翼型绕流湍流流场，结果表明采用非结构重叠网格模拟复杂外形湍流流场不 但网格生成效率高，而且容易推广应用于具有相对运动的非定常问题，具有较高的工程应用价值。 关键词：计算流体力学；N．S方程；非结构重叠网格；SA湍流模型 中图分类号：V211．3 文献标识码：A 网格，这使得复杂外形的网格块数可以大大减少，它 O 引 言 不但可以处理复杂外形问题，而且非常适合于处理具 计算流体力学经过数十年的发展，已经进入了实 有相对大幅运动的非定常问题。 用阶段，可以对工程中复杂外形流动问题进行模拟。 在处理复杂外形流动问题时。非结构重叠网格虽 这主要得益于结构重叠网格技术和非结构网格技术 然在网格生成上比传统重叠网格有较大优势，但在模 拟粘性流动问题时，由于非结构网格本身的无序性， 的发展。结构重叠网格的概念首先由Steger等在文 献[1]中提出，其基本思想是把计算流场划分成多个 使得求解N．S方程所需的储存量相当大，且计算复 具有重叠部分的区域，各区域单独生成结构子网格． 杂，另外粘性层内的网格结点数也极难控制，给具体 流场计算过程中各个子网格之间通过插值进行信息 求解不同雷诺数下的N．S方程带来了很大麻烦。为 传递。重叠网格技术降低了复杂几何外形的结构网 解决该问题，Fumiya等旧J在模拟粘性问题时把非结构 格生成的难度，特别是在处理多体间具有相对运动的 重叠网格方法推广应用到混合网格上。但是，他们在 非定常问题时。由于绕各个运动物体的网格不需要再 宿主单元查找上采用了单元剖分法，即对物面附近的 生，这就大大减少了网格生成的工作量。而非结构网 结构化网格进行剖分，把四边形(二维)和三棱柱、六 格由于在复杂外形网格生成和其自适应方面的优点 面体(三维)分别剖分为三角形和四面体。由于非结 近年来也得到了普遍应用。 构混合网格的剖分本身就是一个复杂的过程，而且粘 尽管结构重叠网格具有很多优点，但是完成一个 性单元剖分得到的单元异常纤细，常常导致高效的相 复杂外形的结构重叠网格的生成并不是一件简单的 邻单元搜索失败。为解决这些问题，本文针对二维情 事，特别是，随着外形复杂度的增加，流场所需要分的 况提出了可以任意类型单元非结构重叠网格的网格 块数会急剧增加，重叠网格系统的装配的难度更是急 间边界定义、宿主搜索和网格问插值算法，使得程序 剧增加。另外，结构重叠网格的生成需要大量经验。 的通用性和鲁棒性得到大大提高。为验证本文方法 这就影响了其普遍应用。为此，Nakahashi等于1999的有效性。在非结构重叠网格上模拟了绕RAE2822 年在文献[2]中提出了非结构重叠网格的概念。非