在包含细微分散的第2相粒子的三维系统中用相场模拟晶粒的生长.docxVIP

在包含细微分散的第二相粒子的三维系统中用相场模拟晶粒的生长aComputationalMaterials Science Center,National Institute for Materials Science,Tsukuba 305-0047,JapanbDepartment of Materials Science and Engineering, Waseda University, Tokyo 169-8555, JapancMaterials Engineering Laboratory, National Institute for Materials Science, Tsukuba 305-0047, JapanReceived 19 December 2005; revised 15 March 2006; accepted 20 March 2006在包含细微的分散第二相粒子的三维多晶系统中用相场模拟晶粒的生长已经有人实验过。各种不同体积分数的第二相粒子的模拟已经演示过。在晶界的粒子分数已经评估。Rlim=1.42r/被用来预测晶体的最大粒径Rlim，r是半径，f是粒子的体积分数。这些模拟结果和Zener类型的有联系的宝贵的模拟结果相对比。2006 ActaMaterialia Inc. Published by Elsevier Ltd. All rights reserved.关键词：电脑模拟；多晶晶粒生长；Zener定位；显微结构；相场模型多晶材料的晶粒尺寸是从显微结构长度尺寸了解他们的机械性能的关键。冶金学家用一种经典的方法控制显微结构的纹理和长度尺寸，这个方法是加入低溶解度的溶质元素，这些溶质元素会在第二相粒子中沉淀以标出晶界。50多年前。Smith和Zener为包含一个随机色散僵化，粒子半径r稳定的晶粒生长系统提出了一个方程式[1]。最大的可得的晶体尺寸已经可有下式得出：R=a（1）式中f是粒子的体积分数，a和b是常数。在Smith-Zene假说里，a=4/3，b=1[2]。后来一些研究者对Smith-Zene方法[3-11]提出了一些修改。在这些修改中a的值大约为1，b的值为1左右。Zener标记的实验数据对于检验不同的假说[12]来说太有限和分散了。因此一些主要是以Monte Carlo（基于波特模型）[13-19]和界面跟踪法（基于一个顶点模型）为主的模拟方法被用来更详细地研究在晶粒生长中粒子的标记效应。最近，Moelans等人为在大小和形状是时间独立[22,23]的第二相粒子下的晶粒生长提出了一个相场模型。他们研究了单个粒子和单个晶界[22]之间的相互作用，还演示了在多晶系统[23]下Zener标记的2D模型。然而，在多晶系统下的3D模型有一个和大量多晶样品的实验数据接近的先决条件。在这篇论文里，我们利用动力学方向性再赋值[24,25]和平行码技术修改了Moelans的运算法则。在包含细微的分散第二相粒子的3D晶粒生长的晶场模型已经演示过了。在各种体积分数的粒子的简单立方晶格4003阵点中演示了电脑模型。利用这些计算结果在关系式（1）的基础下，粒子体积分数下的最大晶粒尺寸的依赖性可以分析了。这里将比较以前模型的结果。在多晶材料[26,27]的晶粒生长相场模型中，多晶材料的显微结构用一组定向场的变量来描述：η1 (r, t),η2(r, t),. . .,ηQ(r, t),这些用来区别不同的晶粒方向；这里Q是可能取向的数量。根据Moelans等人[22]所说，一个取决于空间的函数Φ（r）被用来描述粒子的分布，这个函数在一个粒子是等于1，在基质颗粒中等于0. Φ（r）被认为时间是常量。这个系统的自由能被这个表达式描述：F=.(2)最先的四个术语做出贡献归因于大量自由能。m和ε是正常数。这最后一个术语和所有的正常数是相场变量局部梯度的作用。只有在晶界和粒子矩阵的相互作用下它才不等于0。时间依赖的Ginzburg–Landau方程式描述了方向场变量的时间、空间变化。这个方程式为每个变量ηi (r, t),i = 1,2,. . . ,Q,提供了一个无保守的顺序参数，=-=-式中参数和晶界活动度相关的运动学系数。假设晶界能量和晶界活动度各向同性，这些系数｛｝和｛｝被固定在｛｝=和｛｝=L，式中i = 1,2,. . . ,Q。在波特模型和Fan和Chen（Moelans等人）的相场运算法则中，结晶粗大化的动力学和拓扑结构学被发现依赖于在模拟单胞中独一无二的晶向Q的数量。为了是3D模型适用于动力学方向性再赋值方法Krill和Chen修正了Fan和Chen的相场运算法则。我们扩充Krill和Chen的方案进入各向异性系统。用他们的方法，每次在晶粒取向和最近的晶粒相比较后动力学晶体取向再赋值。然而