本案例演示利用Fluent中的DPM模型模拟固体颗粒物对结构物壁面的冲蚀。

1 问题描述

教程中使用的几何形状最初是作为侵蚀工况下节流阀的减压装置设计的，该装置如下所示。

该装置使用砂液侵蚀试验设备（测试参数见文献 Wallace ^[1]所述）在侵蚀条件下进行了测试，本教程中使用的条件是根据实际测试得出的。在侵蚀试验中，限流器安装在直径为 53.1 mm的管道之间。分别在限流器的上游 106 mm及下游 318 mm位置测量静压。测试过程中，通过定期称量限流器的重量来获得总体质量损失。

2 Fluent设置

2.1 启动Fluent并导入网格

• 以3D、Double Precision方式启动Fluent

• 利用菜单File → Import → Tecplot… 读取网格Erosion.plt

  （文件下载链接在文末）

局部网格如图所示。

2.2 General设置

• 鼠标双击模型树节点General，右侧面板点击按钮Check检查网格

确保计算区域的尺寸符合要求。

2.3 Models设置

1、设置湍流模型

• 双击模型树节点Models → Viscous 打开湍流模型设置对话框，如下图所示选择Realizable k-epsilon湍流模型

2、Discrete Phase模型

• 双击模型树节点Models → Discrete Phase打开模型设置对话框

• 进入Physical Models选项卡，如下图所示激活模型
• 激活选项Virtural Mass Force及Pressure Gradient Force
• 激活选项Erosion/Accretion

注: 对于本案例, 也可以不选择选项 Interaction with Continuous Phase，从而使用单向耦合计算。

”

2.4 设置Injector

• 右键选择模型树节点 Injections，点击弹出菜单项New 打开入射器设置对话框

• 如下图所示设置入射器参数

• 如下图所示设置激活随机游走模型

2.5 Materials设置

• 添加材料介质 water-liquid

• 修改材料 anthracite 的密度

2.6 区域设置

• 指定计算区域内的材料介质为 water-liquid

2.7 边界条件设置

修改边界类型。修改完毕后的边界类型如下图所示，包含1个速度入口，1个压力出口，1个对称边界，其他为壁面边界

1、指定入口边界inflow

• 指定入口速度为 12.901 m/s

• 指定入口边界 DPM 类型为 escape

2、出口边界outflow

• 采用静压为 0 的出口边界

• 指定出口的 DPM 类型为 escape

3、壁面restrictor

案例只考虑 restrictor 边界的冲蚀情况，需要指定该边界上的恢复系数及冲蚀模型。打开边界条件设置对话框的 DPM 选项卡， 如下图所示。

• 进入 Normal 右侧的 Edit... 按钮，如下图所示设置法向恢复系数

• 相同方式设置切向恢复系数

• 设置 Oka 冲蚀模型

2.8 Methods设置

• 计算方法如下图所示

2.9 残差设置

• 将连续方程残差修改为 1e-5

2.9 初始化计算

• 采用 Hybrid 初始化

2.10 进行计算

• 迭代计算 500 步

• 残差如下图所示

3 计算结果

• 对称面上速度分布

• 查看壁面冲蚀

• 壁面冲蚀率如图所示

• 将壁面冲蚀率对面积求积分

计算采用的是四分之一模型，因此整体模型的冲蚀率为：

4 利用DNV模型进行计算

• 打开壁面边界restrictor的参数设置对话框，进入DPM选项卡，激活选项DNV，点击Edit...按钮打开设置对话框

• 采用默认模型参数

重新进行计算。计算完毕后的冲蚀率显示如下图所示。最大冲蚀率为1.95E-3 kg/(m2.s)。

面积积分得到的冲蚀率为4.081656E-8 kg/s，如下图所示。

DNV模型默认参数在固体颗粒物冲蚀碳钢材料时表现非常好。

注：本案例几何模型及参数取自STAR CCM+的Tutorials文件。计算案例：https://pan.baidu.com/s/1Wfo5p4wFUyLP4UpVb_iSQg?pwd=4eek 提取码：4eek

”

参考文献