基于有限体积格式的自适应笛卡尔网格虚拟单元方法及其应用.pdfVIP

第 29卷 第4期 空 气 动 力 学 学 报 Vo1．29， No．4 2011年8月 ACTA AERoDYNAⅣⅡCA SD CA Aug．，2011 文章编号：0258-1825(2011)04-0491-05 基于有限体积格式的自适应笛卡尔网格 虚拟单元方法及其应用 胡 偶，赵 宁，刘剑明，王东红 (南京航空航天大学航空宇航学院，江苏 南京210016) 摘 要：基于 自适应笛卡尔网格，结合虚拟单元方法处理浸入边界，发展了一种有限体积格式实现模式，数值模拟了 二维无粘可压缩流动。本文发展的方法相对于传统的有限差分方法具有执行效率快，易于获得守恒性等优点。最 后，利用一些标准算例和多物体流场的计算问题验证了算法的有效性。 关键词：虚拟单元方法；自适应笛卡尔网格；有限体积格式；多物体；欧拉方程 中图分类号：V211．3 文献标识码 ：A 给出了多种处理方法 ，主要有混合 网格方法 (hyh,4d 0 引 言 d) ，融合单元方法(mergedcellapproach)晗，嵌入 众所周知，计算流体力学中一个持续的障碍是复 单元方法(embeddedcellmethod) ，以及虚拟单元方法 杂几何外形的网格生成。自适应笛卡尔网格作为一种 (gllost—cellmethod) 或浸人边界方法 (immersedbound— 新兴的网格生成方法具有结构网格和非结构网格方法 arymethod) 等。 所无法比拟的优势。 自适应笛卡尔 网格生成简单 ，省 自适应笛卡尔网格与虚拟单元方法结合处理复杂 时，同时网格结构简单，易于实现自动化、具有更强的自 物面边界是一个新的发展方向。最近，Dadone和 Gross— 适应能力，可提高计算精度。尤其对于复杂几何外形问 mail在他们的一系列文章 ’ 中，系统研究了结构 网 题，如果使用传统的结构网格，通常使用嵌套网格和多 格中的虚拟单元方法，并把他们命名为 CCST(Curva- 块对接网格等，虽然这些方法已经得到了成功的使用， tureCorrectedSymmetryTechnique)方法，而且这些有效 但它们需要交换不同网格间的数据信息，会遇到一些 的边界处理方法也被他们推广用于笛卡尔网格 J ¨， 复杂的插值和几何的问题 ；非结构网格虽然易于生成复 并取名为GBCM方法。本文基于有限体积格式，将虚拟 杂外形的网格，但非结构网格生成费时，计算时间存储 单元方法与 自适应笛卡尔网格结合，采用了按边扫描的 量一般都比结构网格高；而使用自适应笛卡尔网格则可 有限体积方法。此种方法不同于通常的叉树结构的自 以避免这些问题，使得网格生成变成一个简单而省时的 适应笛卡尔网格方法，虽然增加了计算存储量，但计算 过程。 效率更高。同时，由于采用了有限体积方法而不是通常 对于复杂几何外形的物体，使用笛卡尔网格方法， 的有限差分方法，能更好地满足守恒性要求。此外，本 在物面边界处网格单元必然与物面相交 ，从而产生不同 文利用基于速度的散度和旋度的判据，对流场进行解的 形状的切割网格单元。切割网格单元存在的一个普遍 自适应，能够更为准确地捕捉激波等流动特征，数值结 问题就是在切割的过程中会产生非常小的切割单元，这 果显示本文所使用的自适应方法具有很好的分辨率， 就导致了方程系统的刚性以及在物面处会产生解的非 且易于实现