Please wait a minute...
Acta Metall Sin  2026, Vol. 62 Issue (3): 523-531    DOI: 10.11900/0412.1961.2024.00138
Research paper Current Issue | Archive | Adv Search |
Phase-Field Simulation of Grain Boundary Groove Formation
SHEN Wenlong1, LIAO Yuxuan1, WU Xuezhi2, JIANG Yanbo1, LIU Wenbo1,3()
1.School of Nuclear Science and Technology, Xi'an Jiaotong University, Xi'an 710049, China
2.Institute of Reactor Engineering and Technology, China Institute of Atomic Energy, Beijing 102413, China
3.Shaanxi Key Laboratory of Advanced Nuclear Energy and Technology, Xi'an Jiaotong University, Xi'an 710049, China
Cite this article: 

SHEN Wenlong, LIAO Yuxuan, WU Xuezhi, JIANG Yanbo, LIU Wenbo. Phase-Field Simulation of Grain Boundary Groove Formation. Acta Metall Sin, 2026, 62(3): 523-531.

Download:  HTML  PDF(1289KB) 
Export:  BibTeX | EndNote (RIS)      
Abstract  

The grain boundary grooves formed on material surface significantly affect the mechanisms of mass transport, altering the driving force for grain boundary migration, and consequently influencing the kinetics of grain growth. Herein, a phase-field model is developed to simulate the coupled evolution of grain boundary grooving and grain growth in UO2 ceramic fuels. The process of developing the model begins by constructing a free energy equation for a polycrystalline system, where the reduction in total free energy drives system evolution. The mobility coefficients are then adjusted to ensure that the formation of grain boundary grooves and grain growth follow different kinetic processes. Results show that groove evolution at individual grain boundary is controlled by surface diffusion, aligning closely with Mullins' classical theoretical solution. Moreover, with the movement of grain boundaries, the groove contours become decreasingly symmetric, with increased material accumulation occurring in the grain in front of the moving direction of the grain boundary. The evolution results of a 3D polycrystalline thin film structure show that all grains evolve toward a columnar crystal structure, with grooves gradually forming at the intersections of grain boundaries and the surface. Moreover, grooves generated by different grain boundaries within the same grain tend to overlap on the grain surface, altering the groove contours. With increasing groove depth, grain boundary movement becomes increasingly slow, ultimately reducing the rate of grain growth.

Key words:  phase-field simulation      grain boundary groove      surface      grain growth      diffusion     
Received:  07 May 2024     
ZTFLH:  TG148  
Fund: National Natural Science Foundation of China(12375258);China Postdoctoral Science Foundation(2019M663738);State Key Laboratory of New Ceramic and Fine Processing Tsinghua University(KF201713);Innovative Scientific Program of China National Nuclear Corporation
Corresponding Authors:  LIU Wenbo, professor, Tel: (029)82668648, E-mail: liuwenbo@xjtu.edu.cn

URL: 

https://www.ams.org.cn/EN/10.11900/0412.1961.2024.00138     OR     https://www.ams.org.cn/EN/Y2026/V62/I3/523

Fig.1  Schematic of phase field variables (ρ—concentration field variable, η1—order parameter of grain 1, η2—order parameter of grain 2)
ParameterSymbolValue
Bulk free energy coefficientsA1.0
B0.33
Z3
Gradient energy coefficient of ρκρ1.0
Gradient energy coefficient of ηiκη0.33
Surface mobility parameterMS1.0
Grain boundary mobility parameterMGB10-4
Bulk mobility parameterMB10-4

Grain boundary relaxation coefficient

Surface relaxation coefficient

LGB

LS

1.0

10LGB

Table 1  Non-dimensional parameters used in the present simulation
Fig.2  Microstructure evolution (a1-d1) and corresponding local magnifications (a2-d2) of a groove of a stationary grain boundary (The blue and yellow parts represent the vapour phase and the grains, respectively)
(a1, a2) initial microstructure (b1, b2) 25 × 105 steps
(c1, c2) 50 × 105 steps (d1, d2) 100 × 105 steps
Fig.3  Temporal evolution of the groove depth and the fitting result (t—time)
Fig.4  Comparison of the normalized groove profile and Mullins' theory[7]
Fig.5  Comparison of the groove profile of a migrating grain boundary and Mullins' theory[8]
Fig.6  Microstructures (a1-c1) and grain size distributions (a2-c2) of three-dimensional polycrystalline structure for t = 0 (a1, a2), t = 400 (b1, b2), and t = 1000 (c1, c2) (The blue parts represent the surface and the grain boundaries, and the gray parts represent the grains, the same in Fig.7)
Fig.7  Temporal evolution of the film surface profile
(a) selected surface location (b) surface profile of the selected area
[1] Akyildiz O, Ogurtani T O. Thermal grooving by surface diffusion: A review of classical thermo-kinetics approach [J]. Hittite J. Sci. Eng., 2017, 4: 7
doi: 10.17350/HJSE19030000000
[2] Coble R L. Thermal grooving in polycrystalline ceramics and initiation of brittle fracture [J]. J. Am. Ceram. Soc., 1971, 54: 59
doi: 10.1111/jace.1971.54.issue-1
[3] Haremski P, Epple L, Wieler M, et al. Grain boundary grooving by surface diffusion in nickel bicrystals [J]. Acta Mater., 2022, 241: 118334
doi: 10.1016/j.actamat.2022.118334
[4] Sun Z Y, Yang C, Liu W B. Phase field simulations of the sintering process of UO2 [J]. Acta Metall. Sin., 2020, 56: 1295
孙正阳, 杨 超, 柳文波. UO2烧结过程的相场模拟 [J]. 金属学报, 2020, 56: 1295
doi: 10.11900/0412.1961.2019.00440
[5] Sun Z Y, Wang Y T, Liu W B. Phase-field simulation of the interaction between pore and grain boundary [J]. Acta Metall. Sin., 2020, 56: 1643
doi: 10.11900/0412.1961.2020.00120
孙正阳, 王昱天, 柳文波. 气孔与晶界相互作用的相场模拟 [J]. 金属学报, 2020, 56: 1643
doi: 10.11900/0412.1961.2020.00120
[6] Jiang Y B, Liu W B, Sun Z P, et al. Phase-field simulation of void evolution in UO2 under applied stress [J]. Acta Phys. Sin., 2022, 71: 026103
姜彦博, 柳文波, 孙志鹏 等. 外加应力作用下UO2中空洞演化过程的相场模拟 [J]. 物理学报, 2022, 71: 026103
[7] Mullins W W. Theory of thermal grooving [J]. J. Appl. Phys., 1957, 28: 333
doi: 10.1063/1.1722742
[8] Mullins W W. The effect of thermal grooving on grain boundary motion [J]. Acta Metall., 1958, 6: 414
doi: 10.1016/0001-6160(58)90020-8
[9] Robertson W M. Grain-boundary grooving by surface diffusion for finite surface slopes [J]. J. Appl. Phys., 1971, 42: 463
doi: 10.1063/1.1659625
[10] Cahn J W, Hilliard J E. Free energy of a nonuniform system. I. Interfacial free energy [J]. J. Chem. Phys., 1958, 28: 258
doi: 10.1063/1.1744102
[11] Chen L Q, Zhao Y H. From classical thermodynamics to phase-field method [J]. Prog. Mater Sci., 2022, 124: 100868
doi: 10.1016/j.pmatsci.2021.100868
[12] Guo S Y, Liu W B, Yang Q C, et al. Phase-field simulation of viscous sintering [J]. Acta Metall. Sin., 2024, 60: 1691
郭松源, 柳文波, 杨庆成 等. 黏性烧结过程的相场模拟 [J]. 金属学报, 2024, 60: 1691
doi: 10.11900/0412.1961.2022.00445
[13] Verma M, Sugathan S, Bhattacharyya S, et al. Effect of concurrent thermal grooving and grain growth on morphological and topological evolution of a polycrystalline thin film: Insights from a 3D phase-field study [J]. Acta Mater., 2023, 261: 119393
doi: 10.1016/j.actamat.2023.119393
[14] Verma M, Sugathan S, Bhattacharya S, et al. A computational analysis of universal behavior of thermal groove in a moving grain boundary [J]. Scr. Mater., 2022, 209: 114383
doi: 10.1016/j.scriptamat.2021.114383
[15] Verma M, Mukherjee R. Grain growth stagnation in solid state thin films: A phase-field study [J]. J. Appl. Phys., 2021, 130: 025305
[16] Mukherjee R, Abinandanan T A. Grain boundary grooving in a polycrystalline thin film: A phase-field study [J]. Comput. Mater. Sci., 2022, 203: 111076
doi: 10.1016/j.commatsci.2021.111076
[17] Chakraborty S, Kumar P, Choudhury A. Phase-field modeling of grain-boundary grooving and migration under electric current and thermal gradient [J]. Acta Mater., 2018, 153: 377
doi: 10.1016/j.actamat.2018.04.019
[18] Zhao Y H. Understanding and design of metallic alloys guided by phase-field simulations [J]. npj Comput. Mater., 2023, 9: 94
doi: 10.1038/s41524-023-01038-z
[19] Mukherjee R, Chakrabarti T, Anumol E A, et al. Thermal stability of spherical nanoporous aggregates and formation of hollow structures by sintering—A phase-field study [J]. ACS Nano, 2011, 5: 2700
doi: 10.1021/nn103036q pmid: 21391705
[20] Ahmed K, Yablinsky C A, Schulte A, et al. Phase field modeling of the effect of porosity on grain growth kinetics in polycrystalline ceramics [J]. Modell. Simul. Mater. Sci. Eng., 2013, 21: 065005
[21] Yang Q C, Gao Y X, Kirshtein A, et al. A free-energy-based and interfacially consistent phase-field model for solid-state sintering without artificial void generation [J]. Comput. Mater. Sci., 2023, 229: 112387
doi: 10.1016/j.commatsci.2023.112387
[22] Moelans N, Blanpain B, Wollants P. Quantitative analysis of grain boundary properties in a generalized phase field model for grain growth in anisotropic systems [J]. Phys. Rev., 2008, 78B: 024113
[23] Ahmed K, Pakarinen J, Allen T, et al. Phase field simulation of grain growth in porous uranium dioxide [J]. J. Nucl. Mater., 2014, 446: 90
doi: 10.1016/j.jnucmat.2013.11.036
[24] Cahn J W. On spinodal decomposition [J]. Acta Metall., 1961, 9: 795
doi: 10.1016/0001-6160(61)90182-1
[25] Bourgeois L, Dehaudt P, Lemaignan C, et al. Pore migration in UO2 and grain growth kinetics [J]. J. Nucl. Mater., 2001, 295: 73
doi: 10.1016/S0022-3115(01)00502-5
[26] Joshi C, Abinandanan T A, Mukherjee R, et al. Destabilisation of nanoporous membranes through GB grooving and grain growth [J]. Comput. Mater. Sci., 2017, 139: 75
doi: 10.1016/j.commatsci.2017.07.026
[27] Vedantam S, Patnaik B S V. Efficient numerical algorithm for multiphase field simulations [J]. Phys. Rev., 2006, 73E: 016703
[28] Gruber J, Ma N, Wang Y, et al. Sparse data structure and algorithm for the phase field method [J]. Modell. Simul. Mater. Sci. Eng., 2006, 14: 1189
doi: 10.1088/0965-0393/14/7/007
[29] Fan D, Chen L Q. Computer simulation of grain growth using a continuum field model [J]. Acta Mater., 1997, 45: 611
doi: 10.1016/S1359-6454(96)00200-5
[30] Herring C. Effect of change of scale on sintering phenomena [J]. J. Appl. Phys., 1950, 21: 301
doi: 10.1063/1.1699658
[31] Wang N, Wen Y H, Chen L Q. Pinning force from multiple second-phase particles in grain growth [J]. Comput. Mater. Sci., 2014, 93: 81
doi: 10.1016/j.commatsci.2014.06.030
[32] Qi X Y, Liu W B, He Z B, et al. Phase-field simulation of the densification process during sintering of UN nuclear fuel [J]. Acta Metall. Sin., 2023, 59: 1513
doi: 10.11900/0412.1961.2022.00182
戚晓勇, 柳文波, 何宗倍 等. UN核燃料烧结致密化过程的相场模拟 [J]. 金属学报, 2023, 59: 1513
doi: 10.11900/0412.1961.2022.00182
[33] Guo Q W, Hou H, Wang K L, et al. Coalescence of Al0.3CoCrFeNi polycrystalline high-entropy alloy in hot-pressed sintering: A molecular dynamics and phase-field study [J]. npj Comput. Mater., 2023, 9: 185
doi: 10.1038/s41524-023-01139-9
[34] Vandermeer R A, Hu H. On the grain growth exponent of pure iron [J]. Acta Metall. Mater., 1994, 42: 3071
doi: 10.1016/0956-7151(94)90404-9
[35] Hackney S A, Ojard G C. Grain boundary grooving at finite grain size [J]. Scr. Metall., 1988, 22: 1731
doi: 10.1016/S0036-9748(88)80274-6
[1] BAO Ergen, WANG Jiantao, XU Wenjing, MA Hui, CHEN Xing-qiu. First-Principles Study on the Influence of Alloying Elements on the Dissolution and Diffusion of Hydrogen Atoms in α-Fe[J]. 金属学报, 2026, 62(3): 489-496.
[2] JIANG Guolong, ZHOU Xia. Molecular Dynamics Simulation on the Interfacial Behavior of Mg/Mg-Zn-Al Brazing with Ultrasonic Assistance[J]. 金属学报, 2026, 62(2): 372-382.
[3] YANG Xu, YANG Zhenwen, WANG Ying, LI Huijun, LI Yongbing. Interfacial Microstructure and Mechanical Properties of Low-Temperature Diffusion-Bonded Zr-4 Alloy with Gradient Nanostructure via Ultrasonic Impact Treatment[J]. 金属学报, 2026, 62(1): 159-172.
[4] YANG Fan, PEI Shichao, LUO Xinrui, CHEN Yuxiang, LI Ningyu, CHANG Yongqin. Microstructure Evolution and Mechanical Properties of 6061 Aluminum Alloy Fabricated by Friction Stir Additive Manufacturing[J]. 金属学报, 2025, 61(8): 1129-1140.
[5] ZHANG Mingchuan, XU Qinsi, LIU Yi, CAI Yusheng, MU Yiqiang, REN Dechun, JI Haibin, LEI Jiafeng. Effect of Hot-Pressing Temperature on the Microstructure and Properties of the Diffusion-Bonded Region of TC4 Alloy[J]. 金属学报, 2025, 61(8): 1183-1192.
[6] XIE Xinliang, ZHOU Liping, YU Jianbo, XUAN Weidong, CHEN Chaoyue, WANG Jiang, REN Zhongming. Effect of Weak Transverse Magnetic Field on the Competitive Grain Growth of Ni-Based Superalloy with Divergent Bi-Crystals[J]. 金属学报, 2025, 61(8): 1203-1216.
[7] TANG Sifan, WEI Jingjing, YUE Yixin, LI Pengyu, YAO Man, WANG Xudong. Field-Variable Diffusion Cellular Automaton Model for Dendritic Growth with Sixfold Symmetry Alloys[J]. 金属学报, 2025, 61(6): 941-952.
[8] GUO Xingxing, SHUAI Meirong, CHU Zhibing, LI Yugui, XIE Guangming. Microstructure Evolution Near Interface and the Element Diffusion Dynamics of the Composite Stainless Steel Rebar[J]. 金属学报, 2025, 61(2): 336-348.
[9] DU Wenli, HOU Chao, LI Yurong, HAN Tielong, SONG Xiaoyan. Effect of Cr and Sc on High-Temperature Stability of Grain Structure in W-Based Alloys[J]. 金属学报, 2025, 61(11): 1664-1672.
[10] WANG Haifeng, PU Zhenxin, ZHANG Jianbao. Non-Equilibrium Interface Dynamics Theory[J]. 金属学报, 2025, 61(1): 29-42.
[11] WANG Binshan, XU Guang, REN Rui, ZHANG Qiang, SHAN Zhaohui, FAN Jianfeng. Effect of Electropulse on Dynamic Precipitation and Microstructure of AZ91 Magnesium Alloy During Warm Extrusion[J]. 金属学报, 2025, 61(1): 129-142.
[12] MA Fang, LU Xingyu, ZHOU Lina, DU Ningyu, LEI Chengshuai, LIU Hongwei, LI Dianzhong. High-Temperature Decomposition Mechanism of M2C Primary Carbide in M50 Steel[J]. 金属学报, 2024, 60(7): 901-914.
[13] XU Zengguang, ZHOU Shiqi, LI Xiao, LIU Zhiquan. Excellent Oxidation Resistance and Solder Wettability of (111)-Oriented Nanotwinned Cu[J]. 金属学报, 2024, 60(7): 957-967.
[14] LI Cong, WANG Meng, TU Hanjun, SHI Liqun. Behavior of Hydrogen Absorption and Desorption in Y-doped Ti Films[J]. 金属学报, 2024, 60(6): 826-836.
[15] XIONG Yi, LUAN Zewei, MA Yunfei, LI Yong, ZHA Xiaoqin. Effect of Surface Nanocrystallization Induced by Supersonic Fine Particles Bombardment on Corrosion Fatigue Behavior of 300M Steel[J]. 金属学报, 2024, 60(5): 627-638.
No Suggested Reading articles found!