自适应高度表格_基于自适应网格的浸漆工艺仿真

1、物理问题

浸漆工艺是工业生产常用工艺流程。为了保证浸漆后表面成膜均匀，需要设计合理的浸润轨迹、通过速度等参数，在STAR-CCM+中非常快速方便的修改相关参数，优化筛选出最优的工艺流程。特别是STAR-CCM+2020版本的自适应网格推出后，极大提升了该类计算的精度。

2、几何及网格

计算域为一油漆池，总长度约为46m，油漆位于在池底部，高度为1.92m，其余区域为开放大气。长4.8m的桁架以0.5m/s的速度沿给定轨迹穿过浸润池，整个工艺过程如图1所示。计算域分为两部分，油漆池作为背景域，桁架为嵌套区域，并对局部模型细节进行加密，均采用Trim网格，如图2、3所示。

图2 背景域网格

图3 嵌套网格

3、物理模型

计算中运用的全部模型如图4所示。

图4 计算物理模型(Adaptive Mesh为2020版本新增功能)

计算域兼有空气和油漆，因此采用VOF多相流模型，新建空气和油漆两相介质。

图5 气液两相定义

STAR-CCM+具备针对船舶行业的专业波浪模型VOF Wave，本算例使用Wave模型中最简单的Flat Wave，定义水平面高度。右键新建Flatwave,定义好液面高度，重力方向，VOF两项密度等参数，本例中油漆池底部Z坐标-4m，定义池内液位自由面Z坐标-1.92m，即液位高度共计1.08m。

图6液位定义

2020版本新增加的Adaptive Mesh模型，简单介绍AdaptiveMesh的可调属性参数，主要包括Properties和Expert两部分。

图7 自适应网格定义

Properties中主要定义自适应的尺寸，包括：

Prism Cell Refinement定义是否对边界层内进行自适应网格划分；Transition Width类似于Volume Growth Rate的概念，定义自适应网格向外增长的速率，效果如下图所示；

图8 Transitionwidth定义

Limit Cell Size主要定义自适应后最小网格尺度，勾选后会自动出现MinAdaption Cell Size选项，设置最小尺度；Expert中主要调整自适应后的网格质量，一般无需做特别设置。展开Adaptive Mesh选项，右键新建可以看到一共有三种AdaptiveMesh方法，注意只有选择VOF模型，具备Overset交界面时才会有选项，一般流动换热计算只有User-Defined一个选项。

图9 自适应网格类型

新建准则后，属性栏里主要有两个可调参数，

图10 自定义网格自适应准则

Max Refinement Level即最大自适应层级，对于Trim网格单次自适应可以剖分出8个新网格单元，对于多面体，单层自适应可剖分12~15个网格单元。默认的2次自适应理论上可以对部分Trim网格加密64倍，对多面体网格加密144~225倍，可以满足大部分应用场景。

图11 Trim和Poly网格自适应示意

最后在Field Function里定义自适应准则，一般推荐特征值的梯度，定义自适应区间，区间内部保持网格，大于上限加密，小于下限加粗。

图12 自适应准则

本算例会用到右键新建Free Surface和Overset两个自适应规则，分别针对VOF交界面和Overset进行网格自适应加密，系统已经做好设置定义，保持默认选项即可。

4、边界条件

本算例边界条件比较简单，底部定义壁面，前后设置对称边界，其他区域为大气环境，全部定义压力出口。

最后需要定义桁架的运动轨迹，本算例采用运动轨迹定义方法。提前在Excel表格中定义好运动轨迹，命名为Track，编辑为csv格式，导入Table表格中。

图13 Trajectory表格定义

在Motion中新建Trajectory，选择Track表格，即为桁架运动轨迹，定义速度函数，本算例中速度为0.5m/s。

图14 运动轨迹定义

最后新建SuperposingMotion选择Constrain Rotation，定义Horizontaloffset为桁架长度4.8m，即运动过程中保持4.8m为直线。

图15 Constrain Rotation运动定义

将定义好的运动赋予Overset区域，即可完成所有运动定义。

图16 Overset区域运动定义

5、后处理

 6、总结

自适应网格功能是STAR-CCM+2020版本新推出的功能，对于VOF和Overset类问题，应用场景比较固定，系统默认的设置对该类问题适用性较好。但是对于普通流动换热问题，自适应标准的选择非常重要，不合理的标准选择会导致网格迅速增加，推荐在使用前先进行测试。