The objectives of this paper are: (1) to conduct the thermal analyses of the disposal cell using COMSOL Multiphysics; (2) to determine whether the design of the disposal cell satisfies the thermal design requirement; and (3) to evaluate the effect of design modifications on the temperature of the disposal cell. Specifically, the analysis incorporated a heterogeneous model of 236 fuel rod heat sources of spent nuclear fuel (SNF) to improve the reality of the modeling. In the reference case, the design, featuring 8 m between deposition holes and 30 m between deposition tunnels for 40 years of the SNF cooling time, did not meet the design requirement. For the first modified case, the designs with 9 m and 10 m between the deposition holes for the cooling time of 40 years and five spacings for 50 and 60 years were found to meet the requirement. For the second modified case, the designs with 35 m and 40 m between the deposition tunnels for 40 years, 25 m to 40 m for 50 years and five spacings for 60 years also met the requirement. This study contributes to the advancement of the thermal analysis technique of a disposal cell.

In order to ensure the long-term safety of a deep geological repository, the performance assessment of the Engineered Barrier System (EBS) considering a thermal process should be performed. The maximum temperature at the side wall of a disposal canister for the technical design requirement should not exceed 100°C. In this study, the thermal modelling was conducted to analyze the effects of the thermal process from a disposal canister to the surrounding near-field host rock using the PFLOTRAN code. The mesh was generated using the LaGriT code and the material properties were assigned by applying the FracMan code. Initial conditions were set as the average geothermal gradient (25.7°C/km) and an average surface temperature (14.7°C) in Korea. The highest temperature was observed at the middle of the canister side wall. The temperature of the buffer was lower than that of the canister, and the temperature increase of the deposition tunnel and the host rock was insignificant due to the lower effect of the heat source. The result of the thermal evolution of the EBS represented the highest thermal effects in the vicinity of the canister. In addition, the thermal effects were largely decreased after 10 years of the entire simulation period. It demonstrated that the model took 3 years to heat up the buffer around the canister. The temperature at the canister side wall increased until 3 years and then decreased after that time. This is because that the radioactive decay heat from the heat source was emitted enough to raise the overall temperature of the EBS by 3 years. However, the decay heat rate of the canister decreased exponentially with the disposal time and then its decay heat was not emitted enough after 3 years. In conclusion, the peak temperature results of the EBS were lower than 70°C to meet the technical design requirement.

The thermal evaluations for the conceptual design of the deep geological repository considering the improved modeling of the spent fuel decay heat were conducted using COMSOL Multiphysics computational program. The maximum temperature at the surface of a disposal canister for the technical design requirement should not exceed 100°C. However, the peak temperature at the canister surface should not exceed 95°C considering the safety margin of 5°C due to several uncertainties. All thermal evaluations were based on the time-dependent simulation from the emplacement time of the canister to 100,000 years later. In particular, the heat source condition was set to the decay heat rate and axial decay heat profile of the PLUS7 fuel with 4.0wt% U-235 and 45 GWD/MTU. The thermal properties of the granitic rock in South Korea were applied to the host rock region. For the reference design case, the cooling time of the SNF was set to 40 years, the distance between the deposition holes 8 meters and that between the deposition tunnels 30 meters. However, the peak temperature at the canister surface at 10 years was 95.979°C greater than 95°C. This design did not meet the thermal safety requirement and needed to be modified. For the first modified case, when the distance between the deposition tunnels was set to 30 meters, three cooling time cases of 40, 50 and 60 years and five distances of 6, 7, 8, 9 and 10 meters between the deposition holes were considered. The design with the distances of 9 and 10 meters between the deposition holes for the cooling time of 40 years and all five distances for 50 and 60 years were less than 95°C. For the second modified case, when the distance between the deposition holes was set to 8 meters, three cooling time cases of 40, 50 and 60 years and five distances of 20, 25, 30, 35 and 40 meters between the deposition tunnels were considered. The design with the distances of 35 and 40 meters between the deposition tunnels for the cooling time of 40 years, the distances of 25, 30, 35 and 40 meters for 50 years and all five distances for 60 years were less than 95°C. As a result, the peak temperature at the canister surface decreased as the cooling time and the distance between the deposition holes and the tunnels increased.

The nuclear criticality analyses considering burnup credit were performed for a spent nuclear fuel (SNF) disposal cell consisting of bentonite buffer and two different types of SNF disposal canister: the KBS-3 canister and small standardized transportation, aging and disposal (STAD) canister. Firstly, the KBS-3 & STAD canister containing four SNFs of the initial enrichment of 4.0wt% 235U and discharge burnup of 45,000 MWD/MTU were modelled. The keff values for the cooling times of 40, 50, and 60 years of SNFs were calculated to be 0.79108, 0.78803, and 0.78484 & 0.76149, 0.75683, and 0.75444, respectively. Secondly, the KBS-3 & STAD canister with four SNFs of 4.5wt% and 55,000 MWD/MTU were modelled. The keff values for the cooling times of 40, 50, and 60 years were 0.78067, 0.77581, and 0.77335 & 0.75024, 0.74647, and 0.74420, respectively. Therefore, all cases met the performance criterion with respect to the keff value, 0.95. The STAD canister had the lower keff values than KBS-3. The neutron absorber plates in the STAD canister significantly affected the reduction in keff values although the distance among the SNFs in the STAD canister was considerably shorter than that in the KBS-3 canister.

In order to dispose of spent nuclear fuel (SNF) in deep geological repository, source term evaluation considering its specification, enrichment, burnup, cooling time should be performed. In this study, the measured values of Takahama-3 pressurized water reactor SNF (WH 17×17) samples were analyzed with SCALE 6.1/ORIGEN-S and TRITON code calculation results for validation. Unlike the ORIGENS code, TRITON code calculations differed from two-dimensional neutron flux distribution by using the multi-group cross-section library. Both calculation results from ORIGEN-S and TRITON code showed higher errors in 234U, 239Pu, and 241Pu compared to other actinide nuclides. In the case of axial locations of fuel rods in fuel assembly, fuel rods located at the edge of the fuel assembly presented increased errors due to nuclear reaction cross-section. Overall, the ORIGEN-S predictions informed more accurate agreement with the measured results compared with TRITON results. Especially to 235U, 239Pu, and 240Pu radionuclides, ORIGEN-S errors were denoted more than twice as low as the TRITON results. Comparing the calculation results with experimental results implied that the ORIGENS code was more accurate code than the TRITON code for source term evaluation.

The criticality analyses considering burnup credit were performed for a spent nuclear fuel (SNF) disposal cell consisting of bentonite buffer and two different types of PWR SNF disposal canister: the KBS-3 type canister and the small standardized transportation, aging and disposal (STAD) canister. The criticality analyses were carried out for four cases as follows: (1) the calculation of isotopic compositions within a SNF using a depletion assessment code and (2) the calculation of the effective multiplication factor (keff) value using a criticality assessment code. Firstly, the KBS-3 type canister containing four SNFs of the initial enrichment of 4.0wt% 235U and discharge burnup of 45,000 MWD/MTU was modelled. The keff values for the cooling times of 40, 50, and 60 years of SNFs were calculated to be 0.74407, 0.74102, and 0.73783, respectively. Secondly, the STAD canister was modelled. The SNFs contained in the STAD canister were assumed to be the enrichment of 4.0wt% and the burnup of 45,000 MWD/MTU. The keff values for the cooling times of 40, 50, and 60 years were estimated to be 0.71448, 0.70982, and 0.70743, respectively. Thirdly, the KBS-3 canister with four SNFs of which the enrichment was 4.5wt% and the burnup was 55,000 MWD/MTU was modelled. The keff values for the cooling times of 40, 50, and 60 years were 0.73366, 0.72880, and 0.72634, respectively. Finally, the calculations were carried out for the STAD canister containing four SNFs of the enrichment of 4.5wt% and the burnup of 55,000 MWD/MTU. The keff values for the cooling times of 40, 50, and 60 years were 0.70323, 0.69946, and 0.69719, respectively. Therefore, all of four cases met the performance target with respect to the keff values, 0.95. The STAD canister showed lower keff values than the KBS-3 canister. This appears to be the neutron absorber plate installed in the STAD canister although the distance among the four SNFs in the STAD canister was shorter than the KBS-3 canister.

Currently, the interim storage pools of spent fuels in South Korea are expected to become saturated from 2024. It is required to prepare an operation plan of a domestic dry storage facility during a long-term period, with the researches on safety evaluation methods. This study modified the FRAPCON code to predict the spent fuel integrity evaluation such as the axial cladding temperature, the hoop stress and hydrogen distribution in dry storage. The cladding temperature in dry storage was calculated using the COBRA-SFS code with the burnup information which was calculated using the FRAPCON code. The hoop stress was calculated using the ideal gas equation with spent fuel information such as rod internal pressure. Numerical analysis method was used to calculate the degree of hydrogen diffusion according to the hydrogen concentration and temperature distribution during a dry storage period. Before 50 years of dry storage, the cladding temperature and hoop stress decreased rapidly. However, after 50 years, they decreased gradually and the cladding temperature was below 400 K. The initial temperature distribution and hydrogen concentration showed a parabolic line, but hydrogen was transferred by the hydrogen concentration and temperature gradient over time.