三维枝晶生长的数值模拟-材料科学研究论文.pdfVIP

三维枝晶生长的数值模拟1 潘诗琰，朱鸣芳 东南大学江苏省先进金属材料高技术研究重点实验室, 南京 (211189) E-mail: zhumf@ 摘 要：本文建立了一个模拟在低Péclet 数条件下三维溶质枝晶生长的数值模型. 本模型采 用Zhu 和Stefanescu (ZS) 提出的溶质平衡方法, 即根据固/液界面的平衡浓度和实际浓度之 差计算固/ 液界面演化的驱动力. 界面的平衡浓度由界面温度和曲率所确定, 实际浓度通过 用有限差分法对溶质扩散控制方程进行数值求解而获得. 该方法能够合理地描述枝晶从初 始的非稳态到稳态的生长过程, 并且具有较高的计算效率. 为了描述具有不同晶体学取向的 三维枝晶生长, 本文提出了一种权值平均曲率算法来计算固/ 液界面的曲率. 在权值平均曲 率的算法中耦合了界面能各向异性的因素. 该算法简单易实现, 并易于从二维推广到三维系 统. 为了对模型进行验证, 将模拟的枝晶尖端稳态生长数据和理论模型的预测结果进行比较. 结果表明: 模拟的 Al-2wt%Cu 合金枝晶尖端稳态生长速度、半径、平衡液相浓度和 Péclet 数随过冷度的变化接近于 Lipton–Glicksman–Kurz (LGK)解析模型的预测结果. 初步比较分 析了二维和三维枝晶生长的模拟结果. 模拟分析了稳态枝晶尖端的形貌. 发现三维枝晶尖端 是非轴对称的, 以四次对称的方式偏离旋转抛物面. 最后, 应用所建立的模型模拟出具有发 达分枝和不同晶体学取向的三维等轴多枝晶和定向凝固中柱状树枝晶的生长形貌. 关键词: 微观组织模拟；溶质枝晶生长；权值平均曲率；三维 中图分类号：TG244, O242 1. 引言 随着凝固理论和数值模拟技术的不断发展, 凝固过程中微观组织形貌的模拟研究已成 为可能并日趋完善. 近年来二维枝晶生长的数值模拟研究已取得了很大的进展. 相比较之下 三维枝晶模拟的难度仍然较大. 数值模拟基于数学物理模型. 在现今描述合金凝固组织的数 学物理模型中, 总体分可为两大类: 相场模型和尖锐界面模型. 相场(phase field, PF)模型使用能量密度泛函统一处理各种复杂的热力学体系问题, 自动 生成弥散的界面. 国内外学者们发展了各种PF模型模拟了纯物质或合金中的三维枝晶生长 [1-5]. 但是, PF模型的界面厚度约束限制了其模拟尺度, 而且为了缩短计算时间，PF模拟常需 采用高于常规铸造工艺条件下的大过冷度. 目前PF方法对三维枝晶的模拟主要针对纯物质 单枝晶的定量模拟[1], 或对合金中的单枝晶和多枝晶的形貌进行定性模拟[2-4]. 尖锐界面模型有元胞自动机方法（cellular automaton, CA ）、水平集方法（level set, LS ）、 界面追踪方法(front tracking, FT). Tan 和 Zabaras 运用LS方法追踪界面, 模拟了纯金属和合 金凝固[6,7] 的三维枝晶形貌. 但是LS方法需要在每一时间步长下重新初始化level set 函数和 延拓界面法向速度, 因此, 算法较复杂, 计算效率较低. CA方法能模拟实验所观测到的各种 凝固微观组织特征并且算法较简单, 具有较高的计算效率. 学者们发展了各种CA模型成功 模拟了实际铸件凝固条件下单枝晶自由生长和多枝晶的竞争生长过程[8-10]. 但至今CA方法 的三维枝晶模拟主要还局限于定性的形貌模拟. 将尖锐界面模型从二维推广到三维的关键问题是合理描述与界面能各向异性相关的三 维枝晶的不同择优取向生长. 根据界面的 Gibbs-Thomson 关系, 晶体生长的曲率过冷度取决 于包含了界面能各向异性作用