研究论文

Al 液诱导凝固过程的分子动力学模拟

展开
  • 上海大学 材料科学与工程学院 省部共建高品质特殊钢冶金与制备国家重点实验室 上海市钢铁冶金新技术开发应用重点实验室, 上海 200444

收稿日期: 2018-04-13

  网络出版日期: 2018-12-21

基金资助

国家自然科学基金资助项目(51374141)

Induced solidification of liquid Al from molecular dynamics simulation

Expand
  • State Key Laboratory of Advanced Special Steel, Shanghai Key Laboratory of Advanced Ferrometallurgy, School of Materials Science and Engineering, Shanghai University, Shanghai 200444, China

Received date: 2018-04-13

  Online published: 2018-12-21

摘要

采用分子动力学(molecular dynamics,MD)方法模拟了嵌入不同半径(0.37~2.4nm)固态 Al 晶核的 Al 液诱导凝固过程. 结果表明,体系临界温度和嵌入晶核半径的倒数呈线性关系,符合吉布斯-汤姆森(Gibbs-Thomson, G-T)效应,得到的 G-T 系数 Γ 为 1.4*10-7 Km,体材料熔点 Tmbulk  为 (985.36±11.25) K.借助 Γ 计算出固液界面能为 (140.35±9.05)mJ/m2, 与通过毛细波动法获得的 149 mJ/m2 比较接近,再次证实了 Turnbull 的实验值(93 mJ/m2)$较低.外推得到诱导凝固失效半径为 0.91 nm,此时体系温度已达到临界温度极限值. 在该极限值下, 体系可以自发形核,孕育时间呈现随机性, 微观结构表现为不稳定的层错互相交织.在有效晶核尺寸范畴内, 孕育时间随嵌入晶核半径的增大而增加,长大速度却随半径的增大而降低, 微观结构是稳定的 lamellar 结构.

本文引用格式

余熔刚, 赖琴梅, 汪昊, 吴永全 . Al 液诱导凝固过程的分子动力学模拟[J]. 上海大学学报(自然科学版), 2020 , 26(2) : 216 -226 . DOI: 10.12066/j.issn.1007-2861.2038

Abstract

Molecular dynamics (MD) simulation was applied to analyze the solidification of liquid Al induced by embedded solid nanoparticles with radii ranging from 0.37 to 2.4 nm. It was found that the critical temperature was proportional to the inverse of the nucleus radius, i.e., Gibbs-Thomson (G-T) effect, through which the G-T coefficient Γ and the bulk melting temperature Tmbulk were obtained, 1.41.4*10-7 Km and (985.36±11.25) K, respectively. The solid-liquid interfacial energy which was estimated (140.35±9.05)mJ/m2 from Γ, was very close to the calculated 149 mJ/m2 by capillary fluctuation method. It proved again that Turnbull underestimated this quantity experimentally (93 mJ/m2). The critical failure radius was extrapolated to be 0.91 nm. Meanwhile, the corresponding system temperature had reached the lower limit for the critical temperature, under which the system could be spontaneously nucleated, and the incubation time was somewhat random. Besides, the microstructure was characterized by the metastable cross stacking faults. When the embedded nucleus could serve as the nucleation core, the incubation time increased with increasing embedded nucleus radius. However, the growth rate decreased separately with the increasing embedded nucleus radius. The microstructure adopted the relatively stable lamellar structure.

参考文献

[1] Turnbull D . Formation of crystal nuclei in liquid metals[J]. Journal of Applied Physics, 1950,21(10):1022-1028.
[2] Espinosa J R, Vega C, Valeriani C , et al. Seeding approach tocrystal nucleation[J]. Journal of Chemical Physics, 2016,144(3):1-10.
[3] Bai X M, Li M . Calculation of solid-liquid interfacial freeenergy: a classical nucleation theory based approach[J]. Journal ofChemical Physics, 2006,124(12):1-12.
[4] Broughton J Q, Gilmer G H . Molecular dynamics investigation ofthe crystal-fluid interface. VI. Excess surface free energies ofcrystal-liquid systems[J]. Journal of Chemical Physics, 1986,84(10):5759-5768.
[5] Hoyt J J, Asta M, Karma A . Method for computing the anisotropyof the solid-liquid interfacial free energy[J]. Physical ReviewLetters, 2001,86(24):5530-5533.
[6] Shibuta Y, Suzuki T . A molecular dynamics study of the phasetransition in BCC metal nanoparticles[J]. Journal of ChemicalPhysics, 2008,129(14):1-10.
[7] Watanbe Y, Shibuta Y, Suzuki T . A molecular dynamics study ofthermodynamic and kinetic properties of solid-liquid interface forBCC iron[J]. ISIJ International, 2010,50(8):1158-1164.
[8] Hashimoto R, Shibuta Y, Suzuki T . Estimation of solid-liquidinterfacial energy from Gibbs-Thomson effect: a molecular dynamicsstudy[J]. ISIJ International, 2011,51(10):1664-1667.
[9] Xia Y, Li C H, Luan Y W , et al. Molecular dynamics studies onthe correlation of undercoolability and thermophysical properties ofliquid Ni-Al alloys[J]. Computational Materials Science, 2016,1112(1):383-394.
[10] Wu Y Q, Shen T, Lu X M , et al. Solidification of liquid Fe withembedded homogeneous solid Fe nanoparticles from molecular dynamicssimulations[J]. Acta Physico-Chimica Sinica, 2013,29(2):245-249.
[11] Plimpton S . Fast parallel algorithms for short-range moleculardynamics[J]. Journal of Computational Physics, 1995,117(1):1-19.
[12] Mendelev M I, Han S, Srojovitz D J , et al. Development of newinteratomic potentials appropriate for crystalline and liquid iron[J]. Philosophical Magazine, 2003,83(35):3977-3994.
[13] Jiang Y W, Luo J, Wu Y Q . The validation and preference amongdifferent EAM potentials to describe the solid-liquid transition ofaluminum[J]. Modelling and Simulation in Materials Science andEngineering, 2017,25(4):1-13.
[14] Li R, Wu Y Q, Xiao J J . The nucleation process and the roles ofstructure and density fluctuations in supercooled liquid Fe[J]. Journal of Chemical Physics, 2014,140(3):1-11.
[15] Zhou H G, Lin X, Wang M , et al. Calculation of crystal-meltinterfacial free energies of FCC metals[J]. Journal of CrystalGrowth, 2013,366(3):82-87.
[16] Liu J, Dong H B . Molecular dynamics calculation ofthermodynamic properties of iron solidification[J]. IOP ConferenceSeries: Materials Science and Engineering, 2012,33(1):1-10.
[17] Chase M W, Curnutt J L, Downey J R , et al. JANAF thermochemicaltables, 1982 supplement[J]. Journal of Physical and ChemicalReference Data, 1982,11(3):695-940.
[18] Kelton K F . Crystal nucleation in liquids and glasses[J]. Solid State Physics, 1991,45(1):75-177.
[19] Morris J R . Complete mapping of the anisotropic free energy ofthe crystal-melt interface in Al[J]. Physical Review B, 2002,66(14):1-7.
[20] Morris J R, Mendelev M I, Srojovitz D J . A comparison ofcrystal-melt interfacial free energies using different Al potentials[J]. Journal of Non-Crystalline Solids, 2007,353(32):3565-3569.
文章导航

/