基于局部特征尺寸场的边界层网格生成算法.pdf

Modeling and Simulation 建模与仿真建模与仿真,2023,12(4),3464-3471 Published Online July 2023 in Hans.https:/www.hanspub.org/journal/mos https:/doi.org/10.12677/mos.2023.124317 文章引用文章引用:陈林祺,周颖.基于局部特征尺寸场的边界层网格生成算法J.建模与仿真,2023,12(4):3464-3471.DOI:10.12677/mos.2023.124317 基于局部特征尺寸场的边界层网格生成算法基于局部特征尺寸场的边界层网格生成算法 陈林祺陈林祺1，周，周 颖颖2*1浙江理工大学，计算机科学与技术学院(人工智能学院)，浙江 杭州 2中国人民解放军92341部队，河南 洛阳 收稿日期：2023年4月24日；录用日期：2023年7月4日；发布日期：2023年7月12日 摘摘 要要 边界层网格的质量是影响计算流体力学中粘性计算精度的关键因素，如何在给定任意复杂模型下生成高边界层网格的质量是影响计算流体力学中粘性计算精度的关键因素，如何在给定任意复杂模型下生成高质量的边界层网格仍是计算流体力学中的瓶颈。针对边界层网格的质量问题，本文提出一种基于局部特质量的边界层网格仍是计算流体力学中的瓶颈。针对边界层网格的质量问题，本文提出一种基于局部特征尺寸场的边界层网格生成算法。该方法将表面网格顶点的最大推进高度作为局部特征尺寸场，并基于征尺寸场的边界层网格生成算法。该方法将表面网格顶点的最大推进高度作为局部特征尺寸场，并基于局部特征尺寸场构造顶点密度场，在密度场引导下进行保质量与保特征的表面三角网格重网格化，最终局部特征尺寸场构造顶点密度场，在密度场引导下进行保质量与保特征的表面三角网格重网格化，最终在新的表面三角网格上生成边界层网格。实验表明本文方法能在给定任意复杂模型下鲁棒地生成高质量在新的表面三角网格上生成边界层网格。实验表明本文方法能在给定任意复杂模型下鲁棒地生成高质量的边界层网格。的边界层网格。关键词关键词 局部特征尺寸场，重网格化，边界层网格局部特征尺寸场，重网格化，边界层网格 Boundary Layer Mesh Optimization Algorithm Based on Local Feature Size Linqi Chen1,Ying Zhou2*1School of Computer Science and Technology(School of Artificial Intelligence),Zhejiang Sci-Tech University,Hangzhou Zhejiang 292341 Troops of the Chinese Peoples Liberation Army,Luoyang Henan Received:Apr.24th,2023;accepted:Jul.4th,2023;published:Jul.12th,2023 Abstract The quality of the boundary layer mesh is a important factor affecting the accuracy of viscous cal-culations in computational fluid dynamics,and how to generate a high quality boundary layer mesh for a given arbitrary complex model is still a bottleneck in computational fluid dynamics.To *通讯作者。陈林祺，周颖 DOI:10.12677/mos.2023.124317 3465 建模与仿真 address the boundary layer mesh quality optimization problem,this paper proposes a boundary layer mesh generation algorithm based on the local feature size.The method takes the maximum advance height of the surface mesh vertices as the local feature size field,and constructs the ver-tex density field based on the local feature size,then remesh the surface triangular mesh under the guidance of the density field to preserve the quality and features,and finally generates the boundary layer mesh on the new surface triangular mesh.Experiments show that the method in this paper can robustly generate high quality boundary layer meshes for a given arbitrary com-plex model.Keywords Local Feature Size,Remesh,Boundary Layer Mesh Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 近年来，随着算力的发展，计算流体力学在航天航空，汽车等各行各业中扮演着重要的角色。计算模拟仿真的核心是用计算机来求解各种物理微分，积分方程，在具有复杂拓扑和几何的曲面或者实体上求解偏微分方程，通常采用有限元法，有限体积法或者等几何分析方法，而这些方法都强烈依赖于曲面实体的网格生成。边界层网格是一种混合网格，因其能更好地兼顾计算精度和网格易用性1，所以在计算流体力学中边界层网格被广泛采用。前沿推进法是应用最为广泛的边界层网格生成方法之一，Lu 等人2提出一种结构多块边界层网格生成方法，该方法首先提取表面网格的几何特征，基于几何特征构造边界层网格的框架，然后在框架内通过 TFI 技术生成边界层网格。曹杰等人3提出一种基于离散中轴面的前沿分块层进算法，首先基于四面体单元的外心构建三角网格内部的离散中轴面，再根据表面与中轴面的距离将边界面分为两类分别生成边界层网格。王硕等人4提出了基于径向基函数的边界层网格生成算法，通过计算表面分区好的网格点的法向构造支撑线，利用径向基函数插值构造网格的轮廓线并生成边界层网格。为解决鲁棒性与网格质量问题，甘洋科等人5通过收缩推进距离的策略来防止局部自交，并且通过建立空间八叉树来提前检测全局相交问题，再对收缩后的推进方向采用拉普拉斯平滑来改善局部网格单元的质量。刘田田等人6提出一种基于双前沿推进的边界层网格生成算法，并结合多种网格技术处理局部几何特征以保证边界层网格的质量。Ye 等人7提出基于三种策略的边界层网格生成算法，该算法进一步提高了边界层网格的质量，并且能够鲁棒地生成拓扑合法的边界层网格。2.方法方法 2.1.问题与方法概述问题与方法概述 本文以外流场仿真应用为例，目的是生成面向外流场流体计算的高质量三棱柱边界层网格。令,ssSV F=表示一个无自交、可定向的流形三角网格，其中sV与sF分别表示网格的顶点与面片集合。本文在三角网格表面与外流场之间生成指定层数 n、初始高度 h 与增长率 r 的三棱柱边界层网格单元 P，在剩余区域内生成四面体网格单元 T。Open AccessOpen Access陈林祺，周颖 DOI:10.12677/mos.2023.124317 3466 建模与仿真 为鲁棒地生成高质量边界层网格，本文提出基于局部特征尺寸场的边界层网格生成算法，该方法基于前沿推进法生成三棱柱网格单元 P，并基于 Delaunay 方法生成四面体网格单元 T。具体而言，为防止产生自交单元，首先计算表面三角网格顶点的最大推进高度，将最大推进高度作为局部特征尺寸场，接着基于局部特征尺寸场与密度函数得到顶点的密度场，对密度场进行拉普拉斯平滑以使密度值分布更为均匀，基于密度场进行表面三角网格的重网格化，最终分别基于前沿推进法与 Delaunay 方法生成三棱柱网格单元 P 与四面体网格单元 T，得到最终的边界层网格。2.2.局部特征尺寸场局部特征尺寸场 为获得局部特征尺寸场，需要计算表面三角网格顶点的推进方向以及最大推进高度。表面三角网格顶点的推进方向可通过平均顶点周围的面法向得到：11fNpiifnnN=(1)其中pn为推进方向，fN为点 p 相邻的总面数，in表示第 i 个面的法向。由指定层数 n、初始高度 h 与增长率 r 可知顶点的总体推进高度 H 为：()011*nHhrrr=+(2)由于在三角网格及外流场边界之间的一些狭窄区域，过高的推进高度会造成自交单元产生，因此为了得到合法单元，需要对不合理的推进高度采取缩减策略，如图 1 所示。图 1 左侧为自交单元，右侧为经过处理后得到的合法单元。Figure 1.Reduction of excessive advance height 图图 1.缩减过大的推进高度 为检测自交单元，可以对单元构建空间八叉树结构，检测单元之间是否发生碰撞。具体而言，从每个表面三角网格顶点发出一条沿着推进反向的射线，并记录下首次于表面三角网格碰撞时所在的三角面片，取三角面片中距离出发点最近的点并计算两点之间的距离l，并以此调整此顶点的最大推进高度maxH：max,33,3HlHHHlH=(3)为防止三棱柱单元局部的自交问题，需要计算三棱柱单元的雅克比行列式，若为负值，则发生了局部自交，同样需要缩减推进高度以得到合法的三棱柱单元。本文对发生局部自交的三棱柱单元的推进高度以 20%的缩减率进行缩减，每缩减一次则检测是否得到合法单元，若依旧不合法，则继续缩减直到单元合法。由此得到表面三角网格顶点的最大推进高度场，此最大推进高度场实际上就是局部特征尺寸场，接下来需要基于局部特征尺寸场构建顶点密度场。陈林祺，周颖 DOI:10.12677/mos.2023.124317 3467 建模与仿真 2.3.顶点密度场顶点密度场 局部特征尺寸场反映了三棱柱单元的尺寸厚度，若尺寸厚度过低，而底面三角形的尺寸过大，则会导致三棱柱单元纵横比过小。纵横比是一种衡量网格质量的评价指标，对于三棱柱单元，其计算方式为三棱柱平均厚度与底面三角形平均边长之比，其值越接近 1 则三棱柱质量越好。因此为了提升棱柱单元的质量，需要通过构建密度场来增加底面三角形的密度，即对表面三角网格进行重网格化。因此，本文计算每个顶点的密度值的方式为：()()()()()maxmaxmaxmax10,010Density1,101,iiiiiLHvLvHvLHvLLHvL (4)其中，()Densityiv代表顶点iv的密度值，()maxiHv代表顶点iv的最大推进高度，即局部特征尺寸。L 代表输入的表面三角网格的目标平均边长，即重网格化后得到的表面三角网格的平均边长。由平均边长可以推导出重网格化后的目标点数，重网格化后的三角网格应是各向同性的，由各向同性的三角网格可知每个三角形为正三角形，那么每个三角形的面积 a 为：234aL=(5)由模型的表面积 A 与每个三角形的面积 a 可得知总三角形个数trin为：triAna=(6)由欧拉公式可知三角形个数为顶点个数的两倍，则triv为：2tritrinv=(7)若某一点与其邻接点的密度值差距过大，则会造成密度值突变，那么可能导致重网格化后的三角形数量过渡不平滑，会对后续边界层网格生成造成影响，所以需要对密度场进行拉普拉斯平滑。那么对于每个顶点iv的密度值，其新的密度值计算方式为：()()()1DensityDensityDensityviiniijjvvvvn=(8)其中，()Densityiv代表每个顶点iv经过平滑后的密度值，()Densityiv代表初始密度值，n 代表顶点iv的1-邻域所有邻接点的数量。由此密度场构建完毕，接下来需要基于密度场对三角网格进行重网格化。2.4.重网格化重网格化 重网格化的目标是对原始三角网格进行重采样以符合密度场要求。本文基于 Shen 等人8提出的重网格化算法进行重网格化，该方法需要输入重网格化后的目标点数，密度场以及三角网格的特征线，由 3.1节中可知本文已计算得到目标点数和密度场，接下来需要提取三角网格的特征线。本文通过检测相邻三角面片的二面角值来进行特征线的提取，即某条边邻接的两个三角面片的法线夹角大于指定阈值，那么将其视为特征边：1212cosn nnn (9)其中，1n和2n分别为两邻接三角面片的单位法向，如图 2 所示。图 2 左侧为输入三角网格，右侧为提取陈林祺，周颖 DOI:10.12677/mos.2023.124317 3468 建模与仿真 出来的特征线。Figure 2.Extract featureline 图图 2.提取特征线 由此，基于目标点数，密度场以及三角网格的特征线对三角网格进行重网格化后，可得到符合密度场要求的三角网格，接下来需要生成边界层网格。2.5.生成边界层网格生成边界层网格 本文基于前沿推进法生成边界层网格的三棱柱单元，即从表面三角网格顶点出发，沿着顶点的推进方向推进得到新的顶点，并根据表面三角网格的拓扑关系连接新的顶点得到三棱柱单元。由层数，初始高度以及增长率可得到多层三棱柱单元，如图 3 所示。Figure 3.Advancing-layer method 图图 3.前沿推进法 本文为得到合法无自交的三棱柱单元，基于局部特征尺寸场进行三棱柱网格的生成，即在保证单元无自交的情况下进行推进高度的调整，并基于局部特征尺寸场得到密度场，在密度场引导下进行表面三角网格的重网格化，最终基于顶点最大推进高度沿着推进方向得到三棱柱单元。生成得到三棱柱单元后，需要在外流场区域生成四面体网格单元。生成四面体网格的方法有很多9 10，本文借助 TetGen 10开源库生成高效鲁棒的四面体网格单元。其生成算法基于 Delaunay 算法，从理论上保证了四面体网格的质量。3.实验结果与分析实验结果与分析 为测试以及评估本文方法，本文通过实验分析得到的结果，实验所用操作系统为 windows11，处理器为 2.9 GHz IntelCorei5-10400 CPU，内存为 16 GB。本文实现该算法为一个单线程的 C+程序，其涉及到矩阵运算库，由 Eigen 11库实现。本文在若干个网格上实现了该算法，并通过一系列质量指标评估结果。陈林祺，周颖 DOI:10.12677/mos.2023.124317 3469 建模与仿真 3.1.边界层网格生成结果边界层网格生成结果 图 4(a)，图 4(b)分别为 airplane 模型，scaffold 模型与 joint 模型生成得到的边界层网格，从图中可以看出，各模型生成的边界层网格单元的过渡性较好，具有良好的正交性，三棱柱单元较为规整，四面体单元大部分接近正四面体。表 1 列出了各个模型的边界层网格的质量评价指标，其中，纵横比指标越接近 1 越好，扭曲度越小越好，雅克比值越接近 1 越好。由表 1 可见，本文生成得到的边界层网格的质量较好。(a)airplane (b)scaffold (c)joint Figure 4.Boundary layer mesh 图图 4.边界层网格 Table 1.Boudary layermesh quality evaluation 表表 1.边界层网格质量评价指标 模型 纵横比 扭曲度 雅克比值 时间(s)airplane 0.374 0.176 0.223 10 scaffold 0.332 0.226 0.281 15 joint 0.328 0.176 0.294 9 3.2.重网格化结果分析重网格化结果分析 本文的重网格化旨在调整表面三角网格的密度以符合密度场要求，目的是为了改善复杂模型狭窄区域的三棱柱网格单元的质量。为了展示重网格化带来的提升，本文对测试模型构造狭窄区域，位于模型与外流场之间。图 5图 7 分别展示了 assy 模型、aries 模型与 bendix 模型重网格化前后得到的边界层网格的切割视 (a)assy 模型与外流场 (b)重网格化前的边界层网格 (c)重网格化后的边界层网格 Figure 5.Boundary layer mesh cutting view before and after remesh of assy model 图图 5.assy 模型重网格化前后生成的边界层网格切割视图 陈林祺，周颖 DOI:10.12677/mos.2023.124317 3470 建模与仿真 (a)aries 模型与外流场 (b)重网格化前的边界层网格 (c)重网格化后的边界层网格 Figure 6.Boundary layer mesh cutting view before and after remesh of aries model 图图 6.aries 模型重网格化前后生成的边界层网格切割视图 (a)bendix 模型与外流场 (b)重网格化前的边界层网格 (c)重网格化后的边界层网格 Figure 7.Boundary layer mesh cutting view before and after remesh of aries model 图图 7.bendix 模型重网格化前后生成的边界层网格切割视图 Table 2.Boudary layermesh quality evaluation 表表 2.边界层网格质量评价指标 模型 三角网格 纵横比 扭曲度 雅克比值 assy 重网格化前 0.273 0.224 0.133 重网格化后 0.306 0.196 0.201 aries 重网格化前 0.266 0.237 0.145 重网格化后 0.315 0.183 0.178 bendix 重网格化前 0.247 0.253 0.130 重网格化后 0.317 0.209 0.219 图。从图中可以看出，重网格化前模型的狭窄区域的三棱柱单元非常扁平，其纵横比过小，经过重网格化后，狭窄区域的表面三角形的密度增加，因而提升了三棱柱单元的纵横比，从而提高了边界层网格的质量。表 2 列出了各个模型重网格化前后的质量指标比较。由表 2 可见，重网格化后能有效提升边界层网格的质量。4.总结总结 本文主要针对计算流体力学中的边界层网格生成问题提出了一种基于局部特征尺寸场的边界层网格生成算法。该算法首先计算三角网格顶点的最大推进高度作为局部特征尺寸场，并基于局部特征尺寸场构造表面三角网格顶点的密度场，为使得密度值过渡平滑，采用了拉普拉斯平滑方法平滑密度场，接着陈林祺，周颖 DOI:10.12677/mos.2023.124317 3471 建模与仿真 基于密度场对表面三角网格进行重网格化，由于重网格化算法需要提取三角网格的特征线，因此基于二面角的方法对三角网格进行了特征线的提取，最后基于前沿推进法得到三棱柱网格单元，并在外流场区域基于 Delaunay 法生成四面体网格单元，得到最终的边界层网格。实验结果表明，本文提出的方法能鲁棒地生成边界层网格，并且边界层网格单元拓扑合法，不会有自交单元的产生，此外，结合重网格化方法对三角网格进行密度调整，在新的三角网格上生成边界层网格，能有效提高复杂模型的狭窄区域处的网格单元质量。参考文献参考文献 1 Park,S.,Jeong,B.,Jin,G.L.,et al.(2013)Hybrid Grid Generation for Viscous Flow Analysis.International Journal for Numerical Methods in Fluids,71,891-909.https:/doi.org/10.1002/fld.3691 2 Lu,F.,Pang,Y.,Jiang,X.,et al.(2017)Automatic Generation of Structured Multiblock Boundary Layer Mesh for Aircrafts.Advances in Engineering Software,115,297-313.https:/doi.org/10.1016/j.advengsoft.2017.10.003 3 曹杰,赵铭伟,余飞,昌继海,关振群.高效可靠的边界层网格分块层进生成算法J.计算机辅助设计与图形学学报,2020,32(8):1199-1207.4 王硕.基于飞行器CFD计算的非结构附面层网格生成技术研究D:硕士学位论文.绵阳:西南科技大学,2021.5 甘洋科,刘剑飞.黏性边界层网格自动生成J.力学学报,2017,49(5):1029-1041.6 刘田田,董庆利,杨洋,冷珏琳,郑澎.一种基于双前沿推进的边界层网格生成算法J.计算力学学报,2022,39(1):42-48.7 沈鑫鑫.流形三角形网格重网格化方法D:硕士学位论文.杭州:浙江大学,2016.8 Ye,H.F,Liu,Y.,Chen,B.,et al.(2020)Hybrid Grid Generation for Viscous Flow Simulations in Complex Geometries.Advances in Aerodynamics,2,Article No.17.https:/doi.org/10.1186/s42774-020-00042-x 9 刘岩.高效可靠的三维约束 Delaunay 四面体有限元网格生成算法D:硕士学位论文.大连:大连理工大学,2010.10 Si,H.(2015)TetGen:A Quality Tetrahedral Mesh Generator.ACM Transactions on Mathematical Software,41,Ar-ticle No.11.https:/doi.org/10.1145/2629697 11 Guennebaud,G.,Jacob,B.,et al.(2010)Eigen.http:/eigen.tuxfamily.org.