数字化研发设计2--(CDF-欧拉多相流、DDPM-密集离散相、VOF-多流体、MixTure-均质模型 再湿法行业应用)
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项目中取得了显著的经济效益。
openEuler 是由开放原子开源基金会孵化的全场景开源操作系统项目,面向数字基础设施四大核心场景(服务器、云计算、边缘计算、嵌入式),全面支持 ARM、x86、RISC-V、loongArch、PowerPC、SW-64 等多样性计算架构
更多推荐
所有评论(0)