COMSOL模拟冰花形成（相场凝固模型）

COMSOL实现过程： 1. 建立一个2D，两个变量的General Form PDE模型(自变量为theta, T)，瞬态求解。 COMSOL模拟冰花形成（相场凝固模型）建模过程 2. Global Definition-> Parameters COMSOL模拟冰花形成（相场凝固模型）建模过程 3. Global Definition->Function->Step COMSOL模拟冰花形成（相场凝固模型）建模过程 4. Model -> Definition -> Variables epsilon EPS*(1+DELTA*cos(ANISO*(angle-ANGLE0))) epsilonp -EPS*ANISO*DELTA*sin(ANISO*(angle-ANGLE0)) epsilon2 epsilon*epsilon m ALPHA/pi*atan(GAMMA*(TEQ-T)) angle if(abs(thetax)>eps,atan2(thetay,thetax)+pi*(atan2(thetay,thetax)<0),pi/2*((thetay>0)+3*(thetay<0))) //注：angle 定义的有点问题有个突变 COMSOL模拟冰花形成（相场凝固模型）建模过程 5. Geometry 建立一个原点对称9×9的矩形 COMSOL模拟冰花形成（相场凝固模型）建模过程 6. PDE设置 COMSOL模拟冰花形成（相场凝固模型）建模过程 Conservative Flux: -thetax*epsilon2+epsilon*epsilonp*thetay -thetay*epsilon2-epsilon*epsilonp*thetax -Tx -Ty Source Term: d(epsilon2,x)*thetax+d(epsilon2,y)*thetay+theta*(1-theta)*(theta-0.5+m) 0 Damping or Mass Coefficient: TAU 0 -K 1 COMSOL模拟冰花形成（相场凝固模型）建模过程 初始值： theta: step1(x^2+y^2) T: 0 COMSOL模拟冰花形成（相场凝固模型）建模过程 7. Mesh 建立Mapped型网格，Maximum element size为0.1(可适当加密，0.06~0.08左右效果比较好) COMSOL模拟冰花形成（相场凝固模型）建模过程 8. Study求解时间为range(0,0.05,0.3) 9. 计算及后处理 mxio 2012.4.4 2013.6.28 补记： 参考文献下载： "A Numerical Approach to Three-Dimensional Dendritic Solidification.pdf" http://vdisk.weibo.com/s/HGiPJ "Modeling and numerical simulations of dendritic crystal growth.pdf" http://vdisk.weibo.com/s/HGhEa 模型下载： http://vdisk.weibo.com/s/avm1xZ8HA1D7f