外场下δ相沉淀早期机制计算机模拟的研究.pdfVIP

_______________________________________________________________________________ 外场下d?相沉淀早期机制的计算机模拟研究 王永欣，陈铮，刘兵，马良，唐丽英，赵宇宏 （西北工业大学 材料科学与工程学院，西安 710072 ） [摘 要]基于离散格点形式的微扩散方程 （Langevin 方程）和非平衡自由能函数，编制了引入原子间相互 作用能变化的Al Li （d?）相沉淀原子层面计算机模拟程序。该程序包容亚稳区到失稳区的全部温度、成分范 3 围，孕育期至粗化的全过程，可以处理与时间相关的过程问题。开展了不同原子间相互作用势下原子图像、 序参数的计算机模拟，进而探讨了最近邻原子相互作用能 （W ）对有序相沉淀的影响机制。发现过渡区合金 1 形核前出现短暂的等成分有序化阶段。探明随W1 的增大，有序相沉淀的孕育期缩短，形核率增加，合金有序 化速度和原子簇聚速度加快，在所研究的时间范围内达到的长程序参数和成分偏离序参数最大值增大。随W1 的增大，有序相的析出呈现出失稳分解的特征。 [关键词] 原子间相互作用势；沉淀；原子图像；序参数；计算机模拟 [中图分类号] TG146 [文献标识吗]A [1] 计算机模拟作为科学研究 的重要手段已经应用于很多方面，分子动力学方法和蒙特卡洛 方法是应用最为广泛的两种方法。分子动力学方法 （MD 方法）又可以分为两类，其一成为经 典分子动力学方法，其二是第一性原理分子动力学方法 （CP 方法）。蒙特卡洛方法是针对某种 事件出现的概率建立一个随机模型，然后制造一系列随机数以模拟这个过程，最后进行统计 [2] 性处理 。在 MD 方法的处理过程中，方程组的建立是通过对物理体系的微观数学描述给出 的。该方法可以处理与时间有关的过程问题，因而可以处理非平衡态问题。但是该方法的计 算程序复杂，计算量很大，内存占用大。蒙特卡洛方法是实现Gibbs 统计力学的途径，在处理 基态、平衡态或近平衡态、均匀体问题时取得了很大成功。随机模拟的方法的计算程序简单， 计算量小，内存占用小，但是该方法难以处理非平衡态问题。 相析出的早期阶段，体系可能出现很大的浓度梯度，可能从母相中产生随机分布的多个 新相胚芽，体系处于高度非线性、高度非平衡态。并且新相胚芽的形成过程恰恰是被关注的 [基金项目]国家自然基金资助项目 ） [作者简介]王永欣（1971－），男，讲师，硕士，西安西北工业大学材料学院（710072 ）电话：（029 e-mail：yongxin@ 1 _______________________________________________________________________________中国科技论文在线 焦点问题。此时应用MD 方法 （包括CP 方法）以及蒙特卡洛方法的应用受到局限。 采用基于离散格点形式的微扩散方程 （Langevin 方程）和非平衡自由能函数的计算可以 比较好的解决上述问题。该计算模型在计算时进行了傅立叶变换和投影，使得原本直接在实 空间求解一个关于三维、非线性问题的偏微分方程问题转变为在倒易空间求解一个关于二维、 -6 -3 -7 线性问题的常微分方程问题，使模拟计算量大为减小、时空间尺度分别达到10 s 到 10 s ，10 m -6 到 10 m 。本文以二元铝锂合金为研究对象，编