Dynamic meshing 动网格设置详细笔记
introduction
- Smoothing Methods:用于三角形或四面体单元。节点移动,但数量和连接性不变。更使用于平移运动。
- dynamic layering:用于棱柱网格(六面体或者棱柱层)。
- remeshing:用于三角形或四面体单元,仅在三维网格可以在边界层网格中重建四面体单元和楔形网格。大变形建议开启。
1. Dynamic Smoothing
1. Diffusion: 比 spring 方法成本更高,但质量通常更好,允许更大边界变形
变形前 | Diffusion 方法变形后 | Spring 方法变形后 |
---|---|---|
- Diffusion function: 可选距离相关/体积相关。设置扩散方程与距离还是体积相关
- Diffusion parameter: 扩散系数。0 为均匀扩散,该系数越高则远处(越大体积)网格吸收运动越多。(通常选用 0-2)
Diffusion parameter = 1 | Diffusion parameter=0 |
---|---|
- AMG stabilization:平滑方法选择
- CG(共轭梯度)法:最快。可能会产生负体积,可通过增加 Maxium number of iterations 到200-500改善;如果用 CG 检测到 divergence (散度),GMRES 方法会自动调用并通知。
- GERMS(广义最小残差)法:速度慢,占用内存高。若用 CG 法时多次调用 GERMS 可选择使用。
- Maximum Number of Iterations : 一般保持默认,出现负体积可以增加。
- Relative Convergence Tolerance: 控制扩散的精度,大部分时间不必要
- Verbosity:输入 1 可以在控制台观察平滑残差
- Generalized Boundary Distance Method:保持使用变形前的边界距离计算扩散函数
- smoothing from reference position(使用参考位置平滑):多周期或准周期运动时,确保始终从同一参考位置进行平滑时,网格质量可能在不同循环之间保持更加一致。但在 layering 和 remeshing 开启时不能用
2. Spring
- 用于边界主要向一个方向移动或运动主要垂直于边界区的情况
3. Linearly Elastic Solid Based Smoothing Method:消耗算力更多,但对于某些网格和网格运动,可以更好保持网格质量。
- 适用于diffusion 可以处理的大部分模型,特别是旋转运动的网格。但在具有旋转边界运动和尖角的情况下,使用 diffusion 可能更好。
- 不能处理多面体网格以及有悬挂节点的网格。这种情况建议使用 diffusion 方法
- 将施加的变形转移到变形网格的内部,就好像网格是一个具有给定材料特性的线性弹性固体一样。可以指定几何体使网格沿其滑动
- 可设置泊松比 Poisson’s Ratio ,允许范围在-1.0和0.5之间。
2. Dynamic Layering
- 与移动边界相邻的单元层分割或合并。用于棱柱网格(六面体或者棱柱层)。
- 除非使用 sliding interfaces,否则单元层必须完全被单边面区域包围.
- 如果移动边界是一个内部区域,则可以考虑将两侧的单元使用不同 h i d e a l h_{ideal} hideal的 Dynamic Layering。
- 在每个移动边界上指定一个 h i d e a l h_{ideal} hideal。
- 当 j 层网格扩张时,网格会扩张,直到
h
m
i
n
>
(
1
+
a
s
)
h
i
d
e
a
l
h_{min} >(1+a_s)h_{ideal}
hmin>(1+as)hideal。使用基于高度的选项,单元格将被分裂,以创建一层具有恒定高度的单元格和一层具有高度的单元格。压缩同理。
-
h
i
d
e
a
l
h_{ideal}
hideal 允许定义为常数值或随时间或比率的函数而变化的值
-
h
i
d
e
a
l
h_{ideal}
hideal 允许定义为常数值或随时间或比率的函数而变化的值
3. Remeshing
- 用于三角形或四面体单元,仅在三维网格可以在边界层网格中重建四面体单元和楔形网格。
- 强烈建议在网格运动范围大,有旋转,变形复杂时与 Dynamic Smoothing 同时开启
- unified remeshing
- 通常建议使用该方法,它结合了各种重新网格方法的各个方面,默认情况下尝试保持初始网格大小分布
- Methods-based remeshing
-
在超过指定偏斜度skewness。大于指定最大尺寸/小于指定最小尺寸,不符合规定的长度比例时重新划分网格
- Methods
- Local cell:默认同时会启用local zone remeshing
- Local face
- Region face
- 2.5D
- Sizing remeshing interval
- 对最大Skewness 的检测与重新划分是每个时间步都进行,而最大/小尺寸的检测则由用户定义Sizing remeshing interval
- Sizing function:根据函数生成网格而不是基于最小和最大长度比例来标记单元格。只能用于local cell 与 2.5D 方法
- Resolution:决定背景箱的大小。默认情况下,2D问题设为3, 3D问题设为1
- Variation: 控制边界内部单元相对于其最近的边界单元的大小。可设置范围为[-1,+∞]。值为0表示远离边界的固定大小分布,值为-0.5表示边界内部单元可以是最近边界单元格的大小的一半。
- Rate:控制变化速度可设置范围为[-0.99, 0.99]。正值表示速度较慢。
- Methods
未使用 Sizing function | 使用 Sizing function |
---|---|