er方程的非结构网格再生成自适应方法.pdfVIP

维普资讯 2)一21 第15卷 第 2期 空 气 动 力 学 学 报 Ve1．15．No．2 1$97年6月 ACTA ^ERODYNAMICA SINICA Jun．-1997 二维Euler方程的非结构网格再生成自 适 应 方 法★ 型 监 马延文 IJ7· 中国科学院力学研究所．北京 100080 摘要 给出了提高=雏Euler方程定常解质量的非结构网格自适应方法和模拟结 果。流场算法为结点型有限体积法和Lax—Wendroff~式。网格用阵面推进法生成，生 成初始网格时使用均匀的直角坐标背景网格。自适应采用再生网格，但不使用背鼍网格， 而以流场解盼涕度信息来构造再生时所需要的网格尺度参数薅e项。计算了无粘擞波在固壁 上的反射，NACA0012翼型踌声速绕流和马赫数为3的葡台阶绕流， 自适应效果较好． 能处理非常复杂的外形且容易实现自适应。本文对现有的二维非结构网格再生成自适应 方法(adaptiveremeshing)作了尝试和改进。在三角形网格上采用了结点型的有限体积 法 “。网格生成采用阵面推进法”，对此法的改进是使用均匀的矩形背景网格。它除了 使网格生成中的搜索运算简单外，还有另一优点，即如果已知流场概况而无需网格自适应 时，可依靠源项加密它们周围的区域，如果需要自适应则可提供质量更好的初始网格， 减少以后自适应的次数。对网格再生成的改进，是取捎背景网格，用源分布代替，源的 网格参数由上一次流场求出。这样做减少了网格再生成的计算量，并使网格疏密更直接 反映流场梯度的影响。 1 非结构网格生成方法 选用阵面推进法来生成非结构网格。首先给出描述求解区域边界的点集和它们之间 的连接信息，构成初始阵面，然后按求解区域内的网格尺度参数分布，产生新的网格 点，将边界点的连接信息向场内推进，使阵面从边界推进到场内区域，直到整个求解区 域被三角形网格覆盖为止 最后用Laplaeian滤波算子对所得到的网格进行光滑化处理。 ·NNSFC和中国科学院科学与工程计算重点实验室资助项目· 本文于1998年4月1日收到，12~1S日收到盛改稿． 维普资讯 第2期 袁 礼等。=维Euler方程的非结构网格再生成白适应方诰 227 具体的方法和步骤见文献E23。 在生成边界结点和三角形网格时，需要知道网格尺度参数，而它可从背景网格或源项 获得。本文在生成初始网格时年Ⅱ用丁均匀的直角坐标背景网格 “。在背景网格区域内分 布一些给定网格尺度参数的源项，而背景网格点上的网格尺度参数分布就类似于热潦扩 散的结果。可用求解带源项的Poisson方程得到，-满足Diriehlet边界条件 VS=G (1) 式中S为网格尺度参数，G为保持某尺度参数的源项，其表达式为 N 。 G。，= ．(sij．r一，) (2) I ‘1L 式中(f，，)代表背景网格节点，Ⅳ是源的总数目， 是第n个源的强度因子。对二维问题 设源影响随距离平方衰减，则函数．r，和L定义为 fS。，r 点源 f1，ri 点源 L击{ia ，．j．击a，线 ’ 用SOR法求解背景网格内结点上的尺度参数分布 k(1)s+∞【( +…+ ( + ．+，)+ A奎 ]／(2+ ’杰 ．) (4) Diriehelet边界条件和迭代的初始条件均为