应 MA婕 同学要求，仿照COMSOL 2006年会Tanai L. Marin的《Solidification of a Liquid Metal Droplet Impinging on a Cold Surface》及《COMSOL Multiphysics有限元多物理场建模与分析》，用COMSOL 4.2A做了一个水平集模型，来模拟存在相变的自由表面流动问题。

参照《建模》中的无量纲参数，结合COMSOL自带的水平集模型及传热模型，制作了本模型。实际计算发现，带入模型原始材料参数，计算不能收敛，这个问题放一边，现将测试参数计算了一遍，效果还算满意。

但是好像结果有待进一步商榷，用扩散场模拟水平集计算模型的图和内置的水平集模型的图有些小出处。

surface为相，contour为等温线，arrow为速度场。

计算参数：

rhog 1
rhol 10
etag 1
etal 1
Ti 1
Tc -0.2
e 0.1
cp 1
kl 0.5
kg 0.01
L 1
C 1600
q 0.001
Tm 0
sigma 1

——————

rho rhol+(rhog-rhol)*phils
eta etal+(etag-etal)*phils
k kl+(kg-kl)*phils
lambd (T-Tm+e)/(2*e)*((T<=(Tm+e))*(T>=(Tm-e)))+(T>(Tm+e))
A -C*(1-lambd)^2/(lambd^3+q)
Sx -A*u
Sy -A*w
cpH L*fldc2hs(T-Tm,e)
cpT cp+cpH

表面张力系数未知，所以先忽略它。

模型基本设置：

 

mxio

2012.4.28

 

模型文件：

http://vdisk.weibo.com/s/avm1xZ8HzG9HW