多极展开为悬浮颗粒周围的溶质浓度提供了远场近似。工程师通过格林函数的各阶矩来表示颗粒的详细表面通量：单极项与颗粒溶质的净释放量或消耗量成正比，并按 (r^{-1}) 衰减；更高阶项则包括偶极项，其按 (r^{-2}) 衰减，并反映通量的方向性变化。当颗粒没有净质量生成时，单极项消失，偶极项通常成为远距离范围内的主导贡献。

核心要点：多极展开用少量具有明确物理意义的源矩，替代了复杂的表面尺度传质计算。单极项描述颗粒与周围介质之间的溶质总交换量，而偶极项描述外加浓度梯度或非均匀表面通量如何使周围场发生畸变。

为什么远场浓度对中试设备很重要

微观问题

在浸出容器、浆料罐或生物反应器中，溶质浓度由颗粒或细胞周围的扩散以及可能存在的对流所决定。每个不规则颗粒表面的边界条件可能涉及局部反应动力学、分配、渗透性或空间变化的质量通量。

在中试规模的悬浮液中，通常不可能解析每个颗粒表面每一点的边界条件。多极分析将由此产生的浓度扰动压缩为一系列项，使其能够在远大于颗粒尺寸的距离处进行计算。

展开如何表示浓度场

从格林函数开始

对于以扩散为主的连续相，浓度扰动可以概略写成涉及格林函数 (G) 的表面积分：

[ c'(\mathbf{x}) \sim \int_{S} G(\mathbf{x},\mathbf{y}),q(\mathbf{y}),dS_{\mathbf{y}}, ]

其中，(q(\mathbf{y})) 是颗粒表面 (S) 上的局部溶质通量，(c') 是相对于未受扰动场的浓度扰动。

对于无界三维流体，主导格林函数的形式为

[ G(\mathbf{x},\mathbf{y}) \propto \frac{1}{|\mathbf{x}-\mathbf{y}|}. ]

将该核函数围绕一个参考点（通常为颗粒中心）展开，即可得到多极级数。

单极项

单极系数由积分表面通量决定：

[ Q=\int_S q(\mathbf{y}),dS. ]

相应的远场浓度表现为

[ c'_{\text{mono}}(\mathbf{x}) \propto \frac{Q}{D r}, ]

其中 (D) 是连续相扩散率。(r^{-1}) 的依赖关系意味着，只要颗粒具有非零的净释放速率或吸收速率，单极项在足够远处就占主导地位。

从物理意义上看，(Q) 是颗粒与周围液体之间的溶质总交换量。在固液浸出过程中，它表示颗粒释放的溶质净量；对于吸收性或反应性颗粒，它表示净消耗量。

偶极项

偶极系数取决于表面通量的一阶空间矩，以及外加浓度梯度和颗粒边界条件。其远场贡献的形式为

[ c'_{\text{dip}}(\mathbf{x}) \propto \frac{\mathbf{p}\cdot \mathbf{r}}{r^3}, ]

其总体幅值与 (r^{-2}) 成正比。

偶极项描述方向性行为。由于外部浓度梯度、局部流体动力学或不对称表面性质，颗粒的一侧释放或吸收的溶质量可能多于另一侧。

单极项为零的情况

如果颗粒没有净质量生成，

[ Q=0, ]

单极贡献就会消失。此时偶极项成为远场中的主导项，尽管颗粒仍可能在局部区域强烈地重新分配或扰动溶质。

这种情况对于不可渗透颗粒和非反应性夹杂物尤其重要。它们不会生成或消耗溶质，但会通过重新引导扩散路径并扰动周围浓度梯度来改变传输过程。

将该方法应用于复杂颗粒

表示任意表面几何形状

颗粒不必是球形，该方法同样适用。其不规则几何形状通过表面通量分布及其产生的多极系数体现出来。

实际计算通常按以下步骤进行：

求解或近似求解局部表面传输问题。
对表面通量进行积分，得到单极项 (Q)。
计算一阶及更高阶空间矩，得到偶极项和更高阶系数。
在远场的观测点处计算截断展开式。

因此，详细几何形状被保留在各个系数中，而不必在每次远场计算中显式携带。

区分局部行为与整体行为

该展开并不能替代近场解。在颗粒附近，单独的表面特征、边界层、尾流以及内部扩散路径都可能十分重要。

相反，多极项为局部颗粒尺度问题向更大尺度悬浮液模型的转换提供了可控衔接。这使其适合将颗粒尺度传质与中试规模的物料衡算相结合。

考虑外加浓度梯度

在悬浮液中，颗粒通常处于不均匀的背景浓度中，而不是完全均匀的背景浓度中。颗粒对该背景梯度的响应会对偶极矩产生贡献。

这样一来，该展开不仅可以描述溶质的净释放，还可以描述颗粒如何改变已有的浓度场。在估算有效传输性质时，这一区分至关重要。

将多极项与悬浮液尺度的扩散率联系起来

从单个颗粒到多个颗粒

对于稀悬浮液，可以将每个颗粒的远场扰动结合起来，以估算系综平均浓度场。单极项控制平均源项或汇项强度，而偶极项描述颗粒对梯度的影响。

在较高浓度下，颗粒之间会通过重叠的扰动场相互作用。此时，简单的单颗粒展开需要针对颗粒排列、流体动力学相互作用以及集体浓度极化进行修正。

麦克斯韦型有效传输模型

对于含有能够分配或扩散溶质的颗粒的稀悬浮液，麦克斯韦型关系为

[ \frac{K_{\mathrm{eff}}D_{\mathrm{eff}}}{D}

1+ 3\left(\frac{\gamma-1}{\gamma+2}\right)\phi + O(\phi^2), ]

其中：

(D) 是连续流体的扩散率；
(D_{\mathrm{eff}}) 是悬浮液的有效扩散率；
(K_{\mathrm{eff}}) 是有效分配系数；
(\phi) 是颗粒体积分数；
(K) 是颗粒-流体溶质分配系数；
(D_s) 是颗粒扩散率；
(\gamma=KD_s/D)。

单个夹杂物的偶极响应是这类有效介质描述背后的关键微观量。它衡量夹杂物对外加大尺度浓度梯度的响应。

不可渗透球形颗粒

对于完全不可渗透的球体，(\gamma=0)。预测的稀释极限结果取决于浓度的定义方式：

[ \frac{D_{\mathrm{eff}}}{D}

1-\frac{3}{2}\phi ]

当以流体相浓度为参照时，或

[ \frac{D_{\mathrm{eff}}}{D}

1-\frac{1}{2}\phi ]

当以体积平均浓度为参照时。

这两个结果并不矛盾。它们采用了不同的平均浓度定义，因此描述的是不同的有效系数。

理解其中的权衡

远场精度与几何细节

截断多极展开之所以高效，是因为它只保留了占主导地位的各阶矩。当观测点接近颗粒，或者颗粒变得高度细长、粗糙，或难以用低阶展开准确表示时，其精度会降低。

更高阶的多极项可以提高局部精度，但也会增加计算和解释的复杂性。

以扩散为主的假设

简单的 (r^{-1}) 和 (r^{-2}) 形式直接适用于有效无界介质中的稳态三维扩散。强对流、容器壁面、瞬态运行或湍流搅拌都可能改变格林函数及其空间衰减规律。

因此，在中试设备中应用无界介质形式之前，需要关注相关的佩克莱数、容器几何形状和混合状态。

稀悬浮液与浓悬浮液

基于孤立颗粒多极项的有效扩散率表达式，在颗粒体积分数较低时最为可靠。在浓浆料中，重叠的扩散场和颗粒接触可能使颗粒对之间以及多体效应变得显著。

如果不检查 (\phi)、颗粒间距和运行条件就直接使用稀悬浮液公式，可能会得到一个表观传质系数，并错误地将其归因于颗粒尺度物理机制。

净通量不等于局部通量

单极项为零并不意味着颗粒对传质没有影响。这只意味着积分表面通量为零。

非零偶极项仍然可以显著改变浓度梯度和有效扩散率，尤其是在不可渗透颗粒周围或颗粒暴露于外部梯度时。

根据目标做出正确选择

应根据需要预测的尺度和物理量来使用该展开：

如果您主要关注溶质的净释放或吸收：根据颗粒表面通量的积分计算单极项，并利用其 (r^{-1}) 场估算长距离浓度扰动。
如果您主要关注非反应性颗粒周围的梯度畸变：将净通量单极项设为零，并评估颗粒对外加浓度梯度的偶极响应。
如果您主要关注中试规模的悬浮液行为：将单颗粒多极系数与体积平均模型或有效介质模型结合，同时检查稀悬浮液假设是否成立。
如果您主要关注不规则颗粒几何形状：使用详细的近场计算来确定各阶矩，然后仅在超出颗粒尺度区域后应用多极展开。
如果您主要关注实验放大：在实际容器几何形状和搅拌状态下，将预测的远场浓度分布和有效扩散率与测量结果进行比较。

多极展开将复杂的颗粒尺度传质物理过程转化为易于解释的源项矩和梯度矩，并可将其与中试规模悬浮液性能联系起来。

总结表：

| 方面 | 关键认识 |\n|---|---|\n| 单极项 | 表示溶质的净释放或吸收；按 (r^{-1}) 衰减。 |\n| 偶极项 | 反映方向性通量和梯度畸变；按 (r^{-2}) 衰减。 |\n| 零净通量 | 单极项消失，使偶极项成为远场中的主导项。 |\n| 适用性 | 当观测距离远大于颗粒尺寸，即 (r \gg a) 时最为有用。 |\n| 复杂颗粒 | 通过表面通量各阶矩表示不规则几何形状。 |\n| 悬浮液建模 | 多极系数支持体积平均模型和有效扩散率模型。 |\n| 主要局限 | 精度取决于扩散主导条件、颗粒浓度、容器几何形状和搅拌状态。 |

需要将颗粒尺度传质与中试规模性能联系起来？联系 LABPARK，讨论适用于化学工程、生物过程与生物技术以及环境与水处理的教育与职业单元操作中试装置。LABPARK 帮助大学、科研机构和企业利用针对实际运行条件定制的实用中试设备，评估溶质释放、吸收、扩散、悬浮液行为和放大过程。