大喷嘴充气器充气性能的仿真研究

2023-09-11孙伟成史帅星董干国段莉莉

有色设备 2023年4期

吴峰 ,孙伟成 ,史帅星 ,董干国 ,段莉莉

(1.矿冶科技集团有限公司矿物加工科学与技术国家重点实验室,北京 100160;2.天津市天河计算机技术有限公司,天津 300457)

当下冶金化工等行业大量的反应设备都需要空气的参与,充气系统设计的优劣是影响反应设备性能的主要因素,因此提高相关反应设备充气系统的性能一直是该领域研究的重点和难点[1-2]。近年来,矿冶科技集团在多个项目的支持下,开发出压力式充气器并将其应用到冶金耗氧反应设备氰化槽和氧化槽充气系统中,在现场的工业试验和实际工程应用中均取得了较好的指标,受到业内广泛关注,应用领域在进一步扩展[3-6]。由于冶金及化工反应设备普遍高度较高,实验室研究无法匹拟现场实际高度,导致无法真实的获得充气器在实际工况下的充气性能数据。本文在实验室前期实验研究的基础上[7],利用仿真手段研究大喷嘴充气器(出气口径10 mm)在真实工作环境状况下的充气特性,即充气压力和浸没深度与充气性能之间的关系,为大喷嘴充气器的实际应用提供理论支持,进一步为其在实际工况应用下时相关参数的选取提供数据支撑。

1 大喷嘴充气器的工作原理

大喷嘴充气器的工作原理示意图如图1 所示。空压机送来的高压气体通过充气器的进气口进入水平管道,进而到达前方的喷嘴处,因为气体压力的作用压缩后面执行机构中的弹簧,弹簧带动运动阀杆向后运动,使得前端喷嘴打开,高压气体通过喷嘴喷出进入反应设备的浆液中,并通过浆液的剪切作用产生大量小气泡,小气泡由于初始速度的作用向喷嘴的前方及左右方向运动,并且在浮力的作用下同步向上运动。

图1 大喷嘴充气器工作原理示意图

表征大喷嘴充气器充气性能的参数主要有单位时间内喷射出的空气体积流量V、气流喷射长度L、气流末端扩散宽度B和喷射气泡直径D等[8],示意图如图2 所示。

图2 大喷嘴充气器充气性能表示示意图

对于大型冶金设备而言,当反应所需气量一定时,单个喷嘴喷射出的气体体积流量决定了充气器的配置数量,是大喷嘴充气性能最重要的表征参数;而气流喷射长度和末端扩散宽度则决定了充气系统中充气器的配置方式。

根据前期的实验研究[5]可知,充气压力和喷嘴结构参数对大喷嘴喷射气泡直径D影响最大,其它因素影响都较小。由于充气压力达到一定值后,喷射的气泡平均直径已经满足工业上的使用要求,同时考虑仿真计算中气泡直径不好评判获取,故本文中不对其进行对比研究。

2 大喷嘴充气器工作仿真模型的建立

由于大喷嘴充气器工作时将高压空气喷射入反应液中,涉及两种流体的变化状态,采用CFD 流体仿真软件对其进行模拟[9-10],反应液简化为清水。

2.1 几何模型

依据前期实验室实验结果[7],大喷嘴喷射气体长度在500 mm 左右,末端扩散宽度也为500 mm 左右,仿真模型采用水槽横截面为1 m ×1 m 的正方形,考虑到喷嘴喷出的气体向下扩散的区域很小,喷嘴中心距水槽底部300 mm,喷嘴前端伸入水槽壁面100 mm,水槽的高度依据大喷嘴真实工作的浸没深度设定,大喷嘴充气器工作时的三维几何模型如图3 所示。

图3 大喷嘴充气器工作三维几何模型示意图

2.2 求解计算模型

由于大喷嘴的充气试验所使用的气量较大,工作时要求具有较好的载气条件,所以采用欧拉-欧拉模型模拟气相与液相的流动情况。

2.3 气泡处理方法

在计算流体力学数值模拟的商业软件中,气泡尺寸分析方法分为单一气泡尺寸模型和多组气泡尺寸模型。单一气泡尺寸适用于气泡流动稳定,粒子与主相的速度大小方向相同,即Storkes 数值较小的情况。多尺寸分组模型(MUSG 模型)是简化的群体平衡模型。群体平衡模型是一种“刻画”多相流中分散相分布的方法,能够考虑到喷气过程中气泡的聚并与破碎现象,且气泡有尺寸范围的分布,对于提高水槽内气含率的准确性具有一定的意义。本文中所使用的MUSIG 方法的破碎模型和聚并模型选用的分别是Luo and Svendsen 模型和Prince and Blanch 模型。

2.4 求解策略与网格无关性验证

对于边界条件的设置,将充气的入口设置在喷嘴靠近槽壁处,顶部出口可以设置为主轴degassing出口边界或Opening 开放边界。为保证容易收敛,顶部出口设置Opening 开放边界,定义仅允许空气流进和流出,压力为0 Pa。除动静流体域的交界面外,其余边界选择为壁面边界。以稳态计算方法为主,求解设定目标残差小于10-4,或模拟至喷嘴入口流量和顶部出口流量守恒为止,数值模拟计算工作使用天津市天河计算机技术有限公司超算平台的商业软件完成计算。

仿真选用四面体的非结构网格,非结构网格尽管计算量较大,但具有适应性强,划分网格单元消耗时间短等显著优点。采用100 万、200 万、400 万网格单元数的模型进行网格无关性验证。如图4 所示,对比了喷嘴压充气力为0.4 MPa,浸没水深2 m的喷气量情况。可以看到相比100 万网格数,在网格200 万的时候喷嘴喷出的气流基本已经连续和网格400 万的情况已经较类似,不影响分析的精度。

图4 不同网格单元数的气含率分布

计算该条件下的喷嘴喷射出的气体流量,结果见表1,表1 数据表明200 万和400 万网格单元数的模型气体流量较为接近。

表1 不同网格单元数时喷嘴喷出的气体流量

综合计算时间和求解精度考虑,最终每2 米深度采用200 万网格的模型方法进行大喷嘴充气器在不同气体压力和不同浸没深度条件下的仿真计算,图5 给出了200 万网格量的模型图。

图5 大喷嘴充气器工作仿真模拟四面体网格模型

2.5 仿真变量参数

利用上面仿真模型的建立策略,本文重点对大喷嘴充气器在实际工作条件下的充气性能进行清水仿真研究,也即对大喷嘴充气器实际浸没深度和充气压力对充气性能的影响,这两个仿真变量参数如表2 所示。

表2 大喷嘴充气器实际工况下仿真变量参数

3 仿真结果及分析

3.1 出流量随浸没深度及充气压力的变化情况

通过后处理软件处理仿真结果,获取到10 mm喷嘴出口截面各点的气体流速,计算出平均气体流速乘以喷嘴截面积即得到不同情况下的喷嘴出流量,利用绘图软件绘制各工况下的气体出流量曲线图,如图6 所示。

图6 浸没深度和充气压力与喷嘴的出流量关系

从图6 中可以看出,在相同的充气压力条件下,随着浸没深度的增大,喷嘴出流量在逐渐减小;在相同的浸没深度条件下,随着充气压力的增大,喷嘴的出流量在增大,这两点变化趋势与实验工况下得到的结论是完全一样的。从表中数据可以得到充气压力为0.3 MPa、0.4 MPa 和0.5 MPa 时,随浸没深度的增大喷嘴出流量的减小幅度分别为39.7%、28.4%、21.6%,可知在喷嘴浸没深度达到实际生产量值时,出流量的减小值非常显著,特别是在充气压力值较低时此影响更加不可忽略。

3.2 喷射气流长度随浸没深度及充气压力的变化情况

通过后处理软件,显示不同浸没深度及充气压力下仿真结果中的含气率云图,截取充气器喷嘴中心竖直截面,可以看到喷嘴喷射气流长度的变化情况,如图7 所示。

图7 浸没深度和充气压力与喷嘴喷射气流长度关系

在图7 中,纯红色代表气相,纯蓝色代表液相,中间的过渡颜色代表含气率的变化,从图中可以直观的看出,在相同的充气压力条件下,随着浸没深度的增大,喷嘴的喷射气流长度在逐渐的减小;在相同的浸没深度条件下,随着充气压力的增大,喷嘴的喷射气流长度在逐渐增大。

通过提取各工况下喷嘴中心线上距喷嘴前端中心点远近不同距离处的含气率数据,绘制曲线图如图8 所示,从图中可以看出,不同浸没深度和充气压力时离开喷嘴的距离小于等于500 mm 的范围内,其水中的含气率都大于60%,可以认为喷嘴喷出气体的气流长度都大于500 mm。

图8 浸没深度和充气压力与离开喷嘴不同距离的含气率关系

3.3 喷射气流宽度随浸没深度及充气压力的变化情况

利用3.2 小节的方法,仍旧在显示不同浸没深度及充气压力下仿真结果中的含气率云图下,截取充气器喷嘴中心水平截面,可以很清楚的看到喷嘴喷射气流宽度的变化情况,如图9 所示。

在图9 中,纯蓝色到纯红色仍旧代表水相到气相的变化情况,纯蓝色表示全是液相,纯红色代表全是气相,中间的过渡颜色代表液相中含气率的大小变化,从图中可以直观的看出,在相同的充气压力条件下,随着浸没深度的增大,喷嘴的喷射气流宽度在逐渐的减小;在相同的浸没深度条件下,随着充气压力的增大,喷嘴的喷射气流宽度在逐渐增大。

通过提取各工况下喷嘴中心线上距喷嘴前端中心点左右不同距离处的含气率数据,绘制曲线图如图10 所示,从图中可以看出不同浸没深度和充气压力时离开喷嘴中心距离等于100 mm 的范围内,其水中的含气率基本都减小为0,可以认为喷嘴喷出气体的末端扩散宽度都小于200 mm。