一种改进的二维平面区域三角形化的前沿推进法.pdfVIP

( ) 2003 年 　第 27 卷 　　　　　　　　　石油大学学报 自然科学版 　　　　　　　　　　Vol . 27 　No . 5 　第 5 期 　　　　　　　　Journal of t he University of Petroleum , China 　　　　　　　　Oct . 2003 　　文章编号 (2003) 一种改进的二维平面区域三角形化的前沿推进法 1 1 1 2 修荣荣 , 徐明海 , 黄善波 , 陶文铨 ( 1. 石油大学储运与建筑工程学院 ,山东东营 25706 1 ; 2 . 西安交通大学能源与动力工程学院 ,西安 7 10049) 　　摘要 :结合两点前沿推进和三点前沿推进法 ,提出了一种改进的二维平面区域前沿推进式三角形网格生成算 法 。交替使用两点前沿的生点连点和三点前沿的补充连点方法 ,避免了两点前沿算法中许多重复无效的操作 ,提高 了算法的计算效率 , 同时网格保留了两点前沿推进法的局部最优特点 。网格剖分实例计算表明 ,当网格单元数很多 时 ,该文中提出的方法较两点前沿推进法省时 50 %以上 ,提高了质量优 良的网格单元比例 ,并可以剔除质量极差的 网格单元 。利用背景网格信息可生成各向异性的、贴体性较好的网格 ,并可以保证第一层网格节点至边界的距离基 本相等 。 关键词 :非结构网格 ;三角形网格生成算法 ;前沿推进法 ; 网格剖分 中图分类号 : T K 124 , TP 39 1. 4 1 　　　文献标识码 :A 法补充连点的方法 。 引　言 1 　算法介绍 近几十年来 , 由于非结构化网格对复杂区域边 界的良好适应性 、易于局部加密以及在构造 自适应 1 . 1 　算法思想与步骤 网格方面的灵活性 ,使其在数值计算中得到了广泛 笔者提出的算法的基本步骤与两点前沿推进法 应用 。为了提高剖分效率和改善网格质量 ,对二维 相似 。首先把待剖分区域按照网格步长要求 ,剖分 区域的非结构化网格的剖分 , 以前沿推进法和 De 成首尾相接的线段 ,作为初始前沿 。内部单元的生 launay 算法为基础 , 提出了许多改进方法 。Delau 成分生点和连点两步 。图 1 是两点前沿法中某一步 nay 方法理论成熟 ,生成的网格系统具有全局最优 的前沿分布 ,按两点法的规则 ,先生成一个候选点 , 的特点 ,但这种算法对实数运算的舍入误差极为敏 然后搜索近距离的前沿点 ,进行排序和相交检验判 感 ,易产生非正常中止问题 。前沿推进法是一种经 断 ,若以 A B ( 图 2 (a) ) 为基边 ,则生成的点位于 C2 验方法 ,又可分为两点前沿推进法和三点前沿推进 点附近, 不满足插入点的要求, 则放弃生成的候选 法 。两点前沿推进法[ 1 ] 是指在选定前沿边的垂直 点 。从图2 中可以看出, 许多相邻前沿的夹角是锐 ( ) 平分线上生成内点 ,生成的三角形接近于正三角形 角 如 C A 与A B 的夹角 , 只要将 B C