基于直角切割点调整的二维网格生成方法.docVIP

基于直角切割点调整的二维网格生成方法 　　摘要：目前网格生成占用了流场计算的大量时间，为了进一步缩短网格生成的计算时间，该文在直角网格上直接进行快速的网格调整，针对多种边界情况进行了分类，保持了对复杂边界良好的适应性，同时在网格质量可以接受的范围内，可以在与网格单元数量成线性增长的时间复杂度下生成网格。最后以NACA0012翼型为例给出计算结果，验证了算法的有效性。 　　关键词：网格生成；直角网格；网格边界调整 　　中图分类号：TP18 文献标识码：A 文章编号：1009-3044（2015）22-0164-03 　　Abstract： Currently， mesh generation for CFD takes up considerable time of the whole calculation. To further reduce the time spent on mesh generation， a new method has been proposed to make adjustment quickly on Cartesian mesh for complex shapes. The mesh is adjusted according to different situations of mesh boundary， so it performs well on complex boundary. The quality of mesh will be in an acceptable range and the time complexity is linear to the number of grid cells. An example of generation on NACA0012 airfoil is given to confirm the effectiveness of the algorithm. 　　Key words： mesh generation； Cartesian mesh； mesh boundary adjustment 　　网格生成技术，是计算流体力学（CFD）的关键技术之一。网格生成的过程，把连续的计算区域划分成许多子区域，以便进行后续的计算。在目前的CFD计算周期中，网格生成占用了很大一部分时间，并且后续的计算精度和计算效率也很大程度上依赖于网格生成的质量和方法。 　　生成的网格从拓扑结构上分，主要分为结构化网格和非结构化网格。结构化网格数据结构简单，生成速度快。结构化网格在网格边界处通常使用适体坐标法[1]，但对于复杂外形的适应性比较差。非结构化网格对于复杂区域的适应性较强，但网格质量不好控制，计算效率相对结构化网格要低。现在非结构网格通常采用带约束的Delaunay三角剖分的方法[2]或四叉树方法[3]。针对不同网格的不同优点及缺点，混合网格应用的也越来越广泛[4]。 　　为了进一步提升网格生成的计算效率，本文所采用了一种在直角网格上调整顶点来适应边界的生成方法。这样做的优点是：与背景的结构化网格衔接时不需要做特殊处理；在边界附近对被切割的网格分情况讨论去调整网格，利用直角网格生成速度快的优势，避免了复杂度较高的Delaunay三角剖分。 　　本文的方法主要分为3个步骤，生成边界网格与计算交点，调整交点位置，生成最终网格，最后以NACA0012翼型为例，给出网格的生成结果。 　　1 光滑边界的网格生成 　　Wenjun Ying和Wei-Cheng Wang的方法[5]可以直接应用到光滑边界的网格生成中。在该方法中，所有边界与网格的交点都被称为控制点。与网格点距离近的称为主要控制点，与网格点距离远的称为次要控制点。 　　将结构网格上顶点 坐标移动到对应的主要控制点，就完成了网格的调整，如图1（b）的箭头方向所示。 　　然而这个方法并没有考虑到复杂外形边界上不光滑的尖点。所以本文在此基础上，针对尖点的边界附近的所有情况给出调整方案，把这个调整边界网格的方法推广到复杂的外形上。 　　2 尖点附近的网格生成 　　这里的“尖点”是指曲线在该点处，两侧的导数不相等。 　　2.1 尖点附近需要处理的问题 　　光滑边界的外形所有的点附近都是平滑的，而尖点两侧曲线的切向不同，如果仅仅用交点去近似，会非常不准确。 　　假设网格足够密，使得外形在3*3的网格范围内最多只有1个尖点。又因为这个顶点两端的切向量不同，所以对于方形结构网格的每个顶点，最多有2条近似的直线与顶点周围的4条边产生交点。 　　针对2.1中的顶点以及图2的6种情况，将会在2.3分别给出调整的方案。在2.2中先给出了计算前的预处理步骤。 　　2.2 网格调整前的预