本算例演示利用STAR CCM+中的网格重构功能模拟仿真摆线泵内部流场。

1 问题描述

摆线泵是一种包含有内转子和外转子的装置，其内转子和外转子具有与其齿数相关的不同转速，当转子旋转时，通过转子间动态变化的容积将流体从入口输送到出口。

摆线泵模拟的最大问题是处理转子之间的狭小间隙。转子旋转会导致小间隙中的网格质量降低。在STAR CCM+中可以使用Remeshing模型处理此类问题。

案例所模拟的摆线泵外转子包含12个齿，内转子包含11个齿，其的转速与齿数满足以下函数关系：

几何模型如下图所示。

2 STAR CCM+设置

• 启动STAR CCM+

2.1 导入文件

• 利用菜单 File > Load导入仿真文件 gerotor_start.sim

几何模型如下图所示。

放大图形可以看到齿间存在小的间隙。

注：计算过程中，两齿不能完全接触。

”

2.2 设置网格参数

导入的仿真文件中已经预置了进出口部分的网格参数，这里只需要指定两个转子的网格参数即可。

当网格重构时，网格生成管道需确保泵处于正确的位置。此外，齿间间隙较小，可以应用Directed Mesh方法以减少轴向方向上的网格数量。为摆线体定义两个运动变换操作，一个布尔减操作及一个Directed Mesh操作。

1、定义两个旋转运动

• 右键选择模型树节点 Tools > Motions，点击弹出菜单项 New → Rotation，修改新增的节点名称为 Rotation_inner

• 选中节点 Rotation_inner，如下图所示设置
• 指定 Axis Direction为**[0,0,1]**
• 指定 Axis Origin为**[0,0,0]**
• 指定 Rotation Specification为 Rotation Rate

• 选中节点 Motions > Rotation_inner > Rotation Rate，指定 Rotation Rate为 1090.90909 rpm

• 相同方式创建另一个运动，将其命名为 Rotation_outer，并指定如下图所示参数

2、定义区域旋转

• 右键选择模型树节点 Geometry > Operations ，点击弹出菜单项 New > Surface Preparation > Transform创建变换

• 在弹出的对话框中选择 Inner并点击OK按钮确认选择

• 修改节点 Transform名称为 Transform_inner，右键选择节点 Transform_inner > Transforms，点击弹出菜单项 New > Motion创建新节点

• 选中节点 Motion，指定参数 Motion为 Rotation_inner

• 相同的方式创建另一个运动，命名为 Transform_outer，指定其Part为 Outer，其运动为 Rotation_outer

3、构造计算区域

• 右键选择节点 Operations，点击弹出菜单项 New → Boolean → Substract打开操作对话框

• 如下图所示，选择 Input Parts为 Inner与Outer，指定 Target Part为 Outer，激活选项 Perform CAD Boolean，点击OK按钮确认选择

• 如下图所示，选中节点 Tessellation Options，指定参数为 Very Fine

• 选中模型树节点 Parts > Subtract ，将节点名称修改为 PumpFluid

4、设置网格参数

• 右键选中模型树节点 Operations，点击弹出菜单项 New → Mesh → Directed Mesh添加扫略网格

• 弹出的对话框中选择部件 PumpFluid，如下图所示

• 指定 Source Surfaces为 PumpFluid.Inner.top

• 指定 Target Surfaces为 PumpFluid.Inner.bottom

• 右键选择节点 Source Meshes，点击弹出菜单项 Automated Source Mesh

• 弹出对话框中选择 PumpFluid，指定网格生成方法为 Quadrilaterial Mesher与 Prism Layer Mesher

• 选中节点 Source Meshes > Auto Mesh > Meshers > Prism Layer Mesher，如下图所示指定参数

• 进入节点 Default Controls，按下表所示指定参数

• 右键选中模型树节点 Mesh Distributions，点击弹出菜单项 New Volume Distribution创建体网格分布

• 如下图所示选择部件 PumpFluid，点击OK按钮确认选择

• 指定 Number of Layers为 10，生成10层边界层网格

2.3 创建部件接触

• 如下图所示，按键盘 CTRL键同时选中节点 Inlet > Surface > interface与 PumpFluid > Outer > top，点击鼠标右键并选择弹出菜单项 Create Weak In-Place Contact创建弱接触

• 相同方式选中节点 Outlet > Surface > interface与 PumpFluid > Outer > top创建另一对弱接触

• 右键选择节点 PumpFluid，点击弹出菜单项 Assign Part to Regions… 分配计算区域

• 如下图所示设置参数

• 右键选择节点 Operations，点击弹出菜单项 Execute All执行所有操作

生成的网格如下图所示。

2.4 设置时间步长

模拟运动时，时间步长应与网格运动相匹配。运动对象周围的流场随对象的运动而变化。因此时间步长应足够小，以捕获运动引起的变化。本案例要求解旋转过程中小间隙内的流动，因此选择时间步长大小以匹配一个旋转角度。当涉及复杂物理时，例如多相流，时间步长还应该进一步减小。

在本教程中，两个转子的配置如下所示：

• 选中模型树节点 Solvers > Implicit Unsteady，设置 Time-Step为 1/6000

2.5 设置网格重构

• 右键选择模型树节点 Physics 1，点击弹出菜单项 Select Models… 打开模型选择对话框

• 如下图所示选择选项 Remeshing

注：其它的选项按常规流动计算设置，这里已经预先选择完毕。

”

• 选中节点 Remeshing，如下图所示设置
• 指定 Mesh Operations为 Directed Mesh
• 选中激活选项 Repeating Meshes

• 选择节点 Remeshing > Repeating Meshes，激活选项 Enable Cyclic Matching

• 选中选项 Remeshing > Repeating Meshes > Time-Step Frequency，指定 Frequency为 30

• 选中节点 Solvers > Remeshing，激活选项 Enable Trigger

• 选中节点 Solvers > Remeshing > Time-Step Frequency的频率为 10

2.6 设置变形运动及边界条件

摆线转子壁面的运动通过变形器来实现。

为了避免小间隙中的网格发生极度变形，需要为PumpFluid区域的内壁面指定为Constraint变形器边界，使其沿内部几何部件的壁面滑动。通过让顶点自由滑动，与强制刚性位移相比，变形网格的强度更小。对于具有小间隙的变形问题，此方法可在不出现负单元格的情况下，在多个时间步上进行变形。

通常，变形过程必须在没有负网格的情况下至少在2-3个时间步内成功，否则网格重构求解器将激活恒定重构，而不会改善网格质量。建议对没有锐边的曲面使用约束方法。否则在旋转期间将无法正确表示曲面。

• 右键选择节点 Tools > Motions ，选择弹出菜单项 New > Morphing创建Morphing

• 选中节点 Regions > PumpFluid > Physics Value > Motion Specification，指定 Motion为 Morphing

• 选择节点 Regions > PumpFluid > Boundaries > Outer.wall > Physics Conditions > Morpher Specification，指定 Specification为 Displacement

• 选择节点 Regions > PumpFluid > Boundaries > Outer.wall > Physics Values > Morpher Rigid Boundary Motion，指定 Rigid Motion为 Rotation_outer

• 选中节点 Regions > PumpFluid > Boundaries > Inner.wall > Physics Conditions > Morpher Specification，指定 Specification为 Constraint

• 选中节点 Physics Conditions > Morpher Constraint Specification，指定 Constraint为 Part Surface

• 选择节点 Physics Values > Part Surface，指定 Part Surface为 Inner.wall

• 选择节点 Physics Values > Morpher Rigid Boundary Motion，指定 Rigid Motion为 Rotation_inner

• 同时选中节点 Outer.bottom, Outer.top, Outer.top[Inlet/PumpFluid], Outer.top[Outlet/PumpFluid]，点击右键菜单项 Edit… 打开设置面板

• 指定 Physics Conditions > Morpher Specification为 Constraint

• 同时选中节点 Outer.bottom及Outer.top，点击右键菜单项 Edit… 打开设置面板

• 指定参数 Physics Conditions > Reference Frame Specification为 Lab Frame

2.7 进行计算

• 选中节点 Stopping Criteria > Maximum Physical Time，取消激活选项 Enabled

• 选中节点 Maximum Steps，设置 Maximum Steps为 720

• 点击菜单 Solution → Run进行计算

3 计算结果

• 网格运动

• 压力随时间变化

本案例来自STAR CCM+随机自带案例，为最近几个版本新增案例。

类似这种狭缝结构在泵类模型中经常会碰到，目前可选的方法包括利用重叠网格、网格重构等。但这些方法都有各自的麻烦之处。

STAR CCM+的Remeshing功能计算效率较低，计算一段时间时间后要停下来重新生成网格，如果重生成的网格与之前的网格差距过大，会造成极大的插值误差，因此在使用此方法的过程中，尽量减小时间步长，以避免网格重构前后的网格节点差距过大。

文章来源：CFD之道