4. 晶体生长速率与粒径演化

成核决定“有多少晶核诞生”，生长则决定“这些晶核能长多大”。在工业结晶器中，最终产品的粒度分布（CSD）是成核速率与生长速率共同作用的结果。若只关注成核而忽略生长动力学，就无法解释为什么同样的过饱和度下，有的批次得到粗大均匀的晶体，有的却得到大量细粉。本节围绕生长速率模型、粒度无关生长（ΔL 定律）、粒度相关生长以及生长分散展开，并说明它们如何与成核速率耦合，共同塑造粒径演化轨迹。

一、生长速率的两种表达：线性生长速率与质量生长速率

晶体生长速率通常用线性生长速率 G 表示，定义为晶体某一特征尺寸随时间的变化率，单位 m/s：

G = dL/dt  ≈  ΔL/Δt

与之对应的是质量生长速率 R_M（kg/s），二者通过形状因子和表面积关联。对球形或近似球形晶体，若特征尺寸为 L、体积形状因子为 k_v、密度为 ρ_c，则单晶质量 m = k_vρ_cL³，质量生长速率可写为：

R_M = dm/dt = 3k_v ρ_c L² · G

这一关系说明：同样的线性生长速率下，大晶体的质量增长远快于小晶体，因为表面积按 L² 增长。工程上常以 G 作为动力学建模的核心变量，而在物料衡算和产量估算时再换算为 R_M。

二、生长速率模型：过饱和度驱动与表面集成控制

生长速率本质上由过饱和度驱动。常用的经验幂律模型为：

G = k_g · (C - C*)ⁿ = k_g · ΔCⁿ

其中 k_g 为生长速率常数，ΔC 为过饱和度（可用浓度差、过冷度或相对过饱和度 σ = (C−C*)/C* 表示），n 为生长级数，通常在 1～2 之间。当 n≈1 时，生长近似为扩散控制；当 n>1 时，表面集成（螺旋位错、二维成核）成为控制步骤。实际建模时，需先通过恒温恒过饱和度实验测定不同 ΔC 下的 G，再对数拟合得到 k_g 与 n。例如某示例体系在 30 ℃、ΔC 从 0.02 到 0.08 kg/kg 溶剂范围内测得 G 从 1.2×10⁻⁸ 增至 6.5×10⁻⁸ m/s，双对数拟合斜率约 1.6，即可判断该体系生长受表面集成显著影响。

三、粒度无关生长与粒度相关生长

粒度无关生长（size-independent growth，ΔL 定律）假设所有粒度的晶体在同一过饱和度下具有相同的线性生长速率，即 G 与 L 无关。这是 MSMPR 理想模型的基础，其粒度分布在对数坐标下呈直线，可由斜率反推 G 与成核速率 B⁰。ΔL 定律在工业上常作为一阶近似，尤其在中粗粒级、过饱和度较低时较为合理。

粒度相关生长（size-dependent growth，SDG）则认为小晶体生长更快、大晶体生长更慢，或反之。常见模型有：

• Abegg 模型：G(L) = G₀(1 + γL)，生长速率随粒度线性增大；
• ASL 模型：G(L) = G₀(1 + γL)ᵇ，引入指数 b 描述非线性；
• MJ2 模型：G(L) = G₀[1 − exp(−aL)]，小晶体生长受抑制、随 L 增大趋于平台。

判断是否采用 SDG，可通过多粒级生长实验：筛分不同粒级晶体，在相同过饱和度下分别测其粒度增长，若各粒级 ΔL/Δt 明显不同，则 ΔL 定律不成立。示例中若 100 μm 与 400 μm 晶体的 G 相差超过 30%，就应引入粒度相关模型修正 CSD 模拟。

四、生长分散及其与成核速率的耦合

生长分散（growth dispersion）指同一粒级、同一时刻的晶体个体之间生长速率存在差异，源于晶体表面缺陷、位错密度、杂质吸附和局部流场不均。它使 CSD 峰变宽、细端拖尾加重。生长分散通常用变异系数 CV 或生长速率分布宽度描述，在群体平衡模型（PBM）中通过生长速率核的分布项体现。

成核与生长的耦合体现在过饱和度分配上：成核消耗过饱和度、产生大量细晶，细晶比表面积大又加速生长消耗过饱和度，形成负反馈。若成核速率 B⁰ 过高，大量细晶来不及长大，CSD 细端抬升；若生长速率常数 k_g 偏大而成核受抑，则晶体粗大但可能团聚。工程上常通过控制过饱和度窗口、加入晶种、调节搅拌与降温速率来平衡 B⁰ 与 G，使粒径分布落在目标区间。

五、工程判断：从模型到粒径调控

在实际结晶工艺开发中，可先通过间歇实验测定 G–ΔC 关系，确定生长级数 n；再用多粒级实验判断是否偏离 ΔL 定律；随后在 PBM 框架下同时拟合 B⁰ 与 G 的动力学参数，并用实测 CSD 校验。若模拟 CSD 细端偏高，说明成核过强或生长分散过大，可降低过饱和度峰值、增加晶种量；若粗端不足，则需提高生长速率或延长养晶时间。需要强调的是，上述参数与数值均为教学示例，实际体系须以实验测定为准。理解生长速率模型与成核速率的耦合，才能把“控制粒径”从经验操作转变为可预测的动力学设计。