稀土钙皂化流程仿真及参数优化研究_邵留宝.pdf

Total No 280December 2022冶金设备METALLURGICAL EQUIPMENT总第 280 期2022 年 12 月第 6 期稀土钙皂化流程仿真及参数优化研究邵留宝唐焱梁冰梁才航(桂林电子科技大学机电工程学院广西桂林 541004)摘要钙皂化是稀土提取过程的重要工艺，而实际工艺生产缺少对其反应机理模型的研究，影响生产质量及效率的提升。考虑钙皂化搅拌釜结构、工艺参数的情况下，构建了钙皂化反应机理模型，实现了对钙皂化生产流程的实时仿真分析。进一步，以提升皂化度、减少皂化时间和功耗为多目标，以搅拌桨类型、搅拌转速、搅拌桨位置和进料位置为设计参数，通过正交实验优化方法优化钙皂化过程，得到优化设计参数数值。与初始状态相比，优化后皂化度提升了 1.8%，皂化时间、功耗分别下降了 2.9%、15.2%。皂化度、皂化时间、功耗仿真结果与实验数据差值分别为 4.96%、6.00%、5.96%，验证了仿真分析的有效性。所建钙皂化反应机理模型准确预测了钙皂化反应过程，其参数优化方法有利于生产质量及效率的提升，为钙皂化工艺过程的数字化设计提供了参考。关键词稀土工艺钙皂化正交试验仿真中图法分类号TF111 3文献标识码ADoi:10.3969/j.issn.1001 1269.2022.06.012Simulation and Parameter Optimization of Rare EarthCalcium Saponification ProcessShao LiubaoTang YanLiang BingLiang Caihang(College of Mechanical and Electrical Engineering，Guilin University of Electronic Technology，Guilin，541004)ABSTRACTCalcium saponification is an important process in the extraction process of rare earths，but theactual process production lacks research on its reaction mechanism model，which affects the improvement ofproduction quality and efficiencyTaking into account the structural process parameters of the calciumsaponification stirred tank，a calcium saponification reaction mechanism model was constructed，and real-timesimulation analysis of the calcium saponification production process was realizedFurthermore，the calciumsaponification process is optimized by the orthogonal experiment optimization method to obtain the optimized designparameter values，the multi-objective of the optimization method is to improve the saponification degree，reduce thesaponification time and power consumption，and the design parameters are the type of stirring blade，the stirringspeed，the position of the stirring blade and the feeding position Compared with the initial state，the optimizedsaponification degree increases by 1.8%，and the saponification time and power consumption decrease by 2.9%and 15.2%respectivelyThe difference between the saponification degree，saponification time，and powerconsumption simulation results and experimental data are 4.96%，6.00%，and 5.96%，respectively，whichverifies the validity of the simulation analysisThe built calcium saponification reaction mechanism modelaccurately predicts the calcium saponification reaction process，and its parameter optimization method is benefit toimprove the production quality and efficiency，it provides a reference for the digital design of the calciumsaponification processKEYWORDSRare earth processCalcium saponificationOrthogonal experimentSimulation55基金项目:国家自然科学基金地区科学基金项目(51566002)作者简介:邵留宝，男，1996 年生，硕士研究生，主要研究方向:稀土钙皂化自动化及配套技术研究，邮箱:1461019056 qq com通信作者:唐焱，邮箱:2283640182 qq com1前言稀土金属是现代制造领域不可取代的战略资源。随着全球科技的持续进步，电子、冶金、军工等领域对稀土金属产量和品质需求急剧增长，加速了生产装备迭代更新速度，推动我国稀土装备自动化、智能化升级已成业内重要问题1。湿法冶金术是各类稀土矿前期处理的主要方法，其中钙皂化是最重要的工序之一，钙皂化质量对稀土的产出有较大影响。近年来，国内外学者对于湿法冶金仿真、优化进行了大量研究，何泳辉2 开发了一款基于网络的湿法冶金全过程监测仿真平台，将湿法冶金生产过程的动态特性通过建立数学模型、编程在计算机上运行，实现了对湿法冶金生产过程的监测，但该仿真平台缺少钙皂化机理模型，因此不能实现钙皂化的过程监测。杨秉宇3 基于粗糙集的数据挖掘方法对积累的现场数据进行优化控制，通过补偿各工序指标及操作变量，实现对湿法冶金全流程的优化，但该优化策略影响钙皂化生产的参数较少，不能实现钙皂化过程的优化。由于钙皂化工艺过程伴随大量复杂的物理、化学反应，涉及原料成分、工艺参数、设备构造、辅材品质等要素，导致难以建立准确的仿真模型，使得对于钙皂化机理模型及工艺参数优化研究较少。因此，目前钙皂化设备关键工艺参数只能依赖人工经验确定，并通过离线检测皂化度即钙皂有机相中 CaA2的含量来重复校正反应结果。基于生产过程影响因素的复杂性，如何建立钙皂化工艺的反应机理模型，成为钙皂化生产质量、效率提升的主要瓶颈。针对缺少钙皂化工艺反应机理模型这一问题，考虑钙皂化搅拌釜的结构工艺参数，构建钙皂化反应机理模型，实现对钙皂化生产流程的实时仿真分析4。并通过正交实验优化方法优化钙皂化过程5，得到搅拌桨类型、搅拌转速、搅拌桨位置和进料位置等参数的优化数值，提升了钙皂化生产的质量及效率。2模型及验证21钙皂化反应机理(1)工艺过程根据企业生产现状在反应釜中加入待皂化原料，即 空 白 有 机 相(HA)，浓 度 2mol/L、总 量2.5m3;按 1L/s 流速逐渐注入浓度为 10mol/L 碳酸钙浆液 0.15m3 0.2m3，同时启动搅拌系统。各相液体在釜内混合发生钙皂化反应，生成钙皂有机相(CaA2);随着釜内皂化反应的持续，钙皂有机相总量不断增加，当钙皂有机相浓度达到0.42mol/L 0.45mol/L 时，钙皂化过程结束。利用钙皂有机相、废液、水三相的密度差异，将钙皂有机相分离并经管道送至下工序车间。钙皂化工艺主要参数如下表 1 所示。表 1钙皂化工艺主要参数碳酸钙浆液浓度 mol/L碳酸钙浆液加入量 m3空白有机相浓度 mol/L100.15 0.22空白有机相总量 m3空白有机相加入速度 L/s钙皂有机相浓度 mol/L2.510.42 0.45mol/L(2)反应机理钙皂化主体反应体系是碳酸钙和空白有机相在常温下的竞争串联反应，该主反应体系又分为主反应和副反应两部分，化学方程式如下:主反应:2HA+CaCO3k1CaA2+H2O+CO2副反应:CaCO3+CaA2k2CaHA322钙皂化流体力学模型CFD(Computational Fluid Dynamics)分析是基于质量守恒方程和动量守恒方程(Navier-Stokes)方程。对于搅拌釜中的湍流流动，需将湍流动量表示为时均值和脉动值的叠加，此时的动量守恒方程称为雷洛方程。因为雷洛应力的存在，方程组不能封闭，所以在方程组中添加标准 湍流模型作为补充方程6。对于搅拌釜中的钙皂化模拟，要求解雷洛方程、质量守恒方程和标准 湍流方程。雷洛方程:ujujxj=pxi+2uixj2+fi uiujx()j(1)式中:i，j=1，2，3;ujujxj对流项;pxi压力梯度;652022 年 12 月第 6 期总第 280 期冶金设备2uixj2扩散项;fi体积力;uiujx()j雷诺应力。质量守恒方程:t+(u)x+(v)y+(w)z=0(2)式中:流体密度;t时间，min;u、v、w在 x、y、z 方向上的速度分量。标准 湍流方程是在湍动能 的基础上，再增加一个湍动能耗散率 的方程。单位质量流体湍流脉动动能:=12u2+v2+w()2(3)单位质量流体湍流脉动动能耗散率:=ux()kuix()k(4)23反应釜三维建模钙皂化设备的核心部件为反应釜，其结构如图 1 所示。其主体呈圆柱形，内壁垂直方向均布一定宽度的四块挡板，釜体中央内置直叶涡轮搅拌桨，搅拌桨底部相对釜底距离 600mm。企业现役反应釜结构工艺参数如表 2 所示。图 1钙皂化搅拌釜模型24数值仿真(1)边界条件运用 meshing 工具对反应釜模型进行网格划分，生成非结构化网格。网格划分时将模型分为动区域和静区域7，划分后网格总数 427034。运用单元质量评定法检测网格划分结果，单元质量指数达到 91.5%，高于一般要求 75%，即划分的网格可用于后续研究。表 2搅拌器结构、工艺参数釜体直径/mm搅拌轴距离釜底距离/mm搅拌桨直径/mm搅拌轴直径/mm150060073380液位高度/mm搅拌桨桨叶高度/mm搅拌桨转速/rpm加料口直径/mm150013370100挡板宽度/mm加料位置釜体高度/mm167釜顶挡板处2000运用滑移网格法对模型进行分析计算，计算域分成动、静两部分，即旋转的搅拌桨和静止的釜体和挡板;两部分的网格在 Fluent 中通过定义非匹配边界条件实现相互独立，外部的网格静止，而内部的网格随搅拌桨转动，两部分网格之间通过通过 interior 面进行数据传递。(2)数据计算方法数值仿真严格采用 2 1 中所述实际生产过程工艺参数。当注入碳酸钙浆液时搅拌系统开始以70rpm 匀速转动，直到规定的皂化反应结束时间终止。通过 Fluent 软件功能，规定以所有计算变量的残差小于 104为收敛条件;运算初期，