基于非结构四边形背景网格的各向异性自适应求解算法研究.pdfVIP

LHD 2012 年度夏季学术研讨会 基于非结构四边形背景网格的 各向异性自适应求解算法研究 郑洪伟 中国科学院力学研究所，高温气体动力学国家重点实验室（筹），北京海淀区 100190 摘要 复杂几何外形问题中多尺度及各向异性的特征的捕捉，对于计算结果的精确性有很大的意义。 本文发展了基于非结构化四边形背景网格的各向异性的自适应网格生成方法，能容易生成高质量的 贴体自适应网格，同时占用内存少、效率高的优点，并且能容易捕捉各向异性的特征。本文详细阐 述了具体的数据结构和自适应算法，最后应用 NACA0012 翼型绕流问题的数值模拟来考察发展的自 适应算法。 关键词 非结构化四边形背景网格, 各向异性, 自适应网格, NACA0012 翼型 本自适应网格方法具有占用内存少的优点。此 引 言 方法一大特点是，所有单元的面是按照有非加 随着计算机技术的迅猛发展, 现在计算流 密分别存储在两个不同列表里，这样就可以用 体力学的能力也越来越强。尽管如此, 对于很多 上基于边的算残量方法，也就是只算没有加密 复杂外形的复杂流动, 其流动现象(如激波、边 的边的通量，然后同时更新这个面的左右单元 界层及多介质界面不稳定等)一般具有多尺度及 的残量。我们知道，有限体积法的计算量主要 各向异性的特征。为了精确模拟和捕捉这些流 集中在面的通量的计算上，因此这样做的结果 动特征，传统的网格生成技术很难胜任。比如， 是极大的减少了计算量，提高了计算效率。此 对于多介质流体的界面的捕捉，如果网格没有 方法另外一大特点是，各个不同层的单元也存 达到足够密的情况下，很难捕捉到界面不稳定 储在不同的列表里，因此求解方程可以只在最 性引起的小涡。这就要求发展能够自适应网格 底层的单元中进行，而上面各个层的单元的解 生成技术。 只要从各自下层单元的解通过守恒律直接得 自适应网格中一大类是笛卡儿网格[1-3]。 到。本文应用 NACA0012 翼型绕流问题的数值 它具有快速自动生成的优点。但是它也有其致 模拟来考察发展的自适应算法, 对压力分布、激 命缺点，那就是边界很难精确拟合，因此边界 波位置及升阻力与其它结果进行了细致的比 处理精度不高。另外一大类是是基于非结构化 较，结果表明，它能捕捉各向异性的流场特征， 背景网格的各向异性的自适应网格生成方法 另外其结果与其它结果很吻合。 [4-6]。比如基于三角形和四面体的自适应网格 1 数据结构及自适应网格生成方法 生成方法[4]。这种方法，很容易产生质量很差 的网格，比如对一个三角形某个顶点陆续加密 本方法是基于非结构化背景网格的各向异 几次后，会导致这个三角形偏斜度特别大。 性的自适应网格生成方法。如图 1 所示，单元 A 、 因此，本文发展了基于非结构化四边形背 B 和 C 构成了非结构化的背景网格。单元 D 和 景网格的各向异性的自适应网格生成方法。由 E 是在单元 A 的基础上沿一个方向的加密。单 于是基于非结构化背景网格下的自适应方法， 元 F 和 G 是在单元 C 的基础上沿另外一个方向 因此可以生成贴体网格，这