使用移动边界热平衡：拟稳态近似将通过现有固体层的热传导视为相对于固液界面缓慢移动的瞬时过程。通过将界面释放的潜热与通过固体、管壁和内部冷却剂膜传递的热量相等，即可获得固体生长速率。随着层厚度增加，其热阻增大，因此生长速率会降低。

实际计算采用基于热阻的斯特凡平衡：确定从相变温度到冷却剂主体温度的热流速率，然后将该速率除以产生额外固体所需的潜热。

拟稳态近似的假设

热分布调整速度快于界面移动速度

尽管固液边界在移动，但固体中的温度场在每一瞬间都被近似为稳态。当通过固体的热扩散速度远快于层厚度发生显著变化所需的时间时，这种假设是合理的。

该近似并不意味着过程是真正的稳态。固体厚度和传热热阻仍随时间变化。

潜热在移动边界平衡中占主导地位

对于冷冻或结晶，与界面运动相关的主要能量需求是熔化潜热或结晶潜热。固体和液体中的显热积累被忽略或作为次要修正处理。

因此，界面平衡基于：

[ \text{移出的热量}=\text{新固体释放的潜热}。 ]

界面保持在相变温度附近

固液界面通常被设定为平衡相变温度 (T_m)，或者在浓度或压力效应显著时设定为有效结晶温度。

这为通过固体和管组件移出热量提供了热驱动力。

建立热阻模型

定义管和固体几何形状

设

(R_i) 为管内半径，
(R_t) 为管外半径，
(R_s(t)) 为移动的固液界面半径，
(k_w) 为管壁导热系数，
(k_s) 为固体层导热系数，
(h_i) 为内部冷却剂侧传热系数，
(T_c) 为冷却剂主体温度，
(\rho_s) 为固体密度，
(L) 为凝固或结晶潜热。

对于从管表面向外生长的沉积物，固体占据 (R_t) 和 (R_s) 之间的区域，其厚度为

[ \delta(t)=R_s(t)-R_t。 ]

表示单位管长的热阻

对于圆柱管，单位长度的冷却剂膜、管壁和固体层热阻分别为

[ \mathcal{R}_i'=\frac{1}{h_i,2\pi R_i}, ]

[ \mathcal{R}_w'=\frac{\ln(R_t/R_i)}{2\pi k_w}, ]

以及

[ \mathcal{R}_s'=\frac{\ln(R_s/R_t)}{2\pi k_s}。 ]

当有正当理由时，可以加入额外的热阻，例如管-固体界面的接触热阻或生长层外部的液体侧热阻。

因此，单位长度的总热阻为

[ \mathcal{R}_{\text{tot}}'

\frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s}。 ]

计算瞬时排热速率

在拟稳态近似下，单位管长的排热速率为

[ q'(R_s)

\frac{T_m-T_c}{\mathcal{R}_{\text{tot}}'}。 ]

该热流速率使用当前界面半径进行评估。因此，该计算是准稳态的，但仍能捕捉到由层生长引起的时间依赖性。

推导固体生长速率

写出潜热平衡

当界面向外移动 (dR_s) 时，单位管长新形成的固体体积为

[ dV'=2\pi R_s,dR_s。 ]

相应的释放潜热为

[ dQ_{\text{latent}}'

\rho_s L,2\pi R_s,dR_s。 ]

除以时间并将其与排出的热量相等，得到

[ q'(R_s)

2\pi R_s\rho_s L\frac{dR_s}{dt}。 ]

代入热阻表达式，得到控制生长方程：

[ \boxed{ \frac{dR_s}{dt}

\frac{T_m-T_c} { 2\pi R_s\rho_s L \left[ \frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s} \right] } } ]

该方程给出了固体层的瞬时径向生长速率。

解释生长速率的下降

随着 (R_s) 增加，固体热阻

[ \frac{\ln(R_s/R_t)}{2\pi k_s} ]

也随之增加。因此，到达相边界的热通量下降，界面推进速度变慢。

潜热项中的因子 (2\pi R_s) 也很重要：随着界面周长的增加，给定的径向增量在单位管长上会产生更多的固体体积。

工程估算的实用简化

对薄层使用平面近似

如果固体厚度相对于管半径很小，即 (\delta \ll R_t)，则圆柱热阻可以近似为平面热阻。生长方程变为

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L \left[ \frac{1}{h_i} + \frac{t_w}{k_w} + \frac{\delta}{k_s} \right]}, ]

其中 (t_w) 为管壁厚度。

如果将冷却剂膜和管壁热阻合并为一个恒定热阻 (R_0)，则

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L\left(R_0+\delta/k_s\right)}。 ]

积分以估算厚度随时间的变化

对于恒定属性和温差，积分得到

[ t

\frac{\rho_s L}{T_m-T_c} \left[ R_0\delta+\frac{\delta^2}{2k_s} \right]。 ]

该结果显示了两种限制机制。

在厚度较小时，冷却剂膜和管壁热阻占主导，因此生长随时间近似呈线性。在厚度较大时，固体导热占主导，出现类似于扩散的行为：

[ \delta \propto \sqrt{t}。 ]

对显著沉积物使用圆柱形式

当层厚度相对于管半径不小时，平面近似可能会低估变化的几何形状和热阻。此时应通过数值积分，或在展开热阻项后进行解析积分来求解圆柱方程。

对于基于管的试验装置，这种区别很重要，因为沉积物可能变得足够厚，以至于局部传热面积和对数导热热阻都会影响结果。

将模型应用于管式单元操作

估算有效传热系数

热阻模型识别了系统中限制生长的部分。内部冷却剂膜由 (1/h_i) 表示，而生长的固体在平面极限下贡献了 (\delta/k_s)。

可以根据总热阻定义一个有效总系数，并基于一致的参考面积：

[ U=\frac{1}{R_{\text{total}}}。 ]

这使得测量的生长速率可以与计算值进行比较，并能揭示主要限制因素是冷却剂侧对流、管壁导热还是固体层导热。

评估冷却剂操作条件

冷却剂温度和流速会影响热驱动力和内部系数 (h_i)。降低 (T_c) 通常会增加驱动力 (T_m-T_c)，而增加冷却剂流速可以降低内部膜热阻。

因此，该模型可用于预测冷却剂供应条件的变化如何影响初始生长速率以及达到目标层厚度所需的时间。

考虑温度和成分效应

在结晶系统中，(T_m) 可能随浓度、压力或界面处的溶质排斥而变化。如果液体逐渐浓缩，相变温度和可用驱动力可能会在操作过程中发生变化。

在这种情况下，平衡应作为位置、时间或局部成分的函数进行更新，而不是使用单一的恒定 (T_m)。

将模型与轴向管行为耦合

真实的传热装置中，冷却剂温度可能沿管长变化。此时局部驱动力为

[ \Delta T(z,t)=T_m(z,t)-T_c(z,t), ]

生长方程必须沿管长和时间方向同时求解。

这将产生局部厚度 (R_s(z,t))，从而识别出可能发生过度沉积、堵塞或非均匀结晶的区域。

理解权衡

近似忽略了瞬态显热

拟稳态模型忽略了冷却新形成固体和附近液体所需的时间。当斯特凡数（Stefan number）较小时，即潜热相对于温度变化相关的显热较大时，这种做法最合理。

当显热效应显著时，使用完整的瞬态移动边界模型更为合适。

恒定属性可能不切实际

导热系数、密度、潜热、粘度和传热系数可能随温度和成分变化。将其视为常数简化了分析，但可能会使预测的厚度和生长速率产生偏差。

当操作范围较广或实验数据表明与简单模型存在系统性偏差时，应考虑属性变化。

表面条件会改变界面行为

该模型通常假设界面定义明确且为一维径向热流。粗糙度、结垢、成核行为、自然对流、剪切和晶体形态可能会违反这些假设。

对于结晶，沉积物也可能是多孔或树枝状的，而不是致密的固体层。其有效导热系数和界面面积可能与理想值有很大差异。

生长会降低有效的过程能力

较厚的固体层会增加热阻，并可能降低装置的有效传热负荷。它还可能收缩流动通道、增加压降，或通过过长的停留时间或不均匀的晶体尺寸导致产品质量变化。

因此，最大化瞬时生长速率并不总是等同于优化整体单元操作。

为您的目标做出正确的选择

将拟稳态模型作为首选工程工具，仅在假设条件实质性影响结果时才增加模型复杂度。

如果您的主要重点是估算随时间变化的固体厚度： 在沉积物较厚时使用圆柱热阻平衡，在薄层估算时使用平面积分方程。
如果您的主要重点是识别传热限制： 分离冷却剂膜、管壁和固体层热阻，以确定哪个操作变量影响最大。
如果您的主要重点是优化冷却剂操作： 同时评估冷却剂温度和流速，因为驱动力和内部对流共同控制潜热移出速率。
如果您的主要重点是预测试验装置中的非均匀生长： 将局部生长方程与轴向冷却剂温度和成分平衡耦合。
如果您的主要重点是高保真预测： 当实验行为无法由拟稳态模型解释时，应包含瞬态显热、可变属性、外部液体热阻和沉积物形态。

拟稳态热阻平衡将移动的凝固边界转化为实用的传热计算，同时清楚地解释了为什么随着固体层变厚，生长速度会变慢。

总结表：

关键要素	工程见解
核心假设	热分布调整速度快于固液界面移动速度。
驱动力	相变温度减去冷却剂主体温度，ΔT = Tm − Tc。
传热路径	内部冷却剂膜、管壁和生长的固体层。
生长平衡	移出的热量等于新形成固体释放的潜热。
圆柱生长	对于厚沉积物，使用当前界面半径和对数固体热阻。
薄层估算	生长初期接近线性，随着固体热阻增加，趋向于 δ ∝ √t。
主要控制因素	冷却剂温度、冷却剂流速、热属性、潜热和沉积物厚度。
模型局限性	瞬态显热、可变属性、成分效应、结垢和多孔沉积物可能需要更高保真的建模。