# Simulation of the Thermal Conditions of Cask with Fuel Assemblies of BN-350 Reactor for Dry Storage

Authors: Ye. T. Koyanbayev; M. K. Skakov; D. A. Ganovichev; Ye. A. Martynenko; A. A. Sitnikov
Journal: Science and Technology of Nuclear Installations
 (2019)

---

## Abstract

The analysis of the thermal condition of spent FA (fuel assembly) of BN-350 reactor in a six-place cask for dry storage is presented. Simulation of the thermal condition of the cask is conducted with finite elements method using ANSYS software. Calculations of fuel temperature, fuel cladding, and assembly structural elements are the part of the safety analysis for storage of spent FA. In conclusion, the results of the thermal calculations in the cases of filling cask with argon and atmospheric air are given when the thickness of the insulation cask with concrete is 0.5 and 1 m. As a result of the calculated studies, the parameters of SNF (spent nuclear fuel) storage are determined, under which the fuel temperatures will have minimum and maximum values.


---

## Body

## 1. Introduction

SNF is most often transferred to a dry storage method after reducing the level of radioactivity and the residual energy release. The general principle of SNF dry storage is that spent fuel is stored in sealed metal baskets, and the baskets are located in the casing of the protective cask [1–3]. Fuel elements or assemblies inside the casks are usually in an inert gas medium, which protects the materials from corrosion and helps to remove residual heat.The main characteristics for the SNF integrity in dry storage are values of fuel temperature and the cladding. There is a need for thermal calculations related to storage conditions to assess the safety of long-term storage of spent assemblies of BN-350 reactor in the cask. The paper shows the results of calculations of the stationary temperature field of a six-place cask with spent FA of the BN-350 reactor. Cases of the cask filling with argon and atmospheric air, when the thickness of the insulation cask concrete is 0.5 and 1 m, are considered; the parameters of the storage in which the fuel temperature will have maximum and minimum values are determined.


## 2. Initial Data for Thermal Calculation

Spent fuel assemblies of the BN-350 reactor are packed in six-place casks. The layout of fuel assemblies in a dry storage cask is shown in Figure1 [4]. Figure 2 shows the layout of the fuel elements in FA, where the temperature distribution calculation area is marked by a circle in the center of FA.Figure 1
The layout of the fuel assembly in six-place cask for dry storage.Figure 2
The layout of fuel elements in FA.The temperature distribution is determined in the area of the central fuel element of one of the fuel assemblies of BN-350 reactor in a six-place cask for dry storage. Materials and geometric dimensions of the cask elements and fuel assemblies are given as follows:Cask: material, austenitic stainless steel 12Х18Н10Т; outer diameter, 406.4 mm; inner diameter, 396.4 mm; height, 3854 mm.Inner basket (six-place): material, austenitic stainless steel 12Х18Н10Т; outer diameter of the cell, 130 mm; inner diameter of the cell, 129 mm.Hexagonal tube: material, austenitic stainless steel 12Х18Н10Т; width across flats, 96 mm; the thickness of hexagonal tube, 2.175 mm; height, 1050 mm.FEs: number of FEs in FA, 127 pcs; material of FE cladding, austenitic stainless steel 0Х16Н15М3B; outer diameter, 6.9 mm; inner diameter, 6.1 mm; height, 1060 mm; the gap among FEs, 1.05 mm.Fuel: material, mixture of plutonium dioxide and uranium dioxide; diameter of fuel pellet, 6 mm with the central hole of ⌀1.7 mm.Calculation is conducted for FA power of 120 and 190 W. Calculations of the profile of temperature distribution are conducted for cases of cask filling with argon and atmospheric air. The cask is surrounded with the concrete of 0.5 or 1 m.


## 3. Calculation

Calculation is conducted with ANSYS software, which implements the method of finite elements [5, 6].For calculations two-dimensional model was used which is provided in Figures1 and 2. The finite element grid of the model consists of 572143 units and 464000 elements; higher grid resolution is concentrated in the FA area.The surface temperature of the concrete is estimated as equal to 293 K. Calculation of stationary distribution of the temperature in the average height of the cross-section of the six-place cask (as most of the heat-stressed section) is conducted for two values of the residual energy in the fuel: ~54.3 kW/m3, which corresponds to the power of FA equal to 190 W, and ~34.3 kW/m3, which corresponds to the power of FA equal to 120 W.The energy release values were determined by the following formula:(1)QvNn×V,whereinN is power of FA (120 and 190 W); n is number of FEs in FA (n = 127 pcs); V is volume of FE (V = 2.756 × 10−5 m3).It is assumed that the energy release in the fuel is constant in height and radius of FA, and the energy release in the construction of FA materials and the cask is equal to zero.The heat flow from the end surfaces of the cask is assumed to be zero.Heat transfer from fuel elements to the rest part of the model is conducted using thermal conductivity and natural convective heat transfer [7]. (2)λэк=εк·λ,wherein εк is coefficient of convection; λ is coefficient of medium heat conductivity (argon, air), W/(m·K).Coefficient of convection is determined by the following formula:(3)εк=α·Gr·Prb,wherein Gr is Grashof number; Pr is Prandtl number.When106<Gr·Pr<1010α=0.4, b=0.2; when 103<Gr·Pr<106α=0.105, b=0.3.Grashof number is equal to(4)Gr=β·g·d3·tmax-tminν2,wherein β is temperature coefficient of volume expansion of the filler, 1/K; g is gravitational acceleration, m/sec2; d is determining size, m (hydraulic diameter, for cask, is 0.191 m, for FA it is equal to 0.0033 m); tmax is maximum temperature of the medium, K; tmin is minimal temperature of the medium, K; ν is coefficient of kinematic viscosity of the medium, m2/sec.Prandtl number is as follows:(5)Pr=μ·Cpλ,wherein μ is dynamic viscosity coefficient of the medium, Pа×sec; Cp is specific heat capacity of the medium under the constant pressure, J/(kg×K).The average calculated value of the convection coefficient when filling the cask cavity with argon under the FA power of 120 W isεк=11.3 and εк=1 for argon, filling the FA cavity. Thus, heat transfer in the FA cavity is conducted only due to thermal conductivity. This is because the value of the Rayleigh number for argon in the FA cavity is less than 103 (is 6.8); in this case, the convection coefficient is assumed to be equal to one.An iterative solution is required to account the effect of natural convection using the above method. The model is calculated with reference values of the thermal conductivity of the medium in zero-order approximation. Then convection coefficients are calculated for the medium that fills the cask cavity. The model is recalculated. The process is repeated until the difference between the calculated values of the temperature of two successive iterations model will be below 0.1 K.Thermophysical properties of materials depending of temperature used for calculations are presented in Tables1, 2, and 3 [8, 9].Table 1
Thermophysical properties of construction material.
Parameter Temperature, K 300 400 500 600 800 Stainless steel 12Х18Н10Т Thermal conductivity, W/(m×°С) 14,5 16,5 17,5 18,5 21,5 Density, kg/m3 7900 7895 7889 7871 7860 Specific heat capacity, J/(kg×°С) 505 520 535 550 600 Fuel Thermal conductivity, W/(m×°С) 11,3 10,1 9,2 8,1 6,1 Density, kg/m3 10697 10687 10674 10661 10627 Specific heat capacity, J/(kg×°С) 245 264 281 292 318 Concrete Thermal conductivity, W/(m×°С) 1,28 Density, kg/m3 2300 Specific heat capacity, J/(kg×°С) 1130Table 2
Thermophysical properties of air.
Parameter Temperature, K 300 400 500 600 800 Thermal conductivity, W/(m×°С) 0,026 0,0336 0,04 0,04354 0,58  Density, kg/m3 1,16 0,871 0,696 0,580 0,435  Specific heat capacity, J/(kg×°С) 1000 1000 1000 1000 1000Table 3
Thermophysical properties of argon.
Parameter Temperature, K 273 373 473 573 773 Thermal conductivity, W/(m×°С) 0,0165 0,0212 0,0256 0,0299 0,0379  Density, kg/m3 1,784 1,305 1,030 0,850 0,627  Specific heat capacity, J/(kg×°С) 519 519 519 519 519According to experimental data, calculation error in ANSYS is from 5 to 10 %.


## 4. Analysis of Obtained Results

Figures3–5 show the calculated temperature field in the average height of the cross-section of the cask filled with argon at FA power of 120 W, while the thickness of concrete is 0.5 m. Figures show that the maximum temperature of the fuel will not exceed 553 K.Figure 3
Temperature field of the cask and the concrete around.Figure 4
Temperature field of the cask.Figure 5
Temperature field of FA.Table4 shows the calculated temperature values determined in the positions indicated in Figures 3–5 for the eight calculation options. Argon and air were used as the filler of the cask and fuel assemblies. The thickness of the cask insulation with concrete is 0.5 and 1 m.Table 4
Value of temperature in the calculation point of the model.
Filler The thickness of the concrete, m Power of FA, W Maximal temperature of FA, K Temperature values in points, K 1 2 3 4 5 6 7 8 9 10 argon 0.5 120 553 293 396 501 450 466 493 503 493 466 552 0.5 190 682 293 460 613 538 562 601 615 601 562 679 1.0 120 594 293 440 546 496 512 538 547 538 512 593 1.0 190 745 293 531 682 610 633 670 684 670 633 743  air 0.5 120 516 293 395 481 434 449 473 482 473 449 513 0.5 190 627 293 455 582 513 534 571 584 571 534 624 1.0 120 557 293 440 525 480 493 518 527 518 493 554 1.0 190 692 293 527 651 585 606 640 653 640 606 689The data obtained as a result of the calculations show that, with an increase of concrete thickness from 0.5 m to 1 m, the fuel temperature increases by 42 K at FA power of 120 W and with filling the cask with argon and by 41 K when cask filling with air.With increasing of the concrete thickness from 0.5 m to 1 m, fuel temperature increases by 66 K at FA power of 190 W and filling of the cask with argon and by 64 K when filling the cavity of the cask with air.When cask filling with air instead of argon, the maximum temperature of the model will be reduced by ~54 K at FA power of 190 W and by ~37 K at FA power of 120 W.It is also worth noting that the maximum value of the temperature is shifted from the center of FA for certain distance. When filling the cask with argon, the maximum temperature will be shifted by ~7 mm, if the gas medium is air, and then displacement will be ~11 mm. Maximal value of the fuel temperature exceeds the temperature of the central fuel element by no more than 3 K.Maximal change in the material temperature in the fuel elements section is not more than 1 K. Since the thermal conductivity of the fuel element repeatedly exceeds the thermal conductivity of the gas medium, the temperature distribution over the cross-section of the cask is uneven.


## 5. Conclusion

The temperature field in the average height of the cross-section of the cask with the fuel assemblies of BN-350 reactor for dry storage was determined as a result of the calculations.Results of thermal calculation allow considering the optimum alternation for storage, where temperature of construction elements will have minimum value. Requirements for BN-350 SNF storage safety assume that as a result of passive cooling temperature of fuel element cladding will not exceed 673 K [10].The calculated maximum fuel temperature is equal to 745 K, obtained in case of cask filling with argon with the concrete thickness of 1 m and FA power of 190 W.The minimum temperature of the fuel among the considered calculation options is 516 K (the maximum in the model for this calculation option) was obtained in the case of cask filling with air, with concrete thickness of 0.5 m and FA power of 120 W.When using air as cask filler instead of argon, the maximum temperature of the model will be reduced by ~54 K at FA power of 190 W and at ~37 K at FA power of 120 W, which is due to the higher conductivity of air compared to argon.It is determined that the maximum temperature in FA is shifted to the center of the cask in the result of calculations. The temperature of the central fuel element is below the maximum temperature by ~3 K. If the cavity of the cask is filled with argon, the displacement of the maximum temperature from the FA axis to the center of the cask will be ~7 mm, if air is used as a gas medium, then displacement will be ~11 mm.



---

*Source: 3045897-2019-04-11.xml*