固定床流体力学

1、固定床流体力学第1页，共29页，2022年，5月20日，9点51分，星期二 固定床中进行催化剂反应时，同时发生传热及传质过程，后两者又与流体在床层内的流动状况密切有关。为了研究固定床中化学反应的宏观反应过程，进行合理的反应器结构设计，必须先讨论固定床的传递过程，即固定床中的流体力学、传热及传质问题。第2页，共29页，2022年，5月20日，9点51分，星期二1.非中空固体颗粒的当量直径及形状系数 常用的非中空颗粒当量直径的表示方法有三种，即等体积圆球直径、等外表面积圆球直径和等比外表面积圆球直径。 若以Sp和Vp表示非中空颗粒的外表面积和体积，按等体积圆球直径计算的当量直径可表示如下： (5-

2、1) 式（5-1）中Vp为与非中空颗粒等体积的圆球体积。一. 固定床的物理特性第3页，共29页，2022年，5月20日，9点51分，星期二 按等外表面积圆球直径计算的当量直径Dp可表示如下： (5-2) 式（5-2）中Sp为非中空颗粒等外表面积的圆球的外表面积。 按等比外表面积圆球直径计算的当量直径ds可表示如下： (5-3) 式（5-3）中Sv 为与非中空颗粒等比外表面积的圆球比外表面积。第4页，共29页，2022年，5月20日，9点51分，星期二 再以SS表示与非中空颗粒等体积的圆球的外表面积，则 (5-4) 因此，引入一个量纲1数 ，称为颗粒的形状系数，其值如下： (5-5)第5页，共2

3、9页，2022年，5月20日，9点51分，星期二 对于球形颗粒， =1；对于非球形颗粒， 小于1。形状系数说明了颗粒形状与圆球的差异程度。 形状系数 可由颗粒的体积及外表面积算得。非中空颗粒的体积可由实验测得，或由其质量及密度计算。形状规则的颗粒，例如圆柱形颗粒，其外表面积可由直径及高求出；形状不规则的颗粒外表面积却难以直接测量，这时可测定由待测颗粒所组成的固定床压力降来计算形状系数。第6页，共29页，2022年，5月20日，9点51分，星期二 上述三种当量直径dv、Dp、ds与形状系数间的相互关系可表示如下 (5-6) 及 (5-7)第7页，共29页，2022年，5月20日，9点51分，星期

4、二2.混合颗粒的平均直径及形状系数 破碎成的碎块形状不规则，大小也不均匀 算术平均直径 为 (5-8) di两相邻筛孔净宽乘积的平方根；xi为直径di颗粒的质量分数。 调和平均直径 为 (5-9) 在固定床及流化床的流体力学中，用调和平均直径较为符合实验数据。大小不等形状各异的混合颗粒 由固定床P计算。第8页，共29页，2022年，5月20日，9点51分，星期二3. 固定床的当量直径 床层中颗粒的比表面积（不计入接触而减少的表面积）： (5-10) 水力半径： 固定床的当量直径： (5-11) 床层由中空颗粒，如单孔环柱体，多通孔环柱体等组成时，不能使用式（5-11）。第9页，共29页，202

5、2年，5月20日，9点51分，星期二4.固定床的空隙率及径向流速分布 与下列因素有关：颗粒形状，颗粒的粒度分布，充填方式，dP /dt之比（壁效应） 图5-13表达了上述关系，当dTdP时，壁效应可不计，一般工程上认为当di /dP8时，可不计壁效应。 图5-14（b）表示不同Re数时流体的径向分布，Re大分布较平坦。第10页，共29页，2022年，5月20日，9点51分，星期二二. 单相流体在固定床颗粒层中的流动及压力降 固定床中，流体在颗粒物料组成的孔道中流动，孔道相互交错联通、弯曲，各孔道的几何形状相差甚大，孔截面积也很不规则且不相等。流体在固定床中流动时，旋涡的数目比空管多，由滞留转入

6、湍流是一个逐渐过渡的过程。1. 流动特性第11页，共29页，2022年，5月20日，9点51分，星期二2. 单相流体通过固定床颗粒层的压力降 单相流体通过固定床时要产生压力损失，主要来自颗粒的黏滞曳力，即流体与颗粒表面间的摩擦，流体流动过程中孔道截面积突然扩大和收缩，及流体对颗粒的撞击及流体的再分布。低流速时，压力降主要是由于表面摩擦而产生，高流速及薄床层中流动时，扩大、收缩则起着主要作用。如果容器直径与颗粒直径之比值较小，还应计入壁效应对压力降的影响。第12页，共29页，2022年，5月20日，9点51分，星期二 计算单相流体通过固定床颗粒层压力降的方法很多，其中有许多都是利用流体在空圆管中

7、流动时的压力降公式，加以合理地修改而成。下面介绍其中的一个使用最广泛的Ergun方法。流体在空圆管中作等温流动，且流体密度的变化可不计时，压力降可表示为 (5-12) 式中：L管长，m；d圆管的内直径，m； 流体的密度，kg/m3；u0e流体的平均流速，m/s；f摩擦系数，量纲1数； 压力降，N/m2。第13页，共29页，2022年，5月20日，9点51分，星期二 上式应用于固定床时，u0e应为流体在床层空隙中的真正平均流速ue，圆管的直径应以固定床的当量直径de代替，而管长则应以流体在固定床中的流动途径来代替。将 和 ， 代入式(5-12)，又考虑到流体在固定床中的流动途径远大于固定床的高度

8、L，并等于L的若干倍，则固定床的压力降可表示为 (5-13) 式中u0以床层不含内构件和所载固体的空截面积Ac和操作状态下气体体积流量计算的流速，或称表观流速，m/s；fM修正摩擦系数。第14页，共29页，2022年，5月20日，9点51分，星期二 经多种颗粒和工况的实验测定，修正摩擦系数fM与修正雷诺数的关系可表示如下 (5-14) 而 (5-15) 式中 流体黏度， ；G流体质量流率或质量通量， 。第15页，共29页，2022年，5月20日，9点51分，星期二 当ReM1.75，即式(5-13)可化简为 (5-16) 当ReM1000时，处于完全湍流状况，式(5-16)中150/ReM10

9、00时， 正比于G2L，即VR一定时， 正比于(W/Ac)2L，或 正比于L3。改为径向流动，即减小了气流在催化床内流动的路程，同时增加了垂直于流向的截面积和减小了气流速度，从而大幅度地减少了床层的压力降，为使用小颗粒催化剂，提高催化剂的总体速率及催化床的空时产率创造了良好的条件。第20页，共29页，2022年，5月20日，9点51分，星期二 颗粒粒度和形状是影响床层压力降的另一重要因素。形状相同的颗粒，减小颗粒的当量直径，会导致固定床压力降增加。当ReM1000时，压力降反比于当量直径。颗粒的筛析范围相同，形状系数小的颗粒，如片状，其当量直径减小，床层的压力降增大。第21页，共29页，202

10、2年，5月20日，9点51分，星期二 床层空隙率对压力降的影响十分显著，当ReM1000时，压力降正比于 ， 由0.4增至0.5时，压力降可降至原来的1/2.3。床层空隙率的大小与颗粒的形状、粒度分布、填充方法、颗粒直径与容器直径之比值等因素有关。混合颗粒的粒度越不均匀，小颗粒填充在大颗粒之间，所组成的床层空隙率越小。催化剂在使用过程中逐渐破碎、粉化，当质量流率不变时，由于空隙率减小，床层压力降相应地逐步增大。催化剂使用后期床层压力降较前期压力降增加的程度随催化剂的机械强度而定。即使不计入破损，操作一段时期后，由于床层中颗粒填实，使床层下沉，空隙率降低而增高压力降。第22页，共29页，2022

11、年，5月20日，9点51分，星期二三. 径向流动反应器中流体的分布径向流动反应器的优点：流体流通截面积大、流速小、流道短，床层压力降小。径向流动反应器的结构如图5-9。第23页，共29页，2022年，5月20日，9点51分，星期二四. 固定床流体的径向及轴向混合 当流体流经固定床时，不断发生分散与汇合，形成了一定程度的径向及轴向混合，尤其当固定床中进行化学反应而又与外界换热时，床层中不同径向位置处流速、温度及反应速率都不相同，也就必然存在着径向浓度分布，更加加剧床层中径向及轴向的混合过程，而其中径向混合比轴向更加显著。1. 固定床径向及轴向混合有效弥散系数 （effective dispers

12、ion coefficient）第24页，共29页，2022年，5月20日，9点51分，星期二 径向混合有效弥散系数Er及轴向混合有效弥散系数Ez一般用Peclet数(Pe)表示，此时 (5-35-A) 或 (5-35-B) 表征径向混合及轴向混合有效弥散系数的Pe数与Re数的关系见图5-17，由图可见，当Re40，即处于湍流状态时，无论对于气体还是液体，径向Per10，几乎不随Re数而变，气体的轴向Pez 2不随雷诺数而变，但液体的轴向值Pez随Re值而有一定程度的变化。第25页，共29页，2022年，5月20日，9点51分，星期二 图5-17 固定床轴向及径向Pe数与Re数的关系第26页，共29页，2022年，5月20日，9点51分，星期二2. 固定床反应器中的轴向返混 第三章已经讨论过可以略去不计逆向混合影响的条件是模型参数Ez/(uL)0.005。在固定床反应器的流动状况下，一般Re40，此时Pez=dsu/Ez=2由此可得 0.5(ds/L)0.005 (5-36) 即固定床反应器中床层高度L超过颗粒直径ds的一百倍时，可以略去轴向返混的影响。第27页，共29页，20