枝晶生长的相场法数值模拟-华侨大学学报自然科学版.PDF

第 3 1 卷 　第 1 期 华 侨 大 学 学 报 ( 自 然 科 学 版 ) Vol . 3 1 　No . 1 　 　 　20 10 年 1 月 J ournal of Huaqiao U niver sit y (Nat ural Science) J an . 20 10 　 　文章编号 : 　100050 13 (20 10) 0 1000606 枝晶生长的相场法数值模拟 刘晶峰 , 赵紫玉 , 蹇崇军 , 江开勇 (华侨大学 机电与自动化学院 , 福建 泉州 36202 1) 摘要 : 　利用一个简单的相场模型 ,对过冷熔体中枝晶的生长形貌进行数值模拟 ,定性分析各向异性强度 、无 量纲潜热及扰动因子等对枝晶形貌的影响 结果表明 ,枝晶主干生长速度随着各向异性强度的增大而增大 , 并且枝晶结构强烈依赖于各向异性强度 ;随着无量纲潜热的增大 ,晶体生长过程中枝晶结构特征愈加明显 ;扰 动对枝晶结构特征有一定的影响 ,随着扰动因子的增加 ,枝晶间竞争生长加剧. 关键词 : 　枝晶生长 ; 相场法 ; 过冷熔体 ; 数值模拟 ; 形貌 中图分类号 : 　T G 111. 4 文献标识码 : 　A 金属凝固过程微观组织数值模拟技术 ,是以金属凝固理论为基础 ,利用数值计算技术和计算机技术 发展的最新成果 ,逼真地再现金属凝固过程中的晶体生长过程. 枝晶生长方式是晶体生长的典型方式. 凝固组织中的柱状晶和等轴晶 ,就是枝晶生长方式的不同表现. 研究枝晶生长规律 ,预测凝固过程复杂 的枝晶形貌 , 已经成为微观组织数值模拟研究领域的热点. 目前 ,凝固过程微观组织数值模拟主要有确 ( ) ( ) ( ) 定性方法 Det er mini stic Met ho d 、随机性方法 Stocha stic Met ho d 和相场方法 Pha seFiel d Met ho d 等. 确定性方法以凝固动力学为基础 ,符合晶体生长物理背景. 但由于它不能处理晶体生长过程中的一 些随机过程 ,所以无法再现凝固时枝晶的生长过程及预测枝晶的具体形貌. 随机性方法考虑了晶体生长 过程中的随机过程 ,利用概率统计理论研究晶粒的形核和长大 ,适用于柱状晶组织的形成及柱状晶与等 轴晶相互转变的模拟. 但该方法只能模拟枝晶的近似形状 ,而不能模拟复杂枝晶的形貌. 相场法通过引 入有序化变量 ———相场变量 ,克服了跟踪复杂固液界面的困难 ,在模拟枝晶生长方面有着其他方法无可 比拟的优势 因此 ,该方法已经成为研究枝晶生长规律的主流方法[ 19 ] . 本文利用相场法 ,对过冷熔体中 的枝晶生长行为进行研究. 1 　计算模型 采用的相场模型由 Kobaya shi[ 1 ] 提出 ,是基于自由能泛函得到的 ,关于纯物质过冷熔体中枝晶生长 的相场模型. 该模型包括相场及温度场两个控制方程. 相场控制方程为 φ 2 2 τ ε φ ( ) = - f φ 1 t ( ) τ (φ ) 式 1 中: 为常数, f φ 为函数f , m 关于相场变量的偏导数. ( ) 在二维直角坐标系中, 方程 1 具有的形式为 φ φ φ 2 1 τ εε εε (