(6.9.1)--5.9-非结构网格与笛卡尔网格生成方法简介讲义.pdfVIP

《计算流体力学基础》讲义 第五章 网格与相应变换 5.9 非结构网格和笛卡尔网格生成方法简介 本节我们简要介绍一下非结构网格和笛卡尔网格的生成方法。之前我们介绍的都是在有 限差分方法框架下的结构网格，需要通过相应的变换将网格从物理平面变换到计算平面。实 际上，物理平面内所看到的非均匀网格在某种意义上都可以看作是一个有限体积单元。由于 有限体积方法并不要求在均匀的正交网格上进行计算，因此有限体积方法可以在物理平面内 的非均匀网格上直接进行计算，不需要任何变换，即，在有限体积方法中，网格生成只涉及 在物理平面上的网格构造。 有限体积方法中并没有要求所采用的网格是结构网格，因此它可以被用于任意形状的网 格单元，这就意味着可以采用非结构网格和笛卡尔网格。非结构网格和笛卡尔网格中没有对 应于ξ,η 的坐标线，没有严格的结构特征，可以最大限度地利用网格的灵活性，将网格布置 在需要的位置上。非结构网格生成方法主要有阵面推进法和 Delaunay 方法。笛卡尔网格的 生成方法主要是叉树法。下面以三角形非结构网格的生成为例，简要介绍阵面推进方法和 Delaunay 方法。最后介绍笛卡尔网格。 阵面推进方法是一种逐个生成网格单元的方法，其基本思路是：首先将物面离散为一系 列初始推进面元，从中选取某一最合适的面元作为当前推进面元，按照一定的准则在计算空 间生成或选取一个节点，这个节点与当前推进面元组合形成网格，同时更新推进面元系列并 重复上述步骤，直到所有流体空间区域都布满网格为止。阵面推进方法具有网格生成原理简 单、自动化程度高、网格疏密容易控制的优点；缺点是每推进一步，仅生成一个单元，效率 较低。 (a)阵面推进方法示意图 (b)阵面推进后的新网格和推进面 图A5.9.1 阵面推进方法 如图所示，运用阵面推进方法生成二维空间网格的过程可以描述如下： 1．确定当前推进边。选择当前推进面元的方法很多，有最小内角法、最小面积法、最 短高方法等，但是不论采用什么方法，其目的都是避免生成较大的网格单元而“覆盖”掉较小 的网格单元。在二维情况下，一般采用最短边原则，例如左图从当前推进面中选出最短边 - 1 - 《计算流体力学基础》讲义 ab 作为当前推进边。 2 ．计算理想生长点。对于当前推进边，先根据所处的位置获取当地的网格背景长度， 然后计算出一个理想生长点。例如左图在 ab 中垂线上生成了理想生长点p 。在没有网格背 景长度约束下，理想生长点是与ab 形成等边三角形的顶点。 3 ．搜寻候选点。从当前阵面中找出所有与理想点的距离小于一定范围，并且位于当前 推进边外法线方向一侧的节点，并将这些点作为候选点。通过综合考虑候选点到当前推进边 的距离、候选点与当前推进边所构成的网格质量等因素来对候选点进行排序。例如左图在推 进面上有一个候选点c ，与理想生长点p 一起组成了候选点序列{p，c} 。 4 ．选择最佳点。按照排好的次序逐个检查候选点与当前阵面推进边构成的三角形单元 的各个边是否与其他推进边相交。如果当前候选点通过了相交性检查，那么就取这一点作为 最佳点，与当前阵面构成下一个三角形网格单元；如果不能通过，那么就取排序列表中的下 一个候选点继续进行相交性检查；如果所有候选点都不能通过相交性检查，那么就应该对局 部的阵面和网格进行剥除，然后重新生长网格。左图中的最佳点为p 点。 5 ．记录网格和更新推进阵面。在生成一个三角形网格单元以后，须将形成三角形的三 个点按照一定的顺序存入相应的数组之中，同时更新推进阵面，即删除被新形成的网格所覆 盖掉的推进阵面单元，加入新形成的网格所带入的新阵面单元。例如左图的推进阵面abc 更 新为右图的apbc 。 6 ．重复上述步骤，直到阵面数目为零，网格生长即告完成。 7 ．最后，进行网格光顺，获得最终网格。 下图是某四段翼型的非结构网格，在翼型表面附近具有较密的网格，可以通过调整背景 网格尺度来调整局部网格疏密程度。可见，非结构网格非常适合此种复杂构型的网格生成。 图A5.9.2 某四段翼型的三角形非结构