相场法模拟枝晶生长karma模型matlab咱们今天来玩点好玩的——用Matlab搞个金属凝固过程的枝晶生长模拟。相场法这玩意儿真是材料模拟界的万金油特别是Karma模型处理枝晶分岔那叫一个丝滑。先整点基础配置% 基础参数设置 N 256; % 网格点数 epsilon 0.02; % 界面厚度参数 tau 0.0003; % 时间松弛因子 delta 0.25; % 过冷度控制 dt 0.005; % 时间步长这几个参数看着简单但调起来能让人头秃。比如epsilon调太小界面容易数值震荡调大了枝晶尖端速度又不准。当年我为了调delta参数硬是盯着屏幕看了两小时动画。核心的相场变量φ和温度场U初始化得有讲究。咱们在晶核位置搞个高斯扰动[x, y] meshgrid(linspace(-1,1,N)); phi exp(-(x.^2 y.^2)*50); % 中心晶核 U -delta*ones(N); % 初始过冷环境这里用指数函数生成初始晶核50这个系数控制晶核尖锐度。实际模拟时会发现初始扰动不够尖的话可能长不出漂亮的分岔结构。时间迭代才是重头戏显式欧拉法虽然简单但胜在直白好改for step 1:2000 % 计算拉普拉斯项 lap_phi del2(phi); % Karma模型动力学方程 phi_new phi dt/tau*(epsilon^2*lap_phi phi.*(1-phi).*(phi - 0.5 U)); % 温度场演化 U U dt*del2(U) phi.*(1-phi)*dt; % 边界处理周期性边界 phi circshift(phi_new, [0,0]); % 每100步画个图 if mod(step,100)0 imagesc(phi); axis equal off; drawnow end end重点看相场演化方程里的(phi - 0.5 U)这一项这就是Karma模型的精髓。相比传统模型这里0.5的偏移量让各向异性控制更灵活。实际跑起来会发现当枝晶尖端U值变化时生长方向会发生微妙调整分岔就是这么来的。相场法模拟枝晶生长karma模型matlab跑个十分钟你就能看到冰花绽放般的生长过程。不过要注意数值稳定性——把时间步长dt调太大分分钟给你颜色看。有次我把dt从0.005改成0.01结果相场值直接爆到NaN整个画面白茫茫一片真干净。想要枝晶长得风骚得在各项异性项上做文章。试试给拉普拉斯算子加点料% 添加各向异性系数 theta atan2(gradient(y),gradient(x)); aniso 1 0.04*cos(4*theta); % 四次对称性 lap_phi aniso.*del2(phi);这个0.04的系数控制枝晶臂的尖锐程度。调到0.06的话能明显看到枝晶臂变出四个棱角特别适合模拟立方晶系的金属凝固。最后说个坑Matlab的del2函数默认是二维拉普拉斯但实际计算时记得除以网格间距的平方。比如咱们的x范围是-1到1真实间距dx2/(N-1)所以严格来说应该dx 2/(N-1); lap del2(phi)/dx^2;不过对于定性演示偷懒直接用del2的结果问题不大。但要是做定量分析这步错了整个动力学方程的量纲都会乱套。跑完程序记得喝杯茶看枝晶慢慢舒展比看ASMR还解压。哪天要是老板问你在干啥就说在研究非平衡态界面动力学逼格瞬间拉满。