增材凝固模拟，八十年代的老算法刚刚再登自然通讯

AM易道分享

Nature Communications最新论文预印本。

他们用的算法叫蒙特卡洛Potts，1984年就有了。

这套东西在材料计算圈算是教科书级别的存在，老得不能再老。

一个四十年前的老算法，怎么还能上Nature Communications？

要解释这件事，得先说说凝固模拟在增材里到底是什么角色。

激光过去那一下，金属熔化，再凝固。

几毫秒里，决定了打印件的微观组织，也就决定了这个零件能不能用、能用在哪里。

性能差异的根，几乎全在这几毫秒里。

但工程师能用的凝固模拟算法不多，就三样。

相场最准，物理上接近完美。

但慢得离谱。

论文里写用相场算一个350微米见方的小熔池要花超过300小时。

350微米比头发丝粗一点点，这种尺寸算两周，整个零件级别根本没法算。

元胞自动机比相场快，能做到工件级别的局部仿真。

问题是结果跟网格大小有关，得小心解读。

而且它只管凝固，晶粒长大还得另写一套代码拼上去。

蒙特卡洛Potts最快，简单粗暴。

但有个致命伤，只能给出大概形貌，凝固速度算不准，跟实验对不上。

工业界一直觉得它就是个示意图工具，不入流。

文章发表的是把第三个工具的致命伤解决掉。

第一，让格子尺度对得上。

Potts用几十微米的格子，真实晶核才几纳米。

这个尺度差距以前是糊弄过去的，他们加了个修正项，让大格子也能反映小尺度的形核物理。

第二是最关键的，他们做出了一个不依赖温度的真实时间换算公式。

原版Potts只有步数没有秒，算一万步都不知道现实过了多久。

新公式里不挂钩温度这一点对增材太重要了，因为熔池附近温差极剧烈，要是时间跟温度挂钩，模拟里不同位置时间流速会不一样，物理上根本说不通。

第三，潜在热释放。

意思是金属凝固会放热反过来影响周围液体。

以前的模拟要么忽略这个过程要么只能把这个过程仅仅硬限1度，现在他们按热量真实扩散的距离来决定节奏。

上面创新的三步，跑了纯铁、铜、镍、Al-4Cu四种材料的凝固速度对比，全部跟实验数据吻合。

模拟一个40乘5乘20毫米的纯铜铸锭，柱状晶平均厚度算出764微米，实验测的753微米，误差不到2%。

这是真实工件尺寸了，相场在这种尺寸面前直接躺平。

对增材的直接价值，论文里给了一个让人眼前一亮的演示。

AA2024铝合金，掺入Zr基金属玻璃粉作形核种子。

0%颗粒长柱状晶，3.7%颗粒部分等轴，8.8%颗粒完全等轴。

模拟图和实验照片几乎重合。

这件事用相场算法来做，工业尺度根本算不动。

粉末厂、打印厂、零件厂，理论上可以用新方法在电脑里先跑一遍，再去烧真粉。

能在仿真里先排除掉一半方案，省的就是真金白银。

畅想几个真实的机会窗口。

把算法做成SaaS或者Agent或者SKILL提供给3D打印服务商，是最直接的一条路。

粉末厂可以提供从配方到性能的预测报告，通过卖优化方案提高溢价。

航空件每改一个工艺参数就要重新认证，仿真先把可行范围圈出来能省下大量盲试。

还有一个容易被忽略的方向，传统铸造、连铸、焊接同样是凝固问题，市场比增材大，这套工具一样能用。

作者提到的局限性

作者提到算法没考虑液相扩散、没精确建模界面各向异性所以不能预测织构、没考虑熔池对流。

代码也开源在了Zenodo上，能下、能尝试复现。

专业读者请自行阅读原论文。

原论文：

Oh, S.-H., Lee, B.-J. A Monte Carlo Potts model for solidification. Nature Communications (2026)

地址：

doi.org/10.5281/zenodo.19426462

更通俗的解读

对于大部分不做底层模拟算法研究的读者来说，如果论文读起来太硬，AM易道做了个微信群聊形式的通俗版解读，让六个不同身份的从业者用聊天把这篇论文聊明白。