混合澄清槽内离子液体体系的三相流体动力学模拟

doi:10.16085/j.issn.1000-6613.2020-2216

摘要/Abstract

摘要：

离子液体是近年来被广泛研究的新兴绿色溶剂，其在乏燃料后处理技术中具有潜在的工业应用前景。但由于缺乏对萃取反应器中离子液体流动特性的研究，制约了离子液体萃取体系的实际应用。本文以萃取工艺中广泛应用的混合澄清槽为对象，以去离子水及1-丁基-3-甲基咪唑双(三氟甲磺酰)亚胺盐（[C₄mim][NTf₂]）分别作为水相及有机相，考虑上层空气对流动行为的影响，对其混合室进行计算流体力学（computational fluid dynamics, CFD）模拟，考察不同转速、流比及温度下的有机相分布、压力场、湍动程度等。结果表明，模拟结果较好地符合实验结果，且最大误差小于6.3%；转速能直观地提升混合性能，但当超过500r/min后，继续提高转速将显著增大出口气量，从而可能对澄清室的分相性能提出更高要求；增大流比、升温均能提升350r/min转速下的有机、水相混合能力，升温还有效减小了桨力矩，但当温度超过303K时，继续升温对于桨力矩、有机相速度的改变不明显。因此，实际工艺条件建议结合升温与转速调节，在实现较好混合性能的同时，减少对澄清室分相性能的要求。本文在建立离子液体三相体系数值模拟方法的同时，为混合澄清槽的工况优化提供合理建议，并为离子液体萃取体系的深入研究提供了参考。

关键词: 离子液体, 计算流体力学, 混合澄清槽, 搅拌容器, 多相流

Abstract:

As a new kind of green solvents, ionic liquids (ILs) have been widely studied in recent years and have potential application prospects in the spent fuel reprocessing technology. However, a lack of study on flow characteristics of ILs in the extraction equipment limits the practical applications of ILs extraction system. Mixer-settler widely used in extraction was taken as the research object in this paper. The deionized water and the 1-butyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ([C₄mim][NTf₂]) were adopted as the aqueous phase and the organic phase, respectively. The computational fluid dynamics (CFD) simulation was carried out to study the organic phase distribution, pressure field and turbulence degree, etc. at different agitation speeds, flow ratios and temperatures. Moreover, the influence of air on the flow behavior was considered. The simulation result accorded with the experiment data and the maximum error was less than 6.3%. The results showed that the mixing performance can be promoted significantly by increasing the agitation speed. The air flow rate at the outlet was increased significantly with increasing agitation speed when it exceeded 500r/min, which may affected the performance of clarifying chamber. Besides, the mixing performance of organic phase and aqueous phase can be promoted with increasing flow ratio and temperature at the agitation speed of 350r/min. The moment of agitator was effectively reduced by heating up. The organic phase velocity and the moment of agitator were not significantly affected by heating when the temperature exceeds 303K. Therefore, the actual process conditions were recommended to combine with heating and the adjustment of agitation speed for achieving better mixing performance while reducing the requirements of the clarifying chamber.A numerical simulation method for the three-phase system of ILs is set up in the study. Moreover, some reasonable suggestions are provided for optimizing the operating conditions of the mixer-settler, and a reference for further study of ILs extraction system is provided.

Key words: ionic liquids (ILs), computational fluid dynamics (CFD), mixer-settler, stirred vessel, multiphase flow

中图分类号:

TB126

汤祺, 鲍迪, 邵少雄, 徐平, 刘联伟, 郑维明. 混合澄清槽内离子液体体系的三相流体动力学模拟[J]. 化工进展, 2021, 40(10): 5468-5479.

TANG Qi, BAO Di, SHAO Shaoxiong, XU Ping, LIU Lianwei, ZHENG Weiming. Three-phase fluid dynamics simulation of ionic liquid system in mixer-settler[J]. Chemical Industry and Engineering Progress, 2021, 40(10): 5468-5479.

图/表 22

图1 研究思路

图2 实验装置

表1 离子液体、水及空气物性参数

（T=293K，p=101.325kPa）

离子液体［C₄mim］［NTf₂］

（纯度99%）

表2 欧拉-欧拉模型控制方程

方程名称	方程式	序号
连续性方程	$? ? t α i ρ i + ? ? a i ρ i U i = 0$	（1）
动量方程	$? ? t a i ρ i U i + ? ? a i ρ i U i U i = - a i ? p + a i ρ i g + ? ? τ i ˉ ˉ + F i j$	（2）
RNG k-ε湍流模型	$? ? t ρ m k + ? ? x i ρ m u m k = ? ? x j α k i μ e f f, m ? k ? x j + G k - ρ m ε$	（3）
	$? ? t ρ m ε + ? ? x i ρ m u m ε = ? ? x j α ε μ e f f, m ? ε ? x j + ε k C 1 ε G k - ε 2 k C 2 ε$	（4）
曳力系数模型		（4）
曳力系数模型	$C D = 24 / R e 1 + 0.15 R e 0.687 ????? R e ≤ 1000 0.44 ???????????????????????????????????????? R e > 1000$	（5）
液滴的Sauter平均直径模型^［30-31］	$d 32 / D = 0.28 1 + 0.92 α j W e i j - 0.6$	（6）
液滴的Sauter平均直径模型^［30-31］	$W e i j = ρ i D 3 N / 60 2 / σ i j$	（7）
N-phase 体积分数方程	$α i + α j + α k = 1$	（8）
水力直径	圆形：D_H=D_i	（9）
水力直径	长方形：D_H=2HW/（H+W）	（10）
湍流强度	$I = 0.16 R e - 1 / 8$	（11）

表2 欧拉-欧拉模型控制方程

方程名称	方程式	序号
连续性方程	$? ? t α i ρ i + ? ? a i ρ i U i = 0$	（1）
动量方程	$? ? t a i ρ i U i + ? ? a i ρ i U i U i = - a i ? p + a i ρ i g + ? ? τ i ˉ ˉ + F i j$	（2）
RNG k-ε湍流模型	$? ? t ρ m k + ? ? x i ρ m u m k = ? ? x j α k i μ e f f, m ? k ? x j + G k - ρ m ε$	（3）
	$? ? t ρ m ε + ? ? x i ρ m u m ε = ? ? x j α ε μ e f f, m ? ε ? x j + ε k C 1 ε G k - ε 2 k C 2 ε$	（4）
曳力系数模型		（4）
曳力系数模型	$C D = 24 / R e 1 + 0.15 R e 0.687 ????? R e ≤ 1000 0.44 ???????????????????????????????????????? R e > 1000$	（5）
液滴的Sauter平均直径模型^［30-31］	$d 32 / D = 0.28 1 + 0.92 α j W e i j - 0.6$	（6）
液滴的Sauter平均直径模型^［30-31］	$W e i j = ρ i D 3 N / 60 2 / σ i j$	（7）
N-phase 体积分数方程	$α i + α j + α k = 1$	（8）
水力直径	圆形：D_H=D_i	（9）
水力直径	长方形：D_H=2HW/（H+W）	（10）
湍流强度	$I = 0.16 R e - 1 / 8$	（11）

图3 流体域的边界条件设置（括号内为实验验证边界条件）

图4 监测线1（x=0，-23≤y≤23，z=5，单位mm）、2（-23≤x≤23，y=0，z=12，单位mm）上物理量与网格数量的关系

图5 最终计算用网格示意图

图6 350r/min下的实验、模拟有机相分布高度对比

图7 y=0截面上不同转速下的有机相、气相体积分数云图

图8 不同转速下搅拌桨中心区域（z=22mm截面）的压力云图

图9 出口处气体流量同转速之间的关系

图10 监测线上（x=0，-23≤y≤23，z=5，单位mm）的湍动能、桨力矩与转速间的关系

表3 4种流比下的有机相及水相流量

水相入口流量/mL·min^-1	有机相入口流量/mL·min^-1	流比（a∶o）
25	25	1∶1
25	20	5∶4
25	15	5∶3
25	10	5∶2

图11 y=0截面上不同流比下的有机相、气相体积分数云图

表4 不同流比下混合室内的水相、有机相相比（体积比）

流比	水相体积/m³	有机相体积/m³	两相相比
1∶1	5.52×10^-5	4.47×10^-5	1.23
5∶4	5.57×10^-5	4.42×10^-5	1.26
5∶3	5.60×10^-5	4.39×10^-5	1.28
5∶2	5.72×10^-5	4.32×10^-5	1.33

表5 不同温度下的离子液体物性参数

温度/K	离子液体密度/kg·m^-3	黏度/mPa·s
293	1440.288	62.674
298	1435.433	50.069
303	1430.595	41.178
308	1425.688	34.779

图12 y=0截面上不同温度下的有机相、气相体积分数云图

表6 不同温度下混合室内水相、有机相相比（体积比）

温度/K	水相体积/m³	有机相体积/m³	两相相比
293	5.52×10^-5	4.47×10^-5	1.23
298	5.53×10^-5	4.47×10^-5	1.24
303	5.56×10^-5	4.45×10^-5	1.25
308	5.58×10^-5	4.44×10^-5	1.26

图13 监测线（x=0，-23≤y≤23，z=5，单位mm）上有机相速度与温度间的关系

图14 不同温度下截面（y=0）上的湍动能云图

图15 桨力矩随温度的变化情况

图16 350r/min时不同温度及293K-400r/min、293K-450r/min条件下监测线（x=0，-23≤y≤23，z=5，单位mm）上的湍动能

参考文献 31

1	ZHANG X, ZHANG X, DONG H, et al. Carbon capture with ionic liquids: overview and progress[J]. Energy & Environmental Science, 2012, 5(5): 6668.
2	DAI S, JU Y H, BARNES C E. Solvent extraction of strontium nitrate by a crown ether using room-temperature ionic liquids[J]. Journal of the Chemical Society Dalton Transactions, 1999(8): 1201-1202.
3	VISSER A E, SWATLOSKI R P, REICHERT W M, et al. Traditional extractants in nontraditional solvents: groups 1 and 2 extraction by crown ethers in room-temperature ionic liquids[J]. Industrial & Engineering Chemistry Research, 2000, 39(10): 3596-3604.
4	DIETZ M L, DZIELAWA J A. Ion-exchange as a mode of cation transfer into room-temperature ionic liquids containing crown ethers: implications for the ‘greenness’ of ionic liquids as diluents in liquid-liquid extraction[J]. Chemical Communications, 2001, 20: 2124-2125.
5	SHI C, JING Y, JIA Y. Tri-n-butyl phosphate-ionic liquid mixtures for Li⁺ extraction from Mg²⁺-containing brines at 303—343K[J]. Russian Journal of Physical Chemistry A, 2017, 91(4): 692-696.
6	DIETZ M, DZIELAWA J, LASZAK I, et al. Influence of solvent structural variations on the mechanism of facilitated ion transfer into room-temperature ionic liquids[J]. Green Chemistry, 2003, 5(6): 682-685.
7	GIRIDHAR P, VENKATESAN K A, SRINIVASAN T G, et al. Extraction of uranium(Ⅵ) from nitric acid medium by 1.1M tri-n-butylphosphate in ionic liquid diluent[J]. Journal of Radioanalytical and Nuclear Chemistry, 2005, 265(1): 31-38.
8	DIETZ M L, STEPINSKI D C. Anion concentration-dependent partitioning mechanism in the extraction of uranium into room-temperature ionic liquids[J]. Talanta, 2008, 75(2): 598-603.
9	路文娟. 微乳液及离子液体萃取金和汞的研究[D]. 济南: 山东大学, 2013.
	LU Wenjuan. Extraction of gold() and mercury(Ⅱ) by microemulition and ionic liquids[D]. Jinan: Shandong University, 2013.
10	敖银勇, 周瀚洋, 彭静, 等. 咪唑离子液体体系辐解产物的鉴定及其对萃取的影响[C]//中国核学会2015年学术年会论文集.绵阳, 2015: 16-22.
	AO Yinyong, ZHOU Hanyang, PENG Jing, et al.Identification of radiolytic products of RTILs-based extraction system and their effect on the removal of metal ions[C]//Progress Report on China Nuclear Science & Technology. Mianyang, 2015: 16-22.
11	束玉珍, 吴继宗, 邓惟勤, 等. 苯并-15-冠-5/离子液体体系液液萃取锂[J]. 核化学与放射化学, 2017, 39(2): 151-157.
	SHU Yuzhen, WU Jizong, DENG Weiqin, et al. Liquid-liquid extraction of lithium by benzo-15-crown-5/ionic liquid system[J]. Journal of Nuclear and Radiochemistry, 2017, 39(2): 151-157.
12	LIU Z, HUANG S, WANG W. A refined force field for molecular simulation of imidazolium-based ionic liquids[J]. The Journal of Physical Chemistry B, 2004, 108(34): 12978-12989.
13	HUANG Y, ZHAO Y, ZENG S, et al. Density prediction of mixtures of ionic liquids and molecular solvents using two new generalized models[J]. Industrial & Engineering Chemistry Research, 2014, 53(39): 15270-15277.
14	HUANG Y, ZHANG X, ZHAO Y, et al. New models for predicting thermophysical properties of ionic liquid mixtures[J]. Physical Chemistry Chemical Physics, 2015, 17(40): 26918-26929.
15	DONG H, WANG X, LIU L, et al. The rise and deformation of a single bubble in ionic liquids[J]. Chemical Engineering Science, 2010, 65(10): 3240-3248.
16	鲍迪, 张香平, 张欣, 等. 非常规介质离子液体中气泡行为研究进展[J]. 工程研究——跨学科视野中的工程, 2015, 7(3): 305-312.
	BAO Di, ZHANG Xiangping, ZHANG Xin, et al. Research progress of bubble behavior in ionic liquid systems[J]. Journal of Engineering Studies, 2015, 7(3): 305-312.
17	汤祺, 白璐, 董海峰, 等. 离子液体体系流体动力学研究现状及发展趋势[J]. 化工进展, 2018, 37(4): 1323-1334.
	TANG Qi, BAI Lu, DONG Haifeng, et al. Research on status and developing trends of hydrodynamics in ionic-liquid system[J]. Chemical Industry and Engineering Progress, 2018, 37(4): 1323-1334.
18	马婷婷. 混合澄清槽内流动的测量与模拟[D]. 天津: 天津大学, 2014.
	MA Tingting. On measurement and simulation of flow in mixer-settlers[D]. Tianjin: Tianjin University, 2014.
19	黄毅, 武斌, 陈葵, 等. 新型大三角桨叶混合槽内流场的数值模拟[J]. 湿法冶金, 2016, 35(4): 361-364.
	HUANG Yi, WU Bin, CHEN Kui, et al. Numerical simulation of flow field of a new triangle impeller in mixer-settler[J]. Hydrometallurgy of China, 2016, 35(4): 361-364.
20	王亮, 聂林. 稀土萃取搅拌槽内两相混合过程的数值计算[J]. 中国有色冶金, 2017, 46(2): 42-44, 76.
	WANG Liang, NIE Lin. Numerical calculation of two-phase mixing process in rare earth extraction agitation tank[J].China Nonferrous Metallurgy, 2017, 46(2): 42-44, 76.
21	唐巧, 叶思施, 王运东. 放大准则对混合澄清槽混合室中混合时间和流动特性的影响[J]. 化工学报, 2016, 67(2): 448-457.
	TANG Qiao, YE Sishi, WANG Yundong. Mixing time and flow characteristic in square pump-mix mixer under different scale-up criteria[J]. CIESC Journal, 2016, 67(2): 448-457.
22	吴富姬, 周雄军. 稀土萃取混合澄清槽的放大设计研究[J]. 湿法冶金, 2014, 33(3): 232-235.
	WU Fuji, ZHOU Xiongjun. Research on enlarge design of mixing settler[J]. Hydrometallurgy of China, 2014, 33(3): 232-235.
23	GU D, LIU Z, XU C, et al. PIV measurement and CFD simulation of liquid-liquid mixing in mixer settler with rigid-flexible impeller[J]. International Journal of Chemical Reactor Engineering, 2019, 17(11): 546-557.
24	郑雄攀. 稀土萃取混合澄清槽内搅拌过程强化[D]. 重庆: 重庆大学, 2015.
	ZHENG Xiongpan. Mixing process intensification in a mixer settler for rare earth extraction[D]. Chongqing: Chongqing University, 2015.
25	许传林. 混合澄清槽内液液两相混沌混合强化规律研究[D]. 重庆: 重庆大学, 2018.
	XU Chuanlin. Study on chaotic mixing intensification of liquid-liquid in mixer settler[D]. Chongqing: Chongqing University, 2018.
26	倪志南, 武斌, 陈葵, 等. 液-液萃取过程中液滴分散特性的数值模拟[J]. 湿法冶金, 2018, 37(5): 402-406.
	NI Zhinan, WU Bin, CHEN Kui, et al. Numerical simulation of droplet dispersion in liquid-liquid mixing chamber[J]. Hydrometallurgy of China, 2018, 37(5): 402-406.
27	TANG Q, YE S, WANG Y, et al. A study on liquid-liquid dispersions in a continuous mixer via computational fluid dynamics (CFD) simulation combined with population balance model (PBM)[J]. The Canadian Journal of Chemical Engineering, 2019, 97(2): 452-464.
28	LUO H, SVENDSEN H F. Theoretical model for drop and bubble breakup in turbulent dispersions[J]. AIChE Journal, 1996, 42(5): 1225-1233.
29	刘秀. 箱式混合澄清萃取槽流场数值模拟研究[D]. 赣州：江西理工大学, 2013.
	LIU Xu. Numerical simulation of flow field of box mixer-settlers in extration tank[D]. Ganzhou: Jiangxi University of Science and Technology, 2013.
30	DESNOYER C, MASBERNAT O, GOURDON C. Experimental study of drop size distributions at high phase ratio in liquid-liquid dispersions[J]. Chemical Engineering Science, 2003, 58(7): 1353-1363.
31	GARDAS R L, GE R, MANAN N AB, et al. Interfacial tensions of imidazolium-based ionic liquids with water and n-alkanes[J]. Fluid Phase Equilibria, 2010, 294(1/2): 139-147.

[1]	王云飞, 秦蕊, 郑利军, 李焱, 李清平. 旋转填充床CFD模拟研究进展[J]. 化工进展, 2023, 42(S1): 1-9.
[2]	黄益平, 李婷, 郑龙云, 戚傲, 陈政霖, 史天昊, 张新宇, 郭凯, 胡猛, 倪泽雨, 刘辉, 夏苗, 主凯, 刘春江. 三级环流反应器中气液流动与传质规律[J]. 化工进展, 2023, 42(S1): 175-188.
[3]	孙继鹏, 韩靖, 唐杨超, 闫汉博, 张杰瑶, 肖苹, 吴峰. 硫黄湿法成型过程数值模拟与操作参数优化[J]. 化工进展, 2023, 42(S1): 189-196.
[4]	罗成, 范晓勇, 朱永红, 田丰, 崔楼伟, 杜崇鹏, 王飞利, 李冬, 郑化安. 中低温煤焦油加氢反应器不同分配器中液体分布的CFD模拟[J]. 化工进展, 2023, 42(9): 4538-4549.
[5]	吴亚, 赵丹, 方荣苗, 李婧瑶, 常娜娜, 杜春保, 王文珍, 史俊. 用于复杂原油乳液的高效破乳剂开发及应用研究进展[J]. 化工进展, 2023, 42(8): 4398-4413.
[6]	卜治丞, 焦波, 林海花, 孙洪源. 脉动热管计算流体力学模型与研究进展[J]. 化工进展, 2023, 42(8): 4167-4181.
[7]	杨许召, 李庆, 袁康康, 张盈盈, 韩敬莉, 吴诗德. 含Gemini离子液体低共熔溶剂热力学性质[J]. 化工进展, 2023, 42(6): 3123-3129.
[8]	吕超, 张习文, 金理健, 杨林军. 新型两相吸收剂-离子液体系统高效捕获CO₂[J]. 化工进展, 2023, 42(6): 3226-3232.
[9]	田启凯, 郑海萍, 张少斌, 张静, 余子夷. 混合增强的微流控通道进展[J]. 化工进展, 2023, 42(4): 1677-1687.
[10]	张成松, 张静, 龚斌, 李明洋, 袁佳新, 李宏业. 自吸射流柔性搅拌桨振动特性[J]. 化工进展, 2023, 42(4): 1728-1738.
[11]	邱沫凡, 蒋琳, 刘荣正, 刘兵, 唐亚平, 刘马林. 气固流化床化学反应数值模拟中颗粒尺度模型研究进展[J]. 化工进展, 2023, 42(10): 5047-5058.
[12]	乔元, 仇畅, 钱锦远, 干瑞彬, 徐春明, 金志江. 笼式调节阀的冲蚀磨损与空蚀[J]. 化工进展, 2023, 42(10): 5111-5120.
[13]	华渠成, 段庆华. 离子液体极压抗磨剂的研究进展[J]. 化工进展, 2022, 41(S1): 331-339.
[14]	陈钰, 刘冲, 邱于荟, 贲梓欣, 牟天成. 离子液体和低共熔溶剂绿色回收废旧锂离子电池的研究进展[J]. 化工进展, 2022, 41(S1): 485-496.
[15]	彭德其, 冯源, 王依然, 谭卓伟, 俞天兰, 吴淑英. 立式缩放管内液固颗粒群浓度分布与汇聚特性[J]. 化工进展, 2022, 41(9): 4662-4672.