COMSOL 相场法模拟各项异性枝晶形貌演变，初学者必备。 带文献
相场法模拟金属凝固过程中的枝晶生长，就像在虚拟实验室里养了一棵会分叉的冰晶树。最近用COMSOL折腾各向异性枝晶的形貌演变，发现这玩意儿对新手来说确实有点劝退。今天咱们就扒开代码看看门道，记得备好咖啡——毕竟相场方程可不是吃素的。
先来点硬核的：相场变量φ从0到1表示固液相变，配合温度场T玩二人转。COMSOL里定义各向异性最骚的操作是在表面能系数γ(θ)里动手脚。举个栗子，四重对称的表达式写成：
gamma_theta = 1 + epsilon*cos(4*(theta - theta0))这里的epsilon控制各向异性强度，theta0是晶体取向。参数设置别瞎填，ε=0.04时能得到漂亮的分叉，再大就等着数值爆炸吧（别问我是怎么知道的）。
模型的核心是这两个耦合方程：
tau*phi_t = W^2*div(grad(phi)) + phi*(1-phi)*(phi - 0.5 + m*(T-T_m)/(L/c_p)) T_t = alpha*laplacian(T) + L/c_p*phi_tW是界面宽度参数，tau是弛豫时间。建议新手先关掉温度场，单独调相场参数找感觉。有次我把W设成0.1，结果界面宽得像银河，枝晶直接变胖头鱼。
各向异性梯度项才是精髓所在。COMSOL内置的d()算子能玩出花：
theta = atan2(d(phi,y), d(phi,x)) // 计算局部法向角 aniso_term = epsilon*cos(4*(theta - theta0)) // 各向异性修正项注意这个theta计算会引发数值震荡，老司机都在后处理加个平滑滤波器。有文献说用五次多项式拟合各向异性函数更稳（Kobayashi, 1993），不过COMSOL自带的解析式已经够新手喝一壶了。
网格划分是另一个大坑。自适应网格必须开，但重划分频率设太高会卡成PPT。建议初始网格用极细化三角形，生长区域局部加密到W的三分之一。见过有人用笛卡尔网格强行跑，结果枝晶边缘全是锯齿，跟狗啃的似的。
收敛性调参堪称玄学。时间步长从1e-5开始试，相场方程用向后差分更稳。非线性求解器建议上自动牛顿法，雅可比矩阵更新频率调到保守。有回我手贱开了自动步长，结果在枝晶尖端处疯狂震荡——活像得了帕金森的程序。
最后安利两篇救命文献：Warren的《相场法入门》把数学物理讲得透透的，Kobayashi那篇PRL给出了最经典的二维枝晶参数。记住，调参时做好实验记录，不然参数乱套时连后悔药都没得吃。毕竟在相场模拟里，每个美丽的枝晶背后，都是成吨的报错日志堆出来的奇迹。
（文献：
- Warren J A, Boettinger W J. Prediction of dendritic growth and microsegregation patterns in a binary alloy[J]. Acta Metallurgica et Materialia, 1995.
- Kobayashi R. Modeling and numerical simulations of dendritic crystal growth[J]. Physica D: Nonlinear Phenomena, 1993.）
