钢液中纳米夹杂物的分子动力学模拟介绍

钢液中纳米夹杂物的分子动力学模拟姓名:韩剑学号：201403047一、宏观动力学的发展历程冶金物理化学是冶金工艺发展的理论基础，而冶金过程宏观动力学是冶金物理化学的一个重要内容,其研究方向和水平同冶金工业的发展密切相关。回溯冶金工业近40年的发展历程，大型设备(诸如高炉、密闭鼓风炉、转炉、精炼炉等)出现后，冶金生产发展到了高效率、高产量、成熟、完整的工艺阶段。冶金科技工作者始终把注意力集中在“强化现有冶炼工艺”上。研究氧与铁液中各元素的反应，从热力学证实降低一氧化碳分压，提高了碳的脱氧能力和固体碳的还原能力。于是,50年代在工艺上相继出现真空冶金、真空处理、真空碳热还原提取金属。元素氧化热力学规律的研究，实现了返回吹氧法冶炼不锈钢及纯氧炼钢，强化了冶金工艺过程。研究液-液分配平衡的成果，导致在湿法冶金中出现液相萃取冶金。60年代，为适应强化了的冶金工艺要求，对熔体反应动力学、冶金电化学以及快速测试方法进行了大量的研究工作。70年代，出现了冶金上三大发明，即炉外精炼、喷吹冶金和浓差电池定氧。此外,高效节能的电极材料相继出现。于是冶炼时间缩短，工艺强化，能耗降低，产品质量提高。80年代科学技术迅速发展，新技术倍出，计算机普及、激光、等离子体、遗传工程等迅速发展并开始向冶金学科渗透。于是开始利用遗传工程研究生物冶金，应用等离子体研究超高温下冶金反应的物理化学规律，利用低能核物理技术研究金属表面与界面的物理化学等等。较为突出的是用计算机模拟冶金反应过程和冶金系统工程，开始探索新的冶金工艺。如熔融态还原与快速凝固配套的冶金生产全流程的设想，在国内外均已开始进行具有相当规模的试验或半工业试验研究。传统冶金工艺的改造要依靠冶金物理化学理论的研究，而冶金新工艺的探素尤其离不开冶金物理化学研究的先行。大量研究和实践证明，分子动力学(MoleculardynamicsMD)模拟是较为成功的研究方法，通过对物质微观结构和运动过程进行数值模拟和统计求和，可获得物质结构、粒子运动规律以及微观结构与宏观性质之间关系等信息。因此，该方法比较广泛地用于原子分子物理、凝聚态物理、材料科学和纳米科技等领域,而目前在冶金领域，特别是钢铁冶金的研究中应用还很少。由于MD研究不受样品制备和测试技术限制，应用于钢铁冶金领域，可弥补高温熔体结构和性质检测的困难，促进对冶金过程微、介观尺度行为的研究。本文概述了MD模拟的基本原理与技术，总结了MD模拟在冶金领域中的应用研究现状，并对采用MD方法模拟研究钢液中夹杂物在纳米尺度下形核与生长的机理进行了探讨。二、分子动力学在冶金研究中的应用进展目前，分子动力学模拟方法在冶金，特别是钢铁冶金领域应用很少。针对硅酸盐炉渣体系，吴永全等⑴釆用MD方法模拟计算了CaO・SiO?、Al2O3-SiO2sCaO-Al2O3和CaO-Al。一SiO?的结构；同时，釆用高温拉曼检测对比分析。研究认为,常压下硅酸盐炉渣中Si以非常稳定的0O44•四面体结构存在，并得到了Si・O键长和键角等相关数据；Ca在渣中起破坏网格作用，A1有助于形成网格，主要以AIOF四面体存在，但稳定性不如SioZ，并模拟得到A1O右和Siof的差异；同时，模拟认为渣中存在桥氧6、非桥氧Onb和自由氧Of,并获得了3类氧之间的平衡关系式。张毅刚卩］利用MD研究计算出Mg2SiO4的粘滞度。谢刚等⑶采用MD方法模拟(LaFBZLiF)"熔盐体系的结构,认为在(LaF3)x(LiF)i.x熔盐中La・F配位数主要为7,且比较稳定，波动范围在7.33〜7.75之间，并计算了F・La・F键角范围；锂离子是电流的主要携带者，其次是氟离子;LaFs浓度增大会降低各离子扩散系数，使正负离子分布更加松散。同时，对LaCls-KCl体系的MD模拟表明,熔体中搁氯络离子很稳定，La-Cl配位数在6.75〜7.75之间。侯怀宇等⑷用MD方法分别模拟了稀土金属La、Gd、Y的氯化物的熔盐结构，认为3种熔盐都以配位数为6的正八面体结构存在，键角分布在8070°之间；LaCh在熔盐中正负离子分布较为混乱松散，配位数较高。在金属和液态金属方面，贾瑜等⑴通过嵌入型原子间相互作用势的MD模拟，计算了理想体心立方结构的金属Fe［001］晶带的表面能随角度变化曲线，得到了更具物理意义的公式，并由此推出表面能与原子的紧密程度成反比。韩秀峰回利用自编的非平衡动力学MD程序，对液态金属Co进行模拟计算，成功得出剪应力使体系的无序度增大的结论，并得到其粘度和激活能的值；同时，采用儿种作用势对液态A1的模拟计算都得到非常接近的粘度值，证明了MD计算熔体粘度的可行性。王荣山等卩］模拟液态Cu的等温凝固过程，得出低温下的等温凝固可得到非晶结构Cu,高温下等温凝固得到过冷液体结构，晶态结构只有在适当温度才能得到，其三维结构取决于Cu原子两个方向的扩散和弛豫。赵毅等［可用MD对金属熔体Ni?Al形核进行了模拟，发现非晶团簇不参与形核过程，晶核为不规则的fee结构和hep结构的混合体。丛红口等叨通过MD证明液态Al8oFe2o中存在较强化学序，导致该熔体呈现一定的有序性。从上述研究可以看出，MD方法是一种较为有效的研究冶金熔体结构的计算机模拟技术。随着洁净钢冶金技术的发展，需要对冶金过程中微、介观尺度的钢液和夹杂物的结构及行为等进行研究，分子动力学方法在钢铁冶金中将有越來越广泛的应用。三、钢液中纳米夹杂物的分子动力学模拟钢中夹杂物形核与长大过程、金属熔体凝固、第二相析出过程等，均要经历纳米尺度的形核与生长阶段，以往的研究基本是基于传统的热力学、界面和动力学理论。关于纳米材料特性的研究表明，纳米尺度材料的性质受材料尺寸影响刨。因此，笔者希望在考虑纳米小尺寸效应的基础上，对夹杂物形核与生长初期过程的热力学进行MD模拟研究，进而研究钢液一纳米夹杂界面性质、夹杂形核与生长的纳米尺度动力学等问题。对钢液中夹杂物形核过程的钢液一纳米夹杂界面能

进行了研究，得到界面能与夹杂尺寸的关系为：5y(2r-3h)式中：Y3为钢液一纳米夹杂界面能，Y为钢液一宏观夹杂界面能，「为夹杂半径，h为夹杂分子直径。可见，y9随着「增大而增大，最大值为10Y/9o以A1?O3夹杂为例，钢液一宏观AI2O3界面能Y=2.30N/m,AI2O3分子直径h二4.78Xl(Tm,根据式计算出钢液一纳米AI2O3界面能与其尺寸的关系，如图1所示。当夹杂半径小于lOnm时，钢液-A12O3夹杂界面能随其尺寸增大而较快增加，尺寸越小，增加越快；当夹杂半径大于lOrnn时，界面能随夹杂尺寸增大缓慢增加，并趋于2.556N/m,与钢液一宏观A12O3界面能2.3ON/m接近。可见,当夹杂长大到大于lOnm后，尺寸对进一步长大的能量影响逐渐减小，当影响很小后即服从宏观界面能规律。2.5562.4 •••*• 2.4342.0・■■1・5・■1.0-/■0.5・0,(b*2 4 6 8 10 12 14 16r/nm图1钢液一纳米Al心界面能对尺寸的依赖关系用MD对钢液中纳米A1?O3夹杂的热容Cp进行了模拟，采用的是文献3中模拟液态AI2O3结构的Matsui型原子间作用势，因为Matsui势比较适合模拟原子作用力较强的体系。+D(B1+BJ+D(B1+BJ)exp(■A+A~ruB’+Bj)-式中：U(£j)是体系势能，等号右边依次为库仑势(Coulombterm)、近程排斥能(Repulsionteam)和范徳华力(vanDerVaalsteam),其中是粒子i和j之间的距离，乙、Z,为粒子的有效电荷，A为表征粒子半径的距离参数，B为排斥常量，C为范徳华力常量，D是标准力常数，数值为41.84KJ-nm^-mor1,A12O3模拟用的其他参数数值见表1。表1AI心模拟所用Matsumiya势函数的参数值Z/eA/nm"AimC/<nm3・kJ"?・mol-12)Al1.41750.078520.003436820O-0.94590.082150.013890610Cp随尺寸的变化规律的模拟结果如图2所示。从图2可以看出，纳米尺度夹杂的热容高于宏观夹杂物的热容，并且夹杂尺寸愈小，热容愈大。在此基础上，可进一步根据基本热力学关系推导其他基本热力学函数G、S、H,研究纳米尺度的平衡热力学、形核和生长动力学等问题，从而为在冶炼和浇注过程中发展夹杂物直接细化控制技术奠定基础；还可尝试釆用MD研究薄板钢中纳米析出物的析出过程及机理。可见，采用MD模拟是一种准确研究微、介观尺度冶金过程行为规律的有效方法。156二14B<144140I36 1—~1 1 ［二T77工■二一［_, ■0 20 40 6() «() 100<Pnin图2 1873K时纳米AI心热容Cp的尺寸依赖关系结语随着社会对高品质钢要求的日趋严格，在冶炼和浇注过程中，对钢液洁净度、夹杂物以及连铸组织进行精确控制已迫在眉睫。要实现对钢液、夹杂物及组织的精确控制，应釆用MD方法对冶金过程熔体的微、介观尺度的性质及行为进行深入研究，考虑微、介观尺度时的小尺寸效应的影响，获得对冶金过程物质变化机理（如夹杂物在纳米尺度时的性质与行为、钢液结晶形核初期过程）的准确认识。对于钢液中内生夹杂物的控制问题，可遵循在冶炼和浇注过程中进行夹杂物百接细化控制的思路，首先用MD模拟钢液一纳米夹杂界面性质规律，进而得出夹杂形核一生长的纳米尺度动力学规律，找出纳米尺度生长时的限制性环节，获得夹杂物尺寸细化控制的工艺技术。参考文献吴永全，尤静林，蒋国昌.铝酸钙熔体结构的分子动力学模拟研究[J].无机材料学报，2003,18(3)：619张毅刚，聂高众.直至核幔边界的Mg2SiO.成分熔体剪切粘滞度的平衡分子动力学研究[J].岩石学报，1998,14(4)：409谢刚，李荣兴.冶金熔体结构和性质的计算机模拟计算[M].北京：冶金工业出版社,2006：44侯怀宇，谢刚，刘国华，等.用分子动力学方法研究棉土金属氯化物熔体的结构[J].有色金属学报，2000,10(2)：270贾瑜，王晓春，姚乾凯，等.金属高指数面表面能的分子动力学模拟[J].郑州人学学报，2004,36(2)：38韩秀峰.金属熔体粘度的分子动力学模拟[D].济南：山东大学，2005王荣山，侯怀宇，陈国良.非晶Cu在晶化过程中的分子动力学模拟[J].金属学报,2009,45(6)：692赵毅，赵九洲，胡壮麟.过冷NijAl熔体形核的分子动力学模拟[J].金属学报，2008,44(10)：1157丛红口，边秀房，李喜珍，等.液态AUF"在快速冷却中的MD模拟[J].物理化学学报，2002,18(5)：414KimEun—Ha,LeeByeong—Jo