1、2022年7月31日星期日 2022年7月31日星期日第四章 非理想流动第一节 停留时间分布1 停留时间分布的定量描述2 停留时间分布的实验方法3 寿命与年龄分布的关系 2022年7月31日星期日第一节 停留时间分布停留时间:停留时间:流体从进入系统时算起,到其离开系统时流体从进入系统时算起,到其离开系统时 为止,在系统内总共经历的时间,即为止,在系统内总共经历的时间,即流体流体 从系统的进口至出口所耗费的时间。从系统的进口至出口所耗费的时间。寿命分布:寿命分布:指流体粒子从进入系统到离开系统的停留指流体粒子从进入系统到离开系统的停留时间。时间。年龄分布:年龄分布:指流体粒子进入系统在系统中停
2、留的时间。指流体粒子进入系统在系统中停留的时间。2022年7月31日星期日区别:区别:寿命分布指的是寿命分布指的是系统出口处系统出口处的流体粒子的停留时间,的流体粒子的停留时间,年龄分布是对年龄分布是对系统中的系统中的流体粒子而言的停留时间。流体粒子而言的停留时间。实际测定得到的且应用价值又较大的是寿命分布。实际测定得到的且应用价值又较大的是寿命分布。通常所说的通常所说的停留时间分布指的是寿命分布停留时间分布指的是寿命分布。2022年7月31日星期日1 1 停留时间分布的定量描述停留时间分布的定量描述1 1)停留时间分布函数)停留时间分布函数出口中红色粒子数t2022年7月31日星期日 tEt
3、 tFt dttEdtt 2022年7月31日星期日 tE:停留时间分布密度函数,或寿命分布密度函数。停留时间分布密度函数,或寿命分布密度函数。0t时:0tE0t时:0tE且 10dttE tdttEtF0 停留时间分布函数停留时间分布函数:停留时间小于停留时间小于t t的流体粒子所的流体粒子所 占的分数占的分数 1,00FF可改写成:可改写成:dttdFtE2022年7月31日星期日可以用可以用年龄分布密度函数年龄分布密度函数tI和和年龄分布函数年龄分布函数ty来描述流体在反应器内的停留时间分布。来描述流体在反应器内的停留时间分布。dttdytI tdttIty0 00dttIy 00 y2
4、022年7月31日星期日)平均停留时间)平均停留时间t 00dttEdtttEt平均停留时间平均停留时间t应是应是 tE曲线的分布中心,即曲线的分布中心,即 tE在所围的面积的重心在在所围的面积的重心在t t坐标上的投影坐标上的投影 曲线曲线 0dtttE 0dtdttdFt 10tFtFttdF)停留时间分布函数的特征值停留时间分布函数的特征值在数学上称在数学上称t t为为 tE曲线对于坐标原点的曲线对于坐标原点的一次矩一次矩,又称,又称 tE的的数学期望数学期望。2022年7月31日星期日)方差)方差2t表示表示停留时间分布的分散程度的量停留时间分布的分散程度的量,在数学上是指对于平,在数
5、学上是指对于平均停留时间的二次矩。均停留时间的二次矩。0022dttEdttEttt 022tdttEt 02dttEtt2022年7月31日星期日令:无因次时间令:无因次时间:tt则:无因次平均停留时间无因次平均停留时间1tt 001dttEdE tEtE tEtEdttEdE)()(tt tFF tdFdF2022年7月31日星期日 0221dE221tt 0221dttEttt 021dttEtt 021dEtt若以若以2表示以表示以为自变量的方差,则它与为自变量的方差,则它与2t的关系为:的关系为:2022年7月31日星期日2 2 停留时间分布的实验方法停留时间分布的实验方法1)1)脉
6、冲示踪法脉冲示踪法主流体0vt=0示踪剂0v2022年7月31日星期日0t tc0tt tc响应曲线2022年7月31日星期日设示踪剂加入量为设示踪剂加入量为m m(根据(根据A A的物料衡算)的物料衡算)00dtcvmAdtdtCCmdtCvAAA00 mcvtEA000dtcvA0dtccAA000dtcvcvAAdttmE)(2022年7月31日星期日为了验证实验数据的可靠性,必须根据三个已知量为了验证实验数据的可靠性,必须根据三个已知量m m,0v和Rv进行一致性检验进行一致性检验,即由实验所获得的,即由实验所获得的 tcA数据或数据或 tcctEA0数据应满足:数据应满足:000dt
7、cvmcA 0cctEAt(空时)0vVR 0dtttEt0vvR是否与是否与 0dtttE相等,若不相等,则须检查原因。相等,若不相等,则须检查原因。2022年7月31日星期日由脉冲法测定某反应器的响应曲线结果如下:由脉冲法测定某反应器的响应曲线结果如下:例:例:t(min)00.511.522.533.544.555.566.5702522272622191510743310Ac计算计算min5t时的时的 tEtF,和t解:解:070dtcdtcAA 0dtcctEAA3333.813333.81Ac2022年7月31日星期日 3333.8145ttE tdttEtF0 505dttEtF
8、t 0dtttEt0492.09467.03333.8177503333.811dtcAmin5369.23333.810dttcA2022年7月31日星期日0v流体0v示踪剂系统0v检测示踪剂0Ac0tAct0Ac tc0t入口出口2)2)阶跃示踪法阶跃示踪法升阶跃升阶跃Ac2022年7月31日星期日ActAct0t0Ac入口出口降阶跃降阶跃2022年7月31日星期日如果如果t t时刻出口物料中时刻出口物料中A A的浓度为的浓度为Ac 0AAcctF对于降阶法:对于降阶法:01AAcctF在在tcA图中应满足:图中应满足:对于升阶跃对于升阶跃2022年7月31日星期日00AcACtdcA或:
9、或:00AcACtdcA应用上式应用上式进行一致性检验进行一致性检验 100ttdFcA 00dtttEcAtcA000vVcRA0Ac2022年7月31日星期日3 寿命与年龄分布的关系寿命与年龄分布的关系 设在定常流动系统中,对于恒容过程,在设在定常流动系统中,对于恒容过程,在0 0 t t时间内对时间内对示踪剂示踪剂A A进行物料衡算:进行物料衡算:输入输入tAdtcv000(0 t 内流入系统中内流入系统中A A 的量)的量)输出输出toAtdtFcv00(0 0t t 内由系统流出内由系统流出A A的量)的量)留在系统中的量留在系统中的量 tARdttICV00)(2022年7月31日
10、星期日 dttIcVdttFcvtARtA000001 tI ttItF1 dttdItE有时采用有时采用 强度函数强度函数 tItEt1表示两者的关系表示两者的关系:tItEt1 dtt表示容器中年龄为表示容器中年龄为t t的粒子在的粒子在t t到到tt系统的粒子所占的分率。系统的粒子所占的分率。时间内离开时间内离开 dttItdI tdytdI2022年7月31日星期日第二节 理想流动及层流的停留时间分布1 平推流2 全混流模型3 流体在圆管内作层流流动时的停留时间分布2022年7月31日星期日第二节第二节 理想流动及层流的停留时间分布理想流动及层流的停留时间分布1 平推流平推流 tEtt
11、 tFtt0.12022年7月31日星期日 tEtt 0tt 即:tttE性质:21ttt tf0 21ttdttttf其他2022年7月31日星期日 tdttEtF0tt 0tt 1tdttt0数字特征为:02dtttt0202tt表明所有的流体粒子在反应器内的停留时间相同。方差越小,说明分布越集中,分布曲线就越窄,停留 时间分布方差等于零这一特征说明系统内不存在返混。2022年7月31日星期日2 2 全混流模型全混流模型AC0AC0v2022年7月31日星期日单位时间内反应器内示踪剂做物料衡算得:单位时间内反应器内示踪剂做物料衡算得:输入输入00Acv输出输出Acv0积累积累dtdcVdt
12、cVdARARAAARcvcvdtdcV000AARAccVvdtdc00AAccv00AAcc02022年7月31日星期日即:001AAAAccdtccd 0AAcctF tFdttdF1积分:积分:tFtdttFtdF0011 ttF11lnttF)(1ln2022年7月31日星期日对恒容:对恒容:ttt tetF1 dttdFtE eEe1 Fte1et12022年7月31日星期日11 tEt tFt数字特征:数字特征:0221dE1102de 021dE2022年7月31日星期日10de可见:可见:返混程度达到最大时,停留时间分布的无因次方差返混程度达到最大时,停留时间分布的无因次方差
13、12平推流时方差平推流时方差 02实际反应器停留时间分布的方差应介于实际反应器停留时间分布的方差应介于0 01 1之间之间,值越大,值越大则停留时间分布越分散,因此,由模型模拟实际反应器时则停留时间分布越分散,因此,由模型模拟实际反应器时应从方差入手。应从方差入手。2022年7月31日星期日设两个反应器进行的反应相同,且平均停留时间相等。设两个反应器进行的反应相同,且平均停留时间相等。对于平推流反应器,对于平推流反应器,所有流体粒子的停留时间相等所有流体粒子的停留时间相等,且都,且都等于平均停留时间。等于平均停留时间。对于全混流反应器,对于全混流反应器,停留时间小于平均停留时间的流体粒停留时间
14、小于平均停留时间的流体粒子占全部流体的分率为:子占全部流体的分率为:632.01etF使停留时间分布集中,可以提高反应器的生产强度。使停留时间分布集中,可以提高反应器的生产强度。2022年7月31日星期日3 流体在圆管内作层流流动时的停留时间分布流体在圆管内作层流流动时的停留时间分布 当流体在管内作层流流动时,管内流速随距管轴心的当流体在管内作层流流动时,管内流速随距管轴心的速度分布为:速度分布为:2max1Rruu平均流速:平均流速:max5.0 uu 在管出口的流体质点在管内的停留时间在管出口的流体质点在管内的停留时间t t2max1RruLuLt2022年7月31日星期日平均停留时间:平
15、均停留时间:AuLAuLt管轴心上的流体粒子的停留时间最小,为管轴心上的流体粒子的停留时间最小,为mint22maxmintuLttt21时时,0tF,只有在只有在tt21物物料中才会有料中才会有A A,根据,根据时,管的流出时,管的流出当当 tF的定义有:的定义有:容积流量整个管截面在出口处的口处)截面内的容积流量(出rtF02022年7月31日星期日 RrrudrrudrtF0022 2max1uLuLtF dttdFtErdrRrrdrRrRr0202112211Rr2max1uu221t221t322t2022年7月31日星期日0.1 tF00t 221ttFt210vVt 两个阴影面
16、积相等2022年7月31日星期日 322tdttdFtE tEt00t212022年7月31日星期日第三节 非理想流动模型1 离析流(离集流)模型2 多级全混流串联模型3 几种常见的单参数模型4 组合模型2022年7月31日星期日第三节第三节 非理想流动模型非理想流动模型建立流动模型的方法是:建立流动模型的方法是:1.通过冷态模型实验测定装置的停留时间分布。通过冷态模型实验测定装置的停留时间分布。2.根据所得的根据所得的 tE和 tF流动模型,并根据停留时间分布的实验数据来确定所流动模型,并根据停留时间分布的实验数据来确定所提出的模型中的参数。提出的模型中的参数。的结果通过合理的简化提出可能的
17、的结果通过合理的简化提出可能的3.结合反应动力学数据通过模拟计算来预测反应结果结合反应动力学数据通过模拟计算来预测反应结果4.通过一定规模的热模实验来验证模型的准确性通过一定规模的热模实验来验证模型的准确性。2022年7月31日星期日1 离析流(离集流)模型离析流(离集流)模型反应器内的流体粒子之间不存在任何形式的物质反应器内的流体粒子之间不存在任何形式的物质交换,或者说它们之间不发生微观混合,流体粒交换,或者说它们之间不发生微观混合,流体粒子就像一个有边界的个体,从反应器的进口向出子就像一个有边界的个体,从反应器的进口向出口运动。口运动。离析流:离析流:设反应器进口的流体中反应物设反应器进口
18、的流体中反应物A A的浓度为的浓度为0Ac当反应时间为当反应时间为t t时,浓度为时,浓度为 tcA2022年7月31日星期日 dttEtcA 停留时间在停留时间在t t到到dtt 间的流体粒子所占的分率为间的流体粒子所占的分率为 dttE则这部分流体对反应器出口流体中则这部分流体对反应器出口流体中A A的浓度的浓度Ac的贡献的贡献:反应器出口处反应器出口处A A的平均浓度的平均浓度Ac 0dttEtccAA根据转化率的定义,也可写成:根据转化率的定义,也可写成:0dttEtxxAA 离析流模型离析流模型 2022年7月31日星期日例:例:等温下在反应体积为等温下在反应体积为355.4m的流动
19、反应器内进行液相反应:的流动反应器内进行液相反应:PRA2该反应为二级反应,反应温度下的反应速率常数该反应为二级反应,反应温度下的反应速率常数min/104.233molmk进料流量进料流量:min,/5.03mvoA A的浓度的浓度:30/6.1mkmolcA停留时间分布为:停留时间分布为:2022年7月31日星期日024681012141618202224014798521.510.60.20t/min 3/mgtcA试计算离析流模型反应器出口处试计算离析流模型反应器出口处A A的转化率的转化率Ax2022年7月31日星期日解:解:A A的转化率可由模型方程求取的转化率可由模型方程求取应先
20、求出应先求出Ac与与t t的关系,积分二级反应速率方程:的关系,积分二级反应速率方程:2AAAcdcr积分:积分:AAcctAAkdtcdc002得:得:tkcctcAAA001 0dttEtccAA2Akc2022年7月31日星期日还应先求出还应先求出 tE 0dttctctE 0001dttEtkcccAAA00242220181614121086420t tE781784787789788785782785.1781786.0782.03/05447.0mkmol2022年7月31日星期日转化率为:转化率为:966.000AAAAcccx若用平推流:若用平推流:0vVR001AAAkcc
21、c00AAAAcccx两者结果相近,原因是该反应器的停留时间分布与平推流两者结果相近,原因是该反应器的停留时间分布与平推流偏离不算太大的缘故。偏离不算太大的缘故。min11.95.055.43/04374.0mkmol9728.02022年7月31日星期日2 多级全混流串联模型多级全混流串联模型0C0v1V1C0v2V2C0vNV0v1NC0vNC12N2022年7月31日星期日1)多釜串联时的停留时间分布多釜串联时的停留时间分布对第对第i i釜作示踪剂的物料衡算:釜作示踪剂的物料衡算:输入输入10Aicv输出输出Aicv0积累积累dtdcvdtcvdAiRAiRdtdcvcvcvAiRAiA
22、i0102022年7月31日星期日iAiAiAiccdtdc10vVRi初始条件:初始条件:t=0=0时,时,).,3,2,1(0nicAi当当i i=1=1时,时,101AAAiccdtdc积分:积分:tAAecc1012i时,21AAAccdtdc0t且时,02Ac2022年7月31日星期日1Ac代入得代入得:tAAAeccdtdc11022积分得:积分得:tAAetcc1102递推求解得:递推求解得:NiitAANitecc110!11多釜串联停留时间分布函数多釜串联停留时间分布函数2022年7月31日星期日若以系统的总平均停留时间代入若以系统的总平均停留时间代入,Nt有有:NiitNt
23、iNtetFt11!1)(1令令tt,写成无因次形式写成无因次形式:NiiNiNeF11!11对对求导,可得多釜串联模型的停留时间分布密度:求导,可得多釜串联模型的停留时间分布密度:ddFENNNeNN1!12022年7月31日星期日2 2)E曲线特征曲线特征 E5.00.15.00.11N10N522022年7月31日星期日1N2510N F0.15.05.00.15.12022年7月31日星期日从从 E曲线可见,曲线可见,1N全混流串联模型的全混流串联模型的 E均出现均出现 maxE峰值,相应此峰值的无因次停留时间以峰值,相应此峰值的无因次停留时间以max记,将记,将 E对求导求导,并令导
24、数并令导数=0,可求得此极值。可求得此极值。曲线曲线NN1max 11max!11NNeNNNE在曲线在曲线 maxE的两侧均有一个拐点:的两侧均有一个拐点:022dEd2022年7月31日星期日021122inf2infNNNNN即:即:111max1infN111max2infN数字特征:数字特征:0221dEN111NN1!1012NNNeNN2022年7月31日星期日1N时,12,与全混流模型一致与全混流模型一致 N时,02,与平推流一致与平推流一致 当当N N为任何正数时,其方差应介于为任何正数时,其方差应介于0 0与与1 1之间,对之间,对N N的不的不同取值可模拟不同的停留时间分
25、布。同取值可模拟不同的停留时间分布。如何应用多釜全混流模型来模拟一个实际反应器的流动状况如何应用多釜全混流模型来模拟一个实际反应器的流动状况?2022年7月31日星期日例:例:按上例,由多级全混流串联模型计算,由已知数据可计算按上例,由多级全混流串联模型计算,由已知数据可计算 0221dE21N3824.052.555.4m故采用故采用5 5个个3824.0m的大釜和一个的大釜和一个343.0m的小釜的小釜1812.052.52022年7月31日星期日21AiAiiAiAikcrcc可得:iAiiAikckc214115.0824.0i5.043.06)5.3,2,1min(648.1imin
26、86.02022年7月31日星期日30/6.1mkmolcA33/1587.0mkmolcA32/2583.0mkmolcA31/522.0mkmolcA34/1104.0mkmolcA35/083.0mkmolcA060AAAAcccx36/0723.0mkmolcA9548.02022年7月31日星期日3 3 几种常见的单参数模型几种常见的单参数模型1 1)有死区的全混流模型)有死区的全混流模型设反应器内的有效容积为设反应器内的有效容积为RV全混流状态的体积为全混流状态的体积为RmVf死区的体积为死区的体积为RmVf1此时,相当于一个实际有效容积为此时,相当于一个实际有效容积为RmVf的全
27、混流反应器的全混流反应器:令:令:mRmmfvVf02022年7月31日星期日由物料衡算可知:由物料衡算可知:mtmetE1 mtetF1相应的平均停留时间:相应的平均停留时间:mt此模型并未改变全混流的特征,只是缩短了反应流体在器此模型并未改变全混流的特征,只是缩短了反应流体在器内的平均停留时间。内的平均停留时间。2 2)具有死区的平推流模型)具有死区的平推流模型死区的存在并不改变平推流的特性,只是减少了反应器的死区的存在并不改变平推流的特性,只是减少了反应器的实际有效容积。实际有效容积。2022年7月31日星期日设平推流体积为设平推流体积为 RpVf死区体积为死区体积为RpVf1所以反应流
28、体在反应器内的停留时间:所以反应流体在反应器内的停留时间:ppt即即 tE曲线比无死区场合提前出峰曲线比无死区场合提前出峰0vVfRp2022年7月31日星期日0vvfttRpp0 tEptt ptt ptt 01 tF022022年7月31日星期日3)3)平推流与全混流串并联模型平推流与全混流串并联模型0v0v0vVfm1Vfm平推流区平推流区全混流区全混流区)按串联方式组合(串联模型)按串联方式组合(串联模型)2022年7月31日星期日0v0v0vVfm1Vfm平推流区平推流区全混流区全混流区设反应器内全混流的体积为RmVf平推流体积为RmVf12022年7月31日星期日该模型的该模型的
29、tE tF形式为:形式为:0vvfRmmptt0mpttme1 tEpttptt0mptte1 tF0)1(vVftRmp2022年7月31日星期日)并联模型)并联模型VfmVfm10v0v0vIm01vIm2022年7月31日星期日停留时间分布为:停留时间分布为:dttEcIvdttEcIvdttEcvmAmpAmA0000001 tEItEItEmmpm 1其中:其中:011vIvftmRmpmtmmpmeIttI110vIvfmRmm2022年7月31日星期日 tFpttpttmtmeI1mtmeI1 tdttE02022年7月31日星期日4)带有短路的全混流模型带有短路的全混流模型0v
30、0vV0vIm0vIm01vIm2022年7月31日星期日平推流停留时间分布为:平推流停留时间分布为:tEItEItEmmpm 1其中0vIvmRm tdttEtF0 带有短路的平推流可看作是两个平推流并联,且其中带有短路的平推流可看作是两个平推流并联,且其中一个平均停留时间为一个平均停留时间为0 0。mtmmmetI101mtmeI12022年7月31日星期日0v0Iv0vVfm1Vfm01vI4)有短路的平推流全混流串联2022年7月31日星期日4 组合模型组合模型1)有短路和死区的全混流模型)有短路和死区的全混流模型0v0vVfm死区Vfm12022年7月31日星期日此时此时0vIvfm
31、Rmm停留时间分布为:停留时间分布为:mtmmmeItItE101 0dtttEftm mtmeItF12022年7月31日星期日2)两个不同容积的全混流并联两个不同容积的全混流并联Vf1Vf110v0v01vI011vI2022年7月31日星期日停留时间分布为:停留时间分布为:22111111tteIeItEm0111vIvfR 011211vIvfR tdttEtF0平均停留时间为:平均停留时间为:0dtttEt2111111tteIeI0vvR21111tItI2022年7月31日星期日3)有短路的平推流全混流并联有短路的平推流全混流并联VfmVfm10v0v0vIm0vIp01vIIp
32、m tEs tEm tEp tE2022年7月31日星期日停留时间分布为:停留时间分布为:mtmmpppmeIttItIItE101其中:其中:01vIvftpRmp 0vIvfmRmm2022年7月31日星期日ptt ptt mtmpeII11mtme1 tF平均停留时间为:平均停留时间为:0dttEt001vIvfvIvfmRmpRm2022年7月31日星期日0v0v0vVfm1Vfm0Iv01vI2022年7月31日星期日停留时间分布为:停留时间分布为:tIEtItE串01ptt ptt mpttme10串E2022年7月31日星期日ptt ptt I1mtIexp1 tF01Ivvft
33、Rmp0IvvfRmm2022年7月31日星期日第四节 轴向分散模型1 模型的建立2 应用2022年7月31日星期日第四节第四节 轴向分散模型轴向分散模型1 模型的建立模型的建立假定:假定:流体一恒定的流速流体一恒定的流速u u通过系统,且为一维流动。通过系统,且为一维流动。垂直于流体运动方向的横截面上径向浓度分布均匀,即垂直于流体运动方向的横截面上径向浓度分布均匀,即径向混合达到最大。径向混合达到最大。扩散混合发生在轴向,且可用费克定律加以描述。扩散混合发生在轴向,且可用费克定律加以描述。2022年7月31日星期日根据假设,可建立轴向扩散的数学模型方程,取根据假设,可建立轴向扩散的数学模型方
34、程,取dede的微元的微元管段,对此微元作示踪剂的物料衡算:管段,对此微元作示踪剂的物料衡算:输入输入=流动流动+扩散扩散lAlAAcAu输出输出=lllcDaucAAA积累积累=tclAtlcAAAllcDaucAAA2022年7月31日星期日lcDaucAAAlllcDaucAAAtclAA0tclllcDauclllcDaucAAAAA0tcllcDaucAAA2022年7月31日星期日即:即:022tclcDalcuAAAlculcDatcAAA22轴向扩散模型方程轴向扩散模型方程 将模型化为无因次量:2022年7月31日星期日令令:tt其中:其中:0vvtRLlz 0AAccf 代入
35、则得轴向扩散模型无因次方程为:代入则得轴向扩散模型无因次方程为:zfzfDauLf221AruArLuL F2022年7月31日星期日令:令:DauLPe 皮克列准数皮克列准数 LcDaucPeAA扩散传质速率对流传质速率PePe表示对流流动和扩散传递的相对大小,反映了返混的程度表示对流流动和扩散传递的相对大小,反映了返混的程度 0Pe对流传递对流传递 速率较之扩散传递速率要慢得多,速率较之扩散传递速率要慢得多,全混流全混流。Pe时 0Da为平推流,为平推流,此时扩散传递可略去不计此时扩散传递可略去不计 PePe越大,返混程度越小,越大,返混程度越小,PePe为轴向扩散模型的模型参数。为轴向扩
36、散模型的模型参数。2022年7月31日星期日代入:代入:zfzfPef221Peff平推流平推流 0PePezff122全混流全混流若采用阶跃示踪法,其初始条件为:若采用阶跃示踪法,其初始条件为:1,0,0zzf0,0,0zzf闭闭-闭,开闭,开-开,闭开,闭-开,开开,开-闭闭 对于开对于开-开边界,边界条件为:开边界,边界条件为:边界条件有四种:边界条件有四种:2022年7月31日星期日z,0z,010zf,0代入初始和边界条件,可得代入初始和边界条件,可得的解为:的解为:fFerferf为为误差函数误差函数,其定义为:,其定义为:yxdxeyerf022 不可积的高次误差函数不可积的高次
37、误差函数 121121Peerf2022年7月31日星期日 PePeE41exp21相应的数学期望值和方差的表达式为相应的数学期望值和方差的表达式为:Pe212282PePe对于闭-闭式的边界条件为:2022年7月31日星期日010,0zfPef 0表示刚进入的状态表示刚进入的状态01zf1表示将出口时表示将出口时 代入代入式可得:式可得:122224444expsin81nPeWnPePePeWnPeWnWneF2022年7月31日星期日平均停留时间与方差为:平均停留时间与方差为:1PeePePe12222对于闭对于闭-闭或开闭或开-闭边界,无解析解,通常做法是根据其数闭边界,无解析解,通常
38、做法是根据其数据来获得据来获得 E和 F曲线,计算出相应的数学期望及方差:曲线,计算出相应的数学期望及方差:Pe112232PePe2022年7月31日星期日可见,不管什么边界条件,当可见,不管什么边界条件,当Pe1较小时较小时01.01Pe都会有相同的都会有相同的和2,即:,即:1Pe22具有加和性具有加和性。,即方差均即方差均即若将反应器分成几个小反应器,每个小反应器的平均停即若将反应器分成几个小反应器,每个小反应器的平均停留时间为留时间为i,方差为方差为2i,则有:则有:N.321总2022年7月31日星期日22322212.N总同理,在对一装置进行示踪剂实验时,只要输入示踪剂的同理,在
39、对一装置进行示踪剂实验时,只要输入示踪剂的浓度分布已知,不必一定要是标准的脉冲或阶跃输入,就浓度分布已知,不必一定要是标准的脉冲或阶跃输入,就可以从进出口示踪剂的方差来确定装置内的方差,若以可以从进出口示踪剂的方差来确定装置内的方差,若以2入和和2出分别表示入口和出口处的方差,根据加和性有:分别表示入口和出口处的方差,根据加和性有:Pe2222入出这样就可以计算出模型的参数这样就可以计算出模型的参数2022年7月31日星期日设计反应器时,若停留时间分布未知,还可根据关联式里设计反应器时,若停留时间分布未知,还可根据关联式里估算估算PePe 若为湍流,则若为湍流,则125.0RePe式中:式中:
40、duRe例如:对于空管反应器,在例如:对于空管反应器,在 100023.0,2000Re1Sc的范围内,可由的范围内,可由192ReRe11ScScPe估算。估算。DuDuSc 施密特准数施密特准数2022年7月31日星期日2 应用应用若将轴向分散模型模型用于定常操作的管式操作的管式反若将轴向分散模型模型用于定常操作的管式操作的管式反应器,则有:应器,则有:022AAArdldcudlcdDa边界条件为:边界条件为:时,0l00dldcDaucucAAA时,ll 0LAdldc2022年7月31日星期日对于一级反应,对于一级反应,AAkcr 则可得解析解:则可得解析解:aPeaaPeaaccA
41、A12exp112exp14220其中:其中:2141Peka2022年7月31日星期日当当Pe时时1a将将a展开成:展开成:2)4(81)4(211PekPeka代入上式代入上式)exp(/0kccAA 平推流对一级反应进行计算的结果平推流对一级反应进行计算的结果说明:轴向扩散模型只不过是在平推流模型的基础上迭加一说明:轴向扩散模型只不过是在平推流模型的基础上迭加一 轴向扩散项轴向扩散项2022年7月31日星期日当当0Pe时时将将2/)1(expaPe作级数展开作级数展开)221()1()221()1(4220PeaPeaPeaPeaaccAAPeaaPeaa344k11全混流进行一级反应的
42、计算式全混流进行一级反应的计算式2022年7月31日星期日 可见,具有闭式边界条件的轴向扩散模型,根据模型参数的可见,具有闭式边界条件的轴向扩散模型,根据模型参数的取值不同,可以体现从平推流到全混流之间的任何返混情况。取值不同,可以体现从平推流到全混流之间的任何返混情况。实际反应器的转化率随实际反应器的转化率随Pe倒数的减小而增加。空时越倒数的减小而增加。空时越大,流动状况偏离理想流动的影响也越大。大,流动状况偏离理想流动的影响也越大。当当n1时,难以求出解析解,可用数值法求解。时,难以求出解析解,可用数值法求解。二级反应的转化率受返混的影响比一级反应大。二级反应的转化率受返混的影响比一级反应
43、大。反应级数越高,返混对反应结果的影响越大。反应级数越高,返混对反应结果的影响越大。2022年7月31日星期日第五节 流体混合及其对反应的影响1 混合程度和流体的混合态2 混合态对反应的影响2022年7月31日星期日第五节 流体混合及其对反应的影响1 混合程度和流体的混合态1)调匀度)调匀度S不同组成的流体之间的混合程度不同组成的流体之间的混合程度AvABvB 0ABAAAcvvvc0BBABBcvvvc定义:定义:AAccS 或BBccS 2022年7月31日星期日若完全混合均匀,若完全混合均匀,S S=1.0=1.0,S S偏离偏离1.01.0,表明混合不均匀,表明混合不均匀,S S偏离偏
44、离1.01.0越大,混合就越不均匀。越大,混合就越不均匀。2)2)流体的混合态流体的混合态(离析流)粒子间不发生任何物质交换,离析流)粒子间不发生任何物质交换,相应的混合态称为相应的混合态称为“完全凝集态完全凝集态”,例如,例如气力输送中的气体。气力输送中的气体。宏观流体宏观流体:微观流体:微观流体:达到分子尺度的混合,非凝聚态达到分子尺度的混合,非凝聚态部分微(宏)观流体:部分微(宏)观流体:介于两者之间介于两者之间2022年7月31日星期日2 混合态对反应的影响1)对反应速率的影响)对反应速率的影响设有两个体积相同,浓度分别为设有两个体积相同,浓度分别为1Ac2Ac级不可逆反应级不可逆反应
45、的流体粒子,进行的流体粒子,进行若为宏观流体,则各自的反应速率为:若为宏观流体,则各自的反应速率为:211AAkcr222AAkcr2022年7月31日星期日平均反应速率为:平均反应速率为:2221212121AAAAAcckrrr若为微观流体,则混合后若为微观流体,则混合后A A的浓度为的浓度为2121AAcc 平均反应速率为:平均反应速率为:AAkcr2212AAcck2022年7月31日星期日两种情况下反应速率的相对大小,与两种情况下反应速率的相对大小,与值有关:值有关:1时时 AArr即微观混合降低了反应速率;即微观混合降低了反应速率;1时时 AArr即微观混合提高了反应速率;即微观混
46、合提高了反应速率;1时时 对于一级反应,宏观流体和微观流体的反应效果对于一级反应,宏观流体和微观流体的反应效果是一样的。是一样的。2022年7月31日星期日2)对反应过程(反应器)的影响对反应过程(反应器)的影响 间歇反应器,微观混合的程度不影响反应器的工作情况。间歇反应器,微观混合的程度不影响反应器的工作情况。平推流反应器,微观混合的程度对平推流反应器工作情况平推流反应器,微观混合的程度对平推流反应器工作情况不产生影响。不产生影响。全混流反应器,除一级反应外,微观混合的程度将影响反全混流反应器,除一级反应外,微观混合的程度将影响反应器的工作情况,这种影响随停留时间分布的不同而不同,应器的工作情况,这种影响随停留时间分布的不同而不同,返混程度越严重,微观混合程度的影响越大。返混程度越严重,微观混合程度的影响越大。