聚丙烯环管反应器的流体力学模拟.docxVIP

聚丙烯环管反应器的流体力学模拟

宣爱国;任苏孟;吴元欣;朱晓明;闫志国浏玮【摘要】采用颗粒动力学为基础的Euler-Euler双流体模型研究反应器内液固相的流动力学特性.通过fluent软件，采用标准k-s模型描述流体的湍流状态、Gidaspow曳力模型描述流体的相互作用力，运用SIMPLE算法求解速度场.结果表明:在5~9m/s的浆液速度范围内，环管反应器直管段的固相体积分数分布均匀，弯管段的固相体积分数分布不均匀;浆液速度增大，其固相体积分数分布的非均匀性增加;当固相体积分数为0.35时，浆液密度为563~571千克/立方米.通过对速度场的分析可知，上升管中段的流速为6~7.5m/s,速度呈对称分布;下降管中段的流速为5~8m/s,速度呈非对称分布.计算结果与工厂实际生产情况接近，表明欧拉双流体模型能有效地描述环管反应器内浆液流动形态.

【期刊名称】《武汉工程大学学报》【年（卷），期】2014（036）004【总页数】5页（P17-21）【关键词】聚丙烯;环管反应器;Euler-Euler双流体模型;计算流体力学【作者】宣爱国;任苏孟;吴元欣;朱晓明;闫志国浏玮

【作者单位】武汉工程大学化工与制药学院，湖北武汉430074;绿色化工过程教育部重点实验室（武汉工程大学），湖北武汉430074;武汉工程大学化工与制药学院，湖北武汉430074;绿色化工过程教育部重点实验室（武汉工程大学），湖北武汉430074;武汉工程大学化工与制药学院，湖北武汉430074;绿色化工过程教育部重点实验室（武汉工程大学），湖北武汉430074;中石化武汉分公司，湖北武汉430082;武汉工程

大学化工与制药学院，湖北武汉430074；绿色化工过程教育部重点实验室（武汉工程

大学），湖北武汉430074;武汉工程大学计算机科学与工程学院，湖北武汉430074【正文语种】中文

【中图分类】TQ021.1

0引言

随着合成工艺技术的不断发展，聚丙烯（PP）树脂以其优异的性能越来越广泛地应用于各个领域.它的聚合反应器主要有环管反应器、本体流化床、气相搅拌釜等多种类型，其生产技术来自于Basell的Spheripol工艺、Borealis的bostar工艺、三井化学的hypol工艺、Basf的Novolen工艺等.其中Spheripol工艺的环管反应器具有结构简单、易于工程放大、能耗低等优点而在均聚反应中倍受青睐［1-2］.然而在丙烯聚合过程中，聚丙烯的产出与操作流型、固体表面的流体边界层结构、压降、多相流体之间的相互作用等量化关系尚不清楚，需要借助于计算流体力学来模拟［3］.

本研究以中石化武汉分公司的环管反应器为研究对象，选用将颗粒相处理为〃拟流体”的Euler-Euler双流体模型［4］,模拟环管反应器中的流动形态，研究反应器内浆液密度、固相体积分数分布以及速度场的变化情况.

1环管反应器模型的建立

1.1丙烯聚合动力学

丙烯聚合反应按链引发、链转移、链终止一般历程可得速率方程：

其中，［M］为丙稀单体浓度，kmol/m3;［C*］为催化剂活性中心浓度，kmol/m3;［H2］为氢气浓度，kmol/m3.

3种反应速率常数服从Arrhenius方程，其指前因子和活化能见表1[5].

表1丙烯聚合动力学参数表Table1Propylenepolymerizationkineticsparameters指前因子反应活化能/(kJ-mol-1kp0/(L-mol-l-s-1)

7.65x106Ep)50ktr0/(L?mol-1?s-1)7.65x106Etr50kd0/s-1

7.65x106Ed50

1.2Euler-Euler双流体模型

模型假设反应器内的浆液为液相丙烯和固相聚丙烯.分别对液相和固相建立连续性

方程和动量守恒方程.

液固相连续性方程：

液固相动量守恒方程：

上述方程中存在固体压力ps、体积粘度入s、固体剪切粘度卜等未知量，选取应

用最多的颗粒动力学理论建立其计算公式[6]:

其中代表颗粒温度，

g0=代表径向分布函数.

1.3物理模型及物性参数

模拟的环管反应器由6根①600x26000的直管构成，管间距为4.2m,其几何结构见图1(a).在FLUENT的前置处理器Gambit软件中构建出其几何模型，反应器直管段采用copper网格，弯管采用tgrid网格划分，直管的径向截面网格及弯管段的网格示意图见图1(b)、1(c).P-1至P-6依次代表上升段及下降段，初始、中间、末端三个截面.

图1环管反应器的物理模型及网格划分Fig.1Physical