平均十分钟,包含了哪些不同经历
把一批饼干同时送进烤箱,十分钟后一起取出,每块饼干的烘烤时间大致相同。若换成连续运行的输送带,途中有的饼干走近路,有的在某处绕了一圈,即使最终算出的平均时间仍是十分钟,也会出现一部分受热不足、另一部分受热过久。流动的聚合物溶液、乳液和熔体也会经历类似差别,只是它们的路线藏在容器内部。对于热敏单体、会继续反应的中间产物或容易随时间降解的材料,了解这些差别,往往比单独报出一个平均数更有用。
停留时间指某一小份物料从规定入口进入,到规定出口离开所经历的时间。为了比较许多份物料,可以想象在同一瞬间给进入的流体做一批记号,然后记录这些记号先后何时到达出口。有些很快出现,有些接近平均时间到达,还有些形成很晚才消失的尾部。把每个时间区间里离开的比例画出来,就得到停留时间分布。入口与出口的位置必须先说清楚;若检测器前还有一段长管,测到的时间就同时包含这段管路,而不只是反应容器内的时间。
先用一个离散例子认识分布。假设有一百份等量示踪物,其中二十份在五分钟左右离开,六十份在十分钟左右离开,另外二十份在十五分钟左右离开。平均值是十分钟,但每一份都刚好停留十分钟的情况,显然与这组结果不同。平均值把每个时间乘上对应的比例,再相加;比例越大的那组,对平均数影响越大。这样的“按份额加权”也用于计算平均粒径和分子量,不过权重可能分别按粒子数量、质量或信号强度定义,使用前需要辨认清楚。
| 符号 | 含义 | 单位 |
|---|---|---|
| t | 从入口标记时刻开始计的经过时间 | min |
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| F(t) | 截至时间t已离开的累计比例 | 1 |
| u;du | 积分中的时间变量及其微小时间宽度 | min |
| dt | 很小的时间间隔 | min |
| ∫;∞ | 积分表示累加小面积,∞表示覆盖所有离开时间 | 积分后为无量纲 |
怎样理解密度乘上时间宽度才是该时间段离开的比例,所以每分钟的密度与分钟相乘后没有单位。第一行把所有时间段的比例加为一;第二行只加到指定时刻t。采用u作为途中变量,是为了区分积分上限t与正在累加的各个时刻。
算一个例子若5至6 min这一段的平均密度是0.08 min⁻¹,则约8%的物料在这一分钟内离开。如果F(5)=0.30、F(6)=0.38,两者之差也给出同样的8%。
使用条件以同一入口标记群体为统计对象,时间非负,物料最终从规定出口被完整计入。普通连续密度的单点高度不是单点概率;理想平推流的集中分布另用δ表示。
连续时间下,不能把曲线某一点的高度直接叫作“恰好在这一时刻离开的概率”。密度曲线的高度要乘上一个小时间宽度,才对应这段时间内离开的比例;整个曲线下的面积加起来为一,表示被统计的物料最终全部离开。纵坐标因此带有时间的倒数单位,例如每分钟。另一方面,累计分布直接回答“到某一时刻为止,已经有多少比例离开”,它从零逐渐走向一。两种曲线描述同一批物料,只是一条看每段新增多少,一条看此前一共积累多少。
平均时间之外,还要量化经历有多分散。可以把每个停留时间减去平均值,得到它偏离平均值的距离;由于正负偏差会抵消,先把偏差平方,再按离开比例求平均,这就是方差。方差的单位是时间的平方,开平方后得到标准差,单位重新变回分钟。刚才五、十、十五分钟的例子,平均值为十分钟,方差为十平方分钟,标准差约为三点一六分钟。这里每一步都有明确的计数含义,不需要把方差当成一个只能由计算器给出的抽象标签。
| 符号 | 含义 | 单位 |
|---|---|---|
| t | 从入口标记时刻开始计的经过时间 | min |
| t̄ | 所考察分布的平均停留时间 | min |
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| σt² | 停留时间方差,平均平方偏差 | min² |
| σt | 停留时间标准差 | min |
| dt;∫;∞ | 小时间宽度、对所有时间的累加及无限上限 | 由被积量确定 |
| √ | 非负平方根,使平方时间单位恢复为时间单位 | 在末行得到min |
怎样理解平均值用各时间段的物料比例作权重。计算方差时,先减平均值,再平方,避免早于平均和晚于平均的偏差相互抵消。最后开平方,把离散程度恢复为可直接与平均时间比较的标准差。
算一个例子20%的物料停留5 min,60%停留10 min,20%停留15 min,平均为10 min。平方偏差加权为0.2×25+0.6×0+0.2×25=10 min²,标准差为√10≈3.16 min。
使用条件积分需要相应的一阶和二阶矩存在。长尾分布可能有有限平均值而无有限方差;有限实验窗口计算的是该窗口内记录数据的统计量。
比较大小不同的设备时,还可以用各自的平均时间作为时间单位。某份物料经历了两倍平均时间,就记为无量纲时间二;经历一半则记为零点五。横轴这样缩放后,密度曲线的高度也必须相应调整,才能让曲线下面积仍为一。方差除以平均时间的平方,则成为相对分散程度。它可以比较两套运行快慢不同的设备,但仍不能代替完整曲线:两个分布即使平均值和方差相同,也可能一个有双峰,另一个有很长的尾部。
| 符号 | 含义 | 单位 |
|---|---|---|
| t | 从入口标记时刻开始计的经过时间 | min |
| t̄ | 所考察分布的平均停留时间 | min |
| θ | 无量纲时间,实际时间除以平均时间 | 1 |
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| E_θ(θ) | 以无量纲时间为横轴的概率密度 | 1 |
| σt² | 有时间单位的方差 | min² |
| σ_θ² | 无量纲时间的方差 | 1 |
怎样理解横轴除以平均时间后,同一小段的宽度变小或变大,纵轴必须乘平均时间,才能保持它代表的面积不变。方差包含时间偏差的平方,因此缩放时除以平均时间的平方。
算一个例子平均时间10 min、标准差5 min时,方差为25 min²,相对方差为25/100=0.25。实际时间20 min对应θ=2;若此处E为0.01 min⁻¹,对应E_θ为0.10。
使用条件平均时间为正且有限;使用方差比值还要求方差有限。无量纲方差并没有普遍小于一的限制。
给流体做记号:从示踪实验得到时间分布
在透明容器里滴入一点颜色,观察颜色何时从出口出现,是示踪实验的直观雏形。真正测量时,记号可以来自可检测的染料、盐、荧光分子或其他合适示踪物,具体选择取决于所研究的流体和检测方法。理想示踪物应该跟随被研究的那一相流动,加入量足够小,不明显改变密度、黏度或流量,也不会在途中反应、吸附或沉积。选择示踪物的任务,就是让“记号的经历”尽可能代表“目标物料的经历”,而不是给设备增加一种新的流动行为。
1953年,丹克沃茨关于连续流动系统的研究,把这些先后离开的经历组织成可以定义、测量和用于计算的分布函数。这个思路的重要之处在于,不必先看清设备内部每一条路径,就能从出口响应检验完全混合和平推流等理想图景。注入一次短脉冲、记录整条出口曲线,能把原来笼统的“混得好不好”转化为具体问题:有多少物料过早离开,平均经历多长,哪些部分停留得特别久。
脉冲法在很短时间内向入口加入已知量示踪物,然后持续测量出口浓度。若体积流率保持不变,某一分钟出口的示踪物浓度乘上这一分钟流出的液体体积,就是这一分钟离开的示踪物量。把所有时间的流出量加起来,可以检查是否回收了投入的总量。再把每个小时间区间的流出量除以总回收量,就得到分布。这样做的核心是物料衡算,曲线归一化只是把实际计量结果换成比例;它并不会自动修复漏测的尾部或示踪物在壁面上的损失。
| 符号 | 含义 | 单位 |
|---|---|---|
| Mrec | 回收到出口的示踪物总质量 | mg |
| q | 恒定出口体积流率 | L/min |
| c(t);c(u) | 扣除背景后的出口示踪物质量浓度 | mg/L |
| t | 从入口标记时刻开始计的经过时间 | min |
| u;dt;du | 积分时间变量与微小时间宽度 | min |
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| ∫;∞ | 从注入后持续累加到完整尾部的积分 | q乘浓度积分得到mg |
怎样理解浓度乘流率是每分钟离开的示踪物质量,再乘时间宽度得到该小段离开的质量。对整条曲线累加得到回收总量,除以它便得到单位时间的离开比例。恒定流率可在分子分母中约掉,因此只用浓度曲线也可归一化。
算一个例子出口曲线面积为20 mg·min/L,流率0.50 L/min,回收示踪物10 mg。某时刻浓度1.0 mg/L,其密度为0.50×1.0/10=0.050 min⁻¹。
使用条件短脉冲、稳定流量,示踪物惰性且代表被研究流体。回收量应与实际注入量比较。若存在未测尾部或不可逆损失,按已回收量归一化得到的是所记录部分的分布,需补足相应实验信息。
测量前应记录背景信号,并在分析时扣除。一个始终存在的小背景,累加很长时间后也可能给曲线尾部带来明显影响。检测器的响应速度、取样管的容积和入口脉冲的宽度同样重要:若注入过程持续了两分钟,出口曲线的宽度就同时包含“入口本来就分散在两分钟内”和“设备造成的时间差”。可以把入口与出口都测下来,利用已知输入响应分析设备,或把注入时间设计得远短于希望分辨的流动时间尺度。
阶跃法换一种输入方式:从某个时刻开始,把入口示踪物浓度由零提高到一个固定值,并维持总流量不变。如果设备最初完全不含这种示踪物,出口浓度会从零逐渐走向入口值。在满足线性示踪条件时,出口与入口浓度之比就是累计分布。它像逐渐被新水替换的水箱,读到零点八意味着与这一阶跃响应相对应的累计比例已达到百分之八十。脉冲法直接给出密度曲线的形状,阶跃法直接给出累计曲线,两者可以相互联系。
| 符号 | 含义 | 单位 |
|---|---|---|
| F(t) | 截至时间t已离开的累计比例 | 1 |
| cout(t) | 阶跃输入后出口的背景校正浓度 | mg/L |
| c₀ | 阶跃后保持不变的入口示踪物浓度 | mg/L |
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| t | 从入口标记时刻开始计的经过时间 | min |
| dF/dt | 累计比例随时间的瞬时变化率 | min⁻¹ |
怎样理解保持入口记号浓度不变后,出口逐渐出现的新记号比例反映累计响应。除以入口浓度把浓度读数变成零到一的比例。累计曲线每分钟增加多少,就是密度,因此密度是累计函数的时间导数。
算一个例子入口阶跃浓度为2.0 mg/L,10 min时出口为1.6 mg/L,则F(10)=0.80。若该时段累计比例每分钟约增加0.03,密度约为0.03 min⁻¹。
使用条件初始无示踪物,输入从零瞬时变到固定值,总流量保持不变,示踪输运满足线性且随时间不变的响应条件。噪声数据求导时需考虑分辨率与平滑方法。
对实验数据,积分和求导的表现也不同。把密度曲线累加为累计曲线,相当于把许多小面积相加,通常会平滑局部噪声;从累计曲线求密度需要比较相邻时刻的微小差异,容易把噪声放大。因此选择输入方式时,既要考虑物理过程,也要考虑检测器分辨率、采样间隔和数据处理。示踪物回收率、流量是否稳定、是否记录到足够长的尾部,应与分布曲线一起保留在实验分析中,才能判断曲线的每一部分反映了什么。
一只釜、许多只釜和一列流体薄片
先看一只体积固定、入口与出口流量相同的理想全混釜。把少量示踪物瞬间加入并混匀,出口每分钟取出的总是釜内混合物的一小部分。起初釜内示踪物多,单位时间流走的量也多;后来剩余示踪物少,流走的量随之下降。这个按剩余量成比例减少的过程,形成指数衰减。物料一进入就与整釜混合,因此模型中可以很快从出口看到少量示踪物,也可以在很多个平均时间之后继续看到少量示踪物。
| 符号 | 含义 | 单位 |
|---|---|---|
| τ | 本模型中等于平均时间的空时 | min |
| V | 全部参与流动与混合的有效体积 | L |
| q | 恒定体积流率 | L/min |
| t̄ | 所考察分布的平均停留时间 | min |
| t | 从入口标记时刻开始计的经过时间 | min |
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| F(t) | 截至时间t已离开的累计比例 | 1 |
| σt² | 停留时间方差 | min² |
| e;exp | 自然指数函数的底,以及指数函数记号 | 1 |
怎样理解充分混合使出口每分钟带走釜内示踪物的固定比例,剩余量因此指数下降。密度的1/τ因子保证总面积为一。对密度积分得到累计曲线,再计算一阶和二阶矩,得到平均时间τ和方差τ²。
算一个例子V=10 L、q=1 L/min时,τ=10 min。10 min时F=1−e⁻¹≈0.632,20 min时约0.865,30 min时约0.950;标准差为10 min。
使用条件恒体积、恒密度、稳态理想全混流,所有有效体积参与交换,示踪物无吸附或反应。这里的稳态指原有流动稳定,示踪物脉冲响应本身随时间改变。
如果平均时间是十分钟,到十分钟为止并不是恰好离开一半,而是约有百分之六十三点二已经离开,还有百分之三十六点八尚未离开。达到百分之九十五累计离开,大约需要三倍平均时间。平均值之所以还能是十分钟,是因为较早离开的许多物料与较晚离开的少量物料共同参与加权。完全混合釜的密度曲线从零时刻开始就有非零高度,是这种理想极限的特征;真实设备还可能有入口到出口的有限输送时间。
把相同总容积分成几只等容釜串联,物料必须先后通过每一只,整段经历不再只有一次随机离开的机会。单只釜的平均时间缩短,总平均时间保持不变,而极短或极长的总经历变得较少。这个过程可以用多个等待时间相加来理解:某一段偶然很短,其他几段仍需要时间;某一段很长,也不意味着每一段都同样长。随着级数增加,整体时间分布向平均时间附近集中,曲线的峰逐渐变窄。
| 符号 | 含义 | 单位 |
|---|---|---|
| N | 等容串联全混釜数量,正整数 | 1 |
| θ | t除以总平均时间所得的无量纲时间 | 1 |
| E_θ | 无量纲时间上的密度 | 1 |
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| t̄ | 所考察分布的平均停留时间 | min |
| t | 从入口标记时刻开始计的经过时间 | min |
| σt²;σ_θ² | 实际时间方差与无量纲方差 | 分别为min²与1 |
| (N−1)! | 从1乘到N−1的阶乘;0!定义为1 | 1 |
| e;exp | 自然指数函数的底,以及指数函数记号 | 1 |
怎样理解总体停留时间由N段全混釜等待时间相加。每段平均为总平均时间的1/N,独立段的方差相加,得到总方差随N减小。第一行给出相应分布形状;当N为1时,它恢复单釜指数分布,θ为零处也按这一极限取值。
算一个例子总平均时间10 min、N=4时,方差100/4=25 min²,标准差5 min。N=2时,在θ=1处E_θ=4e⁻²≈0.541,对应E(10)≈0.0541 min⁻¹。
使用条件恒定总有效体积与流量,各釜等体积、完全混合,釜间没有额外加料、选择性分离或额外滞留。拟合等效釜数时应另说明模型含义。
平均时间不变,分布能变多窄
所有釜完全混合、等容、恒定流率。改变釜数时保持总平均时间不变。密度曲线下的面积代表离开比例,图中总显示时间为平均时间的四倍。
交互图中的总平均时间保持固定,改变的是等容串联釜的数量。单釜曲线较早达到最大值,并拖着长尾;两釜或更多釜的曲线先上升再下降,峰值位置也会移动。曲线变高不表示总物料变多,因为它变窄后面积仍保持为一。对于平均时间十分钟的系统,单釜标准差是十分钟,四只等容釜串联时是五分钟。实际设备若并不完全混合,可以把等效釜数作为某种拟合参数,但它不必等于现场真实设备的数量。
理想平推流是另一种极限:同一时刻进入的所有物料都在同一个时间离开,分布宽度为零。它的累计曲线在此之前保持零,到了规定时间跳到一;密度则需要用“面积为一、宽度趋于零”的集中分布来表示,而不是把普通曲线在某一点画得很高就结束。这个数学对象称为狄拉克δ分布,关键是积分面积,而不是赋予某个普通数值高度。画图时用一根标注面积的箭头,更容易保留这一含义。
| 符号 | 含义 | 单位 |
|---|---|---|
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| t | 从入口标记时刻开始计的经过时间 | min |
| τ | 所有物料共同的停留时间 | min |
| δ(t−τ) | 集中在τ、积分面积为1的狄拉克δ分布 | min⁻¹ |
| F(t) | 截至时间t已离开的累计比例 | 1 |
| t̄ | 所考察分布的平均停留时间 | min |
| σt² | 停留时间方差 | min² |
怎样理解同一批标记物在τ之前没有离开,到了τ全部离开,所以累计曲线发生一步跳变。δ分布保留这次集中离开所代表的单位面积。所有时间都与平均时间相同,偏差为零,方差也为零。
算一个例子τ=10 min时,9.9 min以前累计离开比例为0,达到10 min后为1。它与平均也是10 min的全混釜不同,后者在10 min时累计离开约63.2%。
使用条件理想平推流,所有被追踪物料具有相同路径时间,没有轴向返混或其他时间分散。δ应按积分意义使用,不能当成普通有限高度曲线。
这些模型给出了共同的比较语言。理想全混流的相对方差为一,等容串联釜为釜数的倒数,理想平推流为零。但这几个模型的范围并不是所有流动分布的普遍上限。若一部分物料迅速通过,另一部分长时间滞留,整体相对方差可以大于一;某些理想长尾模型甚至没有有限方差。看到一个超过一的实验结果,应首先检查实际流动、尾部测量与示踪物行为,而不是因为它超出单釜数值就立即判为计算错误。
管道里的快慢车道:层流为什么有长尾
即使管道笔直,里面的物料也未必一起到达终点。黏性液体在圆管内流动时,近壁处受壁面影响较强,中心区域可以流得更快。对于充分发展的牛顿流体层流,速度沿半径形成抛物线分布:中心最快,靠近壁面逐渐降低,在无滑移条件下壁面处速度为零。“牛顿流体”表示剪切应力与剪切速率成正比、黏度在这个描述中为常数;许多高分子熔体会剪切变稀,因此这里先用这个清楚的基准模型认识时间差从何而来。
十九世纪四十年代,泊肃叶对细玻璃管内液体流动的实验,系统改变压力、管长、管径与温度,借助可控制的管路研究流动规律。这类工作与理解细小血管中的流动有关,也为后来讨论黏性输送提供了重要基础。今天用速度剖面重新看这些问题,会发现只测量总流量还不够:同一个流量可以由不同空间位置的速度贡献,而分子从入口到出口花费多久,取决于它在什么位置移动,以及途中是否跨越了不同流线。
| 符号 | 含义 | 单位 |
|---|---|---|
| v(r) | 离轴线r处的轴向速度 | m/min |
| vmax | 管中心轴线处的最大速度 | m/min |
| r;R | 离轴线的径向距离,以及管道内半径 | m |
| v̄ | 按截面积平均的轴向速度 | m/min |
| L | 所研究的充分发展流动段长度 | m |
| τ | 该段体积除以流率的空时,也等于此模型平均停留时间 | min |
怎样理解抛物线速度在中心最大、壁面为零。各个环形区域的面积不同,按面积对速度求平均得到最大速度的一半。体积等于截面积乘长度,流率等于截面积乘平均速度,相除后共同截面积消去,得到空时。
算一个例子L=1 m、平均速度0.10 m/min时,τ=10 min,中心速度0.20 m/min,中心通过时间5 min。在r/R=0.5处,速度是中心的0.75倍,即0.15 m/min。
使用条件直圆管内充分发展的稳态、不可压缩牛顿流体层流,壁面无滑移。这里的时间分布进一步假定不同径向流线间没有交换;非牛顿流体需要相应速度剖面。
在这个模型里,平均速度是中心速度的一半,所以沿中心附近流线通过的物料最早在平均停留时间的一半左右到达。如果总平均时间是十分钟,最早到达约在五分钟;五分钟以前,理想模型没有出口示踪物。接近壁面的物料流得很慢,会在较晚时刻陆续到达,形成长尾。计算某一时刻已经离开的比例时,应按这些流线带走的体积流量加权,不能只计算横截面上对应区域的面积,因为较快区域在同样时间内送出的物料更多。
| 符号 | 含义 | 单位 |
|---|---|---|
| F(t) | 截至时间t已离开的累计比例 | 1 |
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| t | 从入口标记时刻开始计的经过时间 | min |
| τ | 本层流模型的平均时间,等于管段体积除以流率 | min |
| τ/2 | 沿最快中心流线通过的最短时间 | min |
怎样理解给定时间t时,只有速度足以在t内通过管长的流线已经到达。将这些流线的体积流率相加并除以总流率,得到累计比例。对累计曲线求时间导数得到密度。密度中的平方时间除以立方时间,单位仍是时间的倒数。
算一个例子τ=10 min时,5 min前F为零;10 min时F=0.75;20 min时F=0.9375。10 min处E=100/(2×1000)=0.050 min⁻¹。
使用条件牛顿圆管充分发展层流、无径向混合或扩散,入口脉冲对各流线按流量均匀标记,采用出口流量权重。此理想模型平均时间有限,二阶矩和方差发散;真实径向交换会改变尾部。
平均时间十分钟时,二十分钟以前离开的累计比例约为百分之九十三点七五。乍看大部分已经离开,但最后少量物料的尾部对某些统计量仍然很敏感。这个忽略径向交换的理想模型具有有限平均时间,而二阶矩随积分上限继续增长,方差没有有限值。原因是靠近壁面的流线可以提供任意长的停留时间,尽管其流量权重越来越小。现实中的径向扩散、有限入口发展区和实验时间窗会改变最晚部分的表现,因此这也是认识模型假设如何影响尾部的一个例子。
若分子会在流动途中跨越半径方向的位置,它就不再一直待在同一条快车道或慢车道里。分子扩散、湍流涡动或静态混合元件可以改变这种交换,使实际分布与“各流线互不交换”的结果不同。对于非牛顿流体,速度剖面的形状也会改变;对于含粒子的分散体系,粒子还可能出现沉降、迁移或壁面停留。判断一个模型是否合适,要同时考虑流体性质、管道尺度、流动状态和所追踪的对象,而不是只凭“黏度较大”或“管道较长”两个印象。
同样的思路可以帮助理解熔体挤出和换热过程。挤出机内部不只是直圆管,还有螺杆推动、沟槽流动、泄漏和局部回流,不能直接用圆管层流公式代替整个设备。出口示踪曲线仍然可以帮助比较不同转速、进料量或螺杆组合下的流动经历。当某些材料对长时间受热特别敏感时,关注尾部中有多少物料超过某个热暴露时间,会比只比较平均时间更贴近材料变化的原因。
从停留时间推算反应,什么时候还需要混合信息
一旦知道不同物料停留多久,一个自然想法是:分别计算停留五分钟、十分钟、十五分钟的反应结果,再按这些物料的比例平均。对恒温不可逆一级反应,这个思路能够直接给出出口剩余比例,因为每个反应物的存留概率随经历时间按同一指数规律下降。平均之前,先算各个时间对应的指数,再乘分布权重。它一般不等于把平均时间直接代入指数;对一个非线性函数,先平均输入和先计算输出,通常会得到不同结果。
| 符号 | 含义 | 单位 |
|---|---|---|
| CA,out;CA0 | 出口平均A浓度与恒定入口浓度 | mol/L |
| k | 不可逆一级反应的速率常数 | min⁻¹ |
| t | 从入口标记时刻开始计的经过时间 | min |
| E(t) | 出口停留时间的概率密度 | min⁻¹ |
| XA | 出口A转化率 | 1 |
| dt;∫;∞ | 时间小段、对完整分布加权累加及无限积分上限 | 整体剩余比例无量纲 |
| e;exp | 自然指数函数的底,以及指数函数记号 | 1 |
怎样理解某一经历时间对应的一级剩余比例是e的负kt次方。乘该时间段的物料比例,再对所有时间累加,得到平均剩余比例。不同时间的指数值不能通常用平均时间处的一个指数值代替,因为指数函数不是直线。
算一个例子k=0.10 min⁻¹、平均时间10 min。平推流所有时间都为10 min,剩余比例e⁻¹≈0.368,转化率63.2%。全混釜的指数时间分布积分给出剩余0.50,转化率50%。
使用条件恒温、恒密度、线性一级消耗,入口组成恒定,RTD代表反应物所经历的输运。若速率依赖其他浓度、温度沿路径变化,或追踪相与反应相不一致,应补充相应动力学与输运信息。
以平均时间十分钟、一级速率常数每分钟零点一为例,理想平推流的转化率约为百分之六十三点二,理想全混釜则为百分之五十。二者的平均时间一样,前者所有物料都经历十分钟,后者包含许多较早离开和少量较晚离开的物料。较早离开的部分保留了更多原料;已经停留很久的部分再多停一会儿,能额外消耗的原料却越来越少,因此长短经历不能相互完全抵消。这正是分布进入反应预测的一个直观理由。
对于二级或更复杂的反应,还需要知道不同年龄的物料在分子尺度上怎样接触。两份流体都停留十分钟,并不说明它们在这十分钟里一直混在一起:它们可以各自留在小液滴中,也可以很快交换分子而形成均一浓度。完全分隔模型把每个流体小单元看作沿自己年龄演化的小间歇体系,最后才把输出混在一起;充分微观混合的体系则不断让不同年龄的分子接触。两种安排即使具有相同出口时间分布,也可能给出不同的非线性反应结果。
| 符号 | 含义 | 单位 |
|---|---|---|
| C₁;C₂ | 两个等体积单元中A的浓度 | mol/L |
| k | 二级消耗反应的速率常数 | L/(mol·min) |
| r̄_separate | 两个单元分别反应时按总体积平均的瞬时速率 | mol/(L·min) |
| rmixed | 两单元先均匀混合后的瞬时速率 | mol/(L·min) |
怎样理解分隔时,每份液体先按自己的浓度平方反应,再对等体积取平均。混合时,先把浓度平均,再平方计算速率。平方会强调较高浓度,因此两种顺序可以产生不同结果。这是在相同平均浓度下考察微观混合的局部例子。
算一个例子C₁=0、C₂=2 mol/L,k=0.10 L/(mol·min)。分隔时平均速率为0.20 mol/(L·min),先混匀后为0.10 mol/(L·min),二者平均浓度均为1 mol/L。
使用条件等体积、相同温度和二级速率常数,比较当前瞬时消耗,混合不改变总体积。完整出口转化率仍需沿实际混合与反应历程积分。
这个差别可以不用复杂流场就看出来。假设两份等体积液体,一份A浓度为零,另一份为每升二摩尔,合起来的平均浓度是每升一摩尔。若二级反应分别在两份液体内发生,平均的浓度平方是二;若先完全混匀,均匀浓度的平方是一。对同一个二级速率常数,两种瞬时总体消耗速率相差一倍。差异来自分子相遇机会的分布。一级反应对浓度是线性的,在这个比较中先平均和后平均一致,因此时间分布对一级反应具有更直接的预测作用。
示踪实验通常要求示踪物本身不反应,但这不意味着发生化学反应的流体就失去了可描述的停留时间分布。若反应显著改变黏度、密度、相态或流动,应该在能够代表实际工况的条件下测量或建立流动模型;若这些因素保持相近,惰性示踪物仍可以提供有用的流动信息。需要补充的是反应与混合的联系,而不是宣布有反应就不能使用时间分布。对于温度也随路径改变的情况,单独的时间还不足以表示全部经历,需要把温度历程一起计入。
乳液和悬浮体系还要分清连续相与粒子相。水中的溶解示踪物可能准确记录水相的经历,却不一定记录每个聚合物粒子的经历。若粒子在设备中出生、聚并、破碎或被选择性滞留,出口粒子年龄分布与入口注入流体的停留时间分布就需要分别定义。真正有用的预测,应让实验追踪对象与模型中的反应对象对应起来:研究水溶性反应物、聚合物链、胶粒或液滴时,所需要的流动与交换信息并不完全相同。
一条链的年龄,和一份物料的年龄
聚合物给时间分布增加了另一层含义。一份流体在釜里停留十分钟,并不意味着其中每条链都持续生长十分钟。链可能在进入后才被引发,生长很短时间就终止,随后作为不再增长的分子继续随流体移动;也可能在较长时间内保持继续接入单体的能力。对于逐步增长反应,两条已经存在的分子还可以连接成一条更长的链,很难把这条新链简单地看作从某个单体出生时就独立计时。流体年龄、活性链端寿命和分子结构形成的时间,应分别理解。
先看理想线性逐步聚合。反应官能团是分子上能发生连接的局部结构,当两条分子之间发生一次有效的连接,体系中的分子条数就减少一条,组成链的结构单元总数则保持相应守恒。在官能团计量平衡、只形成线性链且没有环化的模型里,可以由连接次数推导数均聚合度与官能团转化率的关系。数均聚合度表示随机选一条分子时,其平均包含多少个计数单元。由于分子会彼此连接,越接近完全反应,剩余分子条数越少,平均每条就越长。
| 符号 | 含义 | 单位 |
|---|---|---|
| P_ℓ | 随机选一条链含ℓ个单元的概率,即按链条数计的比例 | 1 |
| ℓ | 链包含的单元数,正整数 | 1 |
| Xn | 数均聚合度,即每条分子平均包含的计数单元数 | 1 |
| p | 反应官能团的转化率 | 1 |
| Đ | 摩尔质量分散度 | 1 |
| Mn;Mw | 按分子数量平均与按质量加权平均的摩尔质量 | g/mol |
怎样理解第一行把连续存在ℓ−1个连接的概率与随后结束的概率相乘,得到按链数计的长度分布。每次分子间连接使分子条数少一条,因此未反应官能团比例与剩余分子条数相关,得到数均聚合度。随机连接还产生链长分布:按质量加权更强调长链,所以M_w大于M_n。在相同单元质量且端基质量可忽略的理想分布中,二者比值给出1+p。
算一个例子p=0.90时X_n=10,分散度约1.90;p=0.99时X_n=100,分散度约1.99。平均链长显著增加,同时分散度接近2,而不是接近1。
使用条件理想线性逐步聚合,官能团计量平衡、等反应性、无环化与支化。摩尔质量分散度采用结构单元质量相同或等效、端基质量贡献很小的统计近似。该式不代表存在不同年龄或其他反应机制时的通用分布。
弗洛里1936年对线性缩聚产物分子尺寸分布的研究进一步表明,即使官能团反应机会相同、所有物料经历相同操作时间,链长仍然有统计分布。平均增长到一百个单元,并不表示每条都正好一百个。随机连接本身就会产生不同长度,在相应理想模型中,分散度随转化率接近一而趋近二。因此,缩窄停留时间分布能够减少由经历不同带来的差别,却不会自动消除聚合反应自身的链长统计。
对某些活性中心长寿命、引发较同步、链转移和终止较弱的链增长过程,可以建立另一种简化图景:每条活性链以近似相同的平均速率接入单体,生长时间越长,平均链长越大。若单体浓度在所考察阶段近似不变,平均新增单元数便与生长时间成正比。此时活性链不同的生长时长,可以直接转化成额外的链长差异。不过“某一种聚合名称”本身并不能保证中心长寿命,例如配位聚合也可能发生明显链转移;需要根据具体机理与时间尺度判断。
| 符号 | 含义 | 单位 |
|---|---|---|
| ⟨n(tg)⟩ | 生长时间t_g后每条活性链的平均单元数;尖括号表示平均 | 1 |
| n₀ | 开始计时的链所含单元数 | 1 |
| kp | 链增长反应的速率常数 | L/(mol·min) |
| CM | 近似恒定的单体浓度 | mol/L |
| tg | 链保持活性的生长时间 | min |
怎样理解每条活性链的平均接入速率为k_p乘单体浓度,单位是每分钟接入的单元数。速率近似不变时,乘生长时间得到新增单元数,再加上最初单元数。这里的生长时间从链真正开始增长计起,不自动等于流体从设备入口计起的时间。
算一个例子若n₀=1、k_p=2 L/(mol·min)、C_M=0.50 mol/L,平均每分钟增加1个单元。生长10 min后平均为11个,20 min后为21个;这里选用便于理解的设定值。
使用条件单体浓度近似恒定、活性中心长寿命、转移与终止可忽略的简化链增长过程。随机接入仍可产生链间差异;此式描述平均值,不要求每条链严格等长。
常规自由基聚合中,单条活性链的寿命常比设备平均停留时间短得多。链长主要由它活着的那段时间里,增长与终止、转移之间的竞争决定,而不是由这条分子后来在容器里又漂浮了多久决定。但是改变设备停留时间,可能改变单体浓度、引发剂浓度和温度,从而间接改变这种竞争。因此可以说“链长不直接等于流体停留时间乘一个固定增长速度”,却不能据此说自由基聚合产物的分子量在任何条件下都与停留时间无关。
已经终止的链也可能在之后参与新的化学变化,例如向聚合物的链转移会在已有链上产生新的活性位置,进而形成支化。长时间停留提高了某些后续变化发生的机会,但色谱曲线上一个高分子量肩峰并不能单独确定它来自长链支化;不同组分、聚集、校准和其他结构差别也会影响信号。把停留时间分布、动力学信息和结构测量结合起来,才能把“哪一段经历”与“哪一种分子变化”联系得更清楚。
乳胶粒子一起长大,为什么仍有大小差别
乳液聚合把许多反应分散到水中的小聚合物粒子里。可以把每个粒子想象成一间接收单体的微型加工室:单体从周围环境进入,在内部转变成聚合物,粒子体积随之增加。真正的过程还包括粒子成核、自由基进入与离开、单体分配以及稳定剂作用。为了单独看清时间分布带来的影响,先考虑一种明确的种子粒子模型:粒子从入口进入时大小相同,数量不因成核或聚并改变,每个粒子的净体积增长速率在考察阶段保持相同。
这里最容易混淆的是单体进入和聚合物生成。流入粒子的单体可以先作为未反应单体留在其中,也可以被反应消耗;若讨论粒子的总溶胀体积,还需要计入这部分单体的存量与体积变化。下面把模型对象明确规定为新生成聚合物所对应的粒子体积,并假定密度近似不变、单体供应足以维持稳定增长。这样,每分钟增加的聚合物质量除以密度,才对应每分钟增加的体积。先把统计对象说清楚,才可以把它与停留时间直接连接。
| 符号 | 含义 | 单位 |
|---|---|---|
| Vp(t) | 经历时间t后与聚合物质量对应的粒子体积 | μm³ |
| Vp0 | 所有入口种子相同的初始聚合物体积 | μm³ |
| H | 每个粒子恒定的净聚合物体积增长速率 | μm³/min |
| t | 从入口标记时刻开始计的经过时间 | min |
| t̄ | 所考察分布的平均停留时间 | min |
| V̄_p | 按出口粒子数量平均的聚合物体积 | μm³ |
| σV² | 按粒子数量加权的体积方差 | μm⁶ |
| σt² | 这些粒子停留时间的方差 | min² |
怎样理解线性增长把每一分钟变成相同的新增体积。平均时常数初始体积保留,时间变成平均时间。每个体积偏差等于时间偏差乘H,平方后便把时间方差乘H²。方差单位是体积单位的平方,因此写为μm⁶。
算一个例子V_p0=0.004 μm³、H=0.0004 μm³/min、平均时间10 min,平均体积0.008 μm³。单釜时间方差100 min²,对应体积方差0.000016 μm⁶;两只等容釜串联将它减为0.000008 μm⁶。
使用条件相同入口种子,无新增成核、聚并、破碎或选择性滞留,粒子随所用RTD运动,H相同且恒定。体积定义为聚合物质量除以近似恒定密度;包含未反应单体的溶胀体积需额外衡算。
在这个线性增长模型里,停留时间每多一分钟,粒子就多增加同样一份体积。因此,平均粒子体积由平均停留时间决定;体积的波动则把时间波动按增长速率成比例放大。两套设备若平均时间相同,算出的平均粒子体积可以相同,但若一套的时间分布更宽,体积分布也更宽。单只全混釜改成两只等容串联,并保持总有效体积与流量不变,时间方差减半,粒子体积方差也随之减半。这里减少的是方差,标准差减少为原来的约百分之七十点七。
同样的平均体积,不同的体积离散
总平均时间固定为10 min,初始聚合物体积0.004 μm³。种子相同、数量不变、无聚并和新成核,增长速率恒定;累计曲线按粒子数量加权。半径是平均体积对应的等效球半径。
体积不是半径。球形粒子的体积与半径的三次方成正比,所以一个粒子的体积变成八倍,半径才变成两倍。把平均体积换算成一个球的半径,得到的是“具有平均体积的等效球半径”,通常不等于所有粒子半径的算术平均值。对于分布很宽的乳胶,这种差别更值得注意。实验仪器可能按散射强度或质量对不同粒子加权,较大的粒子会获得不同程度的强调;要将模型与测量比较,先应把两者的平均方式和尺寸定义对齐。
| 符号 | 含义 | 单位 |
|---|---|---|
| Vp | 单个球形粒子的体积 | μm³ |
| Rp | 该粒子的半径 | μm |
| π | 圆周率,约3.14159 | 1 |
| (…)^(1/3) | 立方根,寻找其立方等于括号中体积量的长度 | 本式结果为μm |
怎样理解球体体积随三个方向的尺度共同放大,因此与半径立方成比例。反过来从体积求半径,需要开立方。对每个粒子分别换算后再平均,与先平均体积再换算通常不同,因为立方根不是线性函数。
算一个例子体积0.008 μm³对应半径约0.124 μm,即124 nm,直径约248 nm。体积增加到0.064 μm³,变为原来8倍时,半径只变为2倍。
使用条件采用球形几何,尺寸与体积描述同一个粒子边界。若把聚合物干体积换成干球半径,应与实验中的水动力半径或溶胀粒径分别辨认。
同样的反应时间,为什么实验中的粒子还可能大小不一?因为时间分布只是影响之一:种子起始尺寸可能不同,粒子捕获自由基的机会可能不同,单体供应可能随位置变化,体系也可能持续产生新粒子。若两粒子聚并成一个,数量分布会同时改变,简单的“跟随一个入口粒子长大”模型便需要增加相应过程。对于连续乳液聚合,设备组合、加料策略和粒子生成机制相互联系,这也是把流动时间与粒子动力学共同建模的原因。
可以先用交互图回答一个范围清楚的问题:在相同入口种子、相同总平均时间和恒定体积增长速率下,增加串联级数能减少多少由停留差异引起的体积离散?得到这个结果后,再加入实验中确实存在的成核、聚并或增长速率变化,比较它们带来的额外差异。这样的逐步分析让每个模型参数都对应一种机制,也让“更窄的粒径分布”成为可以追溯原因的结果。
把一条出口曲线变成可用的判断
假设一次脉冲实验得到的出口曲线,在入口注入后很快出现一个小峰,随后出现主要峰,最后还有拖得很长的尾部。第一步可以检查入口脉冲和检测管路,确认小峰并非加样或测量造成;随后考虑是否存在较快通道与交换较慢的区域。一个带旁路和滞留区的模型可以产生类似形状,但曲线本身未必能唯一确定内部几何。可以通过改变流量、搅拌或设备连接方式,观察哪一部分随之改变,再结合设备结构判断相应机制。
平均时间与有效体积的关系提供另一条线索。在稳态、恒密度、所追踪物料全部参与流动且示踪物能够代表该相的条件下,平均停留时间等于有效体积除以体积流率。如果几何容积明显更大,可能有一部分空间并未有效参与交换,也可能是示踪物没有进入所有区域或尾部记录不完整。反过来,示踪物在壁面发生吸附会把响应拉长,让测得的时间不再只是流体输送时间。应先检查对应的物理过程,而不是把所有差异都归到“死区”一个词里。
有限采样还会影响计算方式。每个取样点只是某个时刻的浓度,真正需要相加的是它代表的时间区间中的流出量。等间隔采样可以用小梯形面积近似积分;间隔不等时,更要把各自的时间宽度保留下来。随后用这些面积作为权重计算平均时间和方差。尾部信号虽然弱,却乘着较大的时间或时间平方,因此统计量可能对尾部截断很敏感。把“已经回收的示踪物比例随记录时长怎样变化”画出来,常能帮助选择足够的实验持续时间。
对于材料工艺,最终指标应与具体变化相连。若关心未充分反应的部分,可以比较过早离开的比例,并用合适动力学计算它保留多少原料;若关心热降解或颜色变化,可以研究超过某个暴露时间的尾部,必要时加入温度历程;若关心乳胶粒子的离散程度,则结合粒子增长与数量守恒。平均值、方差、分位时间和完整响应曲线各有作用,选择哪一个作为比较指标,应由材料变化的机制来决定。
想一想:平均时间一样,能否认为反应程度一样?
不能只由平均时间得出这一结果。对一级反应,应把各停留时间对应的剩余比例按分布加权。平均十分钟、速率常数0.10 min⁻¹时,全混釜与平推流分别给出50%和约63.2%的转化率。对非线性反应,还需说明不同物料之间的微观混合。可以先在串联釜交互中保持平均时间不变,再观察分布变窄时短停留部分的比例如何变化。
想一想:两只等容釜串联让体积方差减半,半径方差也一定减半吗?
在恒定粒子体积增长模型中,体积对时间是线性关系,因此体积方差直接正比于时间方差。半径却是体积的立方根,是非线性变换。要得到半径方差,需要把整个体积分布变换成半径分布后重新计算,或在分布足够窄时使用说明过条件的近似。相同的平均体积也不自动意味着相同的平均半径。
一条停留时间曲线把设备内部许多看不见的经历压缩成了可以测量的时间信息。理解它时,最重要的几个对应关系是:密度曲线的面积对应物料比例,权重决定平均方式,流动经历与反应机制共同决定产物。沿着这些关系,可以从出口的一次示踪响应,逐步追问原料是否来得及反应、某些物料是否停留过久,以及分子和粒子之间的差别究竟在哪一步产生。