Spent nuclear fuels should be safely stored until being disposed and dry storage system is predominantly used to retain the fuels. Thermal analysis to estimate temperatures of spent nuclear fuel and the storage system should be performed to evaluate whether the temperatures exceed safety limit. Recently, thermal hydraulic analysis with CFD codes is widely used to investigate the temperature of spent nuclear fuel in dry storage. COBRA-SFS is a legacy code based on subchannel analysis code, and its fidelity is verified for evaluating the thermal analysis for licensing a dry cask system. Herein, thermal analysis result based on CFD and COBRA-SFS codes is compared and the Dry Cask Simulator (DCS) is assessed as a benchmark experiment in this study. Extended Storage Collaborating Program (ESCP) led by the Electric Power Research Institute (EPRI) is organized to address the degradation effects of spent nuclear fuel during long-term dry storage, and DCS is the first phase of the program. The dry storage system, containing a single BWR assembly in a canister, was designed to produce validation-quality data for thermal analysis model. ANSYS FLUENT was used to simulate DCS. Simulations were conducted in various decay heat and helium pressure inside the canister. In realistic conditions of decay heat and helium pressure of actual dry cask system, CFD and COBRA-SFS analysis results gave good agreement with experimental measurement. Peak temperatures of channel can, basket, canister and shell predicted by CFD simulation also showed good prediction and the discrepancies were less than 7 K while measurements uncertainty was 7 K. In high decay heat and high pressure condition, however, CFD and COBRA-SFS underestimated peak cladding temperature than experimental results.

The thermal integrity of spent nuclear fuels has to be maintained during their long-term dry storage. The detailed temperature distributions of spent fuel assemblies are essential for evaluating the integrity of their dry storage systems. In this study, a subchannel analysis model was developed for a canister of a single fuel assembly using the COBRA-SFS code. The thermal parameters affecting the peak cladding temperature (PCT) of the spent fuel assembly were identified, and sensitivity analyses were performed based on these parameters. The subchannel analysis results indicated the presence of a recirculation flow, based on natural convection, between the fuel assembly and downcomer region. The sensitivity analysis of the thermal parameters indicated that the PCT was affected by the emissivity of the fuel cladding and basket, convective heat transfer coefficient, and thermal conductivity of the fluid. However, the effects of the wall friction factor of the canister, form loss coefficient of the grid spacers, and thermal conductivities of the solid materials, on the PCT were predominantly ignored.

Spent nuclear fuels should be safely stored until being disposed and dry storage system is predominantly used to retain the fuels. During long-term storage, there are several mechanisms that could result in the degradation of spent nuclear fuels, and the temperature is the most important parameter to predict and estimate the degradation behaviors. Therefore, thermal analysis to estimate temperatures of spent nuclear fuel and the storage system should be performed to evaluate whether the temperatures exceed safety limit. Recently, thermal hydraulic analysis with CFD codes is widely used to investigate the temperature of spent nuclear fuel in dry storage. Herein, Explicit CFD analysis model is introduced and validated by estimating the thermal hydraulic response of the dry storage system that is Dry Cask Simulator (DCS). Extended Storage Collaborating Program (ESCP) led by the Electric Power Research Institute (EPRI) is organized to assess degradation effects of spent nuclear fuel during long-term dry storage, and DCS is the first phase of the program. The dry storage system, containing a single BWR assembly in a canister, was designed to produce validation-quality data for thermal analysis model. ANSYS FLUENT is used to simulate DCS, and the test condition of 0.5 kW decay heat and 100 kPa helium pressure was investigated in this study. In case of peak cladding temperature (PCT), PCT from the experiment was 376 K while that of CFD was 374 K. It implies CFD simulation gives good agreement with experimental measurement. Peak temperatures of channel can, basket, canister and shell predicted by CFD simulation also show good prediction and the discrepancies were less than 7 K while measurements uncertainty was 7 K.

Spent nuclear fuels released from the reactor are stored in cooling pools and then stored in dry storage casks. During the transition from the wet storage to dry storage cask, a vacuum drying process is used to remove residual water in the cask. During the vacuum drying process, gas pressure is reduced to below 400 Pa to promote evaporation and water removal. KAERI is developing a PWR single assembly (PLUS7) test equipment to simulate the thermal flow in spent fuel assembly. In this study, the thermal conductivity of air at low pressure was derived to perform the thermal analysis of the canister in vacuum. In addition, thermal analyses were performed for the canister with backfill gases of helium, air, and a vacuum in the vertical orientation using the COBRA-SFS code. At low pressure, the thermal conductivity of air depends on pressure and temperature. The reduced thermal conductivity, kr (W/m-K) was calculated using the curve fit for air at reduced pressure in thin gaps presented in the General Electric Fluid Flow Handbook. 􀝇􀯥/􀝇􀬴 = 􀬵 􀬵􀬾 􀮼􀯍/􀯉􀰋 Where, k0 is the thermal conductivity at atmospheric pressure (W/m-K), P is the reduced (vacuum) pressure (Pa), δ is the gap size (m), T is the temperature (K), and C is the Lasance constant (7.657E-5 N/m-K). The thermal conductivity of air decreases as the pressure decreases. The reduced thermal conductivity of air at pressures of 400 Pa and 40 Pa was calculated to be 0.97 and 0.77, respectively. For the analysis in vacuum, no enhancement of the convective heat transfer was assumed (Nu=1.0). For the helium backfill, the peak cladding temperature was the lowest and the axial temperature profile was the flattest due to the higher thermal conductivity and lower density of the helium. For the vacuum backfill, the peak cladding temperature was the highest and temperature gradient was the sharpest due to the only radiative heat transfer effect in the fuel assembly.

Detailed temperature distributions of the spent fuel are required to evaluate the long-term integrity of the dry storage system. In this study, a subchannel analysis method was established to obtain the detailed temperatures of a spent fuel using the COBRA-SFS code. The SAHTT (Single Assembly Heat Transfer Test) model was selected as the subchannel analysis. It was developed at the PNL to investigate heat transfer characteristics of spent PWR fuel under dry storage conditions. The SAHTT has a 15×15 rod array with simulated rods 0.42 in. (10.7 mm) in diameter. Control rod thimbles were modeled with unheated rods. The COBRA-SFS input consists a detailed subchannel model with 256 subchannels, 225 rods, and 8 slab nodes. The heat generation rate was axially uniform with total power of 1.0 kW. Subchannel analyses were performed for the vertical orientation under three different backfills of air, helium, and vacuum. For the vacuum backfill, the peak temperature was the highest and temperature gradients the sharpest only due to the radiation heat transfer effect. For the helium backfill, peak temperature was lowest and the axial profiles flattest due to the higher conductivity and lower density of helium. Subchannel analyses were also performed to evaluate the effect of thermal parameters such as surface emissivity, convective heat transfer coefficients, and flow resistance coefficients on the PCT (Peak Cladding Temperature). The PCT was affected by the emissivity of the fuel rod and the basket, and in particular, the basket emissivity had a greater effect. The PCT was affected by the Nusselt number, but the range of the Nusselt number is around 3.66. Therefore, the effect of the Nusselt number on the PCT will not be significant. As a result of the analysis according to the flow resistance coefficients, the PCT was affected by the wall friction factor, but the loss coefficients from the space grid had little effect. Subchannel technique obtained from this work can be used to predict the detailed temperature distributions of spent fuel assembly.

Dry storage is a predominantly used method as a spent nuclear fuel storage system after spent nuclear fuel is cooled in the spent fuel pool. Spent nuclear fuel is highly radioactive and it generates heat called decay heat originated by fission products and radiation. Therefore, temperature of spent nuclear fuel should be predicted whether its cladding temperature is maintained under 400°C, which is the allowable temperature limit of cladding in a dry storage. ANSYS Fluent and COBRA-SFS are predominantly used computational method to investigate the temperature of spent nuclear fuels in a dry storage. Herein, thermal analysis results with the methods were compared based on a Single Assembly Heat Transfer Test, which is a heat test with an electrically heated model of a single PWR fuel assembly in a dry cask performed at the Pacific Northwest Laboratory. Decay heat was 1kW and backfill gas was air. Fix temperature boundary condition is applied to inner shell according to measured temperature. In case of peak cladding temperature (PCT), Fluent predicted 240–284°C, while COBRA-SFS gave 243–292°C. The discrepancy between the codes is under 2.5%. The location where PCT took place was 3.65 m from the bottom of the assembly in both results. However, temperature difference between the upper and lower region of the assembly based on the Fluent was smaller than the temperature difference based on the COBRA-SFS. It means the heat was well transferred in an axial direction with Fluent compared to COBRA-SFS. In lower plenum region where air was naturally circulated, COBRASFS had disadvantages compared to Fluent because it modeled the lower plenum by single node, so it was hard to simulate convection heat transfer by natural circulation. natural circulation speed of air in a center region of the assembly was 0.07–0.1 m·s−1 in both cases.

The temperature of the spent fuel cladding is the basis for the evaluation of integrity. It is almost impossible to directly measure the temperature of spent nuclear fuel. Because spent nuclear fuel is dangerous. We are preparing a test to measure the cladding temperature with an equivalent fuel assembly by simulating the characteristics of spent nuclear fuel. PLUS7 was selected as the test target in consideration of the amount of generation, thermal water retention, residual moisture content, and manufacturability of domestic spent nuclear fuel. The nuclear fuel assembly is planned to be manufactured in two main ways. Except for the cladding tube that simulates decay heat, the structure will be manufactured by KEPCO Nuclear Fuel, and fuel rods and canisters will be manufactured by SUKEGAWA Electric in Japan. The same nuclear fuel assemblies as commercial skeleton will be applied. The temperature of the fuel cladding will be measured by attaching a thermocouple directly to the surface of the cladding tube. The canister is composed of a basket, a basket supporter, a heater and drain tube, a lead, and an observation window. The working fluid is water and helium, and the maximum pressure inside the canister is 1.1 MPa and the minimum pressure is 0.05 kPa. The maximum temperature of the surface of the cladding was designed to be 500°C, and the maximum temperature of the sealing to keep airtightness was designed to be 250°C. To satisfy this condition, we plan to evaluate the leak rate below 10−5 std.cm3·s−1, which is equivalent to helium tightness. The maximum heat of decay per fuel rod is 13 W, and one assembly is up to 3 kW. Production of the test equipment is expected to be completed in the first half of next year, and testing is scheduled to begin in the second half of next year. The test will evaluate all environments that the spent nuclear fuel may experience, such as dry normal conditions, abnormal conditions, wet conditions, and dry conditions. All data will be used for interpretation verification purposes.