CN104732589B - 快速混合网格生成方法 - Google Patents

Abstract

本发明公开了一种快速混合网格生成方法，其包括步骤一，网格生成前处理；步骤二，生成八叉树结构直角网格；步骤三，利用标准分割模版对直角网格进行分割，生成混合网格；步骤四，计算混合网格的网格信息，并输出网格信息。本发明解决了计算流体力学中快速生成混合网格网格的问题，采用标准模版式分割的方法，不需求解代数方程或微分方程，即可实现混合网格无人工干预的自动生成。

Description

快速混合网格生成方法

技术领域

本发明涉及一种快速混合网格生成方法。

背景技术

计算流体力学(CFD)是现代流体力学研究方法中的重要手段，而网格生成技术是CFD中的关键环节。一般认为，在计算流体力学研究中，前处理工作占到总的工作量70％以上，而前处理的主要工作为计算网格的划分。显然，发展高效的网格生成方法可以极大地提高CFD研究的效率。

按照网格类型划分，CFD研究中一般采用的计算网格可分为结构网格和非结构网格。结构网格节点变化有序，求解效率、精度较高，但处理复杂外形时较为繁琐；非结构网格节点和单元形成具有随意性，处理外形时比较灵活，更适合复杂外形，根据单元几何结构的不同，非结构网格又可分为四面体网格、直角网格以及混合网格(包含四面体单元和六面体单元)等。

在生成计算网格过程时，首先需要在物面上布置表面网格。由于结构网格和四面体非结构网格的空间网格质量强烈依赖于表面网格的布置，且有的复杂部位空间网格的生成要求表面网格需要仔细设计，因此这些类型计算网格生成过程中有大量的人力耗费在表面网格设计过程中。一般的混合网格是由结构网格和四面体网格、六面体网格结合而来的，其表面网格生成的难度较结构网格要低，但其工作量在整个网格生成过程中仍占据很大的比重。

直角网格的生成过程与前述各类型网格不同。一般的结构网格、四面体网格等首先需要生成高质量的表面网格，空间网格的生成是通过专门的算法从物面向空间推进生成的。直角网格的空间网格生成过程则与之相反，直角网格生成时虽然也需要输入表面网格，但表面网格仅用于描述外形，空间网格的质量与表面网格的质量不相关，表面网格的生成时间也与物体拓扑结构的复杂度不相关，因而表面网格的生成在直角网格的生成过程中的耗时占比较小。由于直角网格生成时的输入参数少，人工干预需求低，算法鲁棒性高，容易实现网格生成的自动化，因而在很多对网格生成自动化要求较高的方向应用较多。

由于直角网格各网格面、边完全正交，而需要拟合的物形可能为任意形状，原始的直角网格在物面边界处往往不能达到完全贴体，因此如何处理物面边界成为直角网格应用中最为关键的技术。主流边界处理方法可分为两类，一类是非贴体的方式，即在处理物面边界时不改变网格形状，只是通过各种数值手段计算边界附近网格的各物理量；另一类为贴体的方式，即在处理边界时将流场内的网格点光顺、投影，生成贴体的网格。目前主流的计算网格，如结构网格、四面体网格等，都采用贴体的方式，因此贴体的直角网格由于方便继承相关的程序，因而应用较为广泛。

当采用贴体的方式处理直角网格物面边界时，通常需要对物面附近的单元进行光顺和投影处理。但光顺过程会造成附近的立方体单元出现扭转、拉伸、偏斜等现象，造成网格质量的下降。

发明内容

本发明解决的技术问题是：针对现有技术的不足，提供了一种快速混合网格生成方法，实现了对空间网格的快速剖分。方法保留了直角网格系统中空间网格自动生成的特性，且通过原始网格的分割，降低了扭曲单元的比例，提高了计算网格的整体质量。

本发明的快速混合网格生成方法包括：步骤一，网格生成前处理；步骤二，生成八叉树结构直角网格；步骤三，利用标准分割模版对直角网格进行分割，生成混合网格；步骤四，计算混合网格的网格信息，并输出网格信息。

优选所述步骤三包括：设计各种类型网格的标准分割模板的步骤；根据所述八叉树结构直角网格的实际情况，选择不同类型的所述标准分割模板进行分割，生成混合网格的步骤。

优选所述设计各种类型网格的标准分割模板的步骤按照如下方式进行：(1)按照同一基本类型单元内不存在彼此相邻的单元的方式，将所述八叉树结构直角网格的单元分为两种基本类型单元，即，将该单元所属父单元占据的空间按照相对父单元体心位置的卦限分布进行排序，如果该单元处于第一、三、六、八卦限，则此单元归为第一类基本单元，如果该单元位于第二、四、五、七卦限，则此单元归类为第二类基本单元。设计确定该两种所述基本类型单元的标准分割模板；(2)分别利用所述标准分割模板，对两种所述基本类型单元进行分割，将所述网格进一步分为多个次一级类型单元；(3)根据每个所述次一级类型单元上节点的数量和分布情况，再设计确定各类型次一级类型单元的标准分割模板。

优选根据每个所述次一级类型单元上节点的分布情况，所述次一级类型单元能够分为：节点类型单元，仅在网格单元的边节点位上存在网格节点，但面心节点位处无网格节点；面类型单元，仅在面心节点位处存在网格节点的单元；混合类型单元，节点类型单元和面类型单元的混合形式。

优选针对所述次一级类型单元，分别根据节点的数量和分布情况设计出各次一级类型单元的标准分割模板。

优选针对所述步骤二生成的所述八叉树结构直角网格单元：首先判断其分属于两个所述基本类型单元的哪一个所述基本类型，利用相应的所述基本类型标准分割模板对其进行分割，再判断分割后形成的次一级单元分属于哪个所述次一级类型单元，并套用相应的所述次一级类型单元的标准分割模板。

本发明与现有技术相比具有如下优点：

(1)本发明采用模版式分割对网格单元进行分割，分割过程中直接套用模版，不用求解代数方程或微分方程，生成效率高，网格质量高，算法稳定性好；

(2)本发明采用旋转变换的方式，将各类型网格中不同的排列组合变换为标准形式，之后套用标准模版对其进行分割。此方法简化了模版的设计，简化了方法的算法，提高了方法的适用性；

(3)本发明基于直角网格生成混合网格，直角网格生成过程以及网格划分过程都能实现自动化，不需要过多的人工干预；

(4)本发明基于直角网格，直角网格在生成过程中不需要仔细设计表面网格，空间网格生成的算法简单，网格生成效率高，可大大降低网格生成在CFD研究中的时间占比，提高CFD应用的整体效率。

附图说明

图1为直角单元网格27个节点位分布情况的示意图。

图2为两类基本单元的标准分割模版，其中图2(1)为第一类基本单元的分割模版，图2(2)为第二类基本单元的分割模板。

图3为单节点类型的标准分割模版。

图4为双节点类型的标准分割模版。

图5为三节点类型的标准分割模版。

图6为单面类型的标准分割模版。

图7为具有单面类型的直角网格分割后示意图。

图8为双面类型的标准分割模版。

图9为具有双面类型的直角网格分割后示意图。

图10为三面类型的标准分割模版。

图11为具有三面类型的直角网格分割后示意图。

具体实施方式

下面将对本发明做进一步的介绍。

为了下文表述方便，首先以一个简单立方体单元为例，介绍本方法中采用的网格节点编号规则。每个立方体单元有八个顶点，如图1中A-H。对于某些单元，除了八个顶点外，在各边的中点、各面的面心位置还可能存在悬空节点，有的单元其体心位置也会存在网格节点。因而对于任意单元，在其单元内部可能存在的节点数最多为27个(8个顶点+12个边中点+6个面心+1个体心)，图1中黑色圆点即为27个可能存在节点的位置，这里称之为27个节点位。

在下面的介绍中还需要用到卦限，这里采用的卦限定义方式为：在笛卡尔坐标系下，包含X、Y、Z正半轴的卦限为第一卦限，二-四卦限处在XOY面以上，按逆时针方向排列，五-八卦限则分列于一-四卦限之下。以图1中单元为例，若将坐标系原点置于单元的体心，八个顶点A-H排列顺序与1-8卦限排列顺序一致。

在应用基于标准分割模版的快速混合网格生成方法时，主要包括以下几个部分：

步骤1：直角网格生成前处理

网格生成的前处理工作主要包括网格生成参数的设置以及表面网格的生成。直角网格生成时，需要设置的参数有：计算域大小、加密层数、网格尺度、分层加密、局部加密等。直角网格生成中，主要靠表面网格描述物体的外形，此表面网格可以是三角形、四边形或其他任意多边形。

步骤2：生成八叉树结构的直角网格

根据输入的背景网格和设置参数，首先在计算域内划分出初始网格，之后对物体附近的单元进行多层八叉树加密，生成符合设置条件的叉树结构网格。删除物体内部以及与物体相交的网格单元，产生锯齿形内面，并对内面进行光顺，最终获得光顺后的内面。至此，获得了模版划分所需的八叉树结构的直角网格。

步骤3：利用标准分割模板对直角网格进行划分，生成混合网格

此步骤是本发明的主要部分。本发明的核心思想是采用标准分割模版的方式快速分割网格单元，其中重点和难度就在于网格标准分割模版的设计。本发明首先设计基本类型单元的标准分割模版，并在基本类型标准分割模版的基础上，设计不同类型单元的标准分割模版。

此步骤可分解为以下几个过程：

(1)确定基本单元类型。为了保证分割后各网格面能够一一对应，在本发明中将所有网格分类为两种基本网格类型，相邻的单元必须分属于不同的类型。不同基本网格类型采用不同的分割方式，保证相邻网格公共面在分割后仍能保证几何上的一一对应关系。直角网格中所有网格单元都是由父单元一分为八而得到的(第一层单元可以看成是虚拟的第0层单元加密得到)，根据八个兄弟单元所处位置，将处在标准笛卡尔坐标系第一、三、六、八卦限的单元归为第一类基本单元，处在第二、四、五、七卦限的单元归类为第二类基本单元。

(2)采用基本类型单元的标准分割模板分割八叉树结构直角网格。为了保证分割过程的简洁，同时控制分割后混合网格的网格量，必须保证分割后不产生新的网格节点，即所有新网格的网格节点仍为原网格的八个顶点。图2所示即为本发明中采用的两类基本单元的分割模版，采用此模版可以将立方体网格分割为五个四面体网格。由于同一层中第一类与第二类基本单元是交错排列的，因而可以保证分割后四面体单元面是一一对应的。

(3)建立次一级类型单元的标准分割模版。

次一级类型网格的分割模版需要根据节点位网格节点的分布情况确定。27个节点位上各节点的分布多种多样，如果对其进行排列组合，则需要确定成千上万种模版。但由于立方体网格以及基本单元分割模版的对称性，可以确定一系列标准模版，其他类似的情况可以通过单元的旋转转化为标准模版形式，进而按照标准模版进行分割。下面将介绍针对各种情况设计的标准模版。

首先分析节点类型单元。设定次一级单元中边节点为存在网格节点，而面心节点位不存在节点的情况为节点类型单元。根据边节点的数目，节点类型单元又可分为单节点类型、双节点类型和三节点类型。选取第二类基本模版中顶点E所在四面体为例演示节点类型标准模版的切割方法，这里称此四面体为四面体E。单节点类型中四面体E的标准切割模版如图3所示，如果边节点处于另两条边的节点位上，可通过旋转变换至图3所示的情况，从而直接套用单节点类型的标准切割模板。双节点类型和三节点类型的标准分割模板如图4和图5所示。

然后考虑面节点类型。次一级单元面心节点位存在节点的情况称为面类型。当单元面心节点位存在节点时，表明分割前的直角网格在方向的邻居单元为四个处于更密一层网格的单元，因而除了面心节点位存在网格节点，此面对应的各条边的边中点节点位也存在网格节点。对于仅有面心节点位存在网格节点的类型称为单面类型，图6中给出了四面体E的单面类型分割模板。单面类型一般为成对出现，如果四面体E为单面类型单元，则四面体G也必为面类型单元。四面体G可以通过旋转变换，再套用单面类型标准分割。图7所示为经过基本类型模板和单面类型模板分割后的直角网格单元。

对于有两个面的面心节点位存在网格节点的类型称为双面类型，图8中所示为双面类型标准模版。当次一级单元存在双面类型单元时，则分割前的直角网格至少在两个方向上存在加密的邻居单元，其分割后的次一级单元中存在一个双面类型单元和两个单面类型单元，图9中为套用各种标准模板分割后的直角网格示意图。

有三个面的面心节点位存在网格节点的类型称为三面类型，图10为三面类型标准切割模版，由于三面类型具有对称性，因此图10中仅给出了一侧的分割方法，另一侧(点划线表示的部分)的分割方式与之对称。当次一级单元存在三面类型单元时，则分割前的直角网格至少在三个方向上存在加密的邻居单元，其分割后的次一级单元中除了一个三面类型单元外，还存在三个单面类型单元，图11中所示为套用各种标准分割模板切割后的直角网格单元。

除了上述各种情况外，还可能存在一些特殊情况，如直角网格单元相对的两个面面心存在网格节点单元，这表明该网格单元在两个相对的方向上都被加密，为了保证网格质量，应该将此单元也进行八叉树加密。

(4)单元的分割。上面仅针对各类型的典型网格建立了相应的标准分割模版，但实际的网格类型可能多种多样，此时需要应用到网格的旋转变换。由于直角网格中各网格单元为立方体，各节点位的排列也满足对称关系，因此任何网格类型经过旋转变换都能有相应的标准模版与之对应。根据网格的类型，经过适当的旋转变换，可以将各网格变换为与标准模版对应的形式，变换后的网格套用标准分割模版，即可获得以四面体单元、四棱锥单元为主的混合网格。

步骤4：计算和输出网格信息

方法的最后是计算混合网格的网格信息，并输出成网格文件供CFD计算使用。涉及的网格信息主要包括网格点位置，网格面大小、指向、面心位置，网格单元体积、中心位置以及单元点、线、面、体之间的对应关系等。

本发明未详细说明部分属本领域技术人员公知常识。

以上对本发明的优选实施方式进行了说明，但本发明并不限定于上述实施例。对本领域的技术人员来说，在权利要求书所记载的范畴内，显而易见地能够想到各种变更例或者修正例，当然也属于本发明的技术范畴。

Claims (4)

1.一种快速混合网格生成方法，其特征在于，包括：

步骤一，网格生成前处理；

步骤二，生成八叉树结构直角网格；

步骤三，利用标准分割模版对直角网格进行分割，生成混合网格；

步骤四，计算混合网格的网格信息，并输出网格信息，

所述步骤三包括：设计各种类型网格的标准分割模板的步骤；根据所述八叉树结构直角网格的实际情况，选择不同类型的所述标准分割模板进行分割，生成混合网格的步骤，

所述设计各种类型网格的标准分割模板的步骤按照如下方式进行：

(1)按照同一基本类型单元内不存在彼此相邻的单元的方式，将所述八叉树结构直角网格的单元分为两种基本类型单元，即，将该单元所属父单元占据的空间按照相对父单元体心位置的卦限分布进行排序，如果该单元处于第一、三、六、八卦限，则此单元归为第一类基本单元，如果该单元位于第二、四、五、七卦限，则此单元归类为第二类基本单元；设计确定该两种所述基本类型单元的标准分割模板；

(2)分别利用所述标准分割模板，对两种所述基本类型单元进行分割，将所述网格进一步分为多个次一级类型单元；

(3)根据每个所述次一级类型单元上节点的数量和分布情况，再设计确定各类型次一级类型单元的标准分割模板。

2.如权利要求1所述的快速混合网格生成方法，其特征在于：根据每个所述次一级类型单元上节点的分布情况，所述次一级类型单元能够分为：节点类型单元，仅在网格单元的边节点位上存在网格节点，但面心节点位处无网格节点；面类型单元，仅在面心节点位处存在网格节点的单元；混合类型单元，节点类型单元和面类型单元的混合形式。

3.如权利要求2所述的一种快速混合网格生成方法，其特征在于：针对所述次一级类型单元，分别根据节点的数量和分布情况设计出各次一级类型单元的标准分割模板。

4.如权利要求3所述的快速混合网格生成方法，其特征在于：针对所述步骤二生成的所述八叉树结构直角网格单元：首先判断其分属于两个所述基本类型单元的哪一个所述基本类型，利用相应的所述基本类型标准分割模板对其进行分割，再判断分割后形成的次一级单元分属于哪个所述次一级类型单元，并套用相应的所述次一级类型单元的标准分割模板。