一种简单快速的多块网格生成方法.PDFVIP

第 17 卷　第 2 期 应 用 力 学 学 报 . 17　 . 2 V ol N o 2000 年 6 月 CH INESE JOURNAL OF APPL IED M ECHAN ICS Jun. 2000 一种简单快速的多块网格生成方法 谭欣星　　　　席　光 ( ) ( ) 长沙电力学院　长沙　410077 　 西安交通大学　西安　710049 摘　　要 利用有限元中的形函数进行网格插值映射分块生成网格, 构造了一种块网格拓 扑结构, 使得块与相邻块之间自动组合, 并利用该方法对二维、三维复杂区域进行分 块网格生成, 说明该方法生成的网格质量较好, 适合于复杂区域分块求解。 关键词: 块分解; 多块网格; 块组合; 拓扑结构 1　引　言 网格生成是计算流体力学数值分析技术的一个重要组成部分, 是促进CFD 工程实用化的 重要环节。由于流动边界和流场的复杂性, 一种好的网格生成方法应当适用任意形状的边界, 可以控制网格的密度和网格的方向, 运算速度快。在有限差分和有限体积数值分析时, 现有的 网格生成方法中, 基本上有代数方法、保角变换方法和偏微分方程法三大类, 但它们大都局限 于二维或单块区域, 对各种复杂边界通用性不强。本文开发了一种用于有限差分, 有限控制体 积方法数值计算的网格生成方法, 它将复杂的流动区域划分成形状相对简单的多块区域, 在每 块中生成网格, 然后将网格组合形成覆盖流动区域的整体网格。这种网格生成能自动进行, 生 成的网格质量好, 适用任何复杂的二维、三维流动区域。多块网格生成技术还有如下优点: ①它 把大问题化为若干小问题, 缩小计算规模; ②允许单块网格一致, 无需整体网格一致, 甚至各子 域可以用不同离散方法进行计算, 具有较大的灵活性; ③允许在不同的块中选用不同的数学模 型, 以便整体模型更适合于工程实际情况。如气体绕飞行器流动, 在边界层附近为粘性流, 而在 边界层外为无粘流, 二者的数学模型不相同; ④可以对各块并行计算, 还降低对内存的要求。 2　多块网格生成的原理及方法 2 1　网格单元排序 有时复杂的流动区域是由几个相互联系的部分组成的, 它可能具有不规则的边界, 也可能 ( ) 有无流体的区域 如固体障碍引起多连通区域、绕流等 。这种区域宜于使用多块网格, 而在这 ( ) 教育部资助优秀年青教师基金项目 1998[ 9 ] 来稿日期: 修回日期: © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 第 2 期 一种简单快速的多块网格生成方法 95 ( ) 种网格系统中要确定网格节点或单元, 采用间接寻址 排序 比起常用的 i、j 、k 顺序指标更合 适。以二维 T 型管为例, 如图 1, 说明一个复杂区域网格寻址的过程。利用间接寻址, 相邻网格 ( ) 的地址 或网格顺序号, 标记 储存在一个被称着关联系数的矩阵里, 为网格生成作准备。在连 续的网格单