《Separations》：Modelling Adsorption Breakthrough Curves
Xin Shen and
Jules Thibault
编辑推荐：
为优化吸附柱设计，研究者系统评估了五种扩散模型（瞬时/非瞬时吸附、Fickian扩散、异常扩散）预测丁醇在活性炭上吸附突破曲线的性能。研究表明，模型表现强烈依赖于吸附剂类型，分数阶扩散模型在特定情况下提供了更精确的描述。本研究为复杂吸附体系模型的合理选择与优化提供了关键指导，具有重要的工程应用价值。
论文解读
吸附作为一种重要的分离技术，在空气分离、水净化、废水处理等领域扮演着关键角色。它的魅力在于能耗低，但设计一个高效的吸附柱却是一门大学问。想象一下，一股含有目标物质的流体流过装满多孔吸附剂颗粒的柱子，目标物质如何被捕获，何时会“突破”而出，直接决定了吸附柱的尺寸、运行时间和处理效率。描述这种动态行为的核心工具就是“突破曲线”——它描绘了吸附柱出口浓度随时间的变化。然而，如何准确预测这条曲线，是困扰工程师和研究人员的一大难题。传统上，我们依赖基于Fick定律的模型，但对于活性炭等结构复杂的多孔材料，物质在其内部的扩散可能并不遵循经典规律，表现出所谓的“异常扩散”。这时，是否需要引入更复杂的分数阶扩散模型？新模型带来的预测精度提升，能否抵得上其增加的复杂度和计算成本？这成了一个“性价比”未知的抉择，阻碍了吸附工艺的精准设计与优化。
针对这一瓶颈，研究人员在《Separations》上发表了一项研究，系统评估了从简单到复杂的五种扩散模型，在预测丁醇水溶液在两种不同活性炭上吸附突破曲线时的表现，旨在阐明分数阶扩散模型的应用价值与适用条件。研究发现，没有“放之四海而皆准”的最佳模型，模型性能与吸附剂类型强相关，在某些活性炭上分数阶模型优势明显，而在另一些上传统模型已足够精确。这为吸附系统模型的“按需选择”提供了重要依据。
主要技术方法
研究结合了实验室自制数据和文献数据。关键技术包括：(1) 通过循环填充床吸附实验和批量实验测定吸附等温线，并使用高效液相色谱(HPLC)分析浓度，数据拟合至Freundlich等温线。(2) 在实验室规模的固定床吸附柱中进行突破实验，通过精确校准的蠕动泵控制进料，并定时收集流出液用HPLC分析，通过积分突破曲线面积获取吸附容量。(3) 建立了五类数学模型（包括基于瞬时/非瞬时吸附、Fickian扩散和分数阶微分方程的扩散模型）并进行数值求解，以模拟和预测突破曲线，最后通过将模拟结果与实验数据拟合来评估模型性能。
研究结果
1. 引言
吸附柱的设计与优化高度依赖于描述目标物在吸附剂上平衡行为的吸附等温线，以及用于模拟突破曲线的动态吸附模型。尽管存在从分子动力学到神经网络等先进等温线预测方法，实验测定仍是获取平衡数据最可靠的手段。模拟突破曲线则需要求解一组耦合的微分方程，描述柱内的对流传质、颗粒内扩散和吸附剂表面的吸附过程。研究指出，模型“应尽可能简单，也需足够复杂”，强调了平衡模型复杂性与计算效率的重要性。本研究正是通过比较不同复杂度的模型，来评估它们再现丁醇在活性炭上吸附突破曲线的能力，并明确分数阶扩散模型在何种条件下能提供更优的描述。
2. 材料与方法
2.1. 材料
研究使用了两种活性炭：实验室实验采用F-400 (AC F-400)，对比文献数据则采用了Norit ROW 0.8。实验原料为纯丁醇和去离子蒸馏水配制的溶液。
2.2. 方法
等温线实验在25°C下，通过装有20g AC F-400的填充床吸附器（高0.175m，直径0.015m）进行。已知初始浓度的丁醇溶液在床内循环直至平衡，通过物料衡算获取单点等温线数据，多点数据拟合成Freundlich等温线。对比的文献数据（Xue等人）在37°C下通过批量实验获得。突破实验则在装有39g AC F-400的不锈钢柱（高0.195m，直径0.02m）中进行，研究了三种不同流速，通过蠕动泵精确控制流量，并用HPLC分析浓度。
3. 不同模型的描述
3.1. 问题陈述
研究考虑了一个长度为L
C、直径为D
C的填充床，内装直径为D
P的均匀多孔吸附剂颗粒。初始浓度为C
L0的含A物质水溶液以恒定表观流速u
S（对应体积流量F）流入。溶液流经床层空隙（空隙率ε
C，实际流速u
L），物质A通过对流传质到颗粒外表面，然后扩散进入颗粒内部的可移动相（占据颗粒空隙率ε
P，浓度C
P），并在吸附相（浓度C
A）上发生吸附。模拟目标是获得出口浓度随时间变化的突破曲线。
3.2. 不同数学模型
3.2.1. 模型1
该模型描述了占据床层空隙空间（εC）的体相流体中A物质浓度的演化，包含了轴向弥散、柱内对流以及从体相溶液到吸附剂颗粒外表面的传质过程。其核心方程表达了浓度随时间的变化（?CL/?t）等于轴向弥散（εCE ?2CL/?z2）减去对流项（εCuL?CL/?z）再减去外表面传质驱动力[kLa (CL- CL*)]，其中kL是A物质的传质系数，E是轴向扩散系数，a是填充床的比表面积，CL*是颗粒外表面液相浓度。该模型假设颗粒内吸附瞬间达到平衡。
3.2.2. 模型2
在模型1的基础上，模型2进一步考虑了颗粒内部的扩散过程，不再假设吸附瞬时完成。它引入了描述颗粒内部可移动相中A物质扩散的方程。这意味着浓度CP不仅随颗粒径向位置变化，也随时间变化，与体相浓度CL通过边界条件耦合。
3.2.3. 模型3
模型3是对模型2的简化，它忽略了轴向弥散效应。这意味着在描述柱内流体流动和传质时，只考虑了对流和外部传质，而不考虑由浓度梯度引起的逆向扩散。这使得模型的控制方程得以简化，可能降低计算成本，适用于轴向弥散效应不显著的情况。
3.2.4. 模型4
模型4是模型2的另一种简化形式，它假设颗粒内部的吸附过程是瞬时发生的。虽然颗粒内部仍有扩散，但吸附反应速率极快，颗粒内任一点的吸附相浓度CA与可移动相浓度CP之间始终满足平衡等温线关系。这简化了颗粒内部的动力学描述。
3.2.5. 模型5 (分数阶扩散模型)
这是本研究中引入的创新性模型。它放弃了传统Fickian扩散的框架，转而采用分数阶微分方程来描述颗粒内部的异常扩散行为。在经典扩散中，均方位移与时间呈线性关系。而在某些多孔介质中，由于结构异质性、死端孔隙或吸附-解吸迟滞等因素，扩散可能呈现“亚扩散”行为，即均方位移与时间的α次幂成正比（0< />
结论与讨论
本研究对五种不同复杂度的吸附扩散模型进行了系统的开发、评估与比较，核心结论是：没有单一的“最佳”模型适用于所有情况，模型性能强烈依赖于所使用的吸附剂类型。对于某些活性炭，传统的基于Fickian扩散的模型（如模型2或4）已能提供足够精确的突破曲线预测。然而，对于另一些活性炭，引入分数阶扩散项（模型5）的复杂性是必要的，因为它能够更准确地捕捉颗粒内部的异常扩散现象，从而显著提高预测精度。
这一发现具有重要的理论和实践意义。在理论上，它证实了多孔吸附剂结构的复杂性可以直接导致传质行为的偏离，而分数阶微积分为此类“异常”行为提供了一个强有力的数学描述工具。在工程实践上，研究为吸附工艺的设计者和研究者提供了清晰的模型选择指南：在着手建模前，应首先评估吸附剂的性质。如果吸附剂结构可能导致显著的异常扩散，则值得投入资源开发和校准分数阶扩散模型，以获得可靠的设计参数。反之，如果传统Fickian模型已能很好地拟合实验数据，则无需引入不必要的复杂度，从而节省计算成本和开发时间。这完美呼应了“模型应尽可能简单，也需足够复杂”的原则，为实现吸附系统的高效、精准设计与优化迈出了关键一步。
