计算流体力学网格划分方法的现状与发展_李石德.pdfVIP

科技信息 计算机与网络 计算流体力学网格划分方法的现状与发展 广西电力职业技术学院 李石德 ［摘 要］网格划分方法和技术的不断发展，促进了计算流体力学技术的进步。本文分析了国内外网格划分方法和技术的现状，并对 一些有代表性的方法进行评述，最后提出了网格划分方法的发展方向。 ［关键词］计算流体力学 网格划分 分解法 网格模板法 1.引言 这一类方法分两步进行：（1）节点排布；（2）单元生成。节点排布可 计算流体力学的任务是流体力学的数值模拟,建立数学模型以后需 分为随机和非随机技术。Cavendish 的算法在分析域上覆盖栅格，在栅 要解决的问题是寻求高效率、高准确度的计算方法。计算方法不仅包括 格单元内随机生成合乎节点密度分布要求的节点来完成节点排布。在 数学方程的离散化及求解方法，还包括计算网格的建立、边界条件的处 非随机技术中，Lo’s 和Lee 各自提出了自己的方法。在单元生成阶段， 理，网格的划分属于CFD 计算中的预处理部分。网格的生成是将计算 只有Lee 的方法生成四边形单元。其余的方法如Lewis 和Robison 的方 域划分为较小的、不重叠的子域或单元（网格），单元的数目决定了CFD 法，Frederick 的方法，Delaunay 三角剖分方法和Nelson 的方法等均生成 的计算精度，而网格的细密程度决定了解的精度、必需的计算机硬件和 三角形单元。 计算时间，最佳网格多是非均匀的：点到点之间变化快的区间网格较 Delaunay 三角化是计算几何的重要研究内容之一，由其生成的 细，变化相对较慢的区间网格较粗。 Delarnay 三角网格，具有空外接圆、最小角最大等众多特性，是计算几何 网格划分是计算流体力学中的关键技术之一,网格划分是将结构进 的重要几何结构之一。近年来，也被成功地应用于模式识别过程中的模 行离散化，从而得到有限单元及节点。总的来说，网格划分可分为三大 式匹配，比如指纹细节点匹配，受到人们的密切关注。针对任意平面多 部分：分解法、网格模版法以及综合法。 边形域，采用增量思想和均匀网格，可在局部范围内快速生成约束De- 2.分解法 launay 三角形。丁永祥使用“修正的外接圆”准则来实现任意边界条件 2.1 拓扑分解法 下的Delaunay 三角剖分，而卢朝阳则实现了内部也带约束的Delaunay 由Wordenweber 提出的2D 分解技术中，一个对象是一个由顶点和 三角剖分。 边为代表的可能包括曲线的多边形。每个实体具有一个指标其位置或 3.2 基于栅格的技术 形状的几何信息和一个指示与其它实体有关的拓扑信息。顶点的几何 这类技术近年来国内研究的较多，并取得了一些成果。代表性的一 信息是其坐标，拓扑信息由一个其附属的边表组成，边的几何信息是它 些方法有四叉树法、九分法。 所在的曲线，拓扑信息由一个两端点的列表组成。这种技术通过连结对 四叉树法采用树形数据结构，其基本思想是先找出问题域的最小 象顶点形成三角形而将对象分解成一系列粗略的单元。单元的尺寸和 包络矩形，再把矩形细分为四个全等的子矩形，依次判断每个子矩形与 形状不能被真正考虑到，因为对象被分解是由拓扑关系决定的，这些粗 问题域的关系，根据子矩形与问题域的包络关系来决定取舍或是否继 略的单元必须被细化，以达到网格度要求。 续细分下去。四叉树法生成的网格中间单元形状较好，对任意边界适应