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.
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 [
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.
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 Figure
The layout of the fuel assembly in six-place cask for dry storage.
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.
Calculation is conducted with ANSYS software, which implements the method of finite elements [
For calculations two-dimensional model was used which is provided in Figures
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:
wherein
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 [
Coefficient of convection is determined by the following formula:
When
Grashof number is equal to
Prandtl number is as follows:
The average calculated value of the convection coefficient when filling the cask cavity with argon under the FA power of 120 W is
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 Tables
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×°С) | 1130 |
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 | 1000 |
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 | 519 |
According to experimental data, calculation error in ANSYS is from 5 to 10 %.
Figures
Temperature field of the cask and the concrete around.
Temperature field of the cask.
Temperature field of FA.
Table
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 | 689 |
The 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.
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 [
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.
All data of this study are included within the article.
The authors declare that there are no conflicts of interest regarding the publication of this paper.
M.K. Skakov contributed to task definition; Ye. T. Koyanbayev collected baseline data for calculations; D. A. Ganovichev and Ye. A. Martynenko contributed to calculation; M. K. Skakov, А. А. Sitnikov, Ye. T. Koyanbayev, and D. A. Ganovichev contributed to result analysis and conclusions; Ye. T. Koyanbayev and D. A. Ganovichev are responsible for article writing; Ye. T. Koyanbayev, D. A. Ganovichev, and Ye. A. Martynenko provided feedback on reviewers’ comments.
The research was financially supported by the Committee of Science of the Ministry of Education and Science of the Republic of Kazakhstan for 2018