非结构化四面体网格生成方法及改进.pdfVIP

维普资讯 第 38卷 第 1期 西 安 建 筑 科 技 大 学 学 报(自然科学版) Vo1．38 No．1 2006年 2月 J．Xi’anUniv．ofArch．Tech+(NaturalScienceEdition) Feb．2006 非结构化四面体网格生成方法及改进 王昊利 ，王 元 ，岳斌佑 (1．西安交通大学生态环境与现代农业工程中心，陕西 西安 710049； 2．西安建筑科技大学建筑设计研究院，陕西 西安 710055) 摘 耍：对生成三维非结构化网格的Delaunay方法进行了分析，给出四面体网格生成的基本步骤．改进判断 非结构化网格质量的方法。提出采用四面体单元无量纲棱长的标准差对网格质量进行评价．对于已经生成的 初始化网格，在棱长标准差最大的四面体单元内加入新点重新构造网格，反复进行迭代直到所有的标准差小 于设定值．该方法 比较灵活且易于实现．通过不同的算例对网格生成及改进方法进行了验证，获得了理想的结 果． 关键词：非结构化网格；四面体单元；Delaunay方法；标准差 中田分类号：O351 文献标识码 ：A 文章编号：1006—7930(2006)010058—05 网格生成技术是流场数值计算的第一步，是计算流体动力学的一个重要组成部分．目前，提出并得 以应用的网格生成技术有多种方法，总体来讲分为结构化网格和非结构化网格两大类．非结构化网格 比 结构化网格生成方法复杂，研究的充分程度也不如结构化网格．相 比而言，对三维非结构化 网格的研究 又远不如对二维非结构化网格充分．本文针对三维非结构化网格生成中的关键环节，提出了自己的改进 方法． 非结构化网格生成方法主要有前沿推进法、Delaunay三角形(二维平面)或四面体 (三维空间)方法 以及其他一些方法，如八叉树方法等I1．针对三维空间的四面体网格的生成方法文献[2]采用的是前沿 推进法．这种方法须首先在几何体表面生成三角形网格，然后采用前沿推进法逐渐在几何体 内部生成网 格．目前，绝大多数的非结构化网格生成采用另外一种方法，即Delaunay三角形或四面体方法[3]．采用 Delaunay方法对非结构化网格 的研究最早始于 Bowyer[4]和 Watson_l5的工作．不论对二维平面的三角 形网格还是对三维空间的四面体网格，采用Delaunay方法生成网格的基本原理相同，但三维四面体非 结构化网格生成算法远没有二维问题成熟，健壮性也要差得多．对于二维问题，生成三角形网格需要的 信息包括顶点、边、单元面积以及外接圆圆心和半径，而对于三维问题的四面体网格需要的信息包括了 顶点、棱、面、单元体积以及外接球球心和半径，存储的信息量与计算量显著增加，如何加快三维非结构 化 网格 的生成速度是三维非结构化网格研究中的一个重要课题。在非结构化网格生成 中，判断并改善 网 格单元质量是一个主要环节．对此二维问题有较为成熟的作法，但推广于三维中就显得相当复杂．由于 三维非结构化网格生成很大程度上会增加问题的复杂性 ，从文献来看 。对三维非结构化 网格的研究要远 远少于二维非结构化网格．本文采用的生成四面体网格的方法为Delaunay化方法，在区域内部通过构 造 Delaunay四面体进行网格生成，并在改善网格质量的算法上做了改进． 收稿 日期；2004—10—11 修改稿 日期；2005—05—16 基金项目：国家 自然科学基金项 目(1027089)I教育部博士点基金项 目(20030698004) 作者简介 ：王昊利 (1972一)．男．山西话墟 人．憧+研育堆 ．丰尊 “宴晡 T稿音而 珥 峦 维普资讯 第 1期 王吴利等：非结构化四面体网格生成方法及改进 59 1 Delaunay方法 1．1 Delaunay四面体网格生成原理 Delaunay方法建立在严格的计算几何原理上，其基本原理源于 19世纪 5O年代Diric