comsol枝晶生长相场法模拟 二元合金 考虑溶质偏析枝晶生长这玩意儿在材料模拟里算是经典难题了。咱们用相场法搞COMSOL模拟的时候最刺激的就是看那些枝晶分叉怎么从混乱中长出来。这次搞的是二元合金体系重点得盯着溶质偏析这个捣蛋鬼——它能让晶体长得歪七扭八的。先整模型的核心方程。相场变量φ控制固液相变溶质浓度c的分布直接影响生长速度。这里有个坑溶质扩散方程得和相场演化耦合起来。COMSOL的PDE接口里塞进去的方程长这样// 相场演化方程 epsilon^2 * phi_t epsilon^2 * nabla^2 phi phi*(1-phi)*(phi - 0.5 30*lambda*(c - c0)) // 溶质场方程 c_t D * nabla^2 c (c_l - c_s) * phi_t注意看溶质项里那个30倍系数这可不是随便凑的数。实际调参时发现系数低于25会导致界面模糊成浆糊超过35又会让数值爆炸。这个参数背后藏着界面能和各向异性强度的博弈改天专门开一篇讲。网格划分是个技术活。枝晶尖端曲率半径通常只有纳米级但计算域动辄微米级别。局部加密网格时得用这样的脚本mesh1.create(mesh1,geom1); mesh1.feature(size).set(custom, on); mesh1.feature(size).set(hmax, 0.1*epsilon);这里hmax和相场界面厚度epsilon挂钩是关键。有次手滑写成0.5*epsilon结果尖端直接变铲子形状被导师吐槽像马桶刷。comsol枝晶生长相场法模拟 二元合金 考虑溶质偏析溶质再分配现象在动画里特别明显。跑完模拟导出的浓度场数据做后处理时用MATLAB画了个浓度剖面contourf(X,Y,c_data,20,LineColor,none) colormap(jet) caxis([c_min c_max*0.8]) //故意压低上限值让偏析更明显这时候发现固液界面前沿的溶质堆积形成驼峰跟经典文献里的图对上了。不过咱们模拟出的驼峰高度总比理论值高15%后来发现是各向异性函数里的模式参数没调准。收敛性问题最头疼。当枝晶开始分叉时残差曲线经常抽风。解决方案是启用自动时间步长加上这个阻尼项model.param.set(gamma, 0.7*epsilon^2) //数值阻尼系数有个反直觉的现象增大溶质扩散系数D反而会让计算更稳定。因为高扩散速率能快速抹平浓度梯度相当于给系统加了天然稳定器。最后来个炫酷的——用COMSOL的粒子追踪模块模拟枝晶生长时被甩出的溶质原子运动轨迹。虽然物理上有点开脑洞但视觉效果绝对够发朋友圈study.step(time).set(tlist, range(0,0.1,10)); model.component(comp1).physics(pt).feature(init).set(u0, c*1e-6);跑完发现溶质原子在枝晶臂之间玩起了漂移形成类似星际尘埃带的图案。这或许能解释某些合金中观察到的周期性偏析条纹当然要验证还得做TEM。总之相场法这工具用好了比XBox刺激多了。