分享：真空感应气体雾化技术中紧耦合环缝喷嘴 雾化过程的仿真

孙迎建1,2,3,周利杰1,2,3,宫 翔1,2,3,张鹏程1,2,3,郝瑞林1,2,3 (1.河北省工业机械手控制与可靠性技术创新中心,沧州 061001;2.沧州市工业机械手控制与可靠性技术 创新中心,沧州 061001;3.河北水利电力学院机械工程系,沧州 061001) 

摘 要:采用流体体积方法耦合雷诺应力模型与离散相模型结合泰勒类比不稳定性破碎模型, 利用计算流体力学软件 Fluent19.2,对紧耦合环缝喷嘴初次雾化与二次雾化进行全流程仿真,并进 行了试验验证。结果表明:初次雾化过程使导流管底部形成环形液膜结构,液膜前端的初次雾化主 要是气体射流自由边界湍流剪切的结果,并且初次雾化形成的液滴直径满足正态分布;随着分散的 液滴群外侧接触气体射流,从液滴群的外侧开始向心部发生二次雾化过程,但是未接触气体射流的 液滴仍保持较高的过热度。仿真得到紧耦合环缝喷嘴二次雾化后粉末的直径与试验结果吻合较 好,相对误差小于5%,验证了仿真的准确性。 

关键词:真空感应气体雾化;紧耦合环缝喷嘴;初次雾化;二次雾化 中图分类号:TF123 文献标志码:A 文章编号:1000-3738(2022)04-0075-07

0 引 言 

随着金属增材制造与粉末冶金等行业的发展, 对高质量金属粉末的需求日益增加[1]。目前主要的 金属粉末制备方法包括以德国 ALD 公司为代表的 75 孙迎建,等:真空感应气体雾化技术中紧耦合环缝喷嘴雾化过程的仿真 电极 感 应 熔 炼 气 体 雾 化 法 (electrodeinduction meltinggasatomization,EIGA),以加拿大为代表 的等离子雾化 (plasmaatomization,PA)工艺与等 离 子 旋 转 电 极 雾 化 (plasmarotating electrode process,PREP)工艺以及传统的真空感应气雾化 (vacuum induction gas atomization,VIGA)工 艺[2-3]。前3种工艺主要用于制备高活性金属粉末, 如钛及钛合金。对于传统的不锈钢、铝合金以及铜 合金等粉末材料主要采用 VIGA 工艺进行制备,因 此该工艺在国内外得到了广泛的应用。 气雾化喷嘴是 VIGA 设备的核心部件,决定了 制备金属粉末的直径分布、空心率、卫星粉率等关键 技术指标[4]。气雾化喷嘴主要分为紧耦合式喷嘴与 自由落体式喷嘴。紧耦合式喷嘴可大大缩短气体射 流与导流管出口金属液流的距离,显著提高雾化效 率,因此应用广泛。VIGA 工艺制粉是在一个密闭 的炉体内进行的,采用试验手段很难直接对其喷嘴 雾化进行研究。随着计算流体力学的发展,大量研 究人员开始采用流体仿真与试验验证相结合的方式 对喷嘴雾化进行研究。喷嘴雾化过程分为初次雾化 与二次雾化[5-6]。ZEOLI等[7-8]忽略金属液流的初 次雾化过程,采用欧拉-拉格朗日方法对紧耦合环缝 喷嘴的液 滴 二 次 雾 化 以 及 凝 固 过 程 进 行 了 研 究。 MOTAMAN 等[9-10]采用单相气流场仿真与纹影技 术相结合,对紧耦合式喷嘴的雾化气体流场进行了 研究。朱玲玲等[11]采用流体体积(volumeoffluid, VOF)多相流模型对紧耦合式喷嘴初次雾化过程进 行了仿真,研究了喷嘴雾化压力和液体表面张力对 金属液流初次破碎过程的影响。紧耦合式喷嘴气雾 化制粉是一个初次雾化与二次雾化连续进行的整体 过程,但是目前大部分报道都只是对其中某部分过 程进行仿真研究。 金属液流全流程雾化过程的仿真对气雾化制粉 机理的探究与喷嘴结构的优化设计都有着重要的作 用。作 者 采 用 VOF 方 法 耦 合 雷 诺 应 力 模 型 (Reynoldsstressmodel,RSM)与离散相(discrete phasemodel,DPM)模 型 结 合 泰 勒 类 比 (Taylor analogybreak-up,TAB)不稳定性破碎模型,利用 基于 有 限 体 积 方 法 的 商 业 计 算 流 体 力 学 软 件 Fluent19.2对紧耦合环缝喷嘴的 初 次 雾 化 与 二 次 雾化进行全流程仿真,并进行了试验验证,研究喷 嘴的雾化过 程,以 期 为 喷 嘴 结 构 优 化 设 计 提 供 理 论基础。

 1 试验方法 

VIGA 设备原理如图1所示,该设备主要由熔 炼室、熔炼坩埚、中间包、喷嘴、雾化室、粉末收集罐 等 关 键 部 分 组 成。 熔 炼 坩 埚 最 大 熔 炼 质 量 为 100kg,最高加热温度为1973K,中间包最高保温 温度为1573K。制粉原料选择316L不锈钢,雾化 气体选择氩气。VIGA 方法制备金属粉末过程中, 熔炼室内气体压力基本稳定为25kPa,底部雾化室 内气体压力基本为-10kPa;316L 不锈钢熔炼温 度为1700K,中间包保温温度为1550K;雾化氩 气的进气压力为3.5MPa,气体温度为300K。制备 的316L不锈钢粉末采用三点取样获得3组粉末试 样,并利用激光粒度仪对粉末粒径进行测试。 图1 VIGA设备原理示意 Fig 1 SchematicofVIGAequipmentprinciple 

2.1 初次雾化仿真模型 

紧耦合环缝喷嘴的初次雾化仿真采用 VOF多 相流模型耦合 RSM 湍流模型进行。在 VOF 模型 中主要通过体积分数追踪多相流之间的相界面,描 述气/液界面的体积分数方程[12]为 1 ρq ∂ ∂t (αqρq)+ ∂ xi ?? (αqρqui) ?? ???? ?? ?? ???? =0 (1) 式中:ρq 为q 相密度;αq 为q 相体积分数;t为时间; ui 为xi 方向速度;xi 为二维轴对称模型的x 坐标 数值。 紧耦合环缝喷嘴的雾化气体马赫数一般大于 0.3,因此需要考虑气体的压缩性。在气/液两相流 雾化过程中,由可压缩的 N-S 方程控制的连续方 程、动量方程与能量方程[13]分别为 ∂ρ ∂t + ∂ ∂xi (ρui)=0 (2) ∂ ∂t (ρui)+ ∂ ∂xj (ρuiuj)=- ∂p ∂xi + ∂τij ∂xj +ρg +Fs (3) 76 孙迎建,等:真空感应气体雾化技术中紧耦合环缝喷嘴雾化过程的仿真 ∂ ∂t (ρT)+ ∂ ∂xi (ρujT)= ∂ ∂xj K∂T c∂xj -ρujT +ST (4) 式中:uj 为xj 方向的速度;xj 为二维轴对称模型 的y 坐标数值;p 为压力;g 为重力加速度;τij 为雷 诺应力张量;Fs 为表面张力;T 为温度;ρg 为气体 密度;ρ为网格单元内的加权密度;K 为气体传热系 数;c为比热容;ST 为黏性耗散项。 

2.2 二次雾化仿真模型 

具有一定过热度的初次雾化液滴在气动力作用 下变形,在黏性力与表面张力作用下收缩发生二次 破碎过程。为了衡量液滴的二次雾化,提出了一个 重要的无量纲参数韦伯数(We),具体计算公式为 We = ρgu 2 reld σ (5) 式中:urel 为气体与液滴的相对速度;d 为初次雾化 液滴直径;σ为液滴表面张力系数。 韦伯数是扰动的空气动力与表面张力之比,韦 伯数越大,液滴越易破碎[14]。目前典型的液滴破碎 模型主要包括布袋式破碎、复合式破碎与脱落式破 碎[15],不同初始液滴韦伯数范围对应的破碎模型也 不同[16]。TAB模型适合描述布袋式破碎与复合式 破碎,开尔文-亥姆霍兹 KH 模型适合描述脱落式破 碎[17]。估算初次雾化液滴的韦伯数小于80,所以液 滴破碎模型为布袋式破碎、复合式破碎,因此采用 TAB模型仿真喷嘴的二次雾化过程。为了简化仿 真,将二次雾化仿真的初始液滴直径设置为初次雾 化液滴的质量中值直径。同时,利用离散相 DPM 模型,考虑液滴二次雾化过程中会出现的液滴碰撞、 液滴合并与液滴随机轨道等相关问题,选取二次雾 化仿真液滴初始条件如表1所示。 表1 二次雾化仿真液滴初始条件 Table1 Initialconditionsofsecondaryatomization simulationdroplet 参数 x 方向 位置/ mm y 方向 位置/ mm x 方向 速度/ (m·s -1) y 方向 速度/ (m·s -1) 液滴 直径/ μm 液滴 温度/ K 起始点 67 3.5 -0.5 -0.1 276 1700 结束点 67 4.0 -0.5 -0.1 276 1700 液滴在二次雾化过程中,冷却凝固转变为金属 粉末的过程是一个快速凝固的过程,液滴内部的温 度梯度与热辐射基本忽略[4]。冷却凝固过程可以采 用牛顿冷却模型进行描述[4],对应传热方程为 Cpd dTd dt =ΔHd dfs dt - 6h ρdD (Td -Tg) (6) Cpd =(1-fs)cl+csfs (7) ΔHd =ΔHf- (cl-cs)(Tl-Td) (8) h= λg D 2.0+0.6 Re 3 Pr (9) 式中:Cpd 为金属液滴的恒定热容;菪Hd 为液滴的潜 热;cl 为液态的比热容;cs 为固态的比热容;ΔHf 为 单位体积金属的熔化焓;fs 为液滴中固相的体积分 数;Td 为瞬时温度;D 为二次雾化液滴直径;λg 为 雾化气体的热导率;Re 为雷诺 数;ρd 为 液 滴 的 密 度;Pr 为普朗特数;Tl 和 Tg 分别为瞬时金属液相 温度和气相温度。

3 结果与讨论 

3.1 初次雾化过程 

由图4可知,紧耦合环缝喷嘴出口的气体射流 主要由具有较高速度的核心区与低速的自由边界组 成。喷嘴出口位置的气流射流由于与环境压力差的 作用,导致高速气流向下喷射过程中出现压缩波与 膨胀波交替的链状结构[5]。另外,气体射流向喷嘴 底部中间位置汇聚过程中,由于气流逐渐向中间偏 转,造成在导流管底部形成了一个气体回流区。回流 区内的气体流动轨迹主要是,中间气流向上进入气体 回流区,流动一定距离后转变流动方向,逐渐沿着径 向方向流动,最终接触高速气体射流,与其一同回到 起点位置。经过气体回流区的部分气体射流,将与主 射流发生融合,形成一致向下流动的射流结构。 金属液流从导流管底部流出,接触到气体回流 图4 紧耦合环缝喷嘴气体速度云图 Fig 4 Gasvelocityclouddistributionofclose-coupledannularnozzle 区后,其向下的继续流动将受到阻碍;此时金属液流 只能沿着不受阻力的径向方向继续流动,形成环形 液膜结构,实现连续液流至液膜的转变过程,如图5 所示。液膜流动过程中受到气体回流区的气流拖拽 挤压,前端的液膜厚度逐渐变薄。当径向流动的环 形液膜靠近雾化气体射流后,在气体射流拽力作用 下,再次改变流动方向,并且前端液膜将直接发生破 碎或厚度进一步减薄,最终完成初次雾化过程。 图5 紧耦合环缝喷嘴初次雾化原理示意 Fig 5 Schematicofprimaryatomizationofclose-coupledannularnozzle 由图6可以看出:紧耦合环缝喷嘴环形液膜前 端在气体射流拖拽作用下改变了流动方向,并直接 发生了液膜的破碎,形成了一定数量的初次雾化液 滴,液滴在气体射流的作用下,逐渐向喷嘴底部中心 位置汇聚;环形液膜前端仅接触了气体射流外侧的 自由边界,液膜前端在射流边界的湍流剪切力作用 下发生破碎,形成初次雾化液滴,同时在破碎液滴下 落 过 程 中,将 在 液 滴 周 围 形 成 漩 涡 绕 流 结 构。 üNAL等[18]采用高速摄像技术研究了 VIGA 方法 78 孙迎建,等:真空感应气体雾化技术中紧耦合环缝喷嘴雾化过程的仿真 中氩气雾化压力为2.4 MPa时,金属铝液的初次雾 化过程,发现铝液在导流管底部形成了一个环形液 膜结构,与仿真结果相似,基本可以验证初次雾化仿 真的可靠度。  基于德洛奈三角形方法,采用 MATLAB 软件 编写计算程序,对初次雾化形成的液滴进行等效面 积直径的计算,结果如图7所示。由图7可以看出, 紧耦合环缝喷嘴初次雾化后形成的液滴直径基本符 合正态分布,主要分布在250~285μm 之间,中值 直径为268μm。综上可知,紧耦合环缝喷嘴初次雾 化的关键是环形液膜的形成,前端液膜的破碎主要 是气体射流的自由边界湍流剪切造成的,并且初次 雾化形成的液滴直径满足正态分布。 

 4 结 论 

(1)紧耦合环缝喷嘴初次雾化过程将在导流管 底部形成环形液膜结构,实现连续液流至液膜的转 变过程是初次雾化的关键。液膜尖端的破碎主要是 气体射流自由边界湍流剪切的结果,并且初次雾化 液滴直径满足正态分布。 (2)紧耦合环缝喷嘴二次雾化过程分为汇聚、 初始分散与分散等3个阶段。随着液滴群外侧逐渐 靠近气体射流,液滴开始发生二次雾化过程,但是开 始未接触气体射流区域的心部液滴将保持较高的过 热度等待二次雾化。仿真得到紧耦合环缝喷嘴二次 雾化后粉末的直径与试验结果吻合较好,相对误差 小于5%,验证了仿真的准确性。 

来源：材料与测试网