1. Introduction
The operation of the fuel elements in a nuclear reactor core is essential to the safety of the nuclear reactor. Ensuring the safety and integrity of nuclear fuels is of the utmost importance for their development in both normal and emergency situations [
1]. Early research reactors have been utilized to investigate theories and conduct experiments on fission gas behavior [
2]. The variation in grain size has no effect on the release threshold in the LWR. The grain size can only affect the release rate when the grain boundary is networked [
3]. Fission gas will enter the gap between the fuel and cladding in the LWR when the fuel fails in the high enthalpy condition of reactivity-initiated accidents (RIAs) [
4]. The intragranular swelling rate is subjected to fission gas in the PWR fuel rod [
5]. The FGR sharply elevates with increasing temperature at the initial phase and then tends towards moderation in the PWR during irradiation. FGR occurs via diffusion of gas atoms at the opening grain boundary [
6]. The FGR is subjected to gas sweeping on the grain boundary in the PHWR [
7]. The simulation results show that fission products raise with grain diameter in the VVER [
8].
After the Fukushima accident, the international community recognized the need for advanced fuel designs that offer substantially improved performance under severe accidents, leading to the U.S. Department of Energy’s Nuclear Energy Office (DOE-NE) launching the Accident Tolerant Fuel Campaign (ATF Campaign). The safety margin of nuclear power plants (NPPs) can be enhanced by improving the oxidation-resistant performance of fuel cladding. There are primarily three types of advanced fuel cladding. The first is Cr-coated Zr-based cladding, which possesses several merits. It can improve the peak cladding temperature limit during the design of basic LOCA and extend the time of before post-quench ductility failure occurs [
9]. Meanwhile, Chromium coated on ATF is viewed as one of the most promising forms of ATF cladding. In terms of its resistance to high-temperature oxidation and deformation, the retention of fission products, and heat transfer between the coolant and cladding, the performance of ATF cladding is better than conventional Zr-based cladding [
10]. The second type of advanced fuel cladding is FeCrAl alloys that can resist oxidation under high temperature and steam conditions. And the last type of cladding is Si/SiC cladding, which is the best among the three cladding materials in terms of its oxidation resistance and high-temperature strength; therefore, it has been called the ideal ATF cladding material [
9]. Innovative fuel systems must be implemented to enhance fuel, plant work, and safety to meet or exceed design benchmark accidental conditions [
11]. While current reactor designs predict the behavior of fission gases within fuel rods over their operational lifetime, current fuel designs are not able to compensate for over-lifetime gas release, and an extended burnup extent in existing reactors is a significant issue [
12].
The behavior of gaseous fission products under irradiation conditions heavily impacts nuclear fuel rod performance [
13]. The released fission gases access the rod’s free space, which increases the internal pressure of the rod or creates bubbles, leading in fuel swelling and enhanced pellet cladding mechanical interactions (PCMIs). The above processes affect the mechanical behavior of rods [
14], while the FGR reduces the thermal conductivity between the fuel and cladding. Gases within bubbles degrade the thermal conductivity of the fuel, ultimately influencing the temperature distribution of the fuel pellets. As a result [
15], fission gases significantly affect the thermal properties of fuel rods. When fission gases enter bubbles, the fuel swells and mechanically interacts with the fuel cladding, closing the gap between the fuel and cladding. Fission gas release (FGR) causes the pressure in the rod free space to increase as well as causing the thermal conductivity of the rod filling gas to degrade.
The safety of nuclear reactors and ultimate fuel repositories are assessed by utilizing nuclear fuel performance codes, with prediction for FGR being a crucial aspect. In fuel performance codes, FGR typically includes two steps [
14]. First, based on the linear power rating, fission gas atoms are generated through fission reactions, and some of the gas atoms move to the grain boundaries via diffusion, with thermal diffusion during high temperatures being the primary form. In addition to thermally activated defects, irradiation-enhanced defect concentrations drive gas atoms to diffuse, with athermal diffusion under low temperatures being important. The released gas atoms are particularly relevant to the thermal conductivity of fuels with fast neutron factors, like nitride fuels, for which the release fraction at high temperatures may be lower than oxide fuels; however, the thermal conductance of the gap also deteriorates. Second, the concentration and release threshold of gas atoms causes the gas atoms to release and be stored on the grain boundaries. This behavior is affected by temperature and other properties.
When analyzing the thermal–mechanical performance of fuel rods, fission gas behavior is a significant aspect, and this behavior is modeled by fuel performance codes. Modeling this behavior in nuclear fuels under irradiation involves accounting for various interactive phenomena. Gas atoms are produced inside fuel grains and diffuse towards grain boundaries. Bubbles precipitating and growing at the grain faces leads to fuel swelling. After bubble growth and interlink, the gas atoms release from the grain faces, causing thermal FGR. Gas atoms and the vacancy diffusion of bubbles contribute to the development of this process. Except for these diffusion ways, experiments have also justified that micro-cracking is mechanism of grain face separation, and gas atoms may be released from grain boundaries through these cracks, resulting in significant FGR in transient situations. These fast kinetics typically take a burst release into account. Under irradiation experiments, the FGR rate raises due to cracking on the grain boundaries, which is caused by fission gas accumulating and inducing thermal stresses. After transient heating, grain surfaces exhibit the phenomenon of planar separation, and the stored gas atoms release through the cracking along a particular area between two grains. Micro-cracking occurs at grain face causing grains unable to store gas atoms and thus leadin to additional fission gas release during transient situations [
14].
According to Tonks et al. 2018, [
16], fuel performance is affected by the amount and release of fission gases in UO
2 fuel, causing the fuel to swell and the gap pressure to increase. The majority of FGR undergoes three processes: first, gases are generated and move from bulk; second, bubbles in the grain face nucleate, grow, and interlink until they adhere to grain boundaries; third, gas atoms are transported through interlinkage grain boundary passages until they arrive at a free surface and are released. Gas atom diffusion primarily determines the above mechanisms, although knockout, recoil, and burst release may cause some gas atoms to release.
As mentioned earlier, nuclear fuel performance is greatly affected by fission gas behavior. Moreover, fission gas behavior is related to the integrity and safety of nuclear power plants [
17].
Indeed, the behavior of fission gases is of utmost importance in understanding and managing nuclear reactors, particularly in the context of nuclear accidents such as Chernobyl and Fukushima. In the case of Chernobyl, the accident was caused by a combination of design flaws, operator errors, and a rapid power excursion that led to a destructive steam explosion and a subsequent graphite fire. The release of fission products, including fission gases such as xenon and krypton, played a significant role in the dispersion of radioactive material into the atmosphere. Similarly, the Fukushima accident was triggered by a massive earthquake and tsunami that resulted in the loss of power and cooling capabilities at the Fukushima Daiichi nuclear power plant. The subsequent overheating and melting of fuel assemblies led to the release of fission gases and other radioactive materials into the environment. Studying the behavior of fission gases under severe accident conditions helps in assessing the consequences and developing strategies for mitigating their release and dispersal. Research on fission gas behavior encompasses various aspects, including gas release mechanisms, diffusion and transport within fuel materials, retention in fuel rods, and the potential effects on fuel performance and reactor safety. By gaining a deeper understanding of these factors, scientists and engineers can develop improved fuel designs, safety systems, and accident management strategies to enhance the overall safety and reliability of nuclear reactors.
It is worth noting that significant advancements have been made in the field of fission gas behavior research over the years, contributing to the development of more robust and accident-tolerant fuel designs, improved reactor operation and maintenance practices, and enhanced safety measures in the nuclear industry. However, fission gas behaviour is very complicated, as it includes FGR, fission gas diffusion, grain growth, and swelling. On the other hand, different types of nuclear fuels may behave differently. It is relatively easy to simulate single fission gas behavior, such as FGR processes under steady state conditions, but this fails to reflect the thermal fission gas behaviour in the nuclear reactor. A more sophisticated fission gas behaviour model is strongly needed to better understand this phenomenon.
UO2 fuels have been successfully and broadly used as primary fuels in commercial light water reactors. The prediction of fission gas behaviour is important for fuel rod performance when a reactor operates at stable conditions. The fission gas is released when the gas reaches a saturation value in the grain boundaries. The fuel rod pressure and temperature continue to increase with gas release due to the degradation of gap thermal conductance between the pellet and cladding. In this study, we considered both FGR and diffusion for the gas behavior. We mainly focused on investigating the parameter sensitivity of FGR in LWR through simulations. The effect of grain growth was implicitly considered by using FGR with different grain sizes. COMSOL Multiphysics was used for the modeling of UO2 grain geometry and for solutions to relevant equations. The FGR over time, the fission gas diffusion coefficient, and burnup with temperature under three different conditions in a light water reactor (LWR) are discussed in this study. The FGR with unequal grain sizes, diffusion coefficients, temperatures and linear heat generation rates with different working conditions are compared.
2. Methods
The extended model is based on previous theories and empirical formulations [
18,
19,
20,
21,
22,
23,
24,
25,
26]. The thermomechanical properties and their relevant equations are given in reference [
15].
The mechanism of FGR includes two components: athermal and thermal FGR [
27]. Athermal release happens at a high burnup when gas atoms escape from the free surfaces of the fuel through recoil and knockout. When the temperature exceeds the incubation threshold, thermal release drops with increasing temperature [
28]. The exact value of the threshold can be derived from the occurrence of intergranular gas. Gas atoms diffuse from the inner side of the fuel grains to the grain boundaries. The diffusion velocity rises with temperature.
In addition, a part of intergranular atoms is brought back into the solution through fission spikes. This is known as resolution. The resolution flux will prevent the flux diffusion to some extent. Resolution occurs within the grain, where atoms enter and form small bubbles. These bubbles eventually reintegrate into the solution through recombination and resolution, restoring equilibrium.
The mathematical framework of Booth [
29] and Speight [
24] provided a good solution method to the diffusion problem. Booth offered analytical expressions for constant irradiation under the LWR operating conditions, and two cases were used to present steady and transient state conditions in LWR that regarded the constant gas production case as the steady state condition and the zero-gas production case as the transient condition. Diffusion is neglected during transient operation since the power level is high enough for the gas production. The limitations of this study include the assumption of a perfect sink boundary condition, where there is a release as soon as gas atoms reach the grain boundary. In this case, the concentration at the grain boundary is zero in the mathematical expression. Resolution was later introduced by Speight to improve the Booth model [
24]. Furthermore, the concentration at the grain boundary under imperfect sink conditions remains finite. Turnbull derived a closed form analytical solution for the model [
25].
The extended model we present here is based on real physical process. The magnitude order of the physical quantities is important for matching the numerical algorithm to best suit the physical problem. In order to ensure the accuracy of a wide range of algorithm conditions, all the situations encountered in the normal and abnormal operations of the LWR were restricted to Class-II incidents.
The FGR model adopts an idealized spherical UO2 grain. The parameters and expressions are given, such as the grain size and the diffusion coefficient of Xenon inside UO2. These parameters are essential to solve the ensuing mathematical problem efficiently. The details of the parameters are presented below.
2.1. Grain Size
The grain grows with increasing temperature, which introduces negative effects. However, to simplify the model, one can fix the grain size and ignore the effects of grain growth. This approach overlooks the impact of grain growth on the diffusion coefficient and subsequently on gas release. When grain growth occurs, the diffusion coefficient tends to decrease, leading to reduced gas release. This relationship poses a complex problem that requires solving the diffusion equation with a moving boundary. The decreasing diffusion coefficient affects the rate at which gas is released due to grain growth. Additionally, grain growth results in the clustering of gas on the grain boundaries, which enhances gas release. This gas sweeping effect occurs as the grains grow, leading to increased gas release from the material. By comparing the computational results of fixed grain size without gas sweeping and grain growth with gas sweeping, it is evident that the latter exhibits higher gas release.
2.2. Diffusion Coefficient
In Turnbull et al. [
30,
31], the diffusion coefficient consists of the intrinsic diffusion in the absence of irradiation, and the thermal and athermal diffusion produced by irradiation can be derived. Owing to the irradiation-induced resolution of intragranular bubbles, the diffusion coefficient combined with Speight’s theory is revised as follows [
24]:
Part of the parameters are accessed from reference [
30,
31,
32], where
D is diffusion coefficient;
is sum of the intrinsic, thermal, and athermal diffusion coefficients;
g is probability of fission spikes trapping atoms within intragranular bubbles [
27];
b′ is probability of the resolution of these atoms by fission spikes;
D1 is intrinsic diffusion coefficient;
D2 is thermal diffusion coefficient;
D3 is athermal diffusion coefficient;
L is mean free path between two bubbles;
D01,
D02, and
D03 are diffusion constants;
TK is the local temperature;
T01 and
T02 are the parameters of the diffusion coefficient calculation; and
P′ is linear heat generation rate (LHGR). The parameters of these equations are given as:
D01 = 3.9 × 10
−6 m
2/s,
T01 = 45,275 K,
D02 = 1.77 × 10
−15 m
2/s,
T02 = 13,800 K,
D03 = 4 × 10
−21 m
2/s,
L2b′ = 10
−15 m
2/s.
Intrinsic diffusion represents the motion of clusters of gas atoms. Compared with low temperatures, the phenomenon of intrinsic diffusion is more easily observed at high temperatures. Gas atoms diffuse towards grain boundaries and accumulate on the grain boundaries. When the number of gas atoms saturates the grain boundaries and the temperature exceeds the incubation threshold, gas atoms begin to diffuse from the grain boundary towards the outside under the temperature gradient. This is the process of thermal diffusion. Athermal diffusion dominates at low temperatures. Recoil is the primary mechanism by which gas atoms move, resulting from collision between these atoms. This phenomenon leads to a slower saturation.
The above diffusion coefficients vary with temperature and are calculated and presented in
Figure 1. This figure indicates that fission gas diffusion is affected by intrinsic diffusion, thermal diffusion, and athermal diffusion. In particular, intrinsic diffusion is one of the most influential items for the diffusion coefficient. Furthermore, the range of variation with temperature is the largest among all the diffusion coefficients. Without considering the athermal fission gas release, the value of the athermal diffusion coefficient is constant.
Figure 1 shows a summary of the previous models and experimental results. Comparing
Figure 1a,b, the diffusion coefficient in our model is within a reasonable range compared with the other models and experimental results. In addition, the diffusion coefficient of different fuels is compared, such as UO
2, MOX, and UN. Although the results of UO
2 are greater than other fuels, there is a difference between its diffusion coefficient in the calculations and measurements due to different methods being used. For example, Combette and Zacharie used
to calculate the diffusion coefficient. Different temperature and burnup ranges were selected.
2.3. The Incubation Threshold
According to the previous theoretical and experimental results [
18], some of the gas that remained inside the intergranular bubbles is brought back into solution through fission spikes. The equation for resolution flux is as follows:
where
Jrs is the resolution flux,
b is the probability of the intergranular atoms re-entering the solution, and
N is density per unit area at the grain boundaries. Here, we assume that the resolution flux is similar to the diffusion flux in Fick’s law. Then, the mean concentration at the grain boundaries is obtained:
where
cδ is the mean concentration at the grain boundaries and
δ is the resolution depth. The flux exit of the grain is as follows [
24]:
where
J is flux exit of the grain,
J0 is flux without resolution,
β is creation rate, and
t is time. For a spherical grain of a defined radius, the flux is:
where
a is grain radius and
F0 is fission gas release fraction without resolution. The fission gas release fraction (%) without resolution approximately derives from the Booth diffusion model at a low release level:
The surface density at the grain boundaries (1/m
2) is as follows:
Based on the above equation, the surface density is time dependent. The incubation time is also related to the threshold value of the surface density. Its value is composed of two asymptotic solutions at low and high temperatures when the diffusion coefficient is constant over time:
where
Thigh and
Tlow are the high and low temperatures for diffusion, respectively;
ti is incubation time; and
Ns is the threshold value of the surface density. The incubation burnup is then taken as the incubation time:
where
Bu is the burnup. The above two asymptotic expressions can be integrated into a single expression due to the following two important remarks:
Second, many parameters are uncertain, such as
b,
δ, and
Ns. Thus, the solution form of the incubation threshold must be sought:
where
B1,
B2,
B3, and
BMIN are all adjustable parameters of the incubation threshold calculation in the model and
TB is the parameter of the incubation threshold calculation. Its values are taken from previous literature [
44]. In this case, Parameters are given in this calculation:
B1 = 1 MWd/tM,
B2 = 3.3 × 10
−8 K
−1,
B3 = 603 K,
BMIN = 1500 MWd/tM, and
TB =
T02 = 13,800 K. The variation in the thermal threshold is presented in
Figure 2. The burnup decreases with temperature. High temperatures correspond to a low burnup. The burnup of nuclear fuel increases within the reactor over time and is highest at the end of its operational life. The temperature drops with burnup. This corresponds to the reality of engineering. The variations in burnup at different temperatures represents the incubation threshold or thermal threshold, which is related to gas atoms diffusion. At low temperatures, athermal diffusion is dominant. At high temperatures, intrinsic and thermal diffusion are more important than athermal diffusion for gas atoms. The results of the incubation threshold calculated using different equations in the Bernard model are extremely closed. The values of the threshold calculated using different equations in this model are slightly different, and there is a little difference between the Bernard model and this model. However, the threshold value calculated using Equation (19) is different. This difference may be the result of some uncertain fuel parameters. The threshold value calculated using Equation (13) in this model is comparable with the Bernard model.
The value of
Ns is derived by assuming that the intergranular atoms accumulate in bubbles [
26]. The calculated equations are as follows:
by combining the above relations, we obtain:
where
P is the bubble pressure,
VB is the bubble volume,
n is number of atoms in a bubble,
k is Boltzmann constant,
T is the temperature,
nB is the bubble density,
FB is the fraction coverage of the bubbles,
RB is the bubble radius,
θ is the semihedral angle,
PEXT is the pressure,
γ is the free surface energy, and
f(
θ) is the geometry factor of the grain-edge bubble to the sphere. The exact values are given in the previous literature [
44]:
γ = 1 J/m
2,
RB = 5 × 10
−7 m,
FB = 0.25,
θ = 50°, and
f(θ)/sin2
θ = 0.288. The threshold value of the surface density varies from 10
19 m
−2 (like the value of a mono-atomic layer) to 3 × 10
20 m
−2 between the limiting two cases
PEXT = 0 MPa and
PEXT~120 MPa (UO
2 fraction threshold) [
45]. For instance, the value of
P (~90 MPa) [
46] corresponds to the threshold of the value of the surface density (~10
20 m
−2) [
47]. Equation (13) is presented in
Figure 2 as the dashed curve, where
NS = 3 × 10
20 m
−2,
TK = 1273 K, and
b = 2 × 10
−6 s
−1.
Figure 2 shows that some minor differences exist between the solid and the dashed curve, but the two lines approximately match.
The preceding parameter values will change when another gas equation of state is used in the ideal gas law. However, the order of magnitude of the threshold of value of the surface density does not change. For example, when a hard sphere equation of state is used, the value varies from 2 × 10
20 m
−2 to 3 × 10
20 m
−2 [
48].
2.4. The Diffusion Equation
According to the kinetics of diffusion, fission gases diffuse from the inner side of the grain to the grain boundary. This mechanism is used in this model. The diffusion formulation is expressed as follows [
44]:
where
c is local concentration. The local concentration
c (
r,
t) is related to the space and time variable. The boundary condition is given by Equation (9).
4. Model Application
By comparing the results of the above cases between this model and Bernard’s model, we have shown that our model can predict well the above the cases. In this section, the parameters of other models will be incorporated into this model. The validated scenarios in the model primarily consist aforementioned two cases. Then, the results will be compared and analyzed between this study and other models’ processing to validate its effectiveness.
Fission gas release for uranium dioxide fuel was modeled using a new mechanistic and engineering approach [
49]. The parameters
T = 1050 °C,
P′ = 18 kW/m,
D = 5.8510 × 10
−20 m
2/s, and
t = 3000 h were selected, and the irradiation range was approximately 70 GWd/tU. This situation is close to the first case. Namely, the initial concentration of gas atoms both on the grain center and boundary are zero. The gas atoms are generated constantly and then diffuse towards the grain boundaries during LWR steady-state operation. The gas atoms begin to release when they reach the grain boundaries. The release of gas atoms lasts for 3000 h. There are three radiuses for fission gas release: 10 μm, 20 μm and 40 μm.
Figure 6 shows that the results are similar in these two models. Furthermore, this figure demonstrates that the FGR fraction decreases with increasing grain size.
Coupling with the FALCON code, fission gas release was modeled using uranium dioxide fuel [
5]. The conditions were as follows:
T = 500 ℃,
P′ = 20 kW/m,
D = 4.0312 × 10
−21 m
2/s, and
t = 40,000 h, the irradiation range was approximately 75 GWd/tU. The first case is chosen here. The initial concentration of gas atoms both on the grain center and the boundary are zero. Gas atoms are generated constantly and then diffuse towards the grain boundaries during LWR steady-state operation. The gas atoms begin to release when they reach the grain boundaries. The release of gas atoms lasts 40,000 h. There are two radiuses for fission gas release: 9 μm and 51 μm. As shown in
Figure 7, the results of this model and Khvostov’s model are not fitted as well. This may be related to the fact that the FGR fraction is sensitive to the grain size.
Another case of the FGR was reported in oxide fuels for sodium fast reactors [
50]. Its parameters were
T = 1250 K,
P′ = 42 kW/m,
D = 5.0283 × 10
−20 m
2/s, and
t = 40,000 h, and the burnup was approximately 12 GWd/tU. The radius was 10 μm for fission gas release. This situation is similar to the second case. The initial concentration of gas atoms is zero on the grain boundary and remains finite on the grain center. Gas atoms are not generated during LWR transient operation. The gas atoms begin to release when they the reach grain boundaries. The release of gas atoms lasts for 40,000 h.
Figure 8 indicates that the FGR fraction rises with increasing temperature. The results are slightly different between this model and Karahan’s model, and the difference of results is within 15%.
Fission product gas release with microstructure dependence has been modeled for UO
2 fuel [
51]. The main parameters were as follows:
T = 675 K,
P′ = 50 kW/m for
D = 6.3283 × 10
−21 m
2/s,
P′ = 60 kW/m for
D = 6.9323 × 10
−21 m
2/s,
P′ = 70 kW/m for
D = 7.4877 × 10
−21 m
2/s, and
t = 40,000 h, and the burnup was approximately 30 MWd/kgU. There were three radiuses for fission gas release: 10 μm, 20 μm and 50 μm,. This situation corresponds to the second case. The initial concentration of gas atoms is zero on the grain boundary and remains finite on the grain center. Gas atoms are not generated during LWR transient operation. The gas atoms begin to release when they reach the grain boundaries. The release of gas atoms lasts for 40,000 h.
Figure 9 reflects that the FGR fraction drops with increasing linear heat generation rates, and the relationship is same as the above models between the FGR fraction and grain size. A small difference exists between this model and Bernard’s model, but the greatest error was restricted within 5%, and thus, the results are convincing. Most of the error was within 5%, indicating that our model can provide reasonable results.
In order to investigate the sensitivity of fission gas release for diffusion coefficient, all parameters and conditions in Notly’s model were held constant except for the variation in the diffusion coefficient.
Figure 10 demonstrates that the FGR fraction increases with a rising diffusion coefficient. The results of
a = 10 um,
P′ = 70 kW/m under one diffusion coefficient and 30 MWd/kgU,
a = 20 um,
P′ = 60 kW/m under half of diffusion coefficient and 30 MWd/kgU, and
a = 20 um under all diffusion coefficients and liner heat rates in this model are close to Notly’s model. While some differences exist between this model and Notly’s model, the overall trend is consistent. Therefore, the outcomes are effective.
The isothermal fission gas release was modeled using the following parameters [
3]:
T = 1300℃,
P′ = 38 kW/m,
D = 1.6173 × 10
−18 m
2/s, and
t = 5 h, and the irradiation range was approximately 10 GWd/tUO
2. There were three radiuses for fission gas release, 5 μm, 7.5 μm, and 11.5 μm. The second case is applied in this situation. The initial concentration of gas atoms is zero on the grain boundary and remains finite in the grain center. Gas atoms are not generated during LWR transient operation. The gas atoms begin to release when they reach the grain boundaries. The release of gas atoms lasts for 3000 h. As shown in
Figure 11, aside from the 5 μm case, the results in both of these two models are similar. The effect of grain size on the FGR fraction also corresponds to above models’ results.
Fission gas release in oxide fuel was analyzed using BISON and TRANSURANUS with the following parameters [
14]:
T = 1673 K,
P′ = 41.3 kW/m,
D = 7.5595 × 10
−18 m
2/s,
t = 10 h, and the irradiation range was approximately 46 GWd/tU. The radius was 4.68 μm. This condition suits in the second case. The initial concentration of gas atoms is zero on the grain boundary and remains finite on the grain center. Gas atoms are not generated during LWR transient operation. The gas atoms begin to release when they reach the grain boundaries. The release of gas atoms lasts for 80 h. As shown in
Figure 12, the time in Barani’s experiments is 80 h, but in this model, the time 10 h. The final values of the FGR fraction are close in these two models.
5. Conclusions
We presented an extended thermal activated fission gas release model for UO2 fuel. The relevant physical phenomenon addressed in this model involves the diffusion of fission gases from the inner region of a UO2 grain to the grain boundary. It focuses on the resolution between the grain and the grain boundary, as well as the accumulation and saturation of gas at the grain boundary. This study considers three cases under steady-state conditions: a perfect sink with gas production, a perfect sink without gas production, and an imperfect sink. A perfect sink with gas production is viewed as steady-state, and a perfect sink without gas production is regarded are a transient operation of a LWR. In particular, transient operation belongs to typical Class-II incidents in LWRs. In addition, the effect of FGR resolution on the grain boundary was also considered.
This model investigates how the fraction of fission gas release (FGR) varies over time, considering different factors such as the grain size, diffusion coefficient, temperature, and linear heat generation rate. Through the validation of the FGR model, several conclusions can be drawn. Gas atoms are release during transient operation faster than during steady-state operation, and the time of FGR is shorter than later in the model validation process when FGR fraction is equivalent. Fission gas will be released immediately during nuclear accidents. According to the results of the imperfect sink, it is evident that FGR is greatly influenced by resolution. The resolution will delay the onset time of FGR.
Furthermore, it was observed that the FGR fraction increases with higher diffusion coefficients and temperatures, while it decreases with larger grain sizes and higher linear heat generation rates. The diffusion coefficient, which is influenced by the temperature and linear heat generation rate, plays a significant role in driving the variation of the FGR fraction. Therefore, the impact of the diffusion coefficient on the FGR fraction is essentially a result of the effects of temperature and the linear heat generation rate. Importantly, all the results obtained in this model align with existing theories and conclusions. Hence, the effectiveness of this model is deemed acceptable based on its consistency with established knowledge and findings.