基于COMSOL的相变模拟（石蜡、熔盐、金属等）。 材料从完全固态转变到液态（或者液态冷却到固态），考虑液相的自然对流对相变过程的影响。 材料的参数设定与融化或凝固状态相关，如图中所示。
相变材料在储能和热管理领域玩得风生水起，但真要搞明白石蜡融化时液态区域怎么扭着腰对流，COMSOL里可得下点硬功夫。最近被实验室小伙伴抓着调试熔盐储罐模型，发现那些看似简单的相变设置里藏着不少魔鬼细节。
液态区域的华尔兹
传统相变模拟容易把液态当老实人——设定个固定粘度和密度就完事儿。但真实情况是熔化的石蜡会像热锅上的蚂蚁到处乱窜。这时候得在"层流"接口里开启Boussinesq近似，让密度差驱动自然对流。代码里动量方程长这样：
rho * (d(u,t) + u*nabla(u)) = nabla*(-p*I + mu*(nabla(u)+nabla(u)^T)) + rho_ref*g*beta*(T-T_ref)关键是把密度rho改成温度的函数，同时注意参考温度T_ref别设成室温，否则浮力项直接暴走。有次把熔盐的beta值多打了个零，结果液态区域直接上演滚筒洗衣机特效。
材料属性的变脸艺术
相变材料最骚的操作是固态液态属性突变。COMSOL里可以自定义材料属性函数，比如这样定义等效比热：
if (T<T_solidus) cp = 2000; elseif (T>T_liquidus) cp = 2500; else cp = 2000 + 500*(T-T_solidus)/(T_liquidus-T_solidus); end但实测发现这种阶梯式切换容易让求解器抽风。后来改用平滑过渡函数：
alpha = 0.1*(T_liquidus - T_solidus); phase = 0.5*(1 + tanh((T - 0.5*(T_solidus+T_liquidus))/alpha)); cp = cp_solid*(1-phase) + cp_liquid*phase + L*delta/(alpha*sqrt(π));其中delta是相变潜热分布系数。这招让镍基合金的凝固过程收敛速度提升3倍，特别是固液共存区的温度梯度不再跳霹雳舞。
移动网格的障眼法
处理固液相界面时，ALE移动网格和相场法各有拥趸。个人更推荐在"变形几何"里用速度边界条件：
mesh_velocity = phase_transition_rate * normal_vector;配合自适应网格细化，能捕捉到金属液滴凝固时的枝晶生长。有个坑要注意：当液相区出现涡旋时，网格畸变会像被猫抓过的毛线球。这时候别急着砸键盘，把网格运动阻尼系数从默认0.1调到0.3，世界瞬间清净。
实测案例显示，考虑自然对流后，直径10cm的石蜡球完全融化时间缩短23%。而金属相变时，忽略对流会导致固液界面预测误差高达18%。现在遇到最魔幻的情况是熔盐在相变时出现热驱动振荡，像极了烧开水时的卡门涡街，这或许能解释某些储热系统异常温升现象。
玩转这些套路后，终于理解为什么导师总说："相变模拟不是选择题，而是找茬游戏。"下次再看到COMSOL里扭曲的温度云图，不妨检查下液态区的雷诺数是不是偷偷越过了10这个坎——那意味着该从蠕流切到湍流模型了。
