Effects of U-Mo Irradiation Creep Performance on the Thermo-Mechanical Coupling Behavior in U-Mo/Al Monolithic Fuel Assemblies

To obtain the optimized fuel performance, the effects of U-Mo irradiation creep rate coefficient on the thermo-mechanical behavior of a fuel assembly are investigated. In this study, three cases of creep rate coefficient are considered. The distribution and evolution results of temperature, displacement, stress/strain and fuel foil micro-structure are analyzed. The simulation results indicate that with the increase of creep rate coefficient 1) the temperature field in the fuel assembly changes slightly; 2) the maximum out-of-plane displacements in the side plates decrease slightly; the maximum out-of-plane displacements in the fuel plates adjacent to the outside Al plates rise distinctly, induced by the enhanced bending deformation contributions; 3) the peak values of the first principal stresses and skeleton normal stresses in the fuel foils are reduced, while the through-thickness creep strains are enlarged; 4) with an increase of fuel creep rate coefficient from 500 × 10−22 mm3/(fission MPa) to 2,000 × 10−22 mm3/(fission·MPa), the maximum Von Mises stress in the fuel cladding decreases by ∼24% on the 308th day. This work is helpful for advanced fabrication and optimization design of U-Mo/Al fuel assemblies.


INTRODUCTION
To comply with nuclear non-proliferation requirements (Hanson and Diamond, 2011;Woolstenhulme and Nielson, 2011;Robinson et al., 2013), low enriched uranium (LEU) fuels need to be utilized for advanced research and test reactors. Monolithic U-Mo fuels have been demonstrated to be attractive candidates because of their low neutron capture cross section, high uranium density, stable irradiation swelling performance (Burkes et al., 2009;Robinson et al., 2013;Mei et al., 2015;Newell et al., 2017). Rectangular fuel assemblies Wu et al., 2017) can be made by connecting a number of U-Mo/Al monolithic fuel plates (Cheng et al., 2004;Perez et al., 2011;Robinson et al., 2013;Turkoglu et al., 2019) with an Al alloy frame, using the processing method of rolling-swage. In the reactor cores, dozens of fuel assemblies will be placed into specified positions, and control rod channels are retained among the fuel assemblies . The coolant channels exist between the fuel plates to keep the normal coolant flow (Baek et al., 2015;Fan et al., 2016;Guo et al., 2018;Ma et al., 2018), with the fission heat transferred timely (Cheng et al., 2004;Baek et al., 2015;Woolstenhulme et al., 2015). Under the in-service environments, different mechanisms are involved in the deformations of fuel assemblies (Deng et al., 2017), which will result in the variations of control rod channels and coolant channels. In order to achieve the optimization of a reactor design, it is necessary to develop theoretical models and conduct numerical simulation researches on the thermalmechanical behavior of fuel assemblies, combined with experimental researches.
Complicated in-pile thermo-mechanical behavior appears in U-Mo/Al monolithic fuel assemblies (Miller et al., 2010;Kim et al., 2012;Deng et al., 2017), which mainly depends on the following contributions: 1) the deformations (Miller et al., 2010;Kim and Hofman, 2011;Kim et al., 2013;Meyer et al., 2014;Zhao et al., 2015) induced by irradiation swelling and creep of fuel foil (Kim et al., 2013;Kim and Hofman, 2011); 2) the plasticity and thermal creep performances of cladding (Yan et al., 2017;Jian et al., 2019a); 3) formation of porous fuel structure due to the fission-gas-resulted bubbles (Rest, 2010;Meyer et al., 2014); 4) the acceleration of fission gas swelling driven by grain recrystallization (Rest, 2010); 5) the degradations of thermo-mechanical properties and macroscale strength owing to the formed porous structure, pore pressure and the possible creep damage in the fuel foil (Rest, 2010;Robinson et al., 2008;Iltis et al., 2016;Yan et al., 2017;Salvato et al., 2018;Jian et al., 2019a;Schulthess et al., 2019); 6) the complex mechanical interactions across all the parts, including the fuel plates, the outside Al plates and the side plates (Deng et al., 2017). On one hand, the irradiation creep strains will relax the stresses in the fuel foils or the other parts. On the other hand, the high creep strains will result in creep damage and degradation of fuel strength (Iltis et al., 2016;Jian et al., 2019a). It is necessary to examine the effects of fuel creep performance on the thermal-mechanical behavior in U-Mo/Al fuel assemblies, which will supply a basis for advanced fuel design.
The thermo-mechanical behavior analysis for monolithic U-Mo/Al fuel plates can be found in the previous works (Kim et al., 2012;Kim et al., 2013;Jian et al., 2019a). A linear relation of equivalent creep strain rate with fission rate and Von Mises stress was generally adopted, namely _ ε cr ∝ σ · _ f (Dienst, 1977;Ghoshal et al., 2013;Yan et al., 2019). Kim et al. (2013); Yan et al. (2019) and Jian et al., (2019b) made important contributions to the predictions of fuel creep rate effects, and the researches indicated that the fuel creep performance played an important role in the thermo-mechanical behavior of fuel plates. It is noted that these simulations were focused on the mini fuel plates (Miller et al., 2010;. Recently, Mao et al. (2021) simulated the thermo-mechanical coupling behavior of a real-size U-10Mo/Al monolithic fuel assembly with a constant fuel creep rate coefficient of 2,000 × 10 −22 mm 3 /(fission·MPa). On this basis, the effects of fuel creep rate on its thermal-mechanical behavior need to be further investigated.
In this study, simulations of the thermo-mechanical coupling behavior are implemented for three fuel assemblies with respective U-Mo irradiation creep rates, and the corresponding effects on the results of temperature, deformation and stresses are obtained and analyzed, together with the influences on the fuel porosity, pore pressure and skeleton stresses (Jian et al., 2019b).

Finite Element Model
A typical fuel assembly made of sixteen U-Mo/Al monolithic fuel plates and an Al alloy frame is considered (Mao et al., 2021), as  (Mao et al., 2021).
Frontiers in Energy Research | www.frontiersin.org May 2021 | Volume 9 | Article 676881 2 shown in Figure 1A, with the fuel plates assembled well into two grooved side plates through roll-swaging (Kim et al., 2014). The dimensions of fuel plates are 332.0 mm × 68.0 mm × 1.27 mm, and those of U-Mo fuel foil are 279.4 mm × 61.34 mm × 0.216 mm (Baek et al., 2015), as illustrated in Figure 1B. The coolant channels between two neighboring plates (Baek et al., 2015) can be found in Figure 1C, with a width of 2.95 mm. The sizes of grooved side plates are set as 332.0 mm × 75.96 mm × 6.5 mm (Baek et al., 2015).
According to the fission density data  in Figure 2A, a linearly-distributed fission rate is fitted as where _ f depicts the fission rate in the unit of fission/(mm 3 s); x denotes the current coordinates in mm. The heat generation rate can be correspondingly determined (Zhao et al., 2015). It is noted that the fission rate is unvaried during the whole irradiation period.
With the symmetries in geometry, loads and boundary conditions of the considered fuel assembly, 1/4 part is established as the FE model in Figure 2B. Fixed displacement boundary conditions are applied on the two ends. In addition, the continuous conditions of displacement and temperature are considered on the formed interfaces. The temperature of coolant is set as 323 K (Cui et al., 2015), and the convection heat transfer boundary conditions are taken into account, with a heat transfer coefficient of 0.035 W/(mm 2 · K) (Cui et al., 2015). The simulations are performed on the commercial software ABAQUS, and the finite element discretization is carried out with the element C3D8RT. As shown in Figure 2C, 1,098,444 elements are adopted, which has been verified with convergent calculation results.

Governing Equations For the Temperature Field
The temperature field in the fuel assembly is governed by .
where k(T) is the thermal conductivity; _ q is the heat generation rate. It is noted that the porosity-related thermal conductivity models (Jian et al., 2019b) are adopted for the nuclear fuel foils, and the porosity varies with burnup due to the fission-induced bubbles (Jian et al., 2019a). The heat generation rate for nuclear fuel foils is location dependent, which can be obtained with Eq. 1. The thermal conductivity models for Al alloy cladding and frame are the same as that in Jian et al. (2019b).

For the Mechanical Field
In this study, large deformation is considered, and the governing equations of the mechanical field can be expressed as where σ and σ_ are the co-rotational Cauchy stress tensor and its time derivatives; v denotes the co-rotational velocity gradient tensor; d and d e represent the total deformation rate and the elastic deformation rate with respect to the current configuration; ε e depicts the co-rotational elastic logarithm strain tensor, and I refers to the unit tensor; G and λ are Lame constants, and they vary with time.
It should be mentioned that the total deformation rate consists of the deformation contributions of elasticity, irradiation creep and irradiation swelling for the fuel foil. The mechanistic fission gas swelling model has been used. The grain crystallization effect has been involved, and the bubble radius correlates with the grain-scale fission gas atom diffusion, detailed in Jian et al. (2019a); Jian et al., 2019b. For the Al alloy, the total deformation rate is composed of elasticity, creep and plasticity parts, with the models of Lame constants and creep the same as those in Jian et al. (2019b), with the plasticity model adopted from Jian et al. (2019a).
The fission-induced creep rate model is given as Yan et al. (2017).
where, _ ε c denotes the equivalent irradiation creep rate in (/s); σ represents the Von Mises stress in MPa; the fission rate of _ f uses the unit of fission/(mm 3 · s). A UMo is the creep rate coefficient. In this study, three creep rate coefficients with the values of 500 × 10 −22 mm 3 /(fission·MPa), 1,000 × 10 −22 mm 3 /(fission·MPa), and 2,000 × 10 −22 mm 3 /(fission·MPa) are adopted to evaluate the U-Mo creep rate effects on the thermo-mechanical behavior in the fuel assembly, respectively.

Solution Techniques for the Irradiation-Induced Thermal-Mechanical Fields
For the numerical implementation of the thermal-mechanical theoretical models on the commercial software ABAQUS, it is necessary to define the location-irradiation time-temperature dependent thermal-mechanical constitutive relations. The solution techniques in this study are similar to those in our previous works Gong et al., 2013;Gong et al., 2014;Zhao et al., 2014;Cui et al., 2015;Kong et al., 2018). Here, the threedimensional mechanical constitutive relation in an incremental form is briefly introduced as follows.
For a small time increment [t, t + Δt], the incremental constitutive relation for an integration point in a co-rotational coordinate system can be obtained as where Δσ ij and Δε e ij are the incremental components of Cauchy stress tensor and elastic logarithmic strain tensor in the corotational coordinate system; ΔG and Δλ are the increments of Lame constants in a small time increment; ε e(t) ij is the components of elastic logarithmic strain tensor at time t, also in the corotational coordinate system.
In ABAQUS, based on the polar decomposition (Simulia and Fallis, 2013) for the deformation gradient increments, the values of the total logarithmic strain increments are calculated out for stress update. One should develop stress update algorithms and consistent stiffness modulus, with the other strain increment contributions involved, and then the subroutines of UMAT can be programmed to define the irradiation-induced complex mechanical constitutive relations. For the U-Mo fuel foil and Al alloy, the similar solution strategy as those in our previous works Gong et al., 2013;Gong et al., 2014;Zhao et al., 2014;Cui et al., 2015;Zhao et al., 2015;Kong et al., 2018) is adopted in this study, and the developed models, algorithms and user-defined subroutines have been verified. Especially, with the mechanistic gaseous swelling model and irradiation creep model used for U-Mo fuels, the fuel foil thickness increments in a mini plate were calculated to match well with the experimental results (Jian et al., 2019a).

The Skeleton Stress Model
A porous structure of U-Mo fuel foils will be formed due to fission gas bubbles. The skeleton stress model (Jian et al., 2019a;Jian et al., 2019b) is incorporated into the finite element calculations, which depends on the bubble volume fraction ϕ, bubble radius r, surface tension c s and internal pressure of bubbles P, expressed as where σ meso is the skeleton stress in MPa; σ macro in MPa depicts the maximum normal stress on a macroscopic scale; c s equals to 1.0 × 10 −3 N/mm, and r is the bubble radius in mm. Some details can be found in our previous work (Jian et al., 2019a).

RESULTS AND DISCUSSION
In this section, the simulation results of temperature, displacement and stress/strain for three U-Mo/Al plate-type fuel assemblies with different fuel creep rate coefficients are compared and analyzed. The view cuts and output paths are displayed in Figure 3.

Effects of Fuel Creep Rate Coefficient on the Temperature and Displacement Fields
On the Temperatures in the Fuel Plates and Side Plates As the same irradiation conditions are considered, the temperature fields in eight fuel plates are almost the same. Figure 4A gives the contour plot of temperature field in the fuel foil of Plate 8, after irradiation of 308 days for a creep rate coefficient of 500 × 10 −22 mm 3 /(fission·MPa). It can be easily noted that a peak value of ∼377.8 K exists in the heavily irradiated region. To present the effect of U-Mo creep rate coefficient on the temperature field, Path 1 in Figure 3, on the fuel meat top surface of Plate 8, is chosen. Path 2 which goes through the maximum temperature point on the side plate is also selected in Figure 3. As no fission heat is generated in the outside Al plates, a temperature gradient is formed in the width direction of side plate A. As depicted in Figures 4C,D, the temperature results change slightly for the considered three creep rate cases.

On the Displacements of Side Plates
The schematic design of the reactor core in Wu et al. (2017) indicates that narrow control rod channels exist among different fuel assemblies. The irradiation-induced space variations of these channels deserve attention. The out-of-plane displacements U 2 of side plates directly affect the variations of control rod channels, and U 2 depicts the displacement component in the y-direction, as shown in Figure 3. Figure 5A depicts the contour plot of displacement U 2 in side plate A, after irradiation of 308 days for a creep rate coefficient of 500 × 10 −22 mm 3 /(fission·MPa). The maximum displacement of ∼0.084 mm occurs in Path 2 of Figure 3. In Figure 5B, one can observe that the displacements of U 2 show a trend of reduction as a whole, with the increase of creep rate coefficient. After irradiation of Frontiers in Energy Research | www.frontiersin.org May 2021 | Volume 9 | Article 676881 5 308 days, the maximum displacement is ∼0.082 mm for a creep rate coefficient of 2,000 × 10 −22 mm 3 /(fission·MPa), which is ∼2.4% smaller than that for 500 × 10 −22 mm 3 /(fission·MPa). This indicates that the resultant forces from every fuel plate are slightly changed. It should be mentioned that the displacements of side plate mainly result from the bending deformations, driven by the forces from the fuel plates (Mao et al., 2021). Compared to the initial width of 1 cm (Cheng et al., 2004), the effects of out-of-plane displacements in the side plates can be ignored. It is noted that the displacement results are obtained with the two ends of fuel assembly constrained. The effects of creep rate coefficient might be varied under some other constraint conditions.

On the Displacements of Fuel Plates
To avoid flow blockage accident in research reactors (Baek et al., 2015), enough spaces of coolant channels must be guaranteed. As shown in Figure 1, the positions of bottom and top surfaces of fuel plates and outside Al plates can directly reflect the configuration of coolant channels. It should be noted that U 3 denotes the displacement component in the z-direction, as displayed in Figure 3, which can be used to analyze the space variations of coolant channels. As displayed in Mao et al. (2021), the maximum out-of-plane displacements of Plate 8 are much larger than those of the other numbered fuel plates. In Figure 6, the contour plots of displacement U 3 on the 308 th day on the top and bottom surfaces of Plate 8 are given. It can be easily noted that the distributions of out-of-plane displacements on the top and bottom surfaces are asymmetry, especially for the U-Mo creep rate coefficient of 2,000 × 10 −22 mm 3 /(fission·MPa). Besides, one can see from Figures 6A-C that the displacements will be locally enhanced for the creep rate coefficient of 2,000 × 10 −22 mm 3 /(fission·MPa). In Yan et al.    (Mao et al., 2021).
Path 3 on the top surface of a fuel plate is shown in Figure 3. The displacement components of U 3 along the similar paths can be found in Figure 7A, which are the results after irradiation of 308 days for different fuel plates and outside Al plates. The maximum out-of-plane displacement on the top surface of Plate 8 is 0.074 mm, which is larger than those in the other plates. For the creep rate coefficient of 500 × 10 −22 mm 3 /(fission·MPa), the maximum displacements in the fuel plates decrease firstly, and then increase with the distances away from Plate 1. It can be noted that the displacement pattern of Plate 8 differs from those of the other fuel plates. Figure 7B gives the comparison of the Plate 8 displacement results of U 3 along Path 3 on the 308 th day for the three cases. When the creep rate coefficient increases from 500 × 10 −22 mm 3 /(fission·MPa) to 2,000 × 10 −22 mm 3 /(fission·MPa), the maximum out-of-plane displacement increases by ∼124%. In order to interpret the contributions of bending deformations, Path 4 in Plate 8 is chosen (shown in Figure 3), and the bending moments of M 11 along Path 4 are given in Figure 8. The bending moments of M 11 in per unit length are obtained through the integration manipulation, as depicted in Appendix A. It can be noticed that the distribution curve of bending moment along Path 4 correlates well with the displacement distribution curve in Figure 7B. It is demonstrated that the mechanical interactions in the fuel assembly result in the bending deformations of Plate 8, which are responsible for the displacement pattern of Plate 8. In the range of 500 × 10 −22 mm 3 /(fission·MPa) to 2,000 × 10 −22 mm 3 /(fission·MPa), the maximum bending moments of M 11 give a distinct rise from 0.03 to 0.07 N.
From the results in Figures 7A, 9A, one can see that the thickness increments of Plate 1 mainly result from those of the fuel foil. The thickness increments at the foil edges are close to zero, although the irradiation swelling deformations tend to induce the increase of the fuel foil. Besides, a distinct local increase of the foil thickness can be found near the fuel foil edges. This deformation mechanism has been explained in our previous work (Jian et al., 2019a), which is attributed to the mutual interactions of irradiation swelling and creep. In Jian et al. (2019a), the simulation results of the thermal-mechanical behavior in a mini monolithic U-10Mo fuel plate were obtained, with the adopted theoretical models and solution strategies in this study. The calculated results of the fuel foil thickness matched well with the experimental data in Kim et al. (2013), as shown in Figure 9B. The local deformation phenomena in Figures 9A,B are consistent with each other.
For the U-Mo creep rate coefficient of 2,000 × 10 −22 mm 3 / (fission·MPa), it was found the smallest width existed at Channel 7 (Mao et al., 2021). To describe the variations of channel width directly, Figures 10A-D show the width evolution results of Channel 7. It is noted that the channel width of the point for a certain x-coordinate refers to the smallest value in the y-direction. One can see that the channel widths will be narrowed at locations with certain distances away from the two ends. Meanwhile, it can    May 2021 | Volume 9 | Article 676881 9 rate coefficient of 2,000 × 10 −22 mm 3 /(fission·MPa), a position with a local minimum width is observed, with x 80 mm. The point of x 80 mm locates at the plane of view-cut 1 in Figure 3. Nevertheless, the corresponding channel width on the 308 th day is ∼7% smaller than the original value, which is much smaller than the reported maximum increase of 20% in Deng et al. (2017).

On the First Principal Stresses and Irradiation Creep Strains
It was found in some post-irradiation examinations that cracks appeared in the fuel foil near the fuel foil/cladding interface (Meyer et al., 2014), which was explained to have relations with the first principal stresses of fuel foil (Jian et al., 2019a). For the three creep rate cases, the maximum values of the first principal stress all exist in the fuel foil of Plate 1. For 500 × 10 −22 mm 3 / (fission·MPa), the contour plot of the first principal stresses after irradiation of 308 days is given in Figure 11. Relatively large tensile stresses can be found near the corners and edges of fuel foil, and the maximum value becomes ∼54.6 MPa. To analyze the effects of creep rate coefficient on the first principal stresses, Path 5 in the fuel foil of Plate 1 (as shown in Figure 11) is chosen to output results, which passes through the point of the maximum tensile stress. Figures 12A-D give the evolution results of the first principal stress and creep strain component along Path 5. Except for those near the two ends, the tensile stresses can be seen and larger values appear near the locations with ∼2 mm away from the two path ends. With the increase of creep rate coefficient from 500 × 10 −22 mm 3 /(fission·MPa) to 2,000 × 10 −22 mm 3 /(fission·MPa), the maximum tensile stress on the 308 th day decreases by ∼25%, which is similar to that in Yan et al. (2019). After irradiation of 231 days, the maximum tensile stress is ∼58 MPa for 500 × 10 −22 mm 3 /(fission·MPa), which is ∼34% greater than that for 2,000 × 10 −22 mm 3 /(fission·MPa). As mentioned above, creep damages might occur to result in the reduction of the fuel foil strength, and then induce the fuel foil failure (Iltis et al., 2016;Meng and Wang, 2016). Comparing the through-thickness creep strain components in the fuel foil of eight fuel plates, one can obtain that relatively large values also take place in Plate 1. Figures 12A-D also display the evolution results of throughthickness creep strain component ε cr 33 for the considered three creep rate cases. It can be observed that the maximum creep strain component ε cr 33 appears near the two path ends, and increases with the creep rate coefficient. When the coefficient rises from 500 × 10 −22 mm 3 /(fission·MPa) to 2,000 × 10 −22 mm 3 /(fission·MPa), the peak value of ε cr 33 on the 308 th day increases by 63%. Figures 13A,B show the contour plot of the maximum skeleton normal stresses after irradiation of 308 days in the porous fuel  foils of Plate 1 and Plate 8. The maximum skeleton normal stress becomes 97.5 MPa in the fuel foil of Plate 1, which is 9.2% greater than that in Plate 8. As described in Eq. 8, the maximum skeleton normal stress depends on the first principal stress, bubble volume fraction, bubble radius and pressure. As the same irradiation condition is considered in the three cases, the results of bubble volume fraction in % and bubble pressure are almost the same. As depicted in Figure 13C, the maximum bubble volume fraction is up to 12.8% in the heavily irradiated region of fuel foil. In Jian et al. (2019a), it was found that the bubble volume fraction increased with the fission density and was speeded up after grain recrystallization. The values of bubble pressure differ slightly, as described in Figure 13D. After irradiation of 308 days, the maximum bubble pressure is up to 99.7 MPa.

On the Maximum Skeleton Stresses, Bubble Fractions and Bubble Pressures
In Figure 14, the maximum skeleton normal stresses at different irradiation time can be found, which are the results along Path 6 in the fuel foil of Plate 1. It should be noted that Path 6 has the point with the peak value. The maximum values exist at the location with a distance of ∼2 mm away from the two path ends, the same as the peak value location of the first principal stress. With the increase of irradiation time, the tensile stresses are enlarged as a whole. After irradiation of 308 days, a peak value of 97.5 MPa is obtained for the creep rate coefficient of 500 × 10 −22 mm 3 /(fission·MPa), which is ∼23.3% higher than that for 2,000 × 10 −22 mm 3 /(fission·MPa). This indicates that the mechanical interactions between the fuel foil and cladding become weaker with the increase of creep rate coefficient. While, it can't be inferred that a larger creep rate coefficient is

Effects of Fuel Foil Creep Rate Coefficient on the Von Mises Stresses of Cladding
In this section, the Von Mises stresses in the cladding will be discussed, due to the fact that fuel cladding is the first safety barrier. In the considered fuel assembly, the maximum Von Mises stress of cladding takes place in the fuel plates closer to Al outside plates (Mao et al., 2021). Figure 15 gives the contour plot of Von Mises stress on the 308 th day in the cladding of Plate 8, for the creep rate coefficient of 500 × 10 −22 mm 3 / (fission·MPa). The maximum Von Mises stress of 258 MPa appears near the corner of fuel foil. One can find from Marchbanks (1995) that the tensile strength value becomes 320 MPa at a test temperature of 344 K, which is the temperature of the maximum stress point. For 500 × 10 −22 mm 3 /(fission·MPa), it can be obtained that the cladding for all the fuel plates is relatively safe after irradiation of 308 days. To present the effects of creep rate coefficient, Path Seven in the cladding is chosen to output results, which passes through the point with the maximum Von Mises stress. Figure 16 gives the evolution results of Von Mises stress along Path 7 for three cases. The Von Mises stresses can be seen to increase with irradiation time, and decrease obviously with the increase of fuel creep rate. After irradiation of 308 days, the maximum Von Mises stress for 2000 × 10 −22 mm 3 / (fission·MPa) is ∼24% smaller than that for 500 × 10 −22 mm 3 / (fission·MPa). With the increase of creep rate coefficient, the stress relaxation effects are enhanced in the fuel foil, and the mechanical interactions between the fuel foil and the cladding are weakened.

CONCLUSION
In this study, the three-dimensional thermo-mechanical variables in a U-Mo/Al monolithic fuel assembly constrained at the two ends are obtained for the fuel foil creep rate coefficients of 500 × 10 −22 mm 3 /(fission·MPa), 1,000 × 10 −22 mm 3 /(fission·MPa) and 2000 × 10 −22 mm 3 /(fission·MPa). The effects of creep rate coefficient are investigated. With the increase of creep rate coefficient from 500 × 10 −22 to 2,000 × 10 −22 mm 3 / (fission·MPa), it is indicated that: (1) The temperature field in the fuel foil and side plates changes slightly; the maximum out-of-plane displacements in the side plates decrease slightly, with the difference on the 308 th day less than ∼2.4%. After irradiation of 308 days, the peak outof-plane displacement increases by 124% in the fuel plates adjacent to the outside Al plates, compared to that at the initial stage of irradiation, which mainly results from the enhanced bending deformation contributions. Nevertheless, the effects on the widths of coolant channels are relatively small. (2) The first principal stress and the maximum skeleton normal stress in the fuel foil are distinctly reduced, but the throughthickness creep strain component increases. (3) Larger Von Mises stresses in the Al cladding appear near the interface with the fuel foil, and the maximum values exist in the fuel plates closest to the outside Al plates. The Von Mises stresses decrease obviously with the increase of fuel creep rate coefficient, with the maximum value on the 308 th day reduced by ∼24%. From this point of view, a higher creep rate of U-Mo fuel foil needs to be achieved by advanced fabrication.
The failure of U-Mo foil depends possibly on the skeleton stress and the irradiation creep degraded skeleton strength. In our future work, the relation of U-Mo skeleton strength with the through-thickness creep strain component will be developed, and the effects of fuel creep rate on its failure behavior will be further evaluated.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
XM, XJ and JZ contributed to data curation, formal analysis and writting-original draft. SD, YH, HW, LZ and YL contributed to the conception and design of the study. All authors contributed to manuscript revision, read, and approved the submitted version.