快速混合网格生成方法
【技术领域】
[0001]本发明涉及一种快速混合网格生成方法。
【背景技术】
[0002]计算流体力学(CFD)是现代流体力学研宄方法中的重要手段，而网格生成技术是CFD中的关键环节。一般认为，在计算流体力学研宄中，前处理工作占到总的工作量70%以上，而前处理的主要工作为计算网格的划分。显然，发展高效的网格生成方法可以极大地提高CFD研宄的效率。
[0003]按照网格类型划分，CFD研宄中一般采用的计算网格可分为结构网格和非结构网格。结构网格节点变化有序，求解效率、精度较高，但处理复杂外形时较为繁琐；非结构网格节点和单元形成具有随意性，处理外形时比较灵活，更适合复杂外形，根据单元几何结构的不同，非结构网格又可分为四面体网格、直角网格以及混合网格(包含四面体单元和六面体单元)等。
[0004]在生成计算网格过程时，首先需要在物面上布置表面网格。由于结构网格和四面体非结构网格的空间网格质量强烈依赖于表面网格的布置，且有的复杂部位空间网格的生成要求表面网格需要仔细设计，因此这些类型计算网格生成过程中有大量的人力耗费在表面网格设计过程中。一般的混合网格是由结构网格和四面体网格、六面体网格结合而来的，其表面网格生成的难度较结构网格要低，但其工作量在整个网格生成过程中仍占据很大的比重。
[0005]直角网格的生成过程与前述各类型网格不同。一般的结构网格、四面体网格等首先需要生成高质量的表面网格，空间网格的生成是通过专门的算法从物面向空间推进生成的。直角网格的空间网格生成过程则与之相反，直角网格生成时虽然也需要输入表面网格，但表面网格仅用于描述外形，空间网格的质量与表面网格的质量不相关，表面网格的生成时间也与物体拓扑结构的复杂度不相关，因而表面网格的生成在直角网格的生成过程中的耗时占比较小。由于直角网格生成时的输入参数少，人工干预需求低，算法鲁棒性高，容易实现网格生成的自动化，因而在很多对网格生成自动化要求较高的方向应用较多。
[0006]由于直角网格各网格面、边完全正交，而需要拟合的物形可能为任意形状，原始的直角网格在物面边界处往往不能达到完全贴体，因此如何处理物面边界成为直角网格应用中最为关键的技术。主流边界处理方法可分为两类，一类是非贴体的方式，即在处理物面边界时不改变网格形状，只是通过各种数值手段计算边界附近网格的各物理量；另一类为贴体的方式，即在处理边界时将流场内的网格点光顺、投影，生成贴体的网格。目前主流的计算网格，如结构网格、四面体网格等，都采用贴体的方式，因此贴体的直角网格由于方便继承相关的程序，因而应用较为广泛。
[0007]当采用贴体的方式处理直角网格物面边界时，通常需要对物面附近的单元进行光顺和投影处理。但光顺过程会造成附近的立方体单元出现扭转、拉伸、偏斜等现象，造成网格质量的下降。
【发明内容】
[0008]本发明解决的技术问题是:针对现有技术的不足，提供了一种快速混合网格生成方法，实现了对空间网格的快速剖分。方法保留了直角网格系统中空间网格自动生成的特性，且通过原始网格的分割，降低了扭曲单元的比例，提高了计算网格的整体质量。
[0009]本发明的快速混合网格生成方法包括:步骤一，网格生成前处理；步骤二，生成八叉树结构直角网格；步骤三，利用标准分割模版对直角网格进行分割，生成混合网格；步骤四，计算混合网格的网格信息，并输出网格信息。
[0010]优选所述步骤三包括:设计各种类型网格的标准分割模板的步骤；根据所述八叉树结构直角网格的实际情况，选择不同类型的所述标准分割模板进行分割，生成混合网格的步骤。
[0011]优选所述设计各种类型网格的标准分割模板的步骤按照如下方式进行:(I)按照同一基本类型单元内不存在彼此相邻的单元的方式，将所述八叉树结构直角网格的单元分为两种基本类型单元，g卩，将该单元所属父单元占据的空间按照相对父单元体心位置的卦限分布进行排序，如果该单元处于第一、三、六、八卦限，则此单元归为第一类基本单元，如果该单元位于第二、四、五、七卦限，则此单元归类为第二类基本单元。设计确定该两种所述基本类型单元的标准分割模板；(2)分别利用所述标准分割模板，对两种所述基本类型单元进行分割，将所述网格进一步分为多个次一级类型单元；(3)根据每个所述次一级类型单元上节点的数量和分布情况，再设计确定各类型次一级类型单元的标准分割模板。
[0012]优选根据每个所述次一级类型单元上节点的分布情况，所述次一级类型单元能够分为:节点类型单元，仅在网格单元的边节点位上存在网格节点，但面心节点位处无网格节点；面类型单元，仅在面心节点位处存在网格节点的单元；混合类型单元，节点类型单元和面类型单元的混合形式。
[0013]优选针对所述次一级类型单元，分别根据节点的数量和分布情况设计出各次一级类型单元的标准分割模板。
[0014]优选针对所述步骤二生成的所述八叉树结构直角网格单元:首先判断其分属于两个所述基本类型单元的哪一个所述基本类型，利用相应的所述基本类型标准分割模板对其进行分割，再判断分割后形成的次一级单元分属于哪个所述次一级类型单元，并套用相应的所述次一级类型单元的标准分割模板。
[0015]本发明与现有技术相比具有如下优点:
[0016](I)本发明采用模版式分割对网格单元进行分割，分割过程中直接套用模版，不用求解代数方程或微分方程，生成效率高，网格质量高，算法稳定性好；
[0017](2)本发明采用旋转变换的方式，将各类型网格中不同的排列组合变换为标准形式，之后套用标准模版对其进行分割。此方法简化了模版的设计，简化了方法的算法，提高了方法的适用性；
[0018](3)本发明基于直角网格生成混合网格，直角网格生成过程以及网格划分过程都能实现自动化，不需要过多的人工干预；
[0019](4)本发明基于直角网格，直角网格在生成过程中不需要仔细设计表面网格，空间网格生成的算法简单，网格生成效率高，可大大降低网格生成在CFD研宄中的时间占比，提高CFD应用的整体效率。
【附图说明】
[0020]图1为直角单元网格27个节点位分布情况的示意图。
[0021]图2为两类基本单元的标准分割模版，其中图2(1)为第一类基本单元的分割模版，图2 (2)为第二类基本单元的分割模板。
[0022]图3为单节点类型的标准分割模版。
[0023]图4为双节点类型的标准分割模版。
[0024]图5为三节点类型的标准分割模版。
[0025]图6为单面类型的标准分割模版。
[0026]图7为具有单面类型的直角网格分割后示意图。
[0027]图8为双面类型的标准分割模版。
[0028]图9为具有双面类型的直角网格分割后不意图。
[0029]图10为三面类型的标准分割模版。
[0030]图11为具有二面类型的直角网格分割后不意图。
【具体实施方式】
[0031]下面将对本发明做进一步的介绍。
[0032]为了下文表述方便，首先以一个简单立方体单元为例，介绍本方法中采用的网格节点编号规则。每个立方体单元有八个顶点，如图1中A-Η。对于某些单元，除了八个顶点夕卜，在各边的中点、各面的面心位置还可能存在悬空节点，有的单元其体心位置也会存在网格节点。因而对于任意单元，在其单元内部可能存在的节点数最多为27个(8个顶点+12个边中点+6个面心+1个体心)，图1中黑色圆点即为27个可能存在节点的位置，这里称之为27个节点位。
[0033]在下面的介绍中还需要用到卦限，这里采用的卦限定义方式为:在笛卡尔坐标系下，包含Χ、γ、ζ正半轴的卦限为第一卦限，二-四卦限处在XOY面以上，按逆时针方向排列，五-八卦限则分列于一-四卦限之下。以图1中单元为例，若将坐标系原点置于单元的体心，八个顶点A-H排列顺序与1-8卦限排列顺序一致。
[0034]在应用基于标准分割模版的快速混合网格生成方法时，主要包括以下几个部分:
[0035]步骤1:直角网格生成前处理
[0036]网格生成的前处理工作主要包括网格生成参数的设置以及表面网格的生成。直角网格生成时，需要设置的参数有:计算域大小、加密层数、网格尺度、分层加密、局部加密等。直角网格生成中，主要靠表面网格描述物体的外形，此表面网格可以是三角形、四边形或其他任意多边形。
[0037]步骤2:生成八叉树结构的直角网格
[0038]根据输入的背景网格和设置参数，首先在计算域内划分出初始网格，之后对物体附近的单元进行多层八叉树加密，生成符合设置条件的叉树结构网格。删除物体内部以及与物体相交的网格单元，产生锯齿形