黄 聪, 吴可君, 何潮洪,2

不同桨型下含盘管的连续搅拌釜停留时间分布的研究

黄 聪1, 吴可君1, 何潮洪1,2

(1. 浙江省化工高效制造技术重点实验室, 浙江大学 化学工程与生物工程学院, 浙江 杭州 310027;2. 浙江大学衢州研究院, 浙江 衢州 324000)

用于含能化合物制备的连续釜式反应器通常安装有螺旋盘管，釜内流动情况复杂。研究搭建了一套含螺旋盘管的连续搅拌釜装置，对四斜叶桨(PBT)、推进桨(PRO)和翼型桨(CBY)作用下连续搅拌釜的停留时间分布进行了测定，并考察了桨型、进料流量、搅拌转速、表观气速等因素对无因次方差的影响，同时考察了对有气体产生的体系在含盘管的连续搅拌釜中多级混合模型的适用性。研究结果表明：多级混合模型基本适用于对含盘管连续搅拌釜式反应器的非理想流动的描述，可用于含盘管连续搅拌釜式反应器的优化设计。当其他条件相同，随着搅拌转速、表观气速的增大，含盘管的连续搅拌釜的无因次方差逐渐增大。对于强放热反应，推荐使用四斜叶桨或推进桨，可以较好地增强釜内物料的混合。

连续搅拌釜；
螺旋盘管；
停留时间分布；
含能化合物

含能化合物是由单一分子结构构成的含能材料，对于发展先进武器装备有重要作用。含能化合物的制备通常涉及硝化反应、叠氮化反应、缩合反应等强放热反应[1]，通常放出大量热量，而且反应体系含有大量不稳定物质，当反应条件变化引起反应失控时，其中的不稳定物质会分解产生气体。针对上述情况，仅靠普通的夹套无法满足强放热反应的换热需求。为解决强放热反应过程中单一夹套换热能力不足的问题，可在釜内安装螺旋盘管以增大传热面积[2-3]，使反应釜适用于含能化合物的生产。因此有必要测定含盘管反应釜在一定条件下的停留时间分布(RTD)，进而了解釜内流体的实际流动状况，确定其流动模型并计算模型参数[4]，用于反应釜的优化设计。

目前针对连续搅拌釜的停留时间分布研究中，以不含盘管的情况较多。Acquaye[5]通过实验确定了不同的搅拌桨类型(轴向或径向)对混合和平均停留时间的影响。董红星等[6]用内设有单层和双层的直六叶涡轮桨的釜开展了RTD实验，分析了液体流量和转速对流体流动状态的影响。张斌等[7]用装有不同类型搅拌桨的卧式搅拌床反应器，测定了装置的RTD，研究结果表明：安装T型桨的釜返混最小，同平推流最接近；
安装叶片桨的釜返混最大，同全混流最接近。Zhang等[8]在安装有圆盘涡轮桨的搅拌釜进行了RTD实验和计算流体力学(CFD)模拟，研究了反应器几何形状和操作条件等因素对停留时间分布的影响。周国忠等[9]研究了不同的内部组件尺寸和操作条件对搅拌槽内停留时间分布的影响。皇甫功凯等[10]考察了通气条件下不同液位时搅拌釜内物料流动的多釜串联模型参数和液相停留时间分布。李攀[11]和姚占静等[12]研究了不同进料流量大小和不同搅拌转速对连续搅拌釜停留时间分布的影响。Saravanathamizhan等[13]进行了搅拌釜的RTD有关实验，考察了连续搅拌釜式电化学反应器中电解液的流动特性。

关于含盘管反应器的停留时间分布研究，李伟[14]对内设三层45°斜叶桨的搅拌槽开展了RTD实验，考察了桨径、底层桨叶离底高度和单位体积功率输入等因素对搅拌槽内的RTD和流动模型参数的影响。实验结果显示，上述结构和操作条件的改变不影响槽内流体的宏观流动；
该装置的多釜串联模型参数在1.3左右，说明该装置的流动状况与理想的全混流反应釜较为接近，有利于进行相关反应的化工过程开发。

化学反应过程的开发，通常需要将反应过程与传递过程解耦，分别进行热模实验和冷模实验[15]。其中，反应过程的热模实验着重关注过程的反应规律，而包括RTD实验在内的冷模实验则侧重考察过程的传递规律。对于强放热反应而言，所用的反应釜常包含内盘管作为冷却部件。目前关于含盘管的连续釜停留时间分布的研究主要针对单一桨型和不含气体的均相流体，对于推进桨、翼型桨等其他桨型以及通气条件下的停留时间分布研究较少。而含能化合物制备过程中往往有气体产生，故需要了解通气条件下的停留时间分布。基于上述考虑，本研究参考张本贺等[16]、周俊超等[17]关于强放热反应器的设计方法，搭建了一套含盘管的连续搅拌釜装置，釜内主要部件的结构参数基本上按照强放热反应釜的推荐范围进行设计。通过该装置研究桨型、表观气速(表观气速=通气流量/釜内筒体横截面积)、进料流量和搅拌转速等因素对停留时间分布的影响，以期为强放热反应过程连续搅拌釜式反应器的优化设计提供参考。

2.1 实验装置

本实验装置主要包括搅拌釜、示踪剂高位槽、水槽、通气设备和数字采集系统。搅拌釜内设可拆卸的单层搅拌桨，可根据实验需要采用推进式桨、翼型桨和四斜叶桨等形式。实验装置如图1所示。连续搅拌釜、螺旋盘管和气体分布器的具体参数见表1。

图1 含盘管的连续搅拌釜停留时间分布的实验装置示意图

1. water storage tank 2. pump 3. tracer tank 4. stirred tank 5. digital acquisition system 6. waste liquid tank 7. coil 8. gas distributor 9. impeller 10. nitrogen cylinder

表1 连续搅拌釜的结构参数

2.2 实验过程与数据处理

实验采用脉冲示踪法测定液体的停留时间分布，以水为工作介质(当釜内液面静止时水位刚好覆盖盘管，搅拌釜筒体内液面深度约为0.260 m)，通过向釜内通气(表观气速为0～1.70 mm×s-1)以模拟连续搅拌釜内气体的产生情况。本实验以Na2SO4饱和溶液为示踪剂，其电导率在连续搅拌釜出口处测定。数据采集时间间隔为5 s。每组实验条件测定3次，以加强数据的可靠性。

写成离散形式：

示踪剂加入方法如下：开启实验装置，当搅拌釜的流量稳定后，观察釜出口处液体的电导率读数，当电导率稳定后，打开示踪剂高位槽的阀门，向釜内脉冲注入示踪剂。

在本实验装置中，示踪剂由高位槽经管道进入搅拌釜，脉冲示踪剂的加入量主要与脉冲加入时间有关。先通过预实验考察了脉冲示踪剂的加入对测试结果的影响，以确定较为适宜的脉冲示踪剂的加入量和加入时间。当示踪剂的加入时长分别取1、3、5 s时，有以下结果：当其他条件不变，示踪剂加入时间越长，则示踪剂的加入量越多，RTD实验的m曲线峰值越高；
但示踪剂的加入对()曲线的平均停留时间和无因次方差的影响可以忽略。为减少示踪剂对测试结果的影响，每次实验向釜内脉冲加入示踪剂溶液时长为3 s，并将m曲线进行无因次化处理得到()曲线，以便统一比较不同操作条件和桨型对搅拌釜内流体流动的影响。

以其中一组实验数据(四斜叶桨、表观气速为0.34 mm×s-1、进料流量为2.75 L×min-1、搅拌转速300 r×min-1)为例，说明电导率-时间曲线图的绘制过程。通过软件导出时间和电导率m的实验数据，用Excel编辑处理，绘制电导率随时间变化的曲线图，如图2所示。

图2 gm-t曲线图

3.1 通气条件下含盘管的连续搅拌釜的流动模型

多级混合模型是一种常见的描述釜式反应器非理想流动的数学模型，其基本假设为：每个级内为全混流、级间无返混、各级存料量相同，其模型参数表示原装置的返混程度，相当于个全混釜的返混程度，相关公式[19]如下：

不通气情况下的RTD实验研究表明：多级混合模型可用于含盘管的连续搅拌釜[14]。

为考察含盘管的连续搅拌釜，多级混合模型在有气体产生的情况下是否适用。图3给出了几个典型条件下停留时间分布的()实验曲线和多级混合模型拟合的曲线。

图3 在不同操作条件下含盘管的连续搅拌釜的RTD实验曲线和拟合曲线比较结果

图3表明：在实验范围内，含盘管连续搅拌釜的()实验曲线和用多级混合模型拟合的曲线较为接近。因此对于含盘管的连续搅拌釜，在有气体产生的情况下，工程上可采用多级混合模型反映釜内流体的返混情况。

3.2 操作条件对停留时间分布的影响

3.2.1 推进桨作用下不同操作条件对停留时间分布无因次方差的影响

考察在推进桨作用下，表观气速、搅拌转速和进料流量等操作条件对停留时间分布无因次方差和多级混合模型的模型参数的影响，结果如图4所示。其中图4(a)表示当进料流量为2.75 L×min-1，内设推进桨的含盘管搅拌釜在不同搅拌转速和不同表观气速下的无因次方差和模型参数(图5和6同理)。

由图4可见，在推进桨作用下，在实验条件范围内(进料流量为2.25～2.75 L×min-1、搅拌转速为120～360 r×min-1、表观气速为0～1.70 mm×s-1)，当其他条件不变，随着搅拌转速的增大，含盘管的连续搅拌釜的无因次方差逐渐增大，其最大值可达0.763，对应的多级混合模型的模型参数为1.311。以360 r×min-1情况为例，当表观气速从0增大为1.70 mm×s-1，进料流量分别为2.75、2.50、2.25 L×min-1，对应的无因次方差的增幅分别为15.94%、14.60%、12.91%，不同流量下的平均增幅为14.48%。

图4 不同进料流量下，搅拌转速和表观气速对内设推进桨的含盘管的连续搅拌釜无因次方差和多级混合模型参数的影响

3.2.2 四斜叶桨作用下不同操作条件对停留时间分布无因次方差的影响

考察在四斜叶桨作用下，表观气速、搅拌转速和进料流量等操作条件对停留时间分布无因次方差和多级混合模型的模型参数的影响，结果如图5所示。

由图5可见，在四斜叶桨作用下，在实验条件范围内(进料流量为2.25～2.75 L×min-1、搅拌转速为120～360 r×min-1、表观气速为0～1.70 mm×s-1)，当其他条件不变，随着搅拌转速的增大，含盘管的连续搅拌釜的无因次方差逐渐增大，其最大值可达0.759(对应的多级混合模型的模型参数为1.318)。以转速为360 r×min-1的情况为例，当表观气速从0增大为1.70 mm×s-1，进料流量分别为2.75、2.50、2.25 L×min-1，对应的无因次方差的增幅分别为11.92%、11.42%、11.80%，不同流量下的平均增幅为11.71%。

图5 不同进料流量下，搅拌转速和表观气速对内设四斜叶桨的含盘管的连续搅拌釜无因次方差和多级混合模型参数的影响

3.2.3 翼型桨作用下不同操作条件对停留时间分布无因次方差的影响

考察在翼型桨作用下，表观气速、搅拌转速和进料流量等操作条件对停留时间分布无因次方差和多级混合模型的模型参数的影响，结果如图6所示。

图6 不同进料流量下，搅拌转速和表观气速对内设翼型桨的含盘管的连续搅拌釜无因次方差和多级混合模型参数的影响

由图6可见，在翼型桨作用下，在实验条件范围内(进料流量为2.25～2.75 L×min-1、搅拌转速为120～360 r×min-1、表观气速为0～1.70 mm×s-1)，当其他条件不变，随着搅拌转速的增大，含盘管的连续搅拌釜的无因次方差逐渐增大，其最大值可达0.751(对应的多级混合模型的模型参数为1.332)。以转速为360 r×min-1的情况为例，当表观气速从0增大为1.70 mm×s-1，进料流量分别为2.75、2.50、2.25 L×min-1，对应的无因次方差的增幅分别为12.58%、12.55%、12.43%，不同流量下的平均增幅为12.52%。

比较图4～6，当其他条件不变，相比翼型桨，通常四斜叶桨和推进桨的无因次方差较大，对促进物料混合的作用较大。因此对于强放热反应，推荐使用四斜叶桨或推进桨，可较好地增强釜内物料的混合。

本研究搭建了一套含盘管的连续搅拌釜装置，考察了桨型和操作条件(表观气速、进料流量、搅拌转速)等因素对停留时间分布的影响，在实验范围内(桨型为推进桨/四斜叶桨/翼型桨、进料流量为2.25～2.75 L×min-1、表观气速为0.00～1.70 mm×s-1、搅拌转速为120～360 r×min-1)，结论如下：

(a) 通气条件下，含盘管的连续搅拌釜的RTD实验曲线和多级混合模型拟合曲线较为接近，多级混合模型能较好地反映在通气条件下含盘管的连续搅拌釜内流体的返混情况。

(b) 当其他条件相同，随着搅拌转速或表观气速的增大，含盘管的连续搅拌釜的无因次方差逐渐增大。当转速为360 r×min-1、表观气速从0增大为1.70 mm×s-1、进料流量为2.25～2.75 L×min-1，推进桨、四斜叶桨、翼型桨的平均增幅分别为14.48% 、11.71 %、12.52%。

(c) 当其他条件相同，相比翼型桨，四斜叶桨和推进桨对促进物料混合的作用较大。因此对于强放热反应，推荐使用四斜叶桨或推进桨，可以较好地增强釜内物料的混合。

ct— t时刻釜出口处的示踪剂浓度，mol×L-1us— 表观气速，mm×s-1 E(t)— 停留时间分布密度函数g— 釜出口处流体的电导率，μS×cm-1 E(θ)— 无因次化的停留时间分布密度函数gm— 实验直接测量的电导率值，μS×cm-1 N— 多级混合模型的模型参数g∞— 测量时间足够长时的电导率值，μS×cm-1 n— 搅拌转速，r×min-1θ— 对比时间 qV— 进入搅拌釜的液体体积流量，L×min-1n— 示踪剂加入量，mol t— 时间，s— 方差 — 平均停留时间，s— 无因次方差 Δt— 实验记录的时间间隔，s

[1] 覃光明, 葛忠学. 含能化合物合成反应与过程 [M]. 北京: 化学工业出版社, 2011.

QIN G M, GE Z X. Synthesis reaction and process of energetic compounds [M]. Beijing: Chemical Industry Press, 2011.

[2] 毕纪葛, 潘万贵, 周俊超, 等. 四斜叶桨搅拌下釜内盘管非稳态对流传热过程的模拟和实验研究 [J]. 高校化学工程学报, 2015, 29(4): 780-788.

BI J G, PAN W G, ZHOU J C,. Simulation and experimental study on unsteady convective heat transfer process of inner tube in stirred tank with four inclined impeller [J]. Journal of Chemical Engineering of Chinese Universities, 2015, 29(4): 780-788.

[3] 车圆圆, 周俊超, 毕纪葛, 等. 改进CBY桨搅拌釜内单相流体流动与传热特性研究 [J]. 高校化学工程学报, 2014, 28(3): 489-496.

CHE Y Y, ZHOU J C, BI J G,. Study on single phase fluid flow and heat-transfer performance in a stirred tank with an improved CBY hydrofoil impeller [J]. Journal of Chemical Engineering of Chinese Universities, 2014, 28(3): 489-496.

[4] 罗康碧, 罗明河, 李沪萍. 反应工程原理 [M]. 北京: 科学出版社, 2005.

LUO K B, LUO M H, LI H P. Principles of reaction engineering [M]. Beijing: Science Press, 2005.

[5] ACQUAYE F Y. Experimental residence time distribution (RTD) studies on effects of axial and radial flow impellers on hydrodynamic parameters of stirred vessels [D]. Accra, University of Ghana, 2016.

[6] 董红星, 杨晓光, 王兴超, 等. 连续搅拌釜流场数值模拟及停留时间分布 [J]. 石油和化工设备, 2008, 63(3): 19-23.

DONG H X, YANG X G, WANG X C,. Numerical simulation of flow field and residence time distribution in continuous stirred tank [J]. Petroleum and Chemical Equipment, 2008, 63(3): 19-23.

[7] 张斌, 翁惠新. 卧式搅拌床反应器的返混特性 [J]. 华东理工大学学报(自然科学版), 2010, 36(6): 749-754.

ZHANG B, WENG H X. Backmixing characteristics of horizontal stirred bed reactor [J]. Journal of East China University of Science and Technology(Natural Science Edition), 2010, 36(6): 749-754.

[8] ZHANG L, PAN Q, REMPEL G L. Residence time distribution in a multistage agitated contactor with Newtonian fluids: CFD prediction and experimental validation [J]. Industrial & Engineering Chemistry Research, 2007, 46(11): 3538-3546.

[9] 周国忠, 高正明, 王英琛, 等. 多级搅拌槽内流动特性实验及模型研究 [J]. 北京化工大学学报(自然科学版), 2001, 28(1): 18-21.

ZHOU G Z, GAO Z M, WANG Y C,. Experimental and model study on flow characteristics in multistage stirred tank [J]. Journal of Beijing University of Chemical Technology(Natural Science Edition), 2001, 28(1): 18-21.

[10] 皇甫功凯, 李伟, 肖文德, 等. 间二甲苯氧化反应器通气条件下液相停留时间分布的测定 [J]. 石油化工, 2005, 34(6): 565-568.

HUANGFU G K, LI W, XIAO W D,. Determination of liquid residence time distribution in-xylene oxidation reactor under aeration [J]. Petrochemical Industry, 2005, 34(6): 565-568.

[11] 李攀. 生物合成麦芽糖酯反应器的特性研究和优化 [D]. 无锡: 江南大学, 2008.

LI P. Characteristics and optimization of the reactor for biosynthesis of maltose esters [D]. Wuxi: Jiangnan University, 2008.

[12] 姚占静, 郭睿, 张春生, 等. 流动反应器内返混研究 [J]. 化工科技市场, 2008, 31(10): 28-30.

YAO Z J, GUO R, ZHANG C S,. Study on backmixing in flow reactor [J]. Chemical Technology Market, 2008, 31(10): 28-30.

[13] SARAVANATHAMIZHAN R, PARANTHAMAN R, BALASUBRAMANIAN N,Residence time distribution in continuous stirred tank electrochemical reactor [J]. Chemical Engineering Journal, 2008, 142(2): 209-216.

[14] 李伟. 多层斜叶桨混合特性的研究 [D]. 上海：华东理工大学, 2015.

LI W. Research on mixing characteristics of multi-layer inclined blade propeller [D]. Shanghai: East China University of Science and Technology, 2015.

[15] 李希. 化工问题的建模与数学分析方法 [M]. 北京: 化学工业出版社, 2006.

LI X. Modeling and mathematical analysis of chemical engineering problems [M]. Beijing: Chemical Industry Press, 2006.

[16] 张本贺, 何宇晨, 毕纪葛, 等. 基于本质安全的硝化反应釜的概念设计 [J]. 高校化学工程学报, 2015, 29(2): 312-319.

ZHANG B H, HE Y C, BI J G,. Conceptual design of nitrification reactor based on intrinsic safety [J]. Journal of Chemical Engineering of Chinese Universities, 2015, 29(2): 312-319.

[17] 周俊超, 车圆圆, 吴可君, 等. 基于CFD模拟的强放热反应釜盘管空间排布的优化设计 [J]. 高校化学工程学报, 2015, 29(1): 27-34.

ZHOU J C, CHE Y Y, WU K J,. Optimization design of coil pipe space layout of high exothermic reactor based on CFD simulation [J]. Journal of Chemical Engineering of Chinese Universities, 2015, 29(1): 27-34.

[18] 李绍芬. 反应工程 [M]. 3版. 北京: 化学工业出版社, 2013.

LI S F. Reaction engineering [M]. 3rd ed. Beijing: Chemical Industry Press, 2013.

[19] 陈甘棠. 化学反应工程 [M]. 3版. 北京: 化学工业出版社, 2007.

CHEN G T. Chemical reaction engineering [M]. 3rd ed. Beijing: Chemical Industry Press, 2007.

Study on residence time distribution of continuous stirred tanks with coil under different impeller types

HUANG Cong1, WU Ke-jun1, HE Chao-hong1,2

(1. Zhejiang Provincial Key Laboratory of Advanced Chemical Engineering Manufacture Technology,College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310027, China;2. Institute of Zhejiang University-Quzhou, Quzhou 324000, China)

Continuous stirred tank reactors for energetic compound preparation are usually equipped with spiral coils, and the flow behaviors in reactors are complicated. A set of continuous stirred tank system equipped with spiral coil was built. Residence time distribution (RTD) in the continuous stirred tank reactor using four-pitched blade turbine (PBT), propeller (PRO) and hydrofoil impellers (CBY) were studied. Effects of impeller type, feed flow rate, impeller speed and superficial gas velocity on dimensionless variance were determined, and the applicability of the multilevel mixing model for the system with gas generation in the continuous stirred tank reactor with spiral coil was investigated. The results show that the multilevel mixing model is suitable for describing non-ideal flow in a continuous stirred tank reactor equipped with coil, which can be used for design optimization of continuous stirred tank reactors equipped with coil. The dimensionless variance increases with the increase of impeller speed and superficial gas velocity. For a strong exothermic reaction, it is recommended to use PBT or PRO impeller which can strengthen the mixing of the flow in the tank.

continuous stirred tank reactor; spiral coil; residence time distribution(RTD); energetic compounds

1003-9015(2022)04-0554-08

TQ051.7；
TQ217

A

10.3969/j.issn.1003-9015.2022.04.011

2021-05-27；

2022-02-23。

黄聪(1992-)，男，湖南衡阳人，浙江大学硕士生。

何潮洪，E-mail：chhezju@zju.edu.cn

黄聪, 吴可君, 何潮洪. 不同桨型下含盘管的连续搅拌釜停留时间分布的研究 [J]. 高校化学工程学报, 2022, 36(4): 554-561.

：
HUANG Cong, WU Ke-jun, HE Chao-hong. Study on residence time distribution of continuous stirred tanks with coil under different impeller types [J]. Journal of Chemical Engineering of Chinese Universities, 2022, 36(4): 554-561.

猜你喜欢 示踪剂盘管表观 分层示踪剂监测技术在河南油田稠油水驱油藏的研究与应用石油地质与工程(2022年4期)2022-08-06中央空调机组预热盘管防冻控制方法建筑热能通风空调(2022年6期)2022-07-30示踪剂技术在压裂效果评价中的研究进展石油化工应用(2022年3期)2022-04-20基于表观遗传学的针刺治疗慢性疼痛的机制研究世界科学技术-中医药现代化(2021年8期)2021-12-21绿盲蝽为害与赤霞珠葡萄防御互作中的表观响应河北果树(2021年4期)2021-12-02水平井压裂稀土元素示踪监测技术化工管理(2021年23期)2021-08-25缝洞型油藏井间示踪剂分类等效解释模型及其应用西安石油大学学报（自然科学版）(2021年1期)2021-01-27E- GAS气化焦过滤器伴热盘管泄漏分析化工管理(2020年17期)2020-07-17小型洁净室干盘管及其冷凝水管道施工优化建材发展导向(2019年10期)2019-08-24例析对高中表观遗传学的认识福建基础教育研究(2019年10期)2019-05-28