不同多边形网格计算性能的比较研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：073222 不同多边形网格计算性能的比较 王振兴 徐明海 中国石油大学(华东)储运与建筑工程学院 山东 东营 25706l TEL：0546-8393759 摘要：本文首先讨论了利用三角形单元网格转化为多边形网格的方法和结果，进而分析基于多边形 网格求解不可压缩NS方程的特性，从收敛性和精确性方面评价了网格形状的影响，并提出了提高精度 的网格设计措施。 关键词：多边形网格、数值模拟、三角形网格 1、引言 对于计算流体动力学和计算传热学的研究，在某种程度上讲，其实是对N．S方程求 解方法的研究。它们所面对的首要任务就是计算区域的离散化，即网格的生成。目前， 为了复杂区域的流动和传热计算，计算中常采用三角形网格进行，取三角形为控制容积。 然而根据几何关系，在三角形网格系统中三角形的个数约为顶点数目的2倍。因此，从 计算效率的角度讲，使用三角形单元进行计算并不是最经济的。为了克服这一缺点，提 出了生成非结构化四边形网格单元，但是四边形网格单元的自动生成只是一些经验方 法，剖分效率并不高。本文以基本的三角形单元网格为基础，构造与其对偶的多边形为 控制容积，就是将三角形网格转化为任意的多边形网格，相当于文献中的点中心方法。 这样既避开了非最简形体网格单元生成的困难，还保留了非结构化网格的灵活性。由于 三角形网格中节点数目少于单元数目，因此多边形网格单元数目将远远少于三角形单元 的数目，最终使得所求解的线性代数方程组阶数大大降低，将=节省大量计算时间。 目前，把三角形单元转化为多边形单元的方式有：①形心法，将各二角形的形心相 连形成多边形；②外心法，将三角形外心相连形成多边形单元；(3)中点形心法，将三 角形的形心与边的中点相连形成多边形单元。这几种方法生成的多边形单元用于数值计 算各有利弊，形心相连的多边形单元容易生成，对三角形的质量要求不高，但是在离散 扩散项时往往产生交叉扩散项，对于方程的离散和最后求解带来很多不便，而且对计算 精度也有较大影响。外心方法形成的单元具有局部正交的特性，可以避免交叉扩散项的 产生，但要求三角形必须全部为锐角三角形，否则单元顶点可能会位于计算区域外边， 使计算区域偏离物理实际；中点形心法形成的多边形的边数是前两种方法的两倍，但构 造离散方程时相邻单元数与之不匹配，不采取特殊的处理，离散方程与形心法相同。本 文针对前两种单元进行了计算考察，比较了多边形网格在数值计算上的性能。 2、方程离散 在多边形网格离散与结构化网格上离散类似，．也包括对流项、扩散项和源项离散等。 基丁边的数据结构对控制方程进行离散只需针对每条边的两个相关单元计算扩散通量 和对流通量，遍历所有边形成离散方程。采用通用方程是求解NS方程组比较方便的方 法，所需求解的质鼙、能量和动量守恒方程都可以表示成以下通用形式： 212 ·亟璺+讲心砌砌：旃也gra韵+o (1) 甜 式中痧为通用变量，可以代表“，v，w，丁等求解变量，0为广义扩散系数，邑为 广义源项，通过给通用变量和对应的扩导以及源项赋以相应的物理意义就可表示相应的 。 方程。 针对图l所示的控制容积P和其相邻的单元E， J 其公共边IJ。穿过IJ的扩散通量和对流通量： ，， E=一r。((V妒)。·邑)+成(圪‘邑)唬(2) ‘入 耻～!、： 经过选择适当的对流项和扩散