3. 熔池动力学与传热传质
激光或电子束作用于金属粉末/基材后，熔池并不是一潭静止的液体。它内部同时存在由表面张力梯度驱动的马兰戈尼对流、由密度差驱动的浮力对流，以及由蒸发反冲压力驱动的表面凹陷与流动。这三类驱动力共同决定了熔池的形状、元素分布、气孔与裂纹倾向，也决定了熔池内温度梯度与冷却速率的空间分布。理解本节内容，核心是建立“驱动力—流场—温度场—凝固组织”的因果链。
一、三种驱动力的物理来源与量级判断
先明确参照约定：本文以熔池表面为观察面，温度梯度方向由熔池中心指向边缘为正方向；流速方向以从中心向外为“外流”，反之为“内流”。
马兰戈尼对流由表面张力梯度驱动。表面张力 σ 随温度 T 变化，即 dσ/dT。对大多数金属，dσ/dT 为负值，即温度越高表面张力越小。因此熔池中心高温区表面张力低，边缘低温区表面张力高，表面液体被拉向边缘，形成由中心向外的表面流；该表面流到达边缘后沿熔池底部回流，形成典型的外环流。若材料含硫、氧等表面活性元素，dσ/dT 可能变为正值，流动方向反转，变为由边缘向中心的内环流。这一反转会显著改变熔池形状与气孔行为。
浮力对流由温度差引起的密度差驱动。熔池中心温度高、密度低，边缘温度低、密度高，热流体上升、冷流体下沉，形成自然对流。浮力对流的强度可用格拉晓夫数 Gr 与雷诺数 Re 的比值判断，但在激光熔池中，马兰戈尼对流通常占主导，浮力对流在熔池较深、温度梯度较缓时贡献增大。
蒸发反冲压力在功率密度较高时不可忽略。金属蒸发产生反冲压力，使熔池表面凹陷形成匙孔，同时驱动熔池内产生强烈的向下流动。反冲压力与马兰戈尼力、浮力共同作用，决定匙孔的稳定性与气孔形成。
二、温度梯度与冷却速率的空间分布
熔池内温度梯度 G 与冷却速率 Ṫ 是决定凝固组织形态的两个关键参数。温度梯度定义为 G = |∇T|，冷却速率定义为 Ṫ = ∂T/∂t。二者与凝固界面推进速度 V 满足 Ṫ = G·V。在激光熔池中，温度梯度通常沿熔池中心到边缘逐渐降低，而冷却速率在熔池中心最高、边缘最低。这是因为中心区域靠近热源，热量向四周传导的路径短、散热快；边缘区域受基材预热影响，散热较慢。
以 316L 不锈钢激光熔覆为例（示例参数，非实测）：激光功率 800 W，扫描速度 8 mm/s，光斑直径 2 mm，熔池中心温度约 2100 K，边缘温度约 1700 K，熔池宽度约 1.8 mm。此时中心温度梯度可达 10⁶ K/m 量级，边缘降至 10⁵ K/m 量级；冷却速率中心约 10⁵ K/s，边缘约 10⁴ K/s。这种分布直接导致中心区域凝固组织细小、边缘区域组织粗大。
需要特别注意的是，马兰戈尼对流会显著改变温度梯度的分布。当 dσ/dT 为负时，外环流将高温液体从中心输运到边缘，使边缘温度升高、温度梯度变缓；当 dσ/dT 为正时，内环流将低温液体从边缘输运到中心，使中心温度降低、温度梯度变陡。因此，通过控制氧含量或添加表面活性元素，可以主动调节熔池内的温度梯度与冷却速率分布。
三、传热传质对凝固组织的影响
熔池内的传热传质过程直接影响凝固组织的形态与成分分布。传热方面，熔池中心冷却速率高，凝固界面推进速度快，固液界面处的成分过冷度大，容易形成胞状晶或柱状晶；边缘冷却速率低，成分过冷度小，容易形成等轴晶。传质方面，马兰戈尼对流将溶质元素从中心输运到边缘，导致边缘区域溶质富集，可能引发偏析或热裂纹。
在实际工艺中，可以通过以下方式调控熔池动力学与传热传质：
- 调整激光功率与扫描速度：提高扫描速度可增大冷却速率，细化凝固组织，但过高会导致熔池不稳定、气孔增多。
- 控制氧含量：氧含量影响 dσ/dT 的符号，从而改变马兰戈尼对流方向。例如，在氩气保护下氧含量低于 50 ppm 时，316L 不锈钢的 dσ/dT 为负，外环流主导；氧含量升高至 200 ppm 以上时，dσ/dT 可能变为正，内环流主导。
- 预热基材：预热可降低温度梯度与冷却速率，减少裂纹倾向，但会粗化凝固组织。
- 添加表面活性元素：如硫、硒等，可改变表面张力温度系数，主动调控熔池流动方向。
四、实操中的观察与验证方法
在实验或生产中，验证熔池动力学与传热传质状态，可采用以下方法：
- 高速摄像观察熔池表面：记录熔池表面流动纹路与匙孔形态。外环流时表面流线由中心向外辐射；内环流时流线由边缘向中心汇聚。
- 热电偶或红外测温：在熔池边缘不同位置布置热电偶，测量温度随时间变化，计算冷却速率。注意热电偶响应时间需小于冷却时间尺度。
- 金相与 EDS 分析：观察凝固组织形态与元素分布。中心细晶、边缘粗晶且边缘溶质富集，通常对应外环流主导。
- 数值模拟：采用 CFD 软件（如 Fluent、COMSOL）建立熔池模型，输入材料热物性参数（表面张力温度系数、密度、粘度、导热系数），求解动量方程与能量方程，获得流场与温度场分布。模拟结果可与实验观察相互验证。
以 COMSOL 模拟 316L 不锈钢激光熔覆为例（示例设置）：采用二维轴对称模型，熔池区域使用动网格处理自由表面，表面张力梯度作为切向应力边界条件，浮力项采用 Boussinesq 近似。材料参数：密度 7800 kg/m³，导热系数 20 W/(m·K)，比热容 500 J/(kg·K)，表面张力温度系数 -0.0004 N/(m·K)。求解后可得熔池最大流速约 0.5 m/s，中心冷却速率约 8×10⁴ K/s，边缘约 1×10⁴ K/s。这些数值属于合理假设下的模拟结果，实际值需根据具体设备与材料标定。
综上，熔池动力学与传热传质是连接工艺参数与凝固组织的桥梁。掌握马兰戈尼对流、浮力与反冲压力的竞争关系，理解温度梯度与冷却速率的空间分布规律，才能在实际操作中主动调控熔池行为，获得理想的凝固组织与缺陷控制效果。