4. 分子扩散与费克定律在浸取中的应用

固液浸取的传质过程,本质上是溶质从固体内部向液相主体迁移的过程。在浸取剂渗入固体孔隙、溶解目标成分之后,溶质必须穿过孔隙中的滞留液体层,才能进入外部流动的主体溶液。这一"穿过滞留层"的迁移,主要依靠分子扩散完成。因此,理解分子扩散的规律,掌握费克第一定律与第二定律的物理意义,并能用一维扩散模型描述浸取剂或溶质在固体孔隙中的迁移,是分析浸取速率、判断过程控制步骤、选择工艺条件的基础。

一、分子扩散的物理图像与参照约定

分子扩散是指由于浓度差存在,物质分子从高浓度区域向低浓度区域净迁移的现象。它不是宏观流动,而是分子随机热运动在浓度梯度作用下的统计结果。为便于后续公式与图示保持一致,本节统一约定:以固体颗粒表面指向液相主体为 x 轴正方向,浓度 c 随 x 增大而降低,即 dc/dx < 0;扩散通量 J 取正号表示沿 x 正方向(由固体向液相主体)迁移。

在浸取体系中,扩散可发生在两个位置:一是浸取剂从液相主体穿过颗粒外表面液膜进入孔隙;二是溶质从固体内部孔隙向外扩散至颗粒表面。两者都可用同一套扩散方程描述,只是扩散介质(孔隙液体)与扩散路径长度不同。

二、费克第一定律:通量与浓度梯度的关系

费克第一定律描述稳态扩散下通量与浓度梯度的关系,其表达式为:

J = -D · (dc/dx)

式中 J 为扩散通量,单位 mol/(m²·s);D 为扩散系数,单位 m²/s;dc/dx 为浓度梯度。负号表示扩散方向与浓度梯度方向相反,即物质由高浓度向低浓度迁移。物理意义可概括为:扩散通量正比于浓度梯度,比例系数为扩散系数 D。D 越大,同样浓度梯度下迁移越快;D 由溶质性质、介质黏度、温度决定,温度升高、黏度降低时 D 增大。

费克第一定律适用于稳态或准稳态扩散,即浓度分布不随时间变化的情形。在浸取中,当颗粒外液膜很薄、主体溶液充分流动、液膜内浓度分布很快建立稳定梯度时,可用第一定律估算穿过液膜的传质通量。例如,某示例条件下设液膜厚度 δ = 50 μm,溶质在液膜两侧浓度差 Δc = 2 mol/m³,扩散系数 D = 1×10⁻⁹ m²/s,则按一维稳态近似,通量 J ≈ D·Δc/δ = 1×10⁻⁹ × 2 / (50×10⁻⁶) = 4×10⁻⁵ mol/(m²·s)。该数值属于教学示例假设,用于说明量级关系,并非实测结果。

三、费克第二定律:浓度随时间和位置的变化

实际浸取中,固体孔隙内的浓度分布往往随时间变化,属于非稳态扩散。此时需用费克第二定律描述,一维形式为:

∂c/∂t = D · ∂²c/∂x²

其物理意义是:某位置浓度随时间的变化率,正比于该处浓度分布曲线的二阶导数(即浓度曲线的"弯曲程度")。浓度曲线越弯曲,局部浓度变化越快;当浓度分布趋于直线时,二阶导数为零,系统接近稳态。第二定律把"时间"引入扩散过程,是分析浸取初期、颗粒内部浓度随时间演化、以及判断扩散是否达到准稳态的核心工具。

对第二定律的求解需要初始条件和边界条件。浸取中常见的简化情形是:半无限大介质、初始浓度均匀为 c₀、表面浓度突然维持为 c_s,其解可用误差函数表示。工程上更常用的是由第二定律导出的特征时间估算:扩散穿透深度大致与 √(D·t) 同量级。据此可判断:若颗粒半径 R 远大于 √(D·t),则扩散尚未深入颗粒中心,过程受内扩散控制;若 R 远小于 √(D·t),则颗粒内部浓度趋于均匀,过程可能转为外液膜控制。这一判据对工艺选择十分关键。

四、一维扩散模型在孔隙迁移中的应用

把固体孔隙简化为一段长度为 L 的一维通道,是浸取分析中最常用的建模方式。设孔隙内浸取剂或溶质沿 x 方向迁移,两端浓度分别为 c₁(固体侧)与 c₂(液相侧),则稳态下通量可写为 J = D_eff·(c₁ − c₂)/L,其中 D_eff 为有效扩散系数。由于孔隙曲折、截面收缩,D_eff 通常小于纯液体中的 D,可用 D_eff = D·ε/τ 估算,ε 为孔隙率,τ 为曲折度因子。示例中若 ε = 0.3、τ = 3,则 D_eff 约为 D 的 0.1 倍,说明孔隙结构会显著降低迁移速率。

应用一维模型时,可按以下思路操作:先确定扩散路径长度 L(液膜厚度或颗粒半径量级);再估算 D_eff;然后根据过程是否稳态选择第一或第二定律;最后用通量或特征时间判断控制步骤。若外液膜传质系数小、液膜厚,则外扩散控制,应加强搅拌、提高流速;若内扩散阻力大、颗粒粗、D_eff 小,则应减小粒径、提高温度或延长浸取时间。下面的图示把孔隙内的一维浓度分布、通量方向与两定律的适用区间对应起来,便于对照理解。

固体孔隙内一维扩散:浓度分布与通量方向 固体颗粒内部 孔隙液体 c₁ 高 滞留液膜(扩散路径 L) 浓度由 c₁ 降至 c₂ 液相主体 c₂ 低、充分流动 c₁ c(x) c₂ x x=0(固体表面) 液膜厚度 δ J = −D·dc/dx(由高浓度向低浓度) 费克第一定律 J = −D·(dc/dx) 适用:稳态/准稳态扩散 用于估算液膜传质通量 费克第二定律 ∂c/∂t = D·∂²c/∂x² 适用:非稳态扩散 用于判断穿透深度 √(D·t) 注:图中浓度曲线与参数为教学示意,非实测数据

五、从定律到工艺判断

把两定律用于浸取工艺选择,核心是判断"哪一步最慢"。若液膜内浓度梯度大、液膜厚、搅拌不足,则外扩散控制,此时提高搅拌强度、增大流速可显著加快浸取;若颗粒内部孔隙长、D_eff 小、颗粒粗,则内扩散控制,此时减小粒径、提高温度、延长浸取时间更有效。用一维模型估算特征时间 √(D_eff·t) 与颗粒半径 R 比较,是快速判断控制步骤的实用方法。例如示例中 D_eff = 1×10⁻¹⁰ m²/s、t = 3600 s,则 √(D_eff·t) ≈ 6×10⁻⁴ m = 0.6 mm;若颗粒半径 R = 2 mm,明显大于该值,说明一小时内外扩散尚未深入颗粒中心,过程受内扩散控制,应优先考虑减小粒径或升温。

学习建议:遇到具体浸取问题时,先画出"固体—液膜—主体"的一维示意,标出 c₁、c₂、L 与 D_eff,再判断稳态与否,最后用通量或特征时间下结论。这样能把抽象的费克定律转化为可操作的工艺判断依据。

需要强调的是,费克定律描述的是分子扩散这一基本机制,实际浸取中还可能伴随对流、化学反应、颗粒收缩等复杂因素。一维模型是对真实过程的合理简化,其价值在于抓住浓度梯度、扩散系数、路径长度三个关键量,帮助我们在工艺选择时明确主攻方向,而不是追求与实测完全一致的数值。掌握这一点,才算真正把费克定律用到了浸取工程中。