一杯树脂为什么会越放越热
把刚调好的双组分树脂摊成薄薄的一层，与把同样的树脂留在一个较深的杯子里，常常会得到不同的温度经历。薄层有较大的表面与周围交换热量，杯子中心的热却要经过较长的路才能散出去。反应刚开始时，两者都可能没有明显变化；过一会儿，杯内的温度开始上升，固化也随之变快。这里同时发生了两件事：化学反应释放能量，能量又改变反应的快慢。理解这两件事怎样互相推动，就能理解为什么一个在小容器里表现温和的过程，换成大设备后需要重新安排冷却和加料。
温度与热量先要分开。温度描述冷热状态，热量是因温差而传递的能量；同样是升高十摄氏度，一大桶液体需要吸收的能量远多于一小杯。材料吸收能量后温度升高的难易程度由热容量决定。水的比热容较大，意味着相同质量的水升高同样的温度，需要较多能量。反应混合物也有自己的比热容，溶剂、单体、聚合物以及容器壁都可能参与蓄热。温度计给出的是某个位置的温度，不能直接告诉我们整个设备正在产生多少热，或者还有多少化学能尚未释放。
化学反应中的放热也不能理解成“把键打断就会放热”。拆开化学键需要能量，形成新的化学键会释放能量；反应物转变成产物时，全部变化的净结果决定反应吸热还是放热。聚合常把许多小分子连接成长链，某些聚合过程的净放热很明显，另一些则较小，并且会随单体结构、温度和反应路径变化。我们用反应焓表示每发生一定量反应所伴随的热效应，约定放热反应的反应焓为负。为了把“放热功率”写成正数，计算时取反应焓的相反数。
| 符号 | 含义 | 单位 |
|---|---|---|
| Qgen | 反应产生的热功率 | J s⁻¹ = W |
| ΔHr | 按所选反应计量定义的摩尔反应焓，放热时为负 | J mol⁻¹ |
| r | 单位反应混合物体积、单位时间的反应进度速率 | mol L⁻¹ s⁻¹ |
| V | 与速率体积基准一致的反应体积 | L |
怎样理解先用速率乘体积，得到每秒发生多少摩尔反应，再乘每摩尔释放的热。若以单体消耗定义速率，反应焓也按每摩尔该单体换算。负号使放热功率为正。
算一个例子取 ΔH_r = −70 kJ/mol、r = 0.002 mol/(L·s)、V = 10 L，每秒反应 0.020 mol，放热 1.4 kJ/s，即 1.4 kW。
使用条件式中 r 是全体积平均值。存在多个反应时，分别计算各反应的热功率并求和；吸热项按其负的产热贡献保留。
功率中的“每秒”非常关键。假设某一批反应总共释放一百万焦耳，分一个小时释放与集中在十秒内释放，对冷却系统的要求完全不同。前者平均不到三百瓦，后者达到十万瓦；总热量没有增加，峰值功率却相差很大。工业温度控制不仅关心反应最终放出多少热，还关心热量在什么时刻、什么位置释放。自由基聚合的诱导期、自动加速和加料速度，都会改变放热曲线的形状。这也是为什么一张总反应焓表不能代替对反应过程的分析。
现在想象一个混合比较均匀的反应釜。反应向釜内提供热量，夹套把热传给冷却介质，连续流入的较冷物料也需要被加热，出料则带走能量。若液体蒸发，离开的蒸气还能带走汽化所需的能量。把这些过程放进同一笔能量账，剩下的部分才用于改变釜内温度。下面的关系把每项都写成每秒多少焦耳，因此可以直接相加减；它也说明，温度暂时不变只意味着当时收支相等，并不意味着反应停止了。
| 符号 | 含义 | 单位 |
|---|---|---|
| Cth | 有效热容量，表示整体升温 1 K 所需能量 | J K⁻¹ |
| T, Tj, Tf | 分别为均匀釜温、有效夹套温度、进料温度 | K |
| t, dT/dt | 时间与当前温度变化速率；d 表示很小的变化 | s; K s⁻¹ |
| Qgen | 全部反应的净放热功率 | W |
| U | 基于面积 A 的总传热系数 | W m⁻² K⁻¹ |
| A | 有效换热面积 | m² |
| ρ | 液体密度，采用与流量一致的升体积单位 | kg L⁻¹ |
| cp | 定压比热容 | J kg⁻¹ K⁻¹ |
| q | 进出体积流量 | L s⁻¹ |
| Qother | 其余净移热功率，例如蒸发带走的热 | W |
怎样理解温差乘 UA 得到经壁面传出的功率；ρq 是每秒流入的质量，乘比热与进出口温差得到物料显热变化。剩余功率除以热容量，就是每秒升温多少。进料或夹套比釜内热时，相应项自然变成加热。
算一个例子热容量 20 kJ/K，产热 2 kW，壁面移热 1 kW，进出料带走 0.4 kW，其他移热 0.1 kW，剩余 0.5 kW，对应升温 0.025 K/s。若持续相同条件一分钟，局部线性估计升高 1.5 K。
使用条件采用充分混合、近似恒体积与恒物性的液体模型，进出质量近似相等。容器壁只有在能随该温度近似同步变化时才并入热容量；额外加热、搅拌功及相变按需要分别计入。
当放热比移热快时，温度上升；当移热更快时，温度下降。这一步只是能量守恒，还没有回答温度会不会继续偏离。接下来要追问：温度上升一度以后，放热多增加多少，移热又多增加多少？两种增量的比较决定了微小变化是被消除，还是被放大。热稳定性的核心就在这一层反馈关系，而不是某个孤立的“安全温度”数字。
热量离开反应器，要经过哪些路
夹套可以想成包在反应器外面的流动水层，但热量并不会从液体中心直接跳进冷却水。它先从主体物料传到内壁附近，再穿过金属壁，最后传给外侧的冷却介质。如果内壁附着了聚合物沉积层，还要多跨过一道传热障碍。像经过几段相接的窄路一样，整条通道的通行能力受到各段共同影响；只把其中已经很通畅的一段加宽，未必能显著改善全程。这里的“路”是能量传递路径，不是液体一定穿过金属壁。
工程上用总传热系数把这些过程合在一起。对于厚度相对较小、内外面积近似相等的壁面，可以把单位面积上的热阻相加，再取倒数得到总传热系数。热阻越大，在给定温差下通过的热越少。金属导热率高、壁薄，通常有利于传热；不过反应器壁厚还承担压力和结构要求，不能为了散热随意减薄。液体侧的流动、壁面沉积与外侧换热条件，则常是运行中会变化的因素。
| 符号 | 含义 | 单位 |
|---|---|---|
| U | 总传热系数 | W m⁻² K⁻¹ |
| hi, ho | 内侧物料与外侧冷却介质的对流换热系数 | W m⁻² K⁻¹ |
| ℓ_w | 壁厚 | m |
| λw | 壁材导热率 | W m⁻¹ K⁻¹ |
| Rf | 按相同面积基准合计的污垢或沉积层热阻 | m² K W⁻¹ |
| 1/U | 每单位面积的总热阻；倒数表示热阻与传热能力的反向关系 | m² K W⁻¹ |
怎样理解同一稳态热流依次通过各层，各层需要的温差相加，所以热阻相加。金属壁的单位面积导热阻力是厚度除以导热率；它与对流系数的倒数具有相同单位。
算一个例子取 h_i = 500、h_o = 1000 W/(m²·K)，壁厚 0.005 m、导热率 15 W/(m·K)，R_f = 0.0005 m²·K/W。总热阻为 0.003833 m²·K/W，U 约为 261 W/(m²·K)。
使用条件薄壁且内外换热面积近似相等。厚圆筒或不同面积基准需要加入几何面积比，并采用圆筒径向导热阻力。
搅拌的作用不只是把反应物混匀。靠近壁面的液体如果长时间停留，就容易形成温度与主体不同的一层；搅拌让附近液体不断更新，热量更容易被带到壁面。在低黏度液体中，桨叶可能产生较强的循环流动；随着聚合物浓度增加，物料变得黏稠，同样转速未必仍能产生相同的流动。这就把化学变化和传热连接起来：长链增多既会影响反应动力学，也会改变热量离开的能力，两个方向有时会同时促使温度上升。
不同尺寸的设备很难只凭“每分钟转几圈”来比较。小桨叶和大桨叶转速相同，推动的液体范围并不相同；水和浓糖浆密度接近，流动状态却差很多。为整理这些影响，人们把密度、黏度、转速、尺寸、比热容与导热率组合成无量纲数。无量纲意味着单位消掉了，剩下的是两种作用的相对强弱。雷诺数比较惯性与黏性作用，普朗特数比较动量扩散与热扩散，努塞尔数则把对流换热与选定尺度上的导热联系起来。
| 符号 | 含义 | 单位 |
|---|---|---|
| ReM, Pr, Nu | 搅拌雷诺数、普朗特数和努塞尔数 | 1 |
| ρ | 液体密度，此式采用立方米基准 | kg m⁻³ |
| n | 每秒转数，不是以弧度计的角速度 | s⁻¹ |
| d, Dv | 分别为桨叶直径与关联式选用的釜直径 | m |
| η, ηw | 主体温度与壁面温度下的动力黏度 | Pa s |
| cp | 比热容 | J kg⁻¹ K⁻¹ |
| λl | 液体导热率 | W m⁻¹ K⁻¹ |
| hi | 液体侧换热系数 | W m⁻² K⁻¹ |
| B, a, b, c | 在相应设备与流动区间拟合得到的系数和指数；幂表示重复乘积的推广 | 1 |
怎样理解各组中的单位相消，分别比较惯性与黏性、动量扩散与热扩散、实际换热与参考导热。最后一行是实验关联形式，黏度比修正主体与近壁区域的物性差别。
算一个例子ρ = 1000 kg/m³、n = 2 s⁻¹、d = 0.1 m、η = 0.1 Pa·s 时，Re_M = 200。若 c_p = 2000 J/(kg·K)、λ_l = 0.2 W/(m·K)，Pr = 1000；h_i = 100 W/(m²·K)、D_v = 0.3 m 时，Nu = 150。
使用条件关联系数必须对应实际几何和流动范围。非牛顿物料需要匹配剪切条件的黏度或专门关联式；此处算例只计算无量纲组，不用未知系数预测换热。
这类关联式不是普遍适用的自然常数表。桨叶形状、挡板、容器比例、加热还是冷却，以及流动处于什么范围，都会影响系数与指数。聚合物溶液还可能剪切变稀：同一种液体在不同剪切速率下表现出不同的表观黏度。此时需要先说明用哪一种黏度描述搅拌区，而不是把静置时测得的一个数代入所有位置。掌握这些条件，才能明白为什么“增加转速”有时明显提高换热，有时改善很有限，甚至先遇到功率和机械负荷的约束。
冷却介质本身也不是永远保持入口温度的一块冰。它沿夹套流动时会升温，流量不足可能使出口与入口差别很大。把夹套温度近似成一个固定值，适合先建立清楚的能量模型；需要更细致的设备分析时，再把冷却介质的进出口温度与热容量流率列入计算。认识模型中的每个温度究竟对应哪里，比记住一个总传热系数更重要。主体液体、内壁、外壁与冷却水之间的差别，正是热量能够流动的原因。
温度偏离以后，为什么有时回来，有时越走越远
二十世纪初的燃烧研究提出了一个令人困惑的问题：为什么有些反应在同样的环境温度下缓慢进行，改变容器或散热条件后却会突然加速？谢苗诺夫在二十世纪二三十年代发展热爆炸理论，把反应自身产生的热与容器向外散失的热放进同一个问题。他在后来的诺贝尔演讲中回顾了这条思路：自发升温的条件取决于反应动力学、组成和壁面散热，而不是只由物质名称决定。这个认识也适用于许多聚合过程，帮助我们把“突然变快”拆成可以计算的反馈。
先只让温度变化，暂时把反应物浓度等其他条件固定。在许多温度区间内，反应速率常数可以用阿伦尼乌斯关系描述。绝对温度要用开尔文，因为指数里涉及热运动能量；三百开尔文附近升高十开尔文，既是升高十摄氏度，也可能使某个反应的速率明显增加。活化能越高，在同样温度附近对温度变化往往越敏感。这里的“往往”有明确条件：我们正在比较同一个速率表达式，其机理和其他参数在所考察区间内保持适用。
| 符号 | 含义 | 单位 |
|---|---|---|
| k(T), k_* | 温度 T 和参考温度 T_* 下的同一速率常数 | 随速率级数确定；本章一级例子用 s⁻¹ |
| T, T_* | 当前与参考绝对温度 | K |
| Ea | 活化能 | J mol⁻¹ |
| Rg | 气体常数，约 8.314 | J mol⁻¹ K⁻¹ |
| exp, ln | 自然指数和自然对数，互为逆运算；指数须无量纲 | 1 |
| d ln k/dT | 速率常数的相对温度敏感性；严格写法为 ln(k/k_*) 的导数 | K⁻¹ |
怎样理解把两个温度的阿伦尼乌斯式相除，消去前因子，得到第一行。对对数求导得到第二行，表示小幅升温时速率常数的相对变化。
算一个例子若 E_a/R_g = 10000 K，T_* = 350 K、k_* = 0.010 s⁻¹，则 360 K 时 k 约为 0.0221 s⁻¹；350 K 附近的相对敏感性约为每 K 0.0816。
使用条件在选定温区内采用同一机理和近似恒活化能。整个聚合速率还会受到浓度、终止、扩散与相态变化的影响。
把放热功率和移热功率分别画成温度的函数，两条曲线的交点表示热收支相等。站在一个交点上，假想温度略微升高。如果移热增加得更多，额外的热量会被带走，温度便有返回的趋势；如果放热增加得更多，剩余热量又继续升温，最初的偏离就被放大。斜率正是“横向增加一点，纵向变化多少”的数学表达，因此交点附近比较两条曲线的斜率，就能判断这个只有一个温度变量的模型怎样响应小扰动。
| 符号 | 含义 | 单位 |
|---|---|---|
| ε | 相对所研究平衡温度的小偏离，可正可负 | K |
| t | 时间 | s |
| Cth | 热容量 | J K⁻¹ |
| Qgen′, Qrem′ | 在交点处对温度求导的放热、移热曲线斜率，撇号表示导数 | W K⁻¹ |
| τrel | 满足中间不等式时的小扰动衰减时间尺度 | s |
| ≈, > | 近似等于、大于；近似式保留扰动的一次项 | 1 |
怎样理解平衡点上两种功率相等。把温度略偏离后的两种功率各展开为原值加斜率乘偏离，相减后得到第一行。净斜率为负时，正偏离产生降温、负偏离产生升温，均趋向恢复。
算一个例子热容量 20 kJ/K、移热斜率 0.3 kW/K、放热斜率 0.1 kW/K 时，衰减时间为 100 s。在线性近似中，初始 2 K 偏离经过 100 s 减小到约 0.736 K。
使用条件适用于只有温度一个独立动态变量的小扰动模型。若浓度、夹套或控制器也有独立响应，需要分析相应耦合动态。
这一判断还可以变成一幅更容易操纵的图。选一个与温度敏感性有关的温差作刻度，把实际升温除以这个刻度；再用设备自身的冷却时间作时间刻度。单位消去后，模型只留下一个表示初始放热与冷却能力之比的参数。指数曲线代表温度反馈使放热加快，直线代表温差增大使散热增强。这个简化保留了两者竞争的主要形状，同时让不同尺寸和参数的系统可以放在同一张图上比较。
| 符号 | 含义 | 单位 |
|---|---|---|
| Tb | 用于缩放的温差，不是主体温度 | K |
| Rg, Ea | 气体常数和活化能 | J mol⁻¹ K⁻¹; J mol⁻¹ |
| Tj, T | 固定环境或夹套温度与均匀反应温度 | K |
| θ | 相对于环境的无量纲温升 | 1 |
| UA | 总传热导纳，传热系数乘面积 | W K⁻¹ |
| Cth, t | 热容量和实际时间 | J K⁻¹; s |
| s | 无量纲时间；此处字母 s 不表示秒单位 | 1 |
| Q0 | 环境温度下、在固定反应物条件下的放热功率 | W |
| δ | 基准产热与所选温差下散热能力的比值 | 1 |
| e, dθ/ds | 自然常数约 2.718；无量纲温升随无量纲时间的变化率 | 1 |
怎样理解在相对温升较小时把阿伦尼乌斯指数在环境温度附近展开，放热近似为 Q_0 exp(θ)。用 UA T_b 除能量平衡，散热项变成 θ，蓄热项变成对无量纲时间的导数。
算一个例子T_j = 300 K、E_a/R_g = 9000 K 时，T_b = 10 K。若 C_th = 20 kJ/K、UA = 0.2 kW/K、Q_0 = 0.6 kW，则 δ = 0.30，s = 1 对应 100 s，θ = 1 对应温升 10 K。
使用条件均温、固定环境、恒热容量和恒传热能力，反应物消耗在所考察阶段较小。较大升温、相变或明显消耗时，回到未展开的动力学与完整能量平衡。
当放热参数较小时，曲线通常有两个交点。较低温的交点能吸引附近的温度变化，较高温的交点则像一条分界：超过它以后，模型中的放热增长压过散热。随着参数增加，两个交点靠近，最终相切；再增加，就没有可以维持的低温平衡。相切时不仅两条曲线的高度相等，斜率也相等。这两个条件同时成立，给出下面的临界值。它属于当前经过缩放和近似的模型，理解其来历比把这个小数当成所有材料的温度界限更有用。
| 符号 | 含义 | 单位 |
|---|---|---|
| δ, δc | 上一缩放模型的放热参数及其临界值 | 1 |
| θc | 相切点的无量纲温升 | 1 |
| e | 自然指数的底数，约 2.718 | 1 |
| e⁻¹ | e 的倒数 | 1 |
怎样理解第一行要求产热与散热相等，第二行要求两曲线斜率相等。两式左侧相同，因此相切温升为 1，再代回得到参数为 1/e。
算一个例子δ = 0.30 时低温交点 θ 约为 0.4894，高温交点约为 1.7813。δ = 0.40 已超过临界值，此简化模型中不再有交点。
使用条件临界值属于前式定义的单温度近似。它不是温度单位下的数值，也不直接规定不同物料和设备的实际临界条件。
一个小温度偏离，会怎样发展
均匀温度、固定反应物与恒定散热能力的谢苗诺夫近似。上图使用固定坐标，轨迹到 θ = 4 结束；横轴 s 是无量纲时间，不是秒单位。δ 的临界值为 1/e。
在交互图里，先保持初始温升较低，逐渐增大放热参数，观察低温交点如何消失。再选一个仍有两个交点的参数，把初始温度放到较高交点两侧，比较后续轨迹。这样可以看到两种不同的原因：一种是系统条件改变，原来的平衡不再存在；另一种是条件没变，但初始状态跨出了低温平衡能够吸引的范围。曲线快速升高时，实际物料还会经历反应物消耗、沸腾、黏度变化或其他反应，因此后续过程要回到更完整的物料与能量模型继续描述。
连续进料让放热曲线改变了形状
上一部分把反应物浓度暂时固定，适合看清热反馈的起点。连续搅拌釜中，浓度通常会随温度变化：温度升高，反应加快，釜内剩下的反应物反而减少。设想一间不断送入原料的加工室，机器变快后，室内待加工的原料变少；即使机器无限快，总加工量最终也受到每分钟送入多少原料的限制。连续反应器同样存在这样的供料上限，它使稳态放热曲线不再永远沿指数增长。
为把关系算清楚，考虑液体密度和体积近似不变、混合均匀、单一反应物按一级速率消耗的情形。“一级”表示单位体积的消耗速率与该反应物浓度成正比，并不意味着所有聚合都遵循这个规律。平均空间时间等于釜体积除以体积流量。在稳态下，流入、流出与反应消耗恰好相抵，可以先求出每一个温度对应的釜内浓度，再将这个浓度代回放热功率，得到一条沿着稳态物料平衡构造的曲线。
| 符号 | 含义 | 单位 |
|---|---|---|
| τ, V, q | 空间时间、反应体积与体积流量 | s; L; L s⁻¹ |
| Cs(T), Cf | 给定温度下的稳态釜内浓度与进料浓度 | mol L⁻¹ |
| k(T), T | 一级反应速率常数与绝对温度 | s⁻¹; K |
| Qgen,s | 沿稳态物料平衡计算的放热功率 | W |
| ΔHr | 按每摩尔反应物消耗定义的反应焓 | J mol⁻¹ |
| →, ≫ | 趋近、大得多；kτ 是无量纲反应强度 | 1 |
怎样理解稳态时流入减流出等于反应消耗。整理得到浓度，再代入放热式。反应很快时，反应器几乎消耗全部进料，所以产热趋近每秒进料摩尔数乘每摩尔放热。
算一个例子V = 10 L、q = 0.1 L/s、C_f = 1 mol/L、k = 0.01 s⁻¹ 时，τ = 100 s，稳态浓度为 0.5 mol/L。若 −ΔH_r = 200 kJ/mol，放热 10 kW，进料决定的上限为 20 kW。
使用条件恒密度、充分混合的一级不可逆单反应，热效应近似恒定。这里是一条稳态浓度曲线，瞬态浓度另由动态物料平衡决定。
低温时反应慢，进来的原料大多尚未反应就流出，放热较小。中间温度区域，反应对温度比较敏感，转化率迅速增加，放热曲线明显变陡。高温时转化率已经接近一，再加快反应也不能让每秒反应的原料超过每秒进入的量，于是曲线逐渐趋近供料决定的平台。若把冷却和进出料带走的热画成直线，它可能与这条弯曲的放热线相交一次，也可能在合适参数下相交三次。交点数量是模型和参数共同决定的结果，不是每台连续釜固有的三个档位。
范赫尔登在1953年研究自热过程时，系统讨论了反应产生的热怎样维持运行以及怎样出现不同操作状态。所谓自热，是反应释放的能量能够补偿相应热需求，并不表示设备与环境完全隔绝。对工程理解尤其有帮助的是：同一组外部条件可能容许不同的稳态，最终到达哪一个与启动经历有关。由此也能理解点火、熄火与滞后这些术语。点火可以指从低反应速率分支跳向高反应速率分支，未必出现可见火焰；熄火则指反向的状态跳变。
改变夹套温度，通常会平移移热线；改变总传热能力，通常会改变其斜率。改变进料流量则更复杂：流量增加会增强进出料带走的显热，也会缩短停留时间，并改变稳态浓度和最高供料放热。因此不能一边改变流量，一边把原有放热线当成完全不动的背景。把两个作用一起重算，才能判断转化率和温度向哪一边变化。实际设计中常希望把产品要求、可达稳态和可承受扰动范围放在一起考察，而不是一味追求温度最高的交点。
吸热反应提供了一个有益的对照。若温度升高会使吸热反应加快，反应就会多消耗一些热，这通常构成抑制升温的反馈。但一个完整系统还可能有加热控制的延迟、催化剂变化或其他反应。应该根据具体能量与物料关系判断响应，而不能由“吸热”两个字直接推出所有运行状态都稳定。稳定性描述的是某个状态附近的运动趋势，产品是否合格、材料是否降解、压力是否合适，则是需要同时检查的其他性质。
温度和浓度要一起变化，才能看懂过渡过程
稳态图告诉我们可能停在哪里，却没有完整说明怎样走过去。假设连续釜受到一次短暂加热，温度可以立即偏高，但釜内反应物不会同时瞬间跳到新温度对应的稳态浓度。它需要经过反应消耗和流动置换才能改变。于是，在同一个温度下，刚启动的满釜原料与已经运行很久的低浓度物料，可能具有完全不同的放热功率。这一差别解释了为什么沿稳态物料平衡画出的热曲线，不能独自代表全部瞬态运动。
对前面的一级反应模型，可以把浓度和温度各写一条随时间变化的关系。浓度由补入、流出和反应消耗共同决定；温度由进料换热、反应放热和夹套换热共同决定。两条关系中的反应速率常数依赖温度，而温度关系里的放热又依赖浓度，所以必须一起计算。数学上的导数在这里仍然只是“此刻变化得多快”：先由当前状态算出两种变化速度，再沿时间前进，就得到浓度与温度相互影响的轨迹。
| 符号 | 含义 | 单位 |
|---|---|---|
| C, Cf | 釜内与进料反应物浓度 | mol L⁻¹ |
| T, Tf, Tj | 釜内、进料与夹套温度 | K |
| t, τ | 实际时间与空间时间 V/q | s |
| k(T) | 按阿伦尼乌斯关系变化的一级速率常数 | s⁻¹ |
| β | 单位体积消耗单位浓度反应物所对应的温升系数 | K L mol⁻¹ |
| γ | 夹套传热对应的温度响应速率系数 | s⁻¹ |
| ΔHr | 摩尔反应焓 | J mol⁻¹ |
| ρcp | 单位液体体积的热容量 | J L⁻¹ K⁻¹ |
| U, A, V | 总传热系数、换热面积与液体体积 | W m⁻² K⁻¹; m²; L |
| dC/dt, dT/dt | 浓度与温度的瞬时变化率 | mol L⁻¹ s⁻¹; K s⁻¹ |
怎样理解第一式是单位体积物料收支。第二式把能量收支除以液体热容量，所有项都变成温度变化率。温度改变 k，浓度改变产热，因此两个导数必须使用同一时刻的 C 与 T。
算一个例子取 τ = 100 s、β = 200 K·L/mol、γ = 0.01 s⁻¹、C_f = 1 mol/L、T_f = T_j = 300 K。在 T = 350 K、k = 0.01 s⁻¹ 时，若 C = 0.5 mol/L，两导数都为零；若 C = 1 mol/L，温度立刻以 1 K/s 上升，浓度以 0.01 mol/(L·s) 降低。
使用条件恒容、恒物性、单相充分混合液体，固定流量和夹套温度，忽略额外壁体热容量及相变。交互中的数值是这组教学模型参数，不指定真实物料配方。
这种耦合还说明了斜率判据的边界。若浓度能比温度快得多地调整到相应稳态，把它消去只留下一个温度变量，斜率比较就很直接。若两个变量的响应时间相近，则需要检查完整的两变量运动；某些条件下会出现先升后降、衰减振荡，甚至持续振荡。1974年，厄帕尔、雷和普尔对单一一级反应搅拌釜的动态行为进行了分类，展示了温度与浓度轨迹可以比一张稳态交点图丰富得多。这里增加第二个变量，不是为了让模型更复杂，而是为了保留物料真实需要时间才能变化这一事实。
下面的交互使用一组明确给出的教学参数，让反应器体积、流量、热容量和反应速率相互一致。它同时显示稳态的放热与移热关系，以及从选定初始状态出发的温度和浓度轨迹。冷却温度是外部条件，初始温度与初始浓度是起点；改变两类量时，注意分别观察交点是否移动、轨迹是否改变。图中的稳定类型根据同一组浓度—温度方程在交点附近的微小响应计算，而不是只根据图上线条看起来谁更陡来标注。
同一台反应器，不同的起点与冷却条件
模型参数：V = 10 L，q = 0.1 L/s，C_f = 1 mol/L，T_f = 300 K，ρc_p = 1 kJ/(L·K)，−ΔH_r = 200 kJ/mol，UA = 0.1 kW/K。350 K 时 k = 0.01 s⁻¹，E_a/R_g = 10000 K。热与物料方程同时积分；无相变、无反馈控制。
先用较冷、反应物较多的初始状态观察轨迹，再提高初始温度。对于同一冷却条件，两个起点可能到达不同分支；也可能因为开始时储存了较多可反应物料，短时间内出现高于最终稳态的温度峰值。这个峰值尤其值得留意：最终温度合适，不等于整个启动路径都经过相同温度。再把初始浓度降低，比较热峰如何改变，就能把“温度敏感性”和“还有多少反应物可释放热”这两件事分开理解。
实际温度控制又增加了一层动态。传感器测到变化、控制器调整阀门、冷却介质改变流量、热量穿过壁面，都需要时间。加强冷却能力可以改变能量收支，改善测温与控制响应则能更快发现和处理偏离，两者作用不同。如果只看温度设定值，而忽略传热面积、冷却水状态和物料库存，就无法解释同一设定下不同批次为什么产生不同峰值。对过渡过程的理解，最终要把“现在在哪里”“变化有多快”“接下来还有多少热”联系起来。
聚合越黏，为什么有时反而越快
自由基聚合可以从三个动作理解：产生带有活性电子的自由基，使链不断接入单体，以及让活性链失去继续增长的能力。单体是较小的反应分子，活性链随着增长越来越长。两个活性链遇到一起时，可能通过结合或其他终止方式结束增长；这类相遇依赖链在周围物料中的运动。黏度升高以后，长链彼此靠近变困难，并不意味着小单体同样程度地失去运动能力，因此增长和终止的减慢幅度可以很不相同。
1948年，特罗姆斯多夫、克勒与拉加利研究甲基丙烯酸甲酯聚合时，详细讨论了反应在一定转化阶段突然加速的现象。重要的解释是，溶解在单体中的聚合物妨碍增长链相互接近，使终止受到抑制，而单体仍能较快到达活性链端。活性中心于是积累，总增长速率上升。这类现象常叫凝胶效应或自动加速；“凝胶”是沿用的名称，出现它并不要求已经形成由化学交联连接起来的三维网络。
用一个简单的自由基数目平衡可以看清原因。如果单位时间产生的自由基数与双分子终止消耗的自由基数近似相等，自由基浓度由引发速率与终止速率常数共同决定。终止常数减小，在其他量暂时不变时，活性中心浓度就会增加；聚合速率同时取决于活性中心与单体浓度，因此可能上升。关系中的平方根来自终止需要两个自由基相遇，消耗速率与自由基浓度的平方有关。这一解释比“物料变黏所以反应变快”更具体，也能说明为什么后期单体不足或链段运动冻结后，速率还可能再次下降。
| 符号 | 含义 | 单位 |
|---|---|---|
| Ri | 有效自由基生成速率 | mol L⁻¹ s⁻¹ |
| kt | 双分子终止常数，按每次事件消耗两个自由基的约定 | L mol⁻¹ s⁻¹ |
| CR | 活性自由基浓度 | mol L⁻¹ |
| Rp | 单体消耗速率 | mol L⁻¹ s⁻¹ |
| kp | 链增长速率常数 | L mol⁻¹ s⁻¹ |
| CM | 单体浓度 | mol L⁻¹ |
| √, ≈ | 非负平方根、近似相等；平方根来自对二次平衡关系求解 | 1 |
怎样理解准稳态下自由基生成与终止损失近似相等。终止需要一对自由基，所以损失随浓度平方变化；每个活性中心按 k_p C_M 的频率增长，乘活性中心浓度即得总体单体消耗速率。
算一个例子R_i = 2×10⁻⁷ mol/(L·s)、k_t = 10⁷ L/(mol·s) 时，C_R = 10⁻⁷ mol/L。取 k_p = 10³ L/(mol·s)、C_M = 5 mol/L，R_p = 5×10⁻⁴ mol/(L·s)。若只把 k_t 降到原来的百分之一，C_R 和 R_p 都增加十倍。
使用条件均相自由基聚合的简化准稳态模型，采用所写终止常数约定。传质受限、捕获剂、链转移和非稳态活性中心需要相应扩展。
自动加速与热加速可以叠加。终止受限使放热增加，温度上升又改变引发、增长和终止的速率常数；与此同时，黏稠物料可能降低釜内对流换热。要辨认哪种作用先发生，需要比较温度、转化率、黏度和热流随时间的变化。若在严格恒温下仍观察到速率上升，就能把由结构和扩散引起的自动加速，与单纯升温造成的加速区分开。将二者拆开，是建立可迁移的动力学模型的重要一步。
反应开始前还可能有一段诱导期。某些捕获剂会消耗新生成的自由基，使长链暂时难以积累；当捕获能力逐渐耗尽，增长才明显加快。引发剂本身也会随时间分解，所以自由基供应并不一定恒定。下面先把引发剂的一级分解和有效自由基生成写清，再把捕获能力当成一个逐渐减少的库存。一个捕获剂分子可以消耗多少自由基，取决于具体化学过程，用“有效捕获当量”计数可以避免把分子数与可捕获数混在一起。
| 符号 | 含义 | 单位 |
|---|---|---|
| CI(t), CI0 | 当前与初始引发剂浓度 | mol L⁻¹ |
| kd | 一级分解速率常数 | s⁻¹ |
| t | 时间 | s |
| f | 生成自由基中有效引发的比例 | 1 |
| Ri | 有效自由基生成速率 | mol L⁻¹ s⁻¹ |
| e, 2 | 自然指数底数；每个引发剂分解产生两个初生自由基的计量数 | 1 |
怎样理解一级分解意味着每秒分解的比例与剩余引发剂量成正比，积累得到指数衰减。分解事件数乘两个自由基，再乘有效比例，就得到实际进入反应的自由基供应。
算一个例子C_I0 = 0.01 mol/L、k_d = 0.001 s⁻¹、f = 0.5 时，初始 R_i = 10⁻⁵ mol/(L·s)。1000 s 后引发剂剩约 36.8%，R_i 降为约 3.68×10⁻⁶ mol/(L·s)。
使用条件恒温、一级分解，采用每分子产生两个初生自由基的引发剂模型。温度变化时 k_d 随时间变化；其他引发机理使用对应生成规律。
| 符号 | 含义 | 单位 |
|---|---|---|
| B(t), B0 | 当前与初始有效自由基捕获位点浓度，按可捕获自由基当量计 | mol L⁻¹ |
| Ri, u, t | 自由基生成速率、积分用的时间变量、实际时间 | mol L⁻¹ s⁻¹; s; s |
| ∫ | 把各个很短时间内生成的自由基量累加 | 1 |
| CR, CM | 活性自由基浓度与单体浓度 | mol L⁻¹ |
| kB, kp | 有效位点捕获速率常数与链增长速率常数 | L mol⁻¹ s⁻¹ |
| Rp | 单体消耗速率 | mol L⁻¹ s⁻¹ |
| tind | 恒自由基供应近似下的捕获能力耗尽时间 | s |
怎样理解捕获主导时，每产生一个有效自由基近似消耗一个捕获当量。准稳态活性中心由供应与位点捕获相平衡得到。若供应恒定，库存除以消耗速度就是耗尽时间。
算一个例子初始捕获容量 B_0 = 0.001 mol/L，近似恒定 R_i = 10⁻⁵ mol/(L·s)，估计诱导时间为 100 s。50 s 时剩余容量为 0.0005 mol/L；在其他参数相同的捕获主导近似内，活性中心浓度是初始值的两倍。
使用条件有效位点独立捕获、早期增长消耗较小且捕获主导自由基损失。接近耗尽时要加入双分子终止；真实多步骤阻聚体系需要按化学机理定义捕获容量与速率。
当库存接近用尽时，简单的捕获主导公式会预言速率急剧增加。这提示我们应切换到同时包含双分子终止、增长和浓度变化的模型，而不是把分母趋近零理解为真实反应速率无限大。氧与不同阻聚体系的作用还依赖具体化学机制，不能把所有诱导期都归结为“氧逸出所需时间”。从温度曲线看，一段安静期可能是在积累随后的反应条件；把这段时间解释清楚，才知道为什么同样总反应热可以对应很不同的放热峰。
水中的小颗粒，怎样改变反应与传热
乳液聚合把许多反应地点分散在水相中的聚合物小颗粒内。一个直观图景是：较大的单体液滴充当储库，少量单体先溶入水相，再进入正在生长的聚合物颗粒；自由基的进入、离开和终止决定各颗粒有多少活性链。水相有助于带走热量，较低的整体黏度也有利于搅拌。反应场所却并不是所有可见液滴，而主要是满足成核与增长条件的聚合物颗粒。分清储库、运输通道和反应位置，才能理解其速率为何与均相聚合不同。
史密斯和尤尔特在1948年研究乳液聚合动力学时，把活性链分隔在许多小反应场所中，讨论了每个场所平均自由基数不同的情形。某些条件下，一个颗粒常在“没有活性链”和“有一条活性链”之间转换，平均占有数约为二分之一。这个平均值并不是每个颗粒真的装着半条链，而是对许多颗粒或较长时间的统计。隔离使自由基相遇的方式改变，因而仅用总体平均自由基浓度代入均相终止式，可能漏掉关键的颗粒尺度结构。
| 符号 | 含义 | 单位 |
|---|---|---|
| Rp,mix | 按总体反应混合物体积计算的单体消耗速率 | mol L⁻¹ s⁻¹ |
| kp | 颗粒内链增长速率常数 | L mol⁻¹ s⁻¹ |
| CM,p | 按颗粒内部体积定义的单体浓度 | mol L⁻¹ |
| Np | 每升总反应混合物中的反应颗粒数 | L⁻¹ |
| n̄ | 每颗粒平均活性自由基数 | 1 |
| NA | 阿伏伽德罗常数，约 6.022×10²³，用于把个数换成摩尔 | mol⁻¹ |
怎样理解k_p C_M,p 给出单个活性中心每秒接入单体的平均频率。乘每颗粒活性中心数，再乘单位总体积颗粒数，得到分子计数速率，除以阿伏伽德罗常数转换为摩尔。
算一个例子k_p = 1000 L/(mol·s)、C_M,p = 5 mol/L、N_p = 10¹⁷ L⁻¹、n̄ = 0.5 时，R_p,mix 约为 4.15×10⁻⁴ mol/(L·s)。若 −ΔH_r = 70 kJ/mol，对应每升混合物约 29.1 W 产热。
使用条件用代表性颗粒内浓度与平均占有数描述的分散体系。颗粒间差异很大时按群体加和；n̄ = 0.5 是特定动力学条件，不是所有乳液聚合的固定值。
典型间歇乳液聚合常分成成核、近似恒速增长和单体逐渐耗尽三个阶段。在成核阶段，反应颗粒数增加，放热可以随之增加；在一段中间阶段，颗粒数、颗粒内单体浓度以及平均自由基数近似稳定，总速率便可能出现平台。这里的平台并不能单独证明反应由单体扩散控制。即使单体补给足够快，只要决定速率的几个量基本不变，也会得到近似恒速。判断控制步骤，需要比较补给与消耗的能力，并观察改变搅拌、粒径或单体条件后的响应。
若需要估计一个球形颗粒从周围取得单体的能力，可以先采用稳态扩散模型：远处浓度较高，表面浓度较低，单体沿浓度差向颗粒运动。对孤立球体求解扩散方程后，总进入量与扩散系数、球半径和浓度差成正比。注意这里得到的是每个颗粒每秒进入多少摩尔，不是单位面积的通量；两者相差一个表面积。将进入量乘以每摩尔反应热，才能得到供料受限时该颗粒能够产生的热功率。
| 符号 | 含义 | 单位 |
|---|---|---|
| j | 单个球形颗粒获得单体的总摩尔速率 | mol s⁻¹ |
| π | 圆周率，约 3.1416 | 1 |
| D | 外部介质中的单体扩散系数 | m² s⁻¹ |
| rp | 颗粒半径 | m |
| Cb, Cs | 远处主体和颗粒表面的单体浓度，此处均按立方米计 | mol m⁻³ |
| Qp,gen, Qp,rem | 单颗粒产热与向外移热的功率 | W |
| ΔHr | 每摩尔单体消耗的反应焓 | J mol⁻¹ |
| h | 颗粒外表面的换热系数 | W m⁻² K⁻¹ |
| Ts, Tb | 颗粒表面与外部主体温度 | K |
怎样理解球对称稳态扩散给出第一行，它已经对整个表面积积分。第二行假设到达的单体按所述反应消耗。第三行用球面积乘换热系数与温差。
算一个例子D = 10⁻⁹ m²/s、r_p = 10⁻⁷ m、浓度差 1 mol/m³ 时，j 约为 1.26×10⁻¹⁵ mol/s。取每摩尔放热 70 kJ，产热约 8.80×10⁻¹¹ W；h = 10⁶ W/(m²·K) 时，平衡这份功率需要约 0.00070 K 的表面温差。
使用条件孤立球体、无限远近似、恒 D、外部稳态纯扩散。真实流动、多颗粒浓度场、孔内扩散和内部温度梯度需要对应模型；示例中的 h 是独立给定的参数。
同样的思路适用于多孔固体催化剂，但还需要区分外部液膜与颗粒内部孔道。反应物可以先跨过外部边界层，再扩散进孔内与活性位点接触，热量则沿相反方向传出。颗粒表面温度可能高于主体流体，内部还可能存在温度梯度，因此釜内温度计读数平稳并不自动意味着每个微小位置温度相同。把反应速率写成每单位颗粒体积还是每单位反应器体积，也必须始终一致，否则颗粒数量或体积会被重复计算。
在液体中，扩散系数不仅随绝对温度变化，还常常强烈依赖黏度与分子尺度。升温会改变液体黏度，聚合又会改变介质组成，因此不能普遍把扩散控制的放热线画成严格直线。随着单体液滴消失，颗粒内单体可能逐渐下降，速率和热流随之降低；同时，颗粒内部的黏度和玻璃化也会改变增长与终止。水相分散提供了有利的热管理条件，但完整理解仍要把颗粒中的反应、相间运输与整个反应器的冷却放在同一幅图中。
把小试放大，先问还会释放多少热
一个反应在某时刻的放热功率不大，可能因为反应确实接近结束，也可能因为反应暂时受到温度、引发或混合的限制，而系统内仍保留许多可反应物料。区分这两种情况，需要估计剩余反应潜力。绝热温升就是把剩余可释放的热，除以系统吸收热量所需的热容量，回答“如果这部分热暂时没有散出去，温度会增加多少”。它把化学库存与温度尺度联系起来，比只看当前温度更能帮助理解下一阶段可能发生的变化。
| 符号 | 含义 | 单位 |
|---|---|---|
| ΔTad | 剩余反应潜力对应的绝热温升 | K |
| ΔHr | 按所选反应物摩尔基准定义的反应焓 | J mol⁻¹ |
| nremain | 能够继续按该反应消耗的剩余反应物量 | mol |
| Cth | 所计入物料与同步蓄热部分的总热容量 | J K⁻¹ |
| ≈ | 热容量和反应热取常数时的近似 | 1 |
怎样理解剩余摩尔数乘每摩尔放热得到可释放的总能量，再除以每升高 1 K 所需能量。分子、分母的焦耳相消，结果是温差。
算一个例子剩余可反应物为 2 mol，每摩尔放热 70 kJ，有效热容量 10 kJ/K，则绝热温升约 14 K。若在相同反应物库存下热容量增加到 20 kJ/K，该尺度降为 7 K。
使用条件所述反应完成、没有向外散热，采用近似恒反应热和热容量。可逆平衡、其他反应、沸腾、相变或明显变物性时需重新计算能量分配。
绝热温升的计算也说明稀释剂和水为什么会改变温度经历。它们增加需要被加热的质量和热容量，某些情况下还会改变反应物浓度与速率。总放热相同而热容量变大，温升就变小；但材料配方、相态与后续分离也会随之变化，不能只从热容量一个角度决定加入多少液体。若跨越很大温度范围，比热容变化、蒸发和相变需要一起计入。一个用恒定比热得到的数值，应读作明确假设下的能量尺度，而不是所有阶段都会精确达到的终点。
几何放大还会改变产热与散热的比例。形状相似的容器线性尺寸扩大两倍，体积扩大八倍，表面积只扩大四倍。若单位体积反应速率保持相同，产热增加八倍；若单位面积传热能力和温差也保持相同，散热只增加四倍。于是每单位体积可用的换热面积变小。小杯中的温度控制经验之所以不能简单按质量成比例搬进大釜，正是因为体积与面积按不同的幂次增长。
| 符号 | 含义 | 单位 |
|---|---|---|
| V, A, R | 理想球形容器的体积、换热表面积与半径 | m³; m²; m |
| π | 圆周率 | 1 |
| A/V | 每单位体积对应的表面积 | m⁻¹ |
| τth | 无反应、线性换热下的温度松弛时间 | s |
| ρ, cp | 密度和比热容，此处密度按立方米计 | kg m⁻³; J kg⁻¹ K⁻¹ |
| U | 总传热系数 | W m⁻² K⁻¹ |
怎样理解球的面积随半径平方、体积随立方增长，二者相除后留下半径的倒数。热容量除以传热导纳给出冷却响应时间；在物性与 U 相同的相似球中，该时间随半径增大。
算一个例子半径从 0.1 m 增加到 0.2 m，A/V 从 30 m⁻¹ 降到 15 m⁻¹。取 ρc_p = 2×10⁶ J/(m³·K)、U = 200 W/(m²·K)，无反应冷却时间从约 333 s 增到 667 s。
使用条件球体用于展示几何规律；其他相似形状也具有面积平方与体积立方的缩放。实际放大中 U、搅拌、有效面积和反应速率还会变化，需分别计算。
解决放大问题可以沿着这两条尺度关系展开。增加内部盘管、外循环换热或有效换热表面，是增加能量离开的通道；分批或连续地加入单体，是改变反应物库存与放热时间分布；调整溶剂、颗粒分散和搅拌，是同时影响热容量、反应与运输。每一种选择都应与产品要求联系起来。例如过慢加料可能改变共聚组成或链结构，过强剪切可能影响分散颗粒，增加外循环又会改变停留时间与流动经历。热管理不是在既定反应之外简单附加一台冷水机，而是制造过程的一部分。
量热测量提供了把这些问题联系起来的办法。反应量热可以在控制温度的同时记录为维持温度需要移走多少热，据此得到放热功率随时间的变化；绝热或近绝热量热则关注体系在较少散热条件下怎样自行升温。结合单体转化率、进料记录和物料热容量，能够把已经释放的热与仍可能释放的热区分开。对聚合物来说，还应把热流变化与黏度、分子量分布或颗粒结构对应起来，才能判断温度变化是在改变产量、改变材料，还是两者同时发生。
回到最初那杯树脂，现在可以沿一条完整的因果链解释它：反应形成产物并释放净能量，有限的散热使温度上升，温度和黏度又改变反应及传热，剩余物料决定还能继续释放多少热。把薄层改成深杯，改变的是几何和散热；改变初始温度，改变的是起点与速率；改变加料方式，改变的是库存与时间分布。能把这些作用分别辨认，再重新放回同一张能量账中，就能读懂热稳定性图，也能理解它为什么直接影响聚合材料的制造与性能。