三维非结构网格的生成和其应用.pdfVIP

第 32 卷第 1 期 南　京　航　空　航　天　大　学　学　报 V o l. 32 N o. 1 　 2000 年 2 月 Jou rn a l o f N anj in g U n iver sity o f A eron au t ics A st ron au t ics　 　F eb. 2000 三维非结构网格的生成及其应用 夏　健　　伍贻兆 (南京航空航天大学空气动力学系　南京, 2 100 16) 摘要　采用阵面推进法(A dvan cing fron t) 生成了高质量的三维四面体非结构网格。该方法的特 点是在生成单元的同时引入结点, 与其他预布结点的方法相比, 其好处是能更好地控制网格单元 尺度在整个计算域的合理分布并提高网格质量。用D e laun ay 方法连接存储有控制参数的结点生 成了背景网格。网格生成过程中采用了合理的数据结构, 有效地提高了程序运行效率。最后给出 了M 6 机翼的数值算例。 关键词: 计算空气动力学; 数值方法; 有限元法; 非结构网格; 网格自适应 中图分类号: V 2 11. 3; O 24 1; O 242 计算机技术的飞速发展为计算流体力学步入工程实用阶段提供了可能。在实际绕流中, 许多物体的外形是非常复杂的。如何有效处理复杂的物面边界, 生成高质量的计算网格, 是 目前计算流体力学一个重要的研究课题。由于结构网格在拓扑结构上相当于矩形域内的均 匀网格, 其结点定义在每一层的网格线上, 因此对于复杂外形物体要生成贴体的结构网格是 比较困难的。非结构网格结点在空间分布完全是随意的, 没有任何结构特性, 非常适于处理 [ 1 ] [2 ] 复杂边界问题。目前常用的非结构网格生成方法主要有阵面推进法 ,D elaun ay 方法 和改 进的O ct ree 方法。本文用阵面推进法生成了三维非结构网格。在阵面推进过程中自动在流 场内布点并对“最佳点”周围结点分批搜索, 进行比较选择, 既保证了网格的质量又使程序具 有较高的灵活度和适应能力。由于采用了堆栈列表( ) [3 ] 和交替二叉树( H eap list A ltern at in g d ig ita l t ree) [3 ] 等数据结构, 大大提高了程序对点、线、面的搜索效率, 缩短了网格生成时间。 最后给出了M 6 机翼的数值算例, 通过对网格进行自适应处理, 提高了解的分辨率。 1　三维非结构网格的生成 用阵面推进法生成三维非结构网格包括以下步骤: ( 1) 生成背景网格; ( ) 国家自然科学基金 编号: 重点资助项 目。 收稿 日期: 修改稿收到日期: 第一作者: 夏　健, 男, 讲师, 1970 年 3 月生。 100 南　京　航　空　航　天　大　学　学　报 第 32 卷 (2) 生成三维物体的表面网格, 形成初始阵面; (3) 引入新结点, 推进阵面, 生成空间网格。 1. 1　背景网格 在推进阵面形成新网格时, 需要了解所要生成网格的当地特性, 这些特性可以利用背景 网格, 通过线性插值得到。所谓背景网格, 是一组覆盖整个计算域的疏网格, 在它的结点处储 存着生成网格所需的控制参数, 本文用到的控制参数主要是网格的基本尺度。背景网格结点 处控制参数的选择对最终网格结点在整个流场的合理分布是至关重要的。通常背景网格是 靠人工输入数据形成的。当物体形状非常复杂时, 为了得到合理的网格结点分布, 需要输入 较多的