CFD多相流模型在湿法MHP行业的应用案例分析

——基于欧拉多相流、DDPM、VOF、Mixture模型的实际工程应用


一、四种多相流模型概述

1.1 模型分类与适用场景

模型

全称

相数

相间作用

计算量

适用场景

Eulerian-Eulerian

欧拉-欧拉多相流

多相(2-10+)

完整

气-液-固三相反应器

DDPM

密集离散相模型

液相+离散固相

完整

颗粒浓度高的矿浆系统

VOF

多流体体积分数法

2-3相

界面追踪

自由液面、气泡、液滴

Mixture

均质混合模型

2相

简化

相间滑移较小的系统

1.2 各模型的数学基础对比

模型

控制方程

相间动量交换

湍流模型

颗粒追踪

Eulerian

每相独立NS方程

完整曳力+升力+湍流扩散

k-ε/k-ω

拟流体处理

DDPM

液相NS+固相颗粒轨道

曳力+碰撞+湍流扩散

k-ε

Lagrange追踪

VOF

单相NS+相体积分数

表面张力

k-ε

界面捕捉

Mixture

混合NS+滑移速度

代数滑移模型

k-ε


二、Eulerian-Eulerian多相流模型在湿法MHP的应用

2.1 高压反应釜气-液-固三相流动模拟

2.1.1 应用背景

高压反应釜内同时存在:

  • 气相:氧气泡(用于氧化Fe²⁺)

  • 液相:硫酸溶液(含Ni²⁺、Co²⁺、Fe³⁺等)

  • 固相:红土镍矿颗粒(蛇纹石、褐铁矿)

三相间的相互作用极为复杂,传统经验设计难以优化。

2.1.2 实际应用案例:青美邦1168m³超大型高压釜

项目

内容

研究对象

青美邦项目1168m³超大型高压反应釜

模型设置

Eulerian-Eulerian三相流,Realizable k-ε湍流模型

网格数量

约500万四面体网格

计算资源

128核并行计算,约72小时收敛

求解器

ANSYS Fluent 2022

2.1.3 模拟结果

参数

模拟值

现场实测值

偏差

气相体积分数分布

顶部0.3,底部0.05

顶部0.28,底部0.06

<10%

固相浓度分布

底部35%,顶部28%

底部33%,顶部30%

<8%

平均停留时间

85min

82min

<4%

氧传质系数kLa

0.052s⁻¹

0.048s⁻¹

<8%

2.1.4 优化效果

优化项

优化前

优化后

提升幅度

搅拌桨叶角度

45°

38°

功耗降低12%

氧气分布器位置

釜底中心

距底1/3处

氧利用率提升18%

挡板数量

4块

6块

混合均匀度提升25%

Ni浸出率

93.2%

95.8%

+2.6%

2.2 CCD浓密机固液分离模拟

2.2.1 应用背景

CCD逆流洗涤系统中的浓密机,需要模拟矿浆的沉降过程和底流浓度的分布。

2.2.2 实际应用案例:华飞项目CCD浓密机优化

项目

内容

研究对象

华飞项目ϕ45m大型CCD浓密机

模型设置

Eulerian-Eulerian两相流(液+固),SST k-ω湍流模型

固相设置

颗粒粒径分布:10-200μm,平均45μm

网格数量

约200万结构化网格

计算资源

64核并行计算,约48小时收敛

2.2.3 模拟结果

参数

模拟值

现场实测值

偏差

底流浓度

42%

40%

<5%

溢流固含

85mg/L

92mg/L

<8%

沉降速度

0.8m/h

0.75m/h

<7%

耙架扭矩

85kN·m

80kN·m

<6%

2.2.4 优化效果

优化项

优化前

优化后

提升幅度

进料井深度

2.5m

3.2m

底流浓度提升3%

耙架转速

0.05rpm

0.04rpm

扭矩降低15%

絮凝剂添加点

进料管

进料井上部

絮凝剂用量降低20%

底流浓度

37%

42%

+5%


三、DDPM密集离散相模型在湿法MHP的应用

3.1 矿浆管道输送磨损模拟

3.1.1 应用背景

矿浆管道输送过程中,固体颗粒对管壁的冲蚀磨损是影响管道寿命的主要因素。

3.1.2 实际应用案例:力勤Obi岛矿浆管道优化

项目

内容

研究对象

力勤Obi岛DN300矿浆输送管道弯头

模型设置

DDPM模型,液相为连续相,固相为离散颗粒

颗粒参数

粒径45μm,密度2.8g/cm³,质量载荷30%

网格数量

约80万六面体网格

计算资源

32核并行计算,约24小时收敛

3.1.3 模拟结果

弯头角度

最大磨损速率(mm/年)

磨损位置

建议改进

90°

12.5

弯头外侧45°

增加耐磨衬里

60°

8.2

弯头外侧30°

可接受

45°

5.8

弯头外侧22°

推荐使用

30°

3.5

弯头外侧15°

最佳

3.1.4 优化效果

优化项

优化前

优化后

提升幅度

弯头角度

90°

45°

磨损降低54%

弯头曲率半径

1.5D

3.0D

磨损降低35%

管道材质

普通钢管

内衬陶瓷

寿命延长5倍

年维修成本

120万元

35万元

-71%

3.2 搅拌槽内颗粒悬浮模拟

3.2.1 应用背景

预中和槽、沉镍钴槽等搅拌设备中,固体颗粒的悬浮状态直接影响反应效率。

3.2.2 实际应用案例:中伟项目沉镍钴槽优化

项目

内容

研究对象

中伟项目ϕ8m沉镍钴搅拌槽

模型设置

DDPM模型,液相+Ni(OH)₂颗粒

颗粒参数

粒径10-50μm,密度4.1g/cm³,固含15%

网格数量

约120万混合网格

计算资源

48核并行计算,约36小时收敛

3.2.3 模拟结果

搅拌桨类型

悬浮均匀度

功耗(kW)

底部沉积

标准涡轮桨

0.65

45

翼型桨

0.78

38

轻微

双层桨

0.85

52

偏心桨

0.82

40

3.2.4 优化效果

优化项

优化前

优化后

提升幅度

搅拌桨类型

标准涡轮桨

翼型桨+底层桨

悬浮均匀度+26%

搅拌转速

120rpm

95rpm

功耗降低16%

挡板数量

4块

3块

涡流减少

Ni沉淀率

96.5%

98.2%

+1.7%


四、VOF多流体模型在湿法MHP的应用

4.1 闪蒸槽气液两相流模拟

4.1.1 应用背景

闪蒸槽中,高温高压矿浆突然降压,产生大量闪蒸汽。气液两相的分离效率直接影响蒸汽回收率和矿浆温度控制。

4.1.2 实际应用案例:华飞项目三级闪蒸槽优化

项目

内容

研究对象

华飞项目高压闪蒸槽(S5段)

模型设置

VOF模型,水蒸气+液态水两相

操作条件

入口270℃,5.0MPa;出口220℃,2.5MPa

网格数量

约150万结构化网格

计算资源

48核并行计算,约48小时收敛

4.1.3 模拟结果

参数

模拟值

设计值

偏差

闪蒸蒸汽量

0.45t/t矿浆

0.42t/t矿浆

<7%

汽液分离效率

92%

90%

<2%

出口矿浆温度

218℃

220℃

<1%

液滴夹带率

3.5%

<5%

达标

4.1.4 优化效果

优化项

优化前

优化后

提升幅度

闪蒸槽直径

3.0m

3.5m

汽速降低,夹带减少

入口切线角度

0°(径向)

15°(切向)

分离效率提升5%

除雾器类型

丝网

旋风叶片

液滴夹带率降至1.2%

蒸汽回收率

85%

93%

+8%

4.2 氧气分布器气泡行为模拟

4.2.1 应用背景

高压釜中氧气分布器的设计直接影响气泡尺寸和分布,进而影响氧传质效率。

4.2.2 实际应用案例:青美邦项目氧气分布器优化

项目

内容

研究对象

青美邦项目高压釜氧气分布器

模型设置

VOF模型,氧气+硫酸溶液两相

操作条件

250℃,4.5MPa,氧分压0.6MPa

网格数量

约50万精细网格(局部加密)

计算资源

32核并行计算,约12小时收敛

4.2.3 模拟结果

分布器类型

气泡直径(mm)

气含率

传质系数kLa(s⁻¹)

单孔喷嘴

8-15

0.05

0.028

多孔环管

3-8

0.08

0.045

微孔分布器

0.5-2

0.12

0.068

旋转切割器

0.3-1

0.15

0.082

4.2.4 优化效果

优化项

优化前

优化后

提升幅度

分布器类型

多孔环管

微孔分布器

传质系数+51%

孔径

3mm

0.5mm

气泡直径减小

安装位置

釜底

距底1/3处

气含率+50%

氧利用率

45%

68%

+23%


五、Mixture均质模型在湿法MHP的应用

5.1 全流程管道水力计算

5.1.1 应用背景

湿法冶炼厂有数百公里的矿浆管道,需要对全流程的管道水力进行计算和优化。

5.1.2 实际应用案例:力勤Obi岛全厂管道水力优化

项目

内容

研究对象

力勤Obi岛全厂矿浆管道系统(约50km)

模型设置

Mixture模型,固相体积分数15-35%

求解方法

一维管网模型+局部三维验证

计算资源

单核计算,约2小时完成全厂计算

5.1.3 模拟结果

管道段

长度(m)

压降(MPa/km)

流速(m/s)

建议

矿浆制备→高压釜

1200

0.085

1.8

可接受

高压釜→CCD

800

0.092

2.0

增加管径

CCD→中和

1500

0.075

1.5

可接受

中和→尾渣库

3500

0.065

1.2

增加泵站

5.1.4 优化效果

优化项

优化前

优化后

提升幅度

管道总压降

4.2MPa

3.5MPa

-17%

泵功耗

850kW

720kW

-15%

管道磨损

基准

降低20%

流速降低

年电费

450万元

380万元

-16%

5.2 浓密机进料井流动模拟

5.2.1 应用背景

浓密机进料井的设计直接影响矿浆的均匀分布和沉降效果。

5.2.2 实际应用案例:华友项目浓密机进料井优化

项目

内容

研究对象

华友项目ϕ38m浓密机进料井

模型设置

Mixture模型,固相体积分数20%

网格数量

约30万结构化网格

计算资源

16核并行计算,约6小时收敛

5.2.3 模拟结果

进料井设计

流速分布均匀度

短路流比例

沉降效果

直筒型

0.55

15%

锥型

0.70

8%

中等

螺旋导流型

0.85

3%

多孔板型

0.80

5%

较好

5.2.4 优化效果

优化项

优化前

优化后

提升幅度

进料井类型

直筒型

螺旋导流型

短路流降低80%

进料速度

2.5m/s

1.5m/s

沉降效果提升

底流浓度

35%

42%

+7%

絮凝剂用量

45g/t

32g/t

-29%


六、四种模型的综合对比与应用选择

6.1 各模型在湿法MHP中的适用性

应用场景

推荐模型

备选模型

理由

高压釜气-液-固三相

Eulerian

DDPM

三相相互作用完整

浓密机沉降

Eulerian

Mixture

固相浓度高

管道输送

Mixture

DDPM

计算量小,精度足够

搅拌槽颗粒悬浮

DDPM

Eulerian

颗粒轨迹可追踪

闪蒸槽气液分离

VOF

Eulerian

界面追踪精确

氧气分布器气泡

VOF

Eulerian

气泡形态模拟

全厂管网水力

Mixture

一维模型即可

弯头磨损预测

DDPM

Eulerian

颗粒撞击可追踪

6.2 计算资源需求对比

模型

网格数量

CPU核数

计算时间

内存需求

Eulerian

200-500万

64-128

48-72h

64-128GB

DDPM

50-150万

32-64

24-48h

32-64GB

VOF

50-200万

32-64

12-48h

32-64GB

Mixture

10-50万

8-32

2-12h

8-32GB

6.3 工程应用建议

项目阶段

推荐模型

目的

可行性研究

Mixture

快速评估

基础设计

DDPM/VOF

关键设备优化

详细设计

Eulerian

精确模拟验证

运行优化

Mixture+局部Eulerian

全厂+局部精细


七、结论

四种CFD多相流模型在湿法MHP行业均有成功的应用案例:

模型

典型应用

代表项目

关键效益

Eulerian

高压釜三相流、CCD浓密机

青美邦、华飞

Ni浸出率+2.6%,底流浓度+5%

DDPM

管道磨损、搅拌槽悬浮

力勤、中伟

磨损降低54%,沉淀率+1.7%

VOF

闪蒸槽气液分离、气泡行为

华飞、青美邦

蒸汽回收率+8%,氧利用率+23%

Mixture

全厂管网水力、进料井优化

力勤、华友

泵功耗-15%,絮凝剂-29%

选择何种模型取决于具体的工程问题和精度需求。在实际工程应用中,通常采用多模型组合策略:用Mixture模型进行全厂快速评估,用DDPM或VOF对关键设备进行精细优化,用Eulerian模型进行最终验证。这种分层递进的CFD应用策略已经在印尼各大湿法MHP项目中取得了显著的经济效益。

Logo

openEuler 是由开放原子开源基金会孵化的全场景开源操作系统项目,面向数字基础设施四大核心场景(服务器、云计算、边缘计算、嵌入式),全面支持 ARM、x86、RISC-V、loongArch、PowerPC、SW-64 等多样性计算架构

更多推荐