用于分析动态脱气实验的柔性开源质量传输框架

《ACS ES&T Air》：A Flexible Open-Source Mass Transport Framework for Analyzing Dynamic Outgassing Experiments

　　

质量传输建模在表征材料挥发性有机物（VOC）排放中具有关键作用，应用范围涵盖室内空气质量至食品包装。尽管现有模型被广泛使用，但其常因未提供源代码而缺乏可及性，或因依赖拉普拉斯变换方法而无法适应复杂实验条件而缺乏灵活性。研究人员在此提出一种Python语言下的通用开源框架以解决上述局限。该框架采用质量守恒的有限体积法表述方法线（MOL），结合高效常微分方程（ODE）求解器，可模拟板坯、棒材与球状几何结构中具有任意时变温度、流速、进料浓度及化学生成的扩散与质量传输。研究人员通过对解析解验证并复现已有文献结果，证明将模型扩展至包含有限界面传递系数可显著改善对实验吸附数据的拟合。进一步通过综合分析多次顶空萃取（MHE）与程序升温脱附（TPD）实验阐明框架效用，借助计算效率执行高分辨率参数扫描与全局多温度拟合。这些能力共同使该框架成为用于实验设计及解释材料排放与吸附数据的自适应开源工具。

研究背景与问题提出

质量传输模型是表征室内空气质量、食品包装等应用中材料排放的核心工具，用以量化挥发性有机物（VOC）在材料内部的迁移、扩散速率及气相浓度。既有研究大量依赖拉普拉斯变换方法，虽在计算静态顶空平衡等问题上高效，但面对时变温度、时变流速、浓度依赖扩散系数或内部化学反应等复杂边界条件时，需反复推导且难以自然容纳非线性效应；数值模型虽具适应性，却常因源代码不公开、计算效率低而未被广泛复用。为同时突破“不可访问”与“不灵活”双重瓶颈，研究人员开发了统一、高效且开源的动态脱气模拟框架，重点并非提出新控制方程，而是将分散于单一用途代码中的能力整合为单一工具，并借由计算速度赋能传统方法难以完成的参数扫描与多温度全局拟合。

期刊与贡献定位

该论文发表于《ACS ES&T Air》，其贡献在于：（i）提供经验证、免费可用且计算高效的实现，统一多种几何与边界条件；（ii）借助效率揭示新物理洞察，包括顶空萃取中吹扫诱导的分析物损失定量、多温度实验中原位化学生成解析，以及TPD动力学参数中扩散诱导误差的系统量化。

主要关键技术方法

研究人员采用保守有限体积法离散反应-扩散方程的方法线（MOL），引入Kirchhoff势处理浓度依赖扩散系数，表面采用Robin或平衡界面条件，顶空气相作完全混合常微分方程（ODE）；时间积分调用scipy.integrate.solve_ivp中的BDF刚性求解器，并提供稀疏解析Jacobian加速。框架支持板坯（m=0）、圆柱（m=1）、球体（m=2）三种对称几何，允许任意时变温度T(t)、进料浓度c_feed(t)、流速Q(t)，扩散系数D按D=D₀exp(c/c_D?E_a,D/R_gasT)随浓度与温度变化，分配系数K按K=K₀exp(?ΔH_s/R_gasT)变化，并可叠加任意个数一级均相化学源池c_source,i。验证阶段复现Kumar与Little的乙烯基地板酚吸附数据及球粒恒流脱气解析解（拉普拉斯留数反演），数值与解析在N≥50均匀网格下线性重合，早期对数偏移由守恒表面处理引起，经边界加密tanh网格消除；质量平衡误差E_MB维持在10^?5%水平。基准测试显示单台Intel Core i7-12800H笔记本上完整实验模拟耗时O(100 ms)，较COMSOL有限元快约60倍。

研究结果

引言（Introduction）部分阐明拉普拉斯法在线性时不变系统中的优势及其在TPD指数温度耦合、浓度依赖D、有限界面传递下的失效，引出数值MOL框架的动机。

方法（Methods）给出控制方程（式1反应-扩散方程含几何因子x^m，m=0/1/2对应板/柱/球）、源项一级衰减（式2）、原点对称与界面Robin条件（式3-4）、顶空质量平衡（式5）及初始条件（式6）；定义D与K的指数形式（式7-8）及STP溶解度转换（式9）；有限F退化为平衡界面时推得表面浓度ODE（式11）与初始连续条件（式12）；引入参考扩散时间τ_D、停留时间τ_res、萃取佩克莱特数Pe=τ_D/τ_res与传质毕奥数Bi_m=FR/D_ref作诊断无量纲数。

数值实现（Numerical Implementation）详述有限体积面通量用ψ(c)=∫D(c′,T)dc′闭合，守恒离散使内部通量成对抵消，仅表面通量与源项改变总库存；BDF求解器结合用户提供T(t)/Q(t)/c_feed(t)历史，内部时间单位取分钟，D与F做单位换算；提供均匀网格与边界加密tanh网格，GUI与源码发布于GitHub仓库MOLIERE。

验证与概念演示通过图1b球粒脱气、图2时变温度/流量/进料耦合三几何与单源项生成、图3乙烯基地板酚吸附原Laplace拟合复现（R²=0.920）及引入有限F后优化拟合（R²=0.984）三段展示：前者确认数值精度，次者展示框架柔性，末者证明有限界面阻力可修正K与D估计。

多项演示结论：（1）有限F界面模型较无限F平衡假设显著改善历史吸附数据拟合，并修订K、D数值；（2）MHE中顶空吹扫时长过长会引发样品质量损失从而低估分析物回收率，框架可定量该损失；（3）跨温度MHE全局拟合揭示持续VOC排放中存在原位化学生成贡献；（4）扩散限制扭曲MHE曲线与TPD峰形，致Kissinger法提取活化能产生系统误差，框架可量化该误差边界。

讨论与结论翻译

研究人员指出，该框架的新颖性不在控制方程本身，而在可访问、高效、统一的实现与由此启用的物理洞察；保守MOL离散在浓度依赖D下仍保质量守恒，BDF求解使高分辨率参数扫描与多温度全局拟合在秒级可行。框架可模拟任意时变流温料、三几何、多一级源、有限界面传递与指数温浓依赖D/K，弥补拉普拉斯法对TPD、MHE及塑料化系统的覆盖缺口。借由复现文献、量化吹扫损失、分离化学生成、校正TPD动力学误差四类示范，研究人员确立该开源工具在实验设计、排放清单核算及材料低排优化中的通用价值。结论确认：扩展有限界面传递系数可实质改善实验吸附拟合；顶空萃取吹扫诱导损失在特定条件下不可忽略；多温度MHE数据兼容持续化学生成解释；扩散限制对TPD派生动力学参数引入可量化偏差。该框架因而构成适配动态脱气实验设计的开源基础平台。