一种改进的Delaunay三角形化剖分方法-西安交通大学CFD-NHT.PDF

( ) 2001 年 　第 25 卷 　　　　　　　　　石油大学学报 自然科学版 　　　　　　　　　　Vol. 25 　No. 2 　第 2 期 　　　　　　　　Journal of the University of Petroleum , China 　　　　　　　　Apr. 2001 ( ) 　　文章编号 2001 一种改进的 Delaunay 三角形化剖分方法 1 2 3 徐明海 , 　张俨彬 , 　陶文铨 , ( 1. 石油大学储运工程系 , 山东东营257061 ; 2. 胜利石油学校 ,山东东营257061 ; 3. 西安交通大学 ,陕西西安 710049) 　　摘要 :提出了一种基于BowyerWatson 算法的平面区域 Delaunay 三角化剖分的改进方法。它结合了前沿推进法的 内部结点生成技术和Delaunay 联点网格生成技术 ,使得每插入一点所破坏的单元尽可能地少。采用适当的数据结构 , 使 Delaunay 搜索过程限于局部 ,算法大为简化 ,易于编程 ,浮点计算量少 , 同时也避免了使用函数递归调用。采用在基 网格上定义网格步长的办法控制网格的疏密 ,使网格疏密易于控制。几个算例表明 ,该算法是行之有效的。 关键词 :非结构网格 ;三角形单元 ;Delaunay 剖分 ;数值模拟 中图分类号: T K 121 　　　文献标识码 :A 法。 引　言 非结构化网格是相对于结构化网格而言的 ,它具 1 　Delaunay 网格生成算法原理 有分布灵活、适应复杂区域剖分及局部任意加密的优 Delaunay 原理要求每个三角形的外接圆内不包 点。非结构网格需要较多的前期处理工作 ,网格形成 含其他节点 ,这样的三角形系统具有最小角最大的 一般不能用手工操作 ,须采用计算机自动剖分。三角 特点。根据这一原理 ,可以构造出各种生成二维区 形网格自动剖分以 Delaunay 三角化和前沿推进法两 域网格的算法。 类方法为基础。为了提高自动剖分程度、减小剖分的 Shamos 和 Hoey [13 ] 提出了分治法的思想 ,Lee 手工工作量和运行机时 ,近 20 年内发展了许多剖分 和 Schachter[14 ]对其作了改进和完善 , 把某一点集 技巧。Peraire[1 ]给出了基于两点的前沿推进网格生 连接成三角网。Green 和 Sibson 首次实现了 成的基本方法 ,Eymard[2 ]给出了基于三点的前沿推进 Voronoi 图的三角网生长法[15 ] ,其基本步骤是以任 法。Joe[3 ,4 ] 和 Weatherill[5 ] 给出了 Delaunay 剖分的各 意点为起始点找到距其最近的点构成一个基本边 , 种改进方法。两类方法相比较 ,前沿推进法所用到的 然后找第三点形成 Delaunay 三角形 ,依次处理所有 数据结构较为简单 ,但浮点运算工作量较大 ,容易因 边。其找点过程按照 Delaunay 原则进行 , 故生成 硬件舍入误差引起失效。Delaunay 剖分的浮点运算 Delaunay 三角形网。Lawson[16 ] 提出了逐点插入法 主要集中在搜索插入新点所破坏三角形的过程中 ,相 的基本思想 , Bowyer 等人[17～19 ] 对其作了完善 , 对而言 ,其浮点运算较少 ,但一般需要全域枚举或采 Weatherill [5 ]对 Bowyer