CN112613206B - 一种基于各向异性体调和场的附面层网格生成方法 - Google Patents
一种基于各向异性体调和场的附面层网格生成方法 Download PDFInfo
- Publication number
- CN112613206B CN112613206B CN202011470101.7A CN202011470101A CN112613206B CN 112613206 B CN112613206 B CN 112613206B CN 202011470101 A CN202011470101 A CN 202011470101A CN 112613206 B CN112613206 B CN 112613206B
- Authority
- CN
- China
- Prior art keywords
- boundary layer
- mesh
- field
- tetrahedral
- tensor
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000009499 grossing Methods 0.000 claims abstract description 11
- 238000005070 sampling Methods 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 230000006698 induction Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000000007 visual effect Effects 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 6
- 238000011960 computer-aided design Methods 0.000 abstract description 3
- 239000012530 fluid Substances 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 239000002699 waste material Substances 0.000 abstract 1
- 239000010410 layer Substances 0.000 description 68
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000004308 accommodation Effects 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000003102 growth factor Substances 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Image Generation (AREA)
Abstract
本发明公开一种基于各向异性体调和场的附面层网格生成方法，属于计算流体力学、数值模拟仿真、计算机辅助设计、与制造技术领域。首先采用闵科夫斯基和的边界曲面网格构造求解体调和场时所需的四面体背景网格，然后根据需求自动化地添加各向异性张量，并在张量的控制下计算各向异性体调和场，最后结合特殊的加权拉普拉斯平滑计算生成附面层网格所需的前进方向。本发明基于闵科夫斯基和边界曲面网格构造四面体背景网格的策略降低了计算耗时和内存浪费，通过自动化添加各向异性张量能够可控地局部调节附面层网格厚度，结合特殊的加权拉普拉斯光滑对前进方向进行优化显著地提高了附面层网格的生成质量。
Description
技术领域
本发明属于计算流体力学、数值模拟仿真、计算机辅助设计，以及制造技术领域，涉及一种基于各向异性体调和场生成附面层网格的方法，适用于复杂曲面的附面层网格生成。通过张量对体调和场的控制，使其诱导生成的附面层网格更具有灵活性和可控性。
背景技术
附面层是高雷诺数绕流中紧贴物面的粘性力不可忽略的流动薄层，而附面层网格的质量直接决定着数值模拟效果的好坏。在高雷诺数流动的气动模拟中，必须采用垂直于物体的层状各向异性棱柱网格来捕捉粘性壁面附近的附面层。粘性壁面周围的棱柱网格生成一直是计算流体力学、数值模拟仿真、计算机辅助设计与制造技术领域的研究重点。棱柱网格生成的主流方法有两种：前沿结点推进方法，基于求解偏微分方程方法。但是已有的方法一般不考虑各向异性的问题，即每一层增长的厚度基本都一样，面对一些特殊的数值模拟需求时，各向同性很难捕捉一些微小的物理特征。如何使得生成棱柱网格更具有灵活性和可控性，是附面层网格研究的重心。
发明内容
基于上述的问题，本发明提出一种基于各向异性体调和场生成附面层网格的方法。本方法属于基于求解偏微分方程生成附面层网格，包含3个发明内容：
1.基于闵可夫斯基和的边界曲面网格构造，以及附面层空间的四面体背景网格(离散计算域)生成。
2.基于局部张量控制的各向异性体调和场的计算。
3.基于各向异性体调和场计算的前进距离和前进方向的附面层网格(棱柱网格)生成策略。
本发明的技术方案：
一种基于各向异性体调和场的附面层网格生成方法，步骤如下：
(1)基于闵可夫斯基和的边界曲面网格构造，以及附面层空间的四面体背景网格(离散计算域)生成，具体步骤如下：
a)输入原始曲面网格(一般为三角形网格或者四边形网格)和半径为r(一般取值为能够包裹原始曲面网格的最小长方体的对角线长度乘以一个系数c，c一般取值在[0.05,0.3])的球网格，计算闵可夫斯基和的边界曲面网格。
b)对初步得到的闵可夫斯基和的边界曲面网格执行网格优化处理，包含非流形消除和自相交消除，最终得到二维流形的边界曲面网格。
c)定义附面层空间为：计算得到的闵可夫斯基和的边界曲面网格与原始曲面网格(物面)之间的空间。
d)对附面层空间使用软件TetGen进行四面体网格剖分，检测四面体背景网格中是否存在四面体单元的四个点同时位于附面层空间的边界曲面网格上。如果存在则进行细分，直到无上述情况出现。
e)对四面体背景网格的狭缝处进行局部细分。首先，固定原始曲面网格(物面)的体调和能量为常值a(一般取值为1)，固定外包围曲面网格的体调和能量为常值b(一般取值为0)，边权重设置为经典的余切权重(Cotangent Weight)，计算四面体背景网格上的体调和场；然后，从靠近原始曲面网格(物面)的四面体单元开始做广度优先搜索(BFS)，寻找满足每条边上的能量值之差同时小于阈值Tslit(一般取值在[0.01,0.1])的四面体单元，记作集合Rslit；最后，对集合Rslit内的四面体单元进行细分(细分次数一般取值在[2,10])。对局部细分后的四面体背景网格使用delaunay处理进行优化，再使用拉普拉斯平滑辅助优化，最终得到高质量的四面体背景网格。
(2)基于局部张量控制的各向异性体调和场的计算，具体内容如下：
2.1)四面体背景网格顶点上张量的定义
T(υi)＝γ1x1x1 T+γ2x2x2 T+γ3x3x3 T； (1)
其中，υi是四面体背景网格上的顶点，[x1,x2,x3]是标准的三维正交标架，γ1,γ2,γ3分别作为标准正交标架三个方向上的缩放因子。直观上，张量对体调和场的控制可以看作在每一个顶点上放置一个椭球，如图2所示，椭圆的长短轴边长表示着张量沿着长短轴方向对体调和场的控制程度。
2.2)各向异性体调和场在四面体背景网格顶点上的定义
LH＝0； (2)
其中，H是体调和场作用在每一个顶点上的值组合而成的向量；L是权重矩阵，表达式为：
其中，υi是四面体背景网格上的顶点，eij是四面体背景网格上连接υi和υj的边，eij是四面体背景网格上连接υi和υj的边，N(υi)是与υi相邻顶点的集合，W(eij)是求解拉普拉斯方程的边权重：
其中，δ>0，是一个控制因子(一般取值在[0.001，100])，δ越小则边权重W(eij)受到张量T(vi)，T(vj)的影响越大，反之，则边权重W(eij)受到张量T(vi)，T(vj)的影响越小。
2.3)四面体背景网格中的各向异性体调和场的计算
定义各向异性体调和能量K(H)：
其中，E是四面体背景网格上的边集合，H(vi)是在vi处的各向异性体调和能量。
在基于迭代法计算的各向异性体调和场的框架下，边权重可以简化为：
W(eij)＝exp(T(eij)/δ)； (6)
在基于迭代法计算的各向异性体调和场的框架下，对于顶点上体调和能量的优化，表示为：
基于迭代法的各向异性体调和场的计算的算法流程如下：
固定原始曲面网格(物面)的体调和能量为常值a(一般取值为1)，固定外包围曲面网格的体调和能量为常值b(一般取值为0)；根据实际需求，自动化地添加局部各向异性张量控制；设置最大的迭代次数Titer(一般设置为2000)；设置优化体调和能量的截断阈值Tenergy(一般设置为1.0*10-8)；通过公式(7)，迭代地更新H(vi)；每次更新完所有顶点上的体调和能量，计算一次体调和能量K(H)；迭代过程直到满足最大的迭代次数Titer，或者K(H)达到截断阈值Tenergy。
2.4)自动化的各向异性张量控制的构造
2.4.1)一般情况下，限制体调和场沿着特定方向d的梯度变化速率，表示为：
将公式(8)构造的张量作用于体调和场的计算中，则沿着方向d上的梯度变化速率将受到限制；基于公式(8)张量控制的体调和场生成的附面层网格(棱柱网格)，沿着方向d，直观上，附面层网格的总体厚度明显减小。
2.4.2)限制凹边和凹槽处的体调和场梯度变化速率，为了使得计算得到的体调和场等值面更贴合于物面，避免诱导生成的附面层网格在凹边和凹槽处产生较大畸变，自动构造局部张量方式如下：
首先，设置原始曲面网格(物面)的能量为常量a(一般取值为1)，外包围曲面网格的能量常量b(一般取值为0)，要求a>b，边权重设置为经典的余切权重(CotangentWeight)，计算四面体背景网格上的体调和场；然后，从靠近原始曲面网格(物面)的四面体单元开始做广度优先搜索(BFS)，寻找满足每条边上的能量值之差同时小于阈值Tslit(一般取值在[0.01,0.1])的四面体单元，记作集合Rslit；最后，计算张量，表示为：
将公式(9)构造的张量作用于体调和场的计算中，则在凹边和凹槽处的体调和场梯度变化速率将受到限制；基于公式(9)张量控制的体调和场生成的附面层网格(棱柱网格)，在凹边和凹槽处，直观上，附面层网格的畸变明显减小。
2.4.3)限制多连通分支之间狭缝处的体调和场梯度变化速率，延缓体调和场在狭缝处鞍点的产生，提高诱导生成附面层网格在多连通分支之间狭缝处的质量。自动构造局部张量方式如下，设两个相互靠近的体模型P，Q：首先，分别设H1(P)＝a，H1(Q)＝b和H2(P)＝b，H2(Q)＝a作为狄利克雷边界条件(a一般取值为1；b一般取值为0)，边权重设置为经典的余切权重(Cotangent Weight)，计算两个标准体调和场H1，H2；然后，计算张量，表示为：
将公式(10)构造的张量作用于体调和场的计算中，则在多连通分支之间狭缝处的体调和场梯度变化速率将受到限制；基于公式(10)张量控制的体调和场生成的附面层网格(棱柱网格)，直观上，在狭缝处附面层网格的畸变明显减小。
(3)基于各向异性体调和场计算的前进距离和前进方向的附面层网格(棱柱网格)生成策略，具体内容如下：
3.1)前沿结点前进距离的计算
前沿结点的前进距离，通过各向异性体调和场等值面之间的间隙来控制。本发明将用户输入的期望网格厚度转化为采样能量，通过采样能量计算每一层结点的位置。具体实现方式：
首先，根据用户输入的第一层附面层厚度L1，附面层厚度增长速度因子α，附面层层数n，计算出每一层附面层网格的厚度。然后，将物面上的顶点设置为前沿结点，选择一个曲率接近0的前沿结点沿着体调和场的梯度线追溯到外包围曲面网格，并根据计算出的每一层附面层网格厚度，在体调和场中提取n个采样能量；最后，对于离散化在四面体背景网格中的体调和场，每一个四面体单元内都是一个线性空间，在前沿结点前进方向的引导下，前沿结点前进后的位置能够很容易通过采样能量确定。
3.2)前沿结点前进方向的计算
前沿结点的前进方向，通过对体调和场的梯度方向进行加权拉普拉斯光滑化后得到。具体实现方式：
首先，计算前沿结点当前位置的梯度方向，一般为所在四面体单元内等值面的法向量方向；然后，记当前位置为pi，在梯度方向和下一个采样能量的引导下，计算出下一个位置
如图3所示；最后，将拉普拉斯光滑的权重设置为
其中，p(一般取值为4)和q(一般取值为2)是两个控制参数，另外
在这种权重表示下，能够直接检测出在前进方向和前进距离引导下生成的附面层网格是否存在负体积单元，有效保障了生成附面层网格的质量；一般拉普拉斯光滑次数设置为100次。
3.3)附面层网格(棱柱网格)的生成
前沿结点在前进距离和前进方向的引导下，计算得到一族新的前进位置；附面层网格由所有前沿结点的前进位置根据原始曲面网格(物面)拓扑的有向连接得到。
本发明的有益效果：
基于上述的发明内容，本发明提出的基于各向异性体调和场生成附面层网格的方法包含3个有益效果：
(1)传统构造四面体背景网格(离散计算域)的方法采用长方体或者球体作为外包围曲面网格，这样的策略引入大量冗余的四面体单元，带来额外的计算量消耗。基于闵可夫斯基和边界曲面网格的四面体背景网格(离散计算域)生成能够有效地消除冗余四面体单元，从而提高内存的利用率和算法的执行效率。
(2)传统基于求解偏微分方程的附面层网格生成方法只考虑全局信息，欠缺局部上的控制力和灵活性。本发明通过自动化地构造局部各向异性张量来实现对体调和场的控制，一方面，使其能够感知局部的几何信息，加强对局部生成网格的控制力和灵活性(主要是针对附面层网格厚度的控制力)；另一方面，可以根据实际需求，使得生成的附面层网格沿着某个或者某几个方向变得稠密，以至于能够捕捉细微的物理特征。
(3)本发明中对于前沿结点的前进距离和前进方向的计算策略具有以下好处：
a)为了让体调和场的构造与附面层网格的生成更好地结合，本发明将体调和场等值面之间的距离作为引导，控制每一层前沿结点的前进距离，即附面层网格中的每一层结点，在四面体背景网格(离散计算域)中的体调和能量值是相等的，它们都处于同一层体调和场等值面上。基于这种策略，建立体调和场与附面层网格厚度之间的耦合关系，使附面层网格厚度更灵活、更可控，以致于能够通过局部控制体调和场达到局部控制附面层网格厚度的目的，以满足复杂的实际需求。
传统基于偏微分方程生成附面层网格(棱柱网格)的方法，一般直接使用梯度方向作为前进方向，然而容易在凹边和凹槽处引入大量负体积或者零体积的棱柱体单元。除此之外，存在许多工作将加权拉普拉斯平滑应用于梯度方向上，得到更光滑的前进方向。对于加权拉普拉斯平滑策略，重点是权重的设计和选择，这将直接影响附面层网格生成的结果。本发明所采取的策略，利用梯度方向计算初始前进位置，并结合当前位置信息直接计算所生成棱柱体单元的质量，将当前棱柱体单元的质量作为权重设置的根据，一定程度上，有效地避免了负体积棱柱体单元的生成。
附图说明
图1为本发明的算法流程图；
图2为张量在顶点上作用于体调和场的示意图；
图3为优化前进方向的拉普拉斯光滑的权重设计的示意图；
图4为飞机模型基于各向异性体调和场生成附面层网格的示意图，(a)飞机模型原始曲面(物面)网格；(b)飞机模型外包围曲面网格；(c)飞机模型原始(物面)网格；(d)飞机模型标准体调和场等值面的横截面；(e)飞机模型各向异性体调和场等值面的横截面；(f)飞机模型基于各向异性体调和场的附面层网络。
具体实施方式
以下结合附图和技术方案，进一步详述本发明的具体实施方式。
本发明的算法流程如图1所示，总体包含5个步骤：构造闵科夫斯基和的边界曲面网格；对附面层空间生成四面体背景网格(离散计算域)；构造各向异性张量；计算各向异性体调和场；基于各向异性体调和场生成附面层网格(棱柱网格)。本发明的算法的输入包含1份原始曲面网格(物面)和3个参数。其中，原始曲面网格(物面)可以为三角形网格或者四边形网格；3个参数分别为第一层附面层网格厚度L1，附面层网格厚度增长因子α以及附面层网格的层数n。
本实施例以飞机模型基于各向异性体调和场生成附面层网格的具体实施为本发明的示例，如图4所示，具体步骤如下：
1.输入飞机模型(三角形网格)；输入参数L1＝1.0*10-2，α＝1.15，n＝60；根据输入参数L1，α，n，计算每一层附面层网格的期望厚度{L1,L2,...,Ln}，其中，Li+1＝Li*α；计算附面层网格的期望总体厚度
2.设置小球网格半径为Ltotal*1.5，计算原始曲面网格(物面)的闵科夫斯基和的边界曲面网格作为外包围曲面网格；消除外包围曲面网格的非流型区域和自相交区域；
3.为原始曲面网格(物面)和外包围曲面网格之间的附面层空间使用商业化软件TetGen进行四面体网格剖分，得到四面体背景网格(离散计算域)；检测四面体背景网格中是否存在四面体单元的四个点同时位于附面层空间的边界曲面网格上，如果存在则进行细分，直到无上述情况出现；
4.固定原始曲面网格(物面)上的体调和能量为1，固定外包围曲面网格上的体调和能量为0，边权重设置为经典的余切权重(Cotangent Weight)，计算四面体背景网格上的体调和场；以依附着原始曲面网格(物面)的四面体单元为起始位置执行广度优先搜索(BFS)，寻找满足每条边上的能量值之差同时小于阈值0.05的四面体单元，记作集合Rslit；根据公式(9)计算局部张量；
5.固定原始曲面网格(物面)上的体调和能量为1，固定外包围曲面网格上的体调和能量为0；基于各向异性张量，通过公式(6)计算出每一条边上的权重，其中，控制因子δ设置为0.05；设置最大迭代次数为2000，设置能量截断阈值为1.0*10-8；通过公式(7)迭代更新顶点上的体调和能量值，每迭代50次，根据公式根据公式(5)计算新的体调和能量K(H)；若当前体调和能量值与上一次体调和能量值之差小于截断阈值则迭代结束，否则继续迭代，直到达到最大迭代次数；
6.将原始曲面网格(物面)上的顶点设置为前沿结点，选择一个曲率最接近于0的前沿结点沿着体调和场的梯度线追溯到外包围曲面网格，根据计算出的每一层附面层网格厚度{L1,L2,...,Ln}，在追溯的轨迹上提取n个采样能量{K1,K2,...,Kn}；
7.将前沿结点当前所在位置的梯度方向作为初始前进方向，以公式(11)作为加权拉普拉斯光滑的权重对前进方向进行优化，平滑次数设置为100次；
8.基于优化后的前进方向和采样能量Ki，计算出前沿结点推进后的位置；每推进前沿结点一次，重新计算前进方向；根据原始曲面网格(物面)的拓扑对所有前沿结点所分别对应的n个推进位置进行简单的有向连接，得到附面层网格(棱柱网格)。
Claims (1)
1.一种基于各向异性体调和场的附面层网格生成方法，其特征在于，步骤如下：
(1)基于闵可夫斯基和的边界曲面网格构造，以及附面层空间的四面体背景网格生成
a)输入原始曲面网格和半径为r的球网格，计算闵可夫斯基和的边界曲面网格；
b)对初步得到的闵可夫斯基和的边界曲面网格执行网格优化处理，包含非流形消除和自相交消除，最终得到二维流形的边界曲面网格；
c)定义附面层空间为：计算得到的闵可夫斯基和的边界曲面网格与原始曲面网格之间的空间；
d)对附面层空间进行四面体网格剖分，检测四面体背景网格中是否存在四面体单元的四个点同时位于附面层空间的边界曲面网格上；如果存在则进行细分，直到无上述情况出现；
e)对四面体背景网格的狭缝处进行局部细分，并对局部加密后的四面体背景网格进行优化，再使用拉普拉斯平滑辅助优化，最终得到高质量的四面体背景网格；
(2)基于局部张量控制的各向异性体调和场的计算
1)张量在四面体背景网格顶点上的定义：
T(υi)＝γ1x1x1 T+γ2x2x2 T+γ3x3x3 T； (1)
其中，υi是四面体背景网格上的顶点，[x1，x2，x3]是标准的三维正交标架，γ1，γ2，γ3分别作为标准正交标架三个方向上的缩放因子；
2)各向异性体调和场在四面体背景网格顶点上的定义：
LH＝0； (2)
其中，H是各向异性体调和场作用在每一个顶点上的值组合而成的向量；L是权重矩阵，表达式为：
其中，υi是四面体背景网格上的顶点，eij是四面体背景网格上连接υi和υj的边，N(υi)是与υi相邻顶点的集合，W(eij)是求解拉普拉斯方程的边权重：
其中，δ＞0，是一个控制因子，δ越小则边权重W(eij)受到张量T(vi)，T(vj)的影响越大，反之，则边权重W(eij)受到张量T(vi)，T(vj)的影响越小；
3)四面体背景网格中的各向异性体调和场的计算：
定义各向异性体调和能量K(H)：
其中，E是四面体背景网格上的边集合，H(vi)是在vi处的各向异性体调和能量；
在基于迭代法计算的各向异性体调和场的框架下，边权重简化为：
W(eij)＝exp(T(eij)/δ)； (6)
在基于迭代法计算的各向异性体调和场的框架下，对于顶点上体调和能量的优化，表示为：
基于迭代法的各向异性体调和场的计算的算法流程如下：
固定原始曲面网格的体调和能量为常值a，固定外包围曲面网格的体调和能量为常值b；根据实际需求，自动化地添加局部各向异性张量控制；设置最大的迭代次数Titer；设置优化体调和能量的截断阈值Teneergy；通过公式(7)，迭代地更新H(vi)；每次更新完所有顶点上的体调和能量，计算一次体调和能量K(H)；迭代过程直到满足最大的迭代次数Titer，或者K(H)达到截断阈值Tenergy；
4)自动化的各向异性张量控制的构造：
4.1)限制体调和场沿着特定方向d的梯度变化速率，表示为：
将公式(8)构造的张量作用于体调和场的计算中，则沿着方向d上的梯度变化速率将受到限制；基于公式(8)张量控制的体调和场生成的附面层网格，沿着方向d，直观上，附面层网格的总体厚度明显减小；
4.2)限制凹边和凹槽处的体调和场梯度变化速率，为了使得计算得到的体调和场等值面更贴合于物面，避免诱导生成的附面层网格在凹边和凹槽处产生较大畸变，自动构造局部张量方式如下：
首先，设置原始曲面网格的体调和能量为常量a，外包围曲面网格的能量常量b，要求a＞b，计算四面体背景网格上的体调和场；然后，从靠近原始曲面网格的四面体单元开始做广度优先搜索，寻找满足每条边上的能量值之差同时小于阈值Tslit的四面体单元，记作集合Rslit；最后，计算张量，表示为：
将公式(9)构造的张量作用于体调和场的计算中，则在凹边和凹槽处的体调和场梯度变化速率将受到限制；基于公式(9)张量控制的体调和场生成的附面层网格，在凹边和凹槽处，直观上，附面层网格的畸变明显减小；
4.3)限制多连通分支之间狭缝处的体调和场梯度变化速率，延缓体调和场在狭缝处鞍点的产生，提高诱导生成附面层网格在多连通分支之间狭缝处的质量；自动构造局部张量方式如下，设两个相互靠近的体模型P，Q：首先，分别设H1(P)＝a，H1(Q)＝b和H2(P)＝b，H2(Q)＝a作为狄利克雷边界条件，计算两个标准体调和场H1，H2；然后，计算张量，表示为：
将公式(10)构造的张量作用于体调和场的计算中，则在多连通分支之间狭缝处的体调和场梯度变化速率将受到限制；基于公式(10)张量控制的体调和场生成的附面层网格，直观上，在狭缝处附面层网格的畸变明显减小；
(3)基于各向异性体调和场计算的前进距离和前进方向的附面层网格生成策略
前沿结点前进距离的计算：
前沿结点的前进距离，通过各向异性体调和场等值面之间的间隙来控制；将用户输入的期望网格厚度转化为采样能量，通过采样能量计算每一层结点的位置；具体实现方式：
首先，根据用户输入的第一层附面层厚度L1，附面层厚度增长速度因子α，附面层层数n，计算出每一层附面层网格的厚度；然后，将物面上的顶点设置为前沿结点，选择一个曲率接近0的前沿结点沿着体调和场的梯度线追溯到外包围曲面网格，并根据计算出的每一层附面层网格厚度，在体调和场中提取n个采样能量；最后，对于离散化在四面体背景网格中的体调和场，每一个四面体单元内都是一个线性空间，在前沿结点前进方向的引导下，前沿结点前进后的位置能够很容易通过采样能量确定；
a)前沿结点前进方向的计算：
前沿结点的前进方向，通过对体调和场的梯度方向进行加权拉普拉斯光滑化后得到；具体实现方式：
首先，计算前沿结点当前位置的梯度方向，为所在四面体单元内等值面的法向量方向；然后，记当前位置为pi，在梯度方向和下一个采样能量的引导下，计算出下一个位置
最后，将拉普拉斯光滑的权重设置为
其中，p和q是两个控制参数，另外
在这种权重表示下，能直接检测出在前进方向和前进距离引导下生成的附面层网格是否存在负体积单元，有效保障了生成附面层网格的质量；
b)附面层网格的生成：
前沿结点在前进距离和前进方向的引导下，计算得到一族新的前进位置；附面层网格由所有前沿结点的前进位置根据原始曲面网格拓扑的有向连接得到。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011470101.7A CN112613206B (zh) | 2020-12-15 | 2020-12-15 | 一种基于各向异性体调和场的附面层网格生成方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011470101.7A CN112613206B (zh) | 2020-12-15 | 2020-12-15 | 一种基于各向异性体调和场的附面层网格生成方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112613206A CN112613206A (zh) | 2021-04-06 |
| CN112613206B true CN112613206B (zh) | 2022-09-20 |
Family
ID=75233899
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011470101.7A Active CN112613206B (zh) | 2020-12-15 | 2020-12-15 | 一种基于各向异性体调和场的附面层网格生成方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112613206B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022236449A1 (zh) * | 2021-05-08 | 2022-11-17 | 大连理工大学 | 一种基于多重体调和场的三维模型保拓扑变形方法 |
| CN114820991B (zh) * | 2022-06-30 | 2022-09-16 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种非结构附面层网格交叉处理方法及装置 |
| CN115758938B (zh) * | 2022-11-25 | 2023-07-25 | 浙江大学 | 面向粘性边界流场数值模拟的附面层网格生成方法 |
| CN117077578B (zh) * | 2023-10-17 | 2024-01-12 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于桩基点生成附面层网格的方法、装置、设备和介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007264952A (ja) * | 2006-03-28 | 2007-10-11 | Fujitsu Fip Corp | 地盤解析用メッシュ生成方法及び地盤解析用メッシュ生成プログラム |
| CN109410335A (zh) * | 2018-09-21 | 2019-03-01 | 浙江理工大学 | 一种基于可复用拉普拉斯矩阵的高效网格融合方法 |
| KR102097707B1 (ko) * | 2018-10-02 | 2020-04-08 | 주식회사 투썬디지털아이디어 | 특수효과 제작 소프트웨어의 유체 시뮬레이션 방법 |
| CN111079326A (zh) * | 2019-11-20 | 2020-04-28 | 杭州电子科技大学 | 二维各向异性网格单元度量张量场光滑化方法 |
-
2020
- 2020-12-15 CN CN202011470101.7A patent/CN112613206B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007264952A (ja) * | 2006-03-28 | 2007-10-11 | Fujitsu Fip Corp | 地盤解析用メッシュ生成方法及び地盤解析用メッシュ生成プログラム |
| CN109410335A (zh) * | 2018-09-21 | 2019-03-01 | 浙江理工大学 | 一种基于可复用拉普拉斯矩阵的高效网格融合方法 |
| KR102097707B1 (ko) * | 2018-10-02 | 2020-04-08 | 주식회사 투썬디지털아이디어 | 특수효과 제작 소프트웨어의 유체 시뮬레이션 방법 |
| CN111079326A (zh) * | 2019-11-20 | 2020-04-28 | 杭州电子科技大学 | 二维各向异性网格单元度量张量场光滑化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112613206A (zh) | 2021-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112613206B (zh) | 一种基于各向异性体调和场的附面层网格生成方法 | |
| CN110110413B (zh) | 一种基于材料场缩减级数展开的结构拓扑优化方法 | |
| Zingg et al. | A comparative evaluation of genetic and gradient-based algorithms applied to aerodynamic optimization | |
| Toal et al. | Geometric filtration using proper orthogonal decomposition for aerodynamic design optimization | |
| Wang et al. | A review of flux reconstruction or correction procedure via reconstruction method for the Navier-Stokes equations | |
| Yamazaki et al. | Geometry parameterization and computational mesh deformation by physics-based direct manipulation approaches | |
| CN108763658A (zh) | 基于等几何方法的组合薄壁结构固有频率设计方法 | |
| CN109726433B (zh) | 基于曲面边界条件的三维无粘低速绕流的数值模拟方法 | |
| WO2022126348A1 (zh) | 一种基于各向异性体调和场的附面层网格生成方法 | |
| Loppi et al. | Locally adaptive pseudo-time stepping for high-order flux reconstruction | |
| Xu et al. | Adjoint-based optimization for thrust performance of three-dimensional pitching–rolling plate | |
| Arisoy et al. | Design and topology optimization of lattice structures using deformable implicit surfaces for additive manufacturing | |
| Ebrahimi et al. | Shape modeling based on specifying the initial B-spline curve and scaled BFGS optimization method | |
| Balan et al. | Adjoint-based anisotropic mesh adaptation for a stabilized finite-element flow solver | |
| Lavoie et al. | An improved characteristic based volume penalization method for the Euler equations towards icing applications | |
| Toal et al. | Geometric filtration using POD for aerodynamic design optimization | |
| Du et al. | Super Resolution Generative Adversarial Networks for Multi-Fidelity Pressure Distribution Prediction | |
| Fidkowski et al. | Output error estimation and adaptation in computational fluid dynamics: Overview and recent results | |
| Yang et al. | Isogeometric double-objective shape optimization of free-form surface structures with Kirchhoff–Love shell theory | |
| Ito et al. | Efficient computational fluid dynamics evaluation of small-device locations with automatic local remeshing | |
| Shen et al. | Hexahedral mesh adaptation based on posterior-error estimation | |
| Hasanzadeh Lashkajani et al. | Adaptive curvature control grid generation algorithms for complex glaze ice shapes RANS simulations | |
| Chen et al. | Multi-objective topology optimization of heat transfer surface using level-set method and adaptive mesh refinement in OpenFOAM | |
| Smith et al. | A geometry projection method for the design exploration of wing-box structures | |
| Wu et al. | Adjoint-based optimization for blended-wing-body underwater gliders' shape design |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20210406 Assignee: Shanghai Geyu Software Co.,Ltd. Assignor: DALIAN University OF TECHNOLOGY Contract record no.: X2023310000170 Denomination of invention: A boundary layer mesh generation method based on anisotropic volume harmonic field Granted publication date: 20220920 License type: Common License Record date: 20231208 |