基于自适应有限元枝晶自由生长相场模型模拟

兰州理工大学学报 ›› 2020, Vol. 46 ›› Issue (1): 24-31.

基于自适应有限元枝晶自由生长相场模型模拟

朱昶胜^1,2, 徐升¹, 雷鹏¹, 韩丹¹

1.兰州理工大学计算机与通信学院, 甘肃　兰州 730050;
2.兰州理工大学省部共建有色金属先进加工与再利用国家重点实验室, 甘肃兰州 730050

收稿日期:2018-09-27 出版日期:2020-02-28 发布日期:2020-06-23
作者简介:朱昶胜(1972-),男,甘肃秦安人,教授,博导.
基金资助:
国家自然科学基金(51661020,11504149,11364024)

Simulation of dendrite free growth with phase-field model based on adaptive finite element method

ZHU Chang-sheng^1,2, XU Sheng¹, LEI Peng¹, HAN Dan¹

1. College of Computer and Communication, Lanzhou Univ. of Tech. , Lanzhou 730050, China;
2. State Key Laboratory of Advanced Processing and Recycling of Nonferrous Metals, Lanzhou Univ. of Tech. , Lanzhou 730050, China

Received:2018-09-27 Online:2020-02-28 Published:2020-06-23

摘要/Abstract

摘要： 利用相场模型模拟了过冷熔体中单个完整等轴枝晶在凝固过程中的生长和形貌演化,采用自适应有限元方法求解相场模型的控制方程,研究了在计算域较大、界面层厚度较薄的情况下各向异性和过冷度等物理参数对枝晶形状和生长的影响.结果表明:过冷度越大,枝晶间竞争越激烈,生长速度越慢;各向异性系数越大,枝晶沿选定方向生长的趋势越强,生长速度越快,二次枝晶间距越小,侧向分枝越发达. 与有限差分方法(FDM)相比,采用自适应有限元法(AFM)在CPU计算时间和存储空间均降低了一个数量级,并且系统尺寸越大,自适应有限元法优势越明显,便于更大尺度多场耦合相场模型的数值模拟.通过比较均匀网格法和以往的晶体生长数值模拟实验,证明了自适应有限元法的准确性、高效性和鲁棒性.

关键词: 等轴枝晶, 自适应有限元, 数值模拟, 相场模型

Abstract: Phase-field model is used to simulate the growth and morphological evolution of single intact equiaxed dendrite during its solidification in supercooled melt, the governing equations of the phase-field model are solved with adaptive finite element method, and the influence of physical parameters such as anisotropy and supercooling degree on the dendrite shape and growth under condition of larger computational domain and less interfacial layer thickness. The result shows that a competition will occur more intensively among the dendrites and their growth speed will be less when the supercooling degree is greater. However, the greater the anisotropy coefficient is, the stronger the dendrite growth tendency along a chosen direction will be, the faster the growth speed will be, the less the secondary dendrite spacing will be, and the more developed the lateral branch will be. Compared with the finite difference method (FDM), the adaptive finite element method (AFM) will reduce CPU time and storage space by one order of magnitude, and the larger the system size is, the more obvious the superiority of adaptive finite element method will be, which will facilitate the numeric simulation of greater dimensional and multi-field coupled phase field. By means of comparing to FDM and previous numerical simulation of crystal growth, the AFM is proved to be accurate, efficient and robust.

Key words: equiaxial dendrite, adaptive finite element, numerical simulation, phase-field model

中图分类号:

TG244.3

朱昶胜, 徐升, 雷鹏, 韩丹. 基于自适应有限元枝晶自由生长相场模型模拟[J]. 兰州理工大学学报, 2020, 46(1): 24-31.

ZHU Chang-sheng, XU Sheng, LEI Peng, HAN Dan. Simulation of dendrite free growth with phase-field model based on adaptive finite element method[J]. Journal of Lanzhou University of Technology, 2020, 46(1): 24-31.

参考文献

[1] 朱昌盛,王军伟,王智平,等.受迫流动下的枝晶生长相场法模拟 [J].物理学报,2010,59(10):7417-7423.
[2] 朱昶胜,徐升,冯力,等.Multi-GPU加速的二元合金定向凝固三维相场模型 [J].兰州理工大学学报,2018,44(6):24-29.
[3] KOBAYASHI R.Modeling and numerical simulations of dendritic crystal growth [J].Physica D,1993,63(3/4):410-423.
[4] WHEELER A A,MURRAY B T,SCHAEFER R J.Computation of dendrites using a phase field model [C]//Twelfth International Conference of the Center for Nonlinear Studies on Nonlinearity in Materials Science.[S.l.]:Elsevier Science Publishers B.V.1993:243-262.
[5] WANG S L,SEKERKA R F,WHEELER A A,et al.Thermodynamically consistent phase-field models for solidification [J].Physica D-nonlinear Phenomena,1993,69(12):189-200.
[6] PROVATAS N,GOLDENFELD N,DANTZIG J.Adaptive mesh refinement computation of solidification microstructures using dynamic data structures [J].Journal of Computational Physics,1999,148(1):265-290.
[7] PROVATAS N.Efficient computation of dendritic microstructures using adaptive mesh refinement [J].Physical Review Letters,1998,80(15):3308-3311.
[8] 陈云,康秀红,李殿中.自由枝晶生长相场模型的自适应有限元法模拟 [J].物理学报,2009,58(1):390-397.
[9] 朱昶胜,雷鹏,冯力,等.纯物质过冷熔体中自由枝晶生长相场模型的自适应有限元法模拟 [J].兰州理工大学学报,2014,40(4):18-23.
[10] 朱昶胜,雷鹏,肖荣振,等.基于自适应有限元相场模型模拟强制流动条件下的枝晶生长 [J].中国有色金属学报(英文版),2015,10(1):241-248.
[11] 朱昶胜,李椿茂,冯力,等.基于自适应有限元法的纯物质三维枝晶生长模拟 [J].兰州理工大学学报,2015,41(4):1-5.
[12] LI R.On multimesh hadaptive method [J].Journal of Computing,2005,24(3):321-341.
[13] LAN C W,CHANG Y C,SHIH C J.Adaptive phase field simulation of nonisothermal free dendritic growth of a binary alloy [J].Acta Materialia,2003,51(7):1857-1869.
[14] HOYT J J,ASTA M,KARMA A.Atomistic and continuum modeling of dendritic solidification [J].Materials Science & Engineering R,2003,41(6):121-163.
[15] 陈云.金属单相凝固组织演化的定量相场模拟与实验验证[D].沈阳:中国科学院金属研究所,2012.