自适应笛卡尔网格ghostcell方法研究-空气动力学学报.pdfVIP

第28卷第1期 空气动力学学报 V01．28。No．1 ACTA 2010年2月 AERoDYNAMICASINICA Feb．．2010 文章编号；0258-1825(2010J01-0061·05 Cell方法研究 自适应笛卡尔网格Ghost 刘剑明1’2，赵 宁1，胡 偶1，王东红1 ” (1．南京航空航天大学航空宇航学院．江苏南京210016I2．徐州师范大学数学系，江苏徐州221116) 摘要：在笛卡尔网格中．利用GhostCell方法处理浸入边界，模拟二维无粘可压缩流。针对静止物体，比较了各 种不同的GhostCell边界条件下的熵误差．总压误差．以及阻力系数。此外为提高激波分辨率，将GhostCell方法 与基于叉树结构的自适应笛卡尔网格算法相结合，数值结果显示本文的方法是切实可行的。 关键词：GhostCell方法；笛卡尔网格I自适应；欧拉方程 中图分类号：V21I．3 文献标识码：A 或变形等产生的非定常问题，故本文研究笛卡尔网格 0 引 言 方法。 众所周知，计算流体力学中一个持续的障碍是 章睁妇中，系统研究了结构网格中的GhostCell方 复杂几何外形的网格生成。现今存在的流场离散方 法，并把他们命名为ST(Symmetry 法主要包括非结构网格方法、结构网格方法和笛卡尔 以及CCST(CurvatureCorrected Tech- 网格方法[1]。虽然如今网格生成技术得到了进一步 Symmetry nique)方法，而且这些有效的边界处理方法也被他 的发展，但网格生成过程仍然是一个费时费事的工 们推广用于笛卡尔网格[w]，并取名为GBCM方法。 作。非结构网格主要在二维流场使用三角网格，三维 使用四面体或棱柱，其主要优点是易于生成复杂外形 的网格，但非结构网格生成费时，计算时间存储量一 方法也被用于二维非结构网格，得到了很好的结果。 般都比结构网格高。结构网格单独使用单一的贴体 网格很难处理复杂外形，而且会导致网格的高度倾 Cell边界方法求解二维静止与运动物体。本文把ST 斜，因而一般不会使用单一的结构网格处理高度复 方法以及GBCM方法用于自适应叉树结构的笛卡尔 杂的几何外形，通常使用嵌套网格和多块对接网格 等，虽然这些方法已经得到了成功的使用，但它们 GBCM方法，以及其他的一些可用的边界处理方法。 需要交换不同网格间的数据，以及处理其他的一些 数值结果显示GBCM方法熵误差最小，但总压误差 复杂问题。使用笛卡尔网格方法，物体与背景网格 和其他的边界方法基本差不多，而且ST方法与GB— 切割，或者使用物体内部的虚拟控制体利用嵌入边 CM方法的阻力系数相对于其他方法都较小。此外， 界来处理与物体的相交。切割网格需要考虑不同情