Ni3Al合金熔点与比热的分子动力学模拟研究.pdfVIP

中【qT胃热物理学会 工程热力学jj能源利用 学术会议论史 编呼：061119 Ni3Al合金熔点与比热的分子动力学模拟 杨弘 吕勇军 陈民’ (济1}人学航天航空学院，传热‘j能源荦J用北京市章点实验室，北京100084) 摘要：奉义采用嵌入原了法对深过冷条件下Ni3Al合会的熔点和比热进{于J，分了动山学模拟．熔点的 模拟分别采用二叫冶夹层结构法和NVE系综法，两种方法扶f{}的熔点比较一毁。但均略高于实验值。 K 返应该是由：F实验中熔点的钡l量受到普遍存在的表_血熔化过程影响的结果。模拟获得的比热在800 到2000K范嗣内呈现随着温度的升高缓慢线性减小的趋势． 关键词：Ni，Al深过冷熔点比热 I引言 深过冷快速凝矧技术可阻在慢速冷却条件下实现了样品的快速凝固过程，它克服了 急冷凝嘲技术中冷却速度受传热条件限制的缺点，是实现人体积液态金属快速凝固的有 效途径．然而深过冷条件下金属快速凝同过程的定量化研究往往要受到液态金属热物理 性质数据匮乏的限制【¨。目前，对于液态金属，特别是深过冷液态金属的热物理性质的 实验测定主要是借助于电磁悬浮或超声恳浮等方法实现【2—1。这些无容器处理技术能够 有敛地消除容器壁的异质形核作用，有助于获得深过冷状态。但是由于实验成本很高， 目前的实验研究只集中在少数金属和合金系中。而且实验数据也相对有限，不能满足进 一步研究的需要。采用分子模拟方法研究深过冷金属的热物理性质可以不受实验条件的 影响，而且已有的模拟结果表明利用该方法对研究亚稳态金属的热物理性质是一条有效 的途释，可以为深入分析各种宏观热物理性质与微观分子结构之间的本征联系提供可能 性14·5。。 作为最基本的热力学参数，熔点的模拟有助于揭示平衡凝蚓过程中的相结构选择规 律。并准确地确定模拟过冷度。而比热是非平衡液蚓相变过程十分重要的热物理性质之 一。获知过拎熔体的比热值，可以得到在相应过冷度条件下，熔体发生液陶相变时的结 晶清热、熵变和吉布斯自由能改变，从而推算出过冷熔体的形核率，以及开展晶体生长 动力学分析。Ni3AI合金具有高熔点，抗高温氧化性好等优点，是新一代的高温结构材 料，在航天航空等领域有着广泛的应用。本文采J}{1分子动力学方法模拟了Ni3AI合金的 熔点及其在过冷条件下的比热随温度的变化规律。 2势能模型 中每一个原于都看作是嵌入在其它原子组成的墓体中的杂质。系统的能量表示为嵌入熊 和且相作用势能之和，而将多体作用归结,1．土u肌tl，、一匕z中．其基本方程如下嘲：、 一一 nn 巨 Fh 吼k 晖粥 ∑， n= 、●-， (1．2) ‘h ∑川 其中E为原子f的能鼙．Ph，i为箱i个原子处的电子密度，‘表示第，个原子对第f个 原子处电子密度的贡献，r,u为第f个和第，个原子的间距。f(见，i)嵌示第f个原子在 电子密度为Ph』处的嵌入能，qo(ro)表示第f个和第，个原子问的互相作用对势。 势，即vC势【7】’它的敢体对势的表达式如F： 妒(r)=DM{1一exp[-口肼(，．一RM)B2一D肼(1．3) 也子密度的表达式为； ，(，)=，-6e一步+292-2黟) (1．4) 模拟中，对势和电子密度需采用如下方式进行光滑处理： 一附M‰，+(讣圳(垃∽5， 中各参数取值见表1。VC势中嵌入能是通过拟台Rose方程得到。根据Rose方程，利 用合金的对势和电子密度函数并结合(1．1)式可以拟合出台金的嵌入能。