Phase field modeling of irradiation-induced shrinkage fracture in TRISO fuel particle

• Published in: Journal of nuclear materials 2024-04, Vol.592, p.154963, Article 154963
• Main Authors: Tan, Jiatong, Wu, Yingwei, Li, Qun, He, Yanan, Yin, Chunyu
• Format: Article
• Language: English
• Subjects: Crack growth; Fracture mechanics; Nuclear fuel; Phase field; TRISO
• ISSN: 0022-3115; 1873-4820
• DOI: 10.1016/j.jnucmat.2024.154963

•A phase field model of brittle fracture for the shrinkage crack within TRISO fuel particle has been presented, this model could be helpful for the study and optimization of TRISO.•The physical model of TRISO and phase field fracture model have been verified to ensure the reliability of analysis.•Through the analyses with different parameters, crack propagations and stress variations have been obtained. The interaction between stress relaxation and irradiation-induced shrinkage has a significant impact on the crack evolution, which in turn influences the failure of TRISO.•Heat and mass transfer could be coupled into this model, and the variation of temperature and the release of fission products would be studied together with the failure of TRISO. With the features in retaining fission products and maintaining structural stability, TRISO fuel particle has been widely implemented and studied in the development of advanced nuclear fuel element. During service, however, the key silicon carbide (SiC) layer of TRISO may fail due to irradiation-induced shrinkage fracture of inner pyrolytic carbon (IPyC) layer. Hence, it has great significance to analyze the fracture characteristics of TRISO for its application in engineering. This work presents a phase field model of brittle fracture for the shrinkage cracking, which is capable of simulating the propagation of IPyC crack and interfacial crack between IPyC and SiC naturally. At first, the physical model of TRISO and phase field model have been verified respectively. Subsequently, the crack evolution has been investigated under different parameters, including the Bacon Anisotropy Factor (BAF) of IPyC, fuel temperature and density of IPyC. Simulation results reveal that increasing the BAF or decreasing the density and temperature, the fracture of IPyC will be accelerated and the crack deflection may also occur subsequently. In contrast, the fracture of IPyC can be put off and relieve the stress concentration in SiC. The combination effect of shrinkage and stress relaxation from creep ultimately determines the crack evolution, which in turn influences the stress distribution within SiC.