第四章 非理想流动模型 答案
化学反应：第四章 非理想流动
4.1 反应器中的返混现象与停留时间分布
第四章 非理想流动
实际的工业装置
在实际的工业装置中由于物料在反应器内的流 动速率不均匀、或因内部构件的影响造成物料出现 与主体流动方向相反的逆向流动、死角等都会导致 偏离理想流动。
对于所有偏离平推流和全混流的流动模式 统称为非理想流动。
本章将利用停留时间分布定量地对非理想流动 进行讨论，并考察这些非理想流动对反应器性能的 影响。
4.1 反应器中的返混现象与停留时间分布
第四章 非理想流动
假如示踪剂改用红色流体，连续检测出口中红色
流体的浓度，如果将观测的时间间隔缩到非常小，
得到的将是一条连续的停留时间分布曲线。
4.1 反应器中的返混现象与停留时间分布
(1) 停留时间分布密度函数的定义
第四章 非理想流动
定义：在稳定连续流动系统中，同时进入反应器的N 个流体粒子中，其停留时间介于t～t+dt的那部分粒 子dN占总粒子数N的分率记作：
第四章 非理想流动
1. 停留时间分布密度函数E(t)
实验： 在连续反应器内，如果在某一瞬间 (t = 0) 极快地 向入口物流中加入 100 个红色粒子，同时在系 统的出口处记下不同时间间隔流出的红色粒子数， 结果如下表。
4.1 反应器中的返混现象与停留时间分布
第四章 非理想流动
以时间 t 为横坐标，出口流中红色粒子数为 纵坐标，将上表作图如下：
流体的流动速率和方向带有一定的随机性。反应器内的流动状态实际是 随机变化的。
根据概率理论，我们可以借用两种概率分布以定量地描绘物料在流动系统中的 停留时间分布，这两种概率分布就是停留时间分布密度函数 E (t)和停留时间 分布函数 F (t)。
4.1 反应器中的返混现象与停留时间分布
第四章 非理想流动反应器设计
t2
1
最后由模型参数计算非理想反应器的平均结果。
第四章 非理想流动反应器设计
4.1 流体物料粒子的停留时间分布函数和分布密度函数 4.2 停留时间分布（RTD）的数字特征及无量纲化
4.3 理想流动模型
4.4 非理想流动模型 4.5 非理想流动反应器设计 4.6 混合质量对反应的影响
第四章 非理想流动反应器设计
第四章 非理想流动反应器设计
（3）计算分布函数F 根据F (t) 与 E(t) 关系： F (t )
E (t )dt
0
t
（4-4）
当上式用于实验数据计算时，应转换成下列形式：
F (t ) E (t )t
由此就可以计算出不同时间下的 F(t) 值，例如 当 t=5 min 时， F (5) tE (t ) 1 (0.05 0.1 0.16 0.2 0.16) 0.64 当 t=12 min 时，
计算举例：
n0~ 2 F (2) F (0) 0 0 0 N n F (4 ~ 6) 4~ 6 F (6) F (4) 0.38 0.08 0.30 N F (0 ~ 2)
同理可计算其它各时间段的分布函数值，见上表。
第四章 非理想流动反应器设计
（2）停留时间分布密度函数E(t) 定义：停留时间分布函数 F(t) 在某时间段 t→t+dt 内的平
用一定的方法将示踪剂加到反应器进口，然后在反应器出 口物料中检验示踪剂信号，以获得示踪剂在反应器中停留
时间分布的实验数据。
1、不与主流体反应； 2、物理性质相近；
选择示踪剂的原则
3、有别于主流体的可测性； 4、多相检测不发生相转移； 5、易于转变为光、电信号。
化学反应工程第四章答案
4-1 在定态操作反应器的进口物料中脉冲注入示踪物料。
出口处示踪物浓度随时间变化的情况如下。
假设在该过程中物料的密度不发生变化，试求物料的平均停留时间与=100⎰∞=-=-=02222971.8187.625.47)(t dt t E t tσmin 24-2 无量纲方差表达式的推导 （1）推导无量纲方差222/ttσσθ=；（2）推导CSTR 的22tt=σ。
1. τθt=2. ττtet E -=1)(证明： 4-3 设()θF 及()θE 分别为闭式流动反应器的停留时间分布函数及停留时间分布密度函数，θ为对此停留时间。
（1）若该反应器为平推流反应器，试求①F(1); ②E(1);③F(0.8);④E(0.8);⑤F(1.2) （2）若该反应器为全混流反应器，试求①F(1); ②E(1);③F(0.8);④E(0.8);⑤F(1.2) （3）若该反应器为非理想流动反应器，试求 ①F(∞); ②F(0);③E(∞);④E(0);⑤⎰∞0)(θθd E ；⑥⎰∞)(θθθd E解1平推流模型 2 全混流θθ-=e E )( ， θθ--=e F 1)(3非理想流动模型a 多釜串联 θθθN N N e N N E ---=1)!1()(， 0)(C C F N =θ4-4 C(t)t/min4-18图用阶跃法测定某一闭式流动反应器的停留时间分布，得到离开反应器的示踪剂浓度与时间的关系，如图4-18所示。
试求 （1）该反应器的停留时间分布函数)(θF 及分布密度函数)(θE ；（2）数学期望θ 及方差2θσ；（3）若用多釜串联模型来模拟该反应器，则模型参数是多少？ （4）若用轴向扩散模型来模拟该反应器，则模型参数是多少？（5）若在此反应器内进行1级不可逆反应，反应速率常数1min 1-=k ，且无副反应，试求反应器出口转化率。
解（1）.()()θF c t c t F ==0)(⎪⎩⎪⎨⎧-=12)(t t c 3322〉≤≤≤t t t()⎪⎩⎪⎨⎧-==∴120)(0t c t c t F 3322〉≤≤≤t t t ，()θF t F =)(⎪⎩⎪⎨⎧=∆∆=010)(0t c c t E3322〉≤≤≤t t t ，()t t E E =)(θ（2）.1==-ttθ（3）. 多釜串联模型 (4). 轴向扩散模型 试差 Pe=0.001 (5).4-5. 为了测定某一闭式流动反应器的停留时间分布，采用脉冲输入法，反应（1）反应物料在该反应器中的平均停留时间t 及方差2θσ（2）停留时间小于4.0min 的物料所占的分率。
第四章非理想流动反应器(201204)
Vc(t )dt c(t ) F (t ) Vc()dt c()
降 阶 法
(不讲)
c(0)
c(0) c(t)
c0(t) 0
t=0 输入曲线
0 t
响应曲线
t
Vc(t )dt Vc(0) dt
停留时间大于 t 的示踪剂量 示踪剂输入量
t→t +dt
1-F(t)
c( t ) F (t ) 1 c ( 0)
t E (t )t t
2 t 2 0 2
• 取样为等时间间隔时： (t t ) E (t ) t t E (t ) t t E (t )t E (t )t
2 2 2 t
2
t E (t ) t E (t )
2
2
(4)无因次化对比时间及函数
(4) 脉冲法和阶跃法的比较
脉冲法 阶跃法
示踪剂 注入方 法
在原有的流股中加入 示踪剂，不改变原流 股流量
可直接测得
将原有流股换成流量与 其相同的示踪剂流股
dF (t ) dC (t ) E（t ) dt C0 dt
E(t)
F(t)
C (t )dt F (t ) C (t )dt C (t )dt
0
E（t )dt 1
• 停留时间分布函数F(t):停留时间0-t范围内的物料（停留 时间小于t的质点）占进料的分率。 t 有： t=0, F(t)=0, t=∞, F(t)=1,F(t)是单调增函数
F (t ) E（t )dt
0
在某一时间t时， E(t)和F(t)之间的关系为： dF (t ) t E（t ) F (t ) E（t )dt dt 0
第四章 非理想流动模型 答案
0
t20 10tic 2c ii tm 2171 2 1 2 6 .4 61.4272.5275736.374256.31480(s2)
0
无因次方差
θ2tm 2t2217428.5 020.01990.02 ( 接近 PFR )
求平均转化率
（1） 根据一级反应转化率的本征动力学公式： 在该反应器出口处的平均转化率为
F ( t ) 5 0 0 0 .0 0 2 e x p ( 0 .0 0 2 t ) d t e x p ( 0 .0 0 2 t) 5 0 0
2 0 0
2 0 0
11 e 0 .4 e 0 .6 7 0 0 .3 6 8 0 .3 0 2 3 0 .2 %
E(t)
0 .0 0 2 4 0 .0 0 2 2 0 .0 0 2 0 0 .0 0 1 8 0 .0 0 1 6 0 .0 0 1 4 0 .0 0 1 2 0 .0 0 1 0 0 .0 0 0 8 0 .0 0 0 6 0 .0 0 0 4 0 .0 0 0 2 0 .0 0 0 0
第四章 非理想流动模型 答案
② 停留时间在0 ～ tm 的物料所占总料量的百分数为
F (tm )0 tm E (t)d t 0 .0 0 2 0 tm e x p ( 0 .0 0 2 t) d t ( tm= 500 s )
5 0 0 e x p ( 0 .0 0 2 t) d ( 0 .0 0 2 t) e x p ( 0 .0 0 2 t)5 0 0
0
1 0
c i 0 0 0 0 .1 0 .5 1 0 .0 8 .0 4 .0 0 0 2 2 .6
i 1
1 0
t ic i 0 0 0 1 4 .4 9 6 .0 2 4 0 0 2 3 0 4 1 3 4 4 0 0 6 1 5 8 .4
第四章. 非理想流动[1]
![第四章. 非理想流动[1]](https://img.taocdn.com/s3/m/ec0f2f87e53a580216fcfe1c.png)
4.1.2 概率函数的定义与性质
4.1.2.1 RTD密度函数E（t）的定义与性质 • （1）定义：流体以稳定流量进入设备而 不发生化学变化时，任意瞬间（记为t=0） 进入设备的数量为N的流体微元中，从出 口流出的停留时间介于tt+dt的流体微 N d 元数占总数的分率 N ，用统计规律的概 率方法表示时为概率密度与随机变量的 N d 变化值的乘积E（t）dt，即 N =E（t） dt。为时间的函数，称RTD密度函数。
E t dt 1
t 0
t 0
N
即
0
E t dt 1
（4-1）式
15
②E（t）t的图形
图中直方阴影的面 积dS=E(t)dt,密度随 时间的变化曲线与 横坐标围城的面积 是无数个小直方面 积的加和，即
S dS E t dt 1
i
i
13
对离散型数据：（Ni/N）=E（t）t，或：
N i N E t t
因此上述实验（t1=t2==2）得：
N i N
0 100
2 100 6 100 12 100
18 100
22 100
17 100 17 200
12 100
6 100
4 100 4 200
4
（4）. 三种理想反应器RTD的定性说明 • BSTR.间歇操作，无流型可言，故不存在 停留时间分布问题； • PFR. 所有流体微元停留时间相同，集中 在同一时间点=VR/V0，就是说PFR的 RTD是均一的或RTD范围集中在一点。 • CSTR.流体微元停留时间从0应有尽 有，RTD范围最宽。
24
1024 第4章 非理想流动1(1)
E(t)
平行流动
表明反应器内有两股平行 的流体存在
t
接近平推流的 E 曲线形状
10
反 应 工 程 第 四 章 非 理 想 流 动 ·
接近全混流的几种E曲线形状
11
反 应 工 程 第 四 章 非 理 想 流 动 ·
接近全混流的曲线，也同样有着正常形状、出峰太早、 内循环、出峰太晚、由于仪表滞后而造成时间的推迟等。 总之，测定停留时间分布曲线的目的，在于对反应器内 的流动状况作出定性的判断，以确定是否符合工艺要求 或提出相应的改善方案。求取数学期望和方差，以作为 返混的量度，进而求取模型参数N。
1、多级混合流模型（多级全混流串联模型） 反 应 工 程
c0, v0 c1 v0
1 2
c2 v0
cN－1 v0
N
cN v0
第 模型的适用范围：多级混合流模型是描述的轻微偏离理想 四 流动平推流的模型。 当N多个全混流反应器串联，N越大，返混越小，当 章 N→∞时，返混→0，总体相当于平推流。 非 由此：把不同多级全混流反应器串联，来描述不同的返混 理 程度。 想 把实际的反应器，看做若干个等体积的全混流反应区串联， 流 实际反应器，相当多少个等体积的全混流反应区串联，来 动 直观的反映偏离平推流的程度。 ·
c AN N t N 1 N t / t E( t ) ( ) e c0 ( N 1 )! t t t t N ti ,ti N
N
4
反 应 工 程 第 四 章 非 理 想 流 动 ·
无因次化令
t t
上式简化为：
NN E( ) t E( t ) N 1e N ( N 1 )! N 1 max N
因此，对于第N个全混区
第四章非理想流动作业1(1)
形成非理想流动的根源 ① 操作条件：温度、压力、流量、物料组成 ② 流体性质：粘度、重度、扩散系数等 ③ 设备型式及结构
病态流动是指：设计、制造不良，即反应设备缺陷 引起的非理想流动，如死角、短路等
2. 什么是E(t) 、F(t)、 、σt2、σθ2 ？ E(t)dt 和F(t) 物理意义？ 答： E(t) 是停留时间密度分布函数：定常态下的连
4. 示踪剂有何要求（选取原则）？ 答：采用何种示踪剂，要根据物料的物态、
相系及反应器的类型等情况而定。
1.不与主流体发生反应（无化学反应活性）。 2.与所研究的流体完全互溶，除了显著区别
于主流体的某一可检测性质外，二者应具 有尽可能相同的物理性质
便于检测：应具有或易于转变为电信号 或光信号的特点，且浓度很低时也能 检测。
第四章
第四组 组长：钟美弟 组员：韩平、杨丹、 马梦华、郭慧春、曾召名、 曾凡林、土孙江、包生璞
1.非理想流动的起因是什么？其根源有哪些？ 什么是病态流动?
答：起因有两个
1.设备内不均匀的速度分布。包括：死角、沟流、 短路、层流流动、截面突变引起的收缩膨胀等。
2.与物料主体流动方向相反的流动。例如：管式反 应器：扩散、局部循环流动、压差、流体与固体 间的摩擦等引起；釜式反应器：搅拌引起流体循 环运动等。
或：流过反应器的物料中停留时间小于 t (介 于0 ~ t 之间)的质点所占的分率。
σt2停留时间分布分散程度的量度。数学上指 对于平均停留时间的二次矩。
σθ2 是无因此散度 E(t)dt：表示同时进入反应器的N个流体质点
中，停留时间介于t 与t+dt 间的质点所占分 率dN/N。
F(t)：表示出口流体中停留时间小于t 的物料 （ 0~t 范围内的质点）占进料的分率。
化学反应工程第四章习题答案
化学反应工程第四章习题答案work Information Technology Company.2020YEAR第四章 非理想流动1.停留时间分布的密度函数在t ＜0时，E （t ）=_______。
（0） 2.停留时间分布的密度函数在t ≥0时，E （t ）_______。
（＞0） 3.当t=0时，停留时间分布函数F （t ）=_______。
（0） 4.当t=∞时，停留时间分布函数F （t ）=_______。
（1） 5.停留时间分布的密度函数E （θ）=_______E （t ）。
（t ）6.表示停留时间分布的分散程度的量=2θσ_______2tσ。
（21t ）7.反应器物料的停留时间的分布曲线是通过物理示踪法来测定的，根据示踪剂的输入方式不同分为_______、_______、_______、_______。
（脉冲法、阶跃法、周期示踪法、随机输入示踪法） 8.平推流管式反应器t t =时，E （t ）=_______。
（∞） 9.平推流管式反应器t t ≠时，E （t ）=_______。
（0） 10.平推流管式反应器t t ≥时，F （t ）=_______。
（1） 11.平推流管式反应器t ＜t 时，F （t ）=_______。
（0）12.平推流管式反应器其E （θ）曲线的方差=2θσ_______。
（0） 13.平推流管式反应器其E （t ）曲线的方差=2t σ_______。
（0） 14.全混流反应器t=0时E （t ）=_______。
（tte t -1）15.全混流反应器其E （θ）曲线的方差=2θσ_______。
（1） 16.全混流反应器其E （t ）曲线的方差=2t σ_______。
（2t ） 17.偏离全混流、平推流这两种理想流动的非理想流动，E （θ）曲线的方差2θσ为_______。
（0～1）18.当流体在半径为R 的管内作层流流动时，在径向存在流速分布，轴心处的流速以0u 记，则距轴心处距离为r 的流速=r u _______。
化学反应工程 第四章 非理想流动
今用分散模型关联，求
数。
化学反应工程
4.2.1 常见的几种流动模型
解：
换算为无量纲时标，
则得下表数据。
将实验数据标绘成曲线，然后读取
等间隔时的诸E值
见下表。
化学反应工程
4.2.1 常见的几反应工程
4.2.1 常见的几种流动模型
③化学反应的计算 定态情况下平推流管式反应器的物料衡算式为：
流， ；对一般实际流况， 。
；对平推
所以，用
来评价分布的分散程度比较方便。
化学反应工程
4.1.4 用对比时间θ表示的概率函数
例4-1 今有某一均相反应器中测定的下列一组数据（见 ，示踪加入 下表第一栏和第二栏），实验采用
量Q=4.95g，实验完毕时测得反应器内存料量V=1785mL，求 解：
（详见教材P92）
对定态系统的非理想流动，同样可作微元段的物料衡算而得：
若用无量纲参数表示并注意到：
这样式（4-32）便变为：
化学反应工程
4.2.1 常见的几种流动模型
对一级反应可得解析解：
对于二级反应，用数值法求得的结果，表示在图（4-17）
和图（4-18）中。
化学反应工程
4.2.1 常见的几种流动模型
（４）组合模型
化学反应工程
4.1.1 非理想流动与停留时间分布
在一个稳定的连续流动系统中，当在某一瞬间同时进 入系统的一定量流体，其中各流体粒子将经历不同的停留 时间后依次自系统中流出。如果把函数 用曲线表示，
则图4-2（a）中所示阴影部分的面积值也就是停留时间介 于t和t＋dt之间的流体分率。
化学反应工程
4.1.1 非理想流动与停留时间分布
化学反应工程
化学反应工程课程第四章非理想流动模型详解
48
用多釜串联模型来模拟一个实际反应器的步骤
①.测定该反应器的停留时间分布； ②.求出该分布的方差； ③.将方差代入式(4-28)求模型参数N； ④.从第一釜开始，逐釜计算。
采用上述方法来估计模型参数N的值时，可能 出现N为非整数的情况，用四舍五入的办法圆整 成整数是一个粗略的近似处理方法，精确些的办 法是把小数部分视作一个体积较小的反应器。
tE(t)dt
t
0
E(t)dt
tE(t)dt t
0
0
2 t
(t t )2 E(t)dt
0
E(t)dt
(t t )2 E(t)dt 0
0
0
2
0
1 1
34
2. 全混流模型
考察有效体积为Vr、进料体积流量为Q的全混流 反应器，若在某一瞬间t=0，将流体切换成流量相 同的含有示踪剂的流体，同时检测流出物料中示踪 剂浓度变化。
含示踪剂流体Q
C0
流体 Q
C0
切换
系统 V
C0
示踪剂检测
Q
t=0
t
0
阶跃法测定停留时间分布
t
21
在切换成含示踪剂的流体后，t-dt～t时间间
隔内示踪剂流出系统量为vc(t)dt ，这部分示踪剂
在系统内的停留时间必定小于或等于t，任意的dt
时间间隔内流入系统的示踪剂量为vc0dt ，由F(t)
定义可得
i0
25
2.方差(对均值的二次矩)
散度：停留时间分布分散程度的量度
2 t
(t t )2 E(t)dt
0
E(t)dt
(t t )2 E(t)dt
0
n
《化学反应工程》第三版(陈甘堂著)课后习题答案
《化学反应工程》第三版(陈甘堂著)课后习题答案第二章均相反应动力学基础2－4三级气相反应2NO+O22NO2，在30℃及1kgf/cm2下反应，已知反应速率常数2kC=2.65×104L2/(mol2 s)，若以rA=kppApB表示，反应速率常数kp应为何值？解：原速率方程rA=dcA2cB=2.65×104cAdt由气体状态方程有cA=代入式（1）2－5考虑反应A课所以kp=2.65×104×(0.08477×303) 3=1.564后当压力单位为kgf/cm2时，R=0.08477，T=303K。
答p p 2rA=2.65×10 A B =2.65×104(RT) 3pApBRT RTp表示的动力学方程。
解：.因，wwnAp=A，微分得RTVdaw案24网pAp，cB=BRTRT3P，其动力学方程为( rA)=dnAn=kA。
试推导：在恒容下以总压VdtVδA=3 1=21dnA1dpA=VdtRTdt代入原动力学方程整理得wdpA=kpAdt设初始原料为纯A，yA0=1，总量为n0=nA0。
反应过程中总摩尔数根据膨胀因子定义δA=n n0nA0 nA若侵犯了您的版权利益，敬请来信通知我们！Y http://.cn.co（1）mol/[L s (kgf/cm2) 3]m（1）则nA=nA01(n n0)δA1(P P0)δA（2）恒容下上式可转换为pA=P0所以将式（2）和式（3）代入式（1）整理得2－6在700℃及3kgf/cm2恒压下发生下列反应：C4H10发生变化，试求下列各项的变化速率。
（1）乙烯分压；（2）H2的物质的量，mol；（3）丁烷的摩尔分数。
解：P=3kgf/cm2，（1）课MC4H10=58，（2）w.krC2H4=2( rC4H10)=2×2.4=4.8kgf/(cm2 s)PC4H10=PyC4H101 dpC4H10= P dt2.4-1==0.8 s 3w（3）nC4H10=nyC4H10=n0(1+δC4H10yC4H10,0xC4H10)yC4H10dnH2dtdnH2dt=hdaw后n0=nC4H10,0=δC4H10rC4H10=反应开始时，系统中含C4H*****kg，当反应完成50%时，丁烷分压以2.4kgf/(cm2 s)的速率dyC4H10dt答1rCH=2.4224wdnC4H10dt案116×1000=2000mol582+1 1==21网dyC4H10=n0(1+δC4H10yC4H10,0xC4H10) dt=2000×(1+2×1×0.5)×0.8=3200 mol/s若侵犯了您的版权利益，敬请来信通知我们！Y http://.cno2C2H4+H2，dP=k[(δA+1)P0 P]=k(3P0 P)dtm（3）dpA1dP= dtδAdt2－9反应APS，( r1)=k1cA , ( r2)=k2cp，已知t=0时，cA=cA0 ，cp0=cS0=0, k1/k2=0.2。
第4章 非理想流动
对于升阶跃 如果t时刻出口物料中A的浓度为 c A
cA F t c A0
对于降阶法：
cA 1 F t c A0
在 c A ~ t 图中应满足：
2019年4月1日星期一
AC tdcA c A0 tdF t c A0 tEt dt c A0 t
1 ) 脉冲示踪法
t=0示踪剂
主流体 v0
v0
2019年4月1日星期一
c0 t
c t
t0
t
响应曲线
t
2019年4月1日星期一
设示踪剂加入量为m（根据A的物料衡算）
m v0c A dt v0 c A dt
0
0
CA v0C Adt m dt mE(t )dt 0 C Adt
2019年4月1日星期一
区别：
寿命分布指的是系统出口处的流体粒子的停留时间，
年龄分布是对系统中的流体粒子而言的停留时间。 实际测定得到的且应用价值又较大的是寿命分布。 通常所说的停留时间分布指的是寿命分布。
2019年4月1日星期一
1 停留时间分布的定量描述
1）停留时间分布函数
出 口 中 红 色 粒 子 数
2019年4月1日星期一
E t t 5
t
4 0.0492 81.3333
F t E t dt
0
F t t 5 E t
0
5
77 1 c dt A dt 81.3333 81.3333 0.9467 0
5
t tEt dt
0
tc dt
A 0
81.3333
1017 第4章 非理想流动1(1)
填 的，非常必要对实际的流型进行逼近模拟。 料 塔
第 四 章 非 理 想 流 动 ·
所谓模型法，就是通过对复杂的实际过 程的分析，进行合理的简化，然后用一 定的数学方法予以描述，使其符合实际 过程的规律性，此即所谓的数学模型， 然后加以求解。
20
反 应 工 程
1 )方程组的庞大，如裂解反应。多个反应同时发生，因
E( t )
,( t t ) 0 ,( t t )
第 四 章 非 理 想 流 动 ·
E( t )dt 1.0
0
E(t)是典型的
t2
f ( t )E( t )dt
t1
Dirac 函数，该函数有如下性质： f ( t ),t1 t t2
0 , t t1 ,t或 t t1 ,t2 2
提出了一些简化了的定性模型。经这样的简化，边界条件
就可以确定了，当然，在简化模型时所引入的参数，还要 靠实验测定。有了定性模型，然后再结合有关的数学方程， 从而使复杂的过程得以简化，并能予以定量计算。在模型 化中，最重要的是合理的“简化”。把复杂的实际过程简 化为较简单清晰的物理图形，即具有用数学方式能加以表 达的物理模型，然后再变为数学模型。建立数学模型就是 在准确性与简单性之间，寻求妥善的解决方案。
4.模型参数少，便于测定 简化模型中都含有模型参数。模型参数是简化模型偏离 真实过程的归并结果，都要通过实验确定，所以模型参 数越少越好，而且要便于测定。
26
反 应 工 程
数学模拟方法是化学反应工程中主要的研究方法，是可 行的，但并不是一种万能方法。 对于物料的流动状况，可以用流动模型描述。平推流和
均有
非理想流动模型
• 例3-3某非理想流动反应器，其停留时间分 布规律同例3-2。在该反应器内进行一级反 应，动力学方程为-rA=3.33×10-3cA，请确 定该反应器的出口转化率（反应物A的化学 计量系数为1）。
• 解：采用凝集流模型进行计算。
• 对于一级反应，在间歇反应器中转化率 与反应时间关系如下：
t
xA
1 exp
0
kt c
c
1
0.3347
0.6653
多级串联槽模型
• 物理模型： • 反应器是由若干大小相等的全混流反应器
串联而成。 • 这些全混流反应器之间没有返混，没有反
应。 • 定常态操作。
平推流与全混流的比较
1 (-rA )
平推流
= VR
V0
=cA0
dx xAf
A
0 (-rA )
NVRi
V0
t
则：
由Vci1
Vci
VRi
dci dt
得：
dci
d
Nci
Nci1
BC : 0t 0, ci 0
对第一槽 解得：
dc1
d
Nc1
Nc0
c1 c0 1 eN
对第二槽
dc2
d
Nc2
Nc1
Nc0
1
e N
解得： c2 c0 1 eN 1 N
对第三槽
dc3
d
t 1 2 1 2 1.164 Pe 12.186
t 374.4 321.65
s
1.164 1.164
• 计算β值
1 4k 1 4 3.33103 321.65 1.163
Pe
12.186
• 计算出口转化率
第四章非理想流动习题
1 / 1
第四章非理想流动习题
1．脉冲示踪法测得如下数据：
求E 函数和F 函数
2．有一个反应器测得其离散准数τ
0.2D
uL
=。
对此反应进行示踪实验时符合闭式容器条件。
今在同样条件下用多釜串联模型描述其模型参数N 为多少？
3．有一个长12m 的管子，其中左端1m ，装有2mm 的固体颗粒，中间9m ，装有1cm 的固体颗粒，右端2m ，填有4mm 的固体颗粒。
假设空隙率相同
2D
udp
=，流体通过的时间为2min ，试估计出口方差：【提示：1()D udp =2()D udp =3()2D udp =，
()D uLi =()D
udpi ()dpi Li
，22t ti σσ=∑（有加合性）
】
求：（1）平均停留时间和方差；
（2）若在该反应器中进行一级不可逆反应A R →，速率常数0.045k =min -1 求A 的平均转化率？
求：（1）平均停留时间和无因次方差； （2）若用于一级反应，且该反应在单个全混釜中相同的τ下限定组分A 的转化率为0.82，在此真实反应器中转化率为多少？
（3）若进料量为20L/s,该容器的体积多大？
请预测一下相当于几个全混釜的串联效果。
停留时间分布
V v0
化
学
工
程
系
死区 模型：
e
( v ) vm
化
学
工
程
系
平推和全 混串联组 和模型
化
学
工
程
系
4.停留时间分布曲线的定性应用
1.出峰太早 E(t)
短 路
沟流
t
原因：反应器存在沟流、短路现象，使出峰提前。
化
学
工
程
系
化
学
工
程
系
化
学
工
程
系
2.出现多峰，且递降 E（t）
tˆ
原因：反应器内有循环流
dCA J A EZ dz
程
系
uAC
C Ez A l
dl
C uA(C dl ) l C Ez A (C + dl ) l l
化
2
学
工
程
系
C C C = Ez -u 2 t l l
将方程无因次化：
C t l Pe表示对流流动和扩散传递程度的相对大小 ; C= ,θ = ,Z = C0 tˆ L 其数值越大，轴向返混程度越小。
化
学
工
程
系
化
学
工
程
系
对于开开边界
1 1 2 2( ) 8( ) Pe Pe
2
1 ˆ =1+2 Pe
对于闭-闭边界：
2 2 2 (1 e Pe ) Pe Pe
2
ˆ 1
对于开闭（闭开）边界：
1 1 2 3 Pe Pe
塔式反应器内流体流动的模拟。
04-1第四章 非理想流动
10～ 11 4 0.04
11～ 12～ 12 14 1 0.01 0 0
0 0
2
6
12
18
22
17
12
6 0.06
子数
△N∕N
0.02 0.06 0.12 0.18 0.22 0.17 0.12
100粒子：0~10s 95个 5~6s 18个
N
0
10
N
95 95% 100
N 18 18% N 100
第四章 非理想流动
前面介绍了：
这两种反应器，在相同的操作条件下，两者的 结果有很大的差别，究其原因是反应物料在反应 器内的流动状况不同，亦即停留时间不同。
PFR：所有流体质点在反应器内的停留时间
都相等，严格按照先进先出的规律循序而
进，相邻两截面之间没有混合（返混为0）。 CSTR：刚进入的新鲜物料立即和釜内原有 物料充分混合均匀，各物料颗粒在器内具 有一定的停留时间分布（返混为最大） 。 实际反应器：凡是偏离平推流和全混流的所 有流动状态均称为非理想流动。
dF ( ) dF (t ) dF (t ) E ( ) t tE (t ) d dt d t t
2 ( 1) 2 E ( )d 1 E 0 0
t t
(t t ) 2 t
2
0
E (t )dt
t2
t
2
10
0~10s
F (t )
0
dN
N
N
0
10
N
95 95% 100
5~6s
dN N 18 E (t )dt 18% N N 100
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1 1 1 0.135 0.865 1 0.0075 (272.5 / 50)
可见方法（1）与方法（2）得到相同的结果。
（3） 选用PFR流动模型，平均停留时间tm 为 272.5 s
xA 1 exp(ktm ) 1 exp(0.0075 272.5) 1 0.13 0.870
1 exp(0.002t ) 500 s 0.002 0
② 停留时间在0 ～ tm 的物料所占总料量的百分数为
F (tm ) E (t )dt 0.002 exp(0.002t )dt
0 0
tm
tm
( tm= 500 s )
500
500
0
exp(0.002t )d(0.002t ) exp(0.002t ) 0
10
0 0 0
0 0 0
0 0 0
0.004 14.4 2073.6
0.027 96.0 18432
0.469 2400 576000
0.823 2304 663552
1.00 1344 451584
1.00 0 0
1.00 0 0
c
i 1
i
0 0 0 0.1 0.5 10.0 8.0 4.0 0 0 22.6
10
i
C
i 1
2 tm
i
将相关数据列于下表
t/s
c (g / m3)
0
0 0
48
0 0
96
0 0
144
0.1 0.1
192
0.5 0.6
240
10.0 10.6
288
8.0 18.6
336
4.0 22.6
384
0 22.6
432
0 22.6
c
i 1
10
i
F( ti ) t ic i t i2 c i
θ2
t2
2 tm
1480 0.0199 0.02 ( 接近 PFR ) 2 272.5
求平均转化率
（1） 根据一级反应转化率的本征动力学公式： 在该反应器出口处的平均转化率为
xA 1 exp(kt )
xA
0
xA E (t )dt
0
E (t )dt
1 exp(kt )E (t )dt E (t )dt
200
0.002exp(0.002t )dt exp(0.002t ) 200
1 1 0.4 0.670 0.368 0.302 30.2% e e
0.0024 0.0022 0.0020 0.0018
63.2%
E(t)
0.0016 0.0014 0.0012 0.0010 0.0008 0.0006 0.0004 0.0002 0.0000 0 200 400 600 800 1000 1200 1400 1600 1800 2000
30.2%
E(t)=0.002exp(-0.002t)
tm
t (s)
4-3 t/s c (g / m3)
用脉冲法测得某反应器的示踪剂浓度为 0 0 48 0 96 0 144 0.1 192 0.5 240 10.0 288 8.0 336 4.0 384 0 432 0
计算平均停留时间。在反应器中进行一级反应 A → P, k = 7.5×10-3 s-1, 求平均转 化率。若此反应分别在CSTR和PFR中进行，具有相同的平均停留时间，反应结果 如何？
t c
i 1 10 i 1
10
i i
0 0 0 14.4 96.0 2400 2304 1344 0 0 6158.4
2 t i ci 0 2073.6 18432 576000 663552 451584 0 1711641.6
解： ①
tm
由求平均停留时间公式
0 0
tE (t )dt E (t )dt
tE (t )dt 0.002t exp(0.002t )dt
0 0
t exp(0.002t ) 0 exp(0.002t )dt
0
t exp(0.002t ) 0
i 1 i
10
i
c
i 1
10
i
192 0.5
240 10.0
288 8.0
336 4.0
384 0
432 0
xi ci
10
0
0
0
0.066
0.382
8.347
7.077
3.678
0
0
x c
i 1
A i
0 0 0 0.066 0.382 8.347 7.077 3.678 0 0 19.55
0 0
因为 E (t )
V0 c(t ), M
将其代入后，得
xA
1 exp(kt ) c(t )dt c(t )dt
0 0
对于离散型的实验数据，上式又写成
1 exp(kti ) ci
t c
0 0 48 0 96 0 144 0.1
xA
1 exp(kt ) c
解： 根据计算离散型停留时间分布函数公式
F (t )
Ci ti
1 10 i 1 i i
i
C t C
i 1
Ci
1 10 i
i
; tm
tiCi ti
i 1 10
10
C t
i 1 i
tiCi
i 1 10
10
i
C
i 1
; t2
2 t i Ci i 1 10
平1 exp(kt ) c
i 1 10 i i
10
10
c
i 1
i 1
i
i
i
c
i 1
10
19.55 0.865 22.6
i
（2）选用 m-CSTR 非理想流动模型，模型参数 m 为
m
1
θ2
m
1 50 0.02
m
1 1 xA 1 1 1 k 1 k ( t / m ) m
第四章 非理想流动模型
4-1 有一反应器的停留时间密度函数为
E (t ) 0.002exp(0.002t )
（时间单位为s）
试求：① 平均停留时间（即数学期望）tm是多少？
② 停留时间在0 ～ tm 的物料所占总料量的百分数？ ③ 停留时间在 200s 到 500s 的物料所占总料量的百分数？
（4） 选用CSTR流动模型，平均停留时间tm为272.5 s
1 xA 1 1 k
1 1 0.671 1 0.0075 272.5
结论： 1. 根据实验数据计算出的方差为0.02，流型很接近于PFR；
2. 无论采用何种求算反应器出口平均转化率，均有相同的结果， 出口平均转化率为0.865； 3. 相同反应和相同平均停留时间， 在PFR 的出口转化率为0.870，
e1 1 0.632 63.2%
③ 停留时间在200 s ～ 500 s的物料所占总料量的百分数为
F (t ) F (500) F (200) F (t )
500
500
0
E (t )dt
200
0
E (t )dt
500
500
200
E (t )dt
平均停留时间
tm
方差
t c c
0 0 10 i
10
i i
6158.4 272.5 (s) 22.6
t2
t c c
0 0 10 i
10
2 i i
2 tm
1711641.4 272.52 75736.3 74256.3 1480 (s2 ) 22.6
无因次方差
说明该反应器的确接近于PFR；
4. 相同反应和相同平均停留时间， 在CSTR的出口转化率为0.671。