Journal Search Engine

View PDF Download PDF Export Citation Korean Bibliography PMC Previewer
ISSN : 1738-1894(Print)
ISSN : 2288-5471(Online)
Journal of Nuclear Fuel Cycle and Waste Technology Vol.20 No.4 pp.411-428
DOI : https://doi.org/10.7733/jnfcwt.2022.043

Development of Model to Evaluate Thermal Fluid Flow Around a Submerged Transportation Cask of Spent Nuclear Fuel in the Deep Sea

Guhyeon Jeong, Sungyeon Kim, Sanghoon Lee*
Keimyung University, 1095, Dalgubeol-daero, Dalseo-gu, Daegu 42601, Republic of Korea
* Corresponding Author.
Sanghoon Lee, Keimyung University, E-mail: shlee1222@kmu.ac.kr, Tel: +82-53-580-5264

November 10, 2022 ; November 28, 2022 ; December 6, 2022

Abstract


Given the domestic situation, all nuclear power plants are located at the seaside, where interim storage sites are also likely to be located and maritime transportation is considered inevitable. Currently, Korea does not have an independently developed maritime transportation risk assessment code, and no research has been conducted to evaluate the release rate of radioactive waste from a submerged transportation cask in the sea. Therefore, secure technology is necessary to assess the impact of immersion accidents and establish a regulatory framework to assess, mitigate, and prevent maritime transportation accidents causing serious radiological consequences. The flow rate through a gap in a containment boundary should be calculated to determine the accurate release rate of radionuclides. The fluid flow through the micro-scale gap can be evaluated by combining the flow inside and outside the transportation cask. In this study, detailed computational fluid dynamic and simplified models are constructed to evaluate the internal flow in a transportation cask and to capture the flow and heat transfer around the transportation cask in the sea, respectively. In the future, fluid flow through the gap will be evaluated by coupling the models developed in this study.



초록


    1. Introduction

    According to data released by Korea Hydro & Nuclear Power in June 2022 as Fig. 1, the SNF (Spent Nuclear Fuel) storage capacity of the temporary storage facility at the nuclear power plant is expected to be saturated soon. Therefore, transportation from the current temporary storage facility to the intermediate storage facility is required.

    JNFCWT-20-4-411_F1.gif
    Fig. 1

    Status of SNF storage [1] (LWR: light water reactor, HWR: heavy water reactor).

    Such transportation is largely divided into land transportation using vehicles and maritime transportation on the sea. According to the current status of domestic nuclear power plants provided by the Korea Atomic Industry Forum in October 2022 [2], all nuclear power plants are located on seaside. Considering this domestic environment, maritime transportation of radioactive waste for interim storage and disposal is an unavoidable management option. Evaluating the risk caused by accidents during the transportation of radioactive waste is a major item of radiation safety regulation technology [3-4].

    In particular, risk assessment for transportation accidents of high-level radioactive waste such as SNF is a key technology necessary to secure public safety from the radiological hazard of radioactive waste. One of the dangerous situations that can occur during maritime transportation is that the transportation cask loaded with radioactive waste is lost in the deep sea due to a ship accident, etc., and the contents of the transportation cask are released into the ocean. Evaluating how the released radioactive material will diffuse in the ocean and the extent to which the diffused nuclide will cause exposure to the public through various pathways including the food chain of ocean ecosystems are the key factors in risk assessment of the maritime transportation of radioactive waste. In foreign countries such as the United States, France, and Japan, studies to evaluate the risk of radioactive material transportation by sea and the sinking of the transportation cask have been conducted using codes such as MARINRAD, POSEIDON, and Barrier Effect Model [5-6].

    There has been no reported case of self-developed codes in Korea similar to the codes mentioned above, and related technology development is necessary. The abovementioned studies focus on the assessment of the ocean diffusion of radioactive materials and their effects. In MARINRAD and POSEIDON, the effect of barrier effect due to the engineered barrier existing in the transportation cask was ignored. In CRIEPI’s Barrier Effect Model, it was simplified by conservative assumptions and reflected in the calculation.

    A type-B transportation cask is a very robust system and its containment function might not be completely lost even in the deep sea. Therefore, the risk assessment performed without properly considering the barrier effect is likely to give very conservative results. The goal of the second stage of this study is to reflect the barrier effect of the containment system of the transportation cask in the calculation of release rate of its radioactive contents with a scientific and engineering method using CFD. Through this, it would be possible to more realistically determine the release rate of radioactive material from the lost transportation cask in the deep sea.

    As part of the model development, the goal of the first step is to develop detailed flow models inside and around the transportation cask to understand the detailed flow and heat transfer in and around the cask. The results will be coded and linked with the maritime diffusion assessment code and risk assessment code. Thus, it is intended to contribute to the development of an integrated maritime transportation risk assessment code of radioactive waste reflecting the domestic ocean environment.

    2. CRIEPI’s Barrier Effect Model

    2.1 CRIEPI’s Scenario of Radioactive Material Release

    CRIEPI’s Barrier Effect Model is a model to calculate the release rate of radioactive material from a submerged transportation cask where the release rate (Ci·yr−1) is kept smaller than the leaching rate of the contents due to the barrier effect of the transportation cask. The release scenario considered by CRIEPI is shown in Fig. 2.

    JNFCWT-20-4-411_F2.gif
    Fig. 2

    CRIEPI’s dose assessment flowchart [7-9].

    It was assumed that all transportation casks could be salvaged in the near sea below 200 m depth and that there is no barrier effect by the transportation cask submerged in the sea deeper than 200 m. In addition, the barrier effect by the transportation cask was calculated assuming that the transportation cask was intact and only the sealing material was damaged in the sea of 200 m depth. In this case, the process of release rate calculation is as follows.

    JNFCWT-20-4-411_F3.gif
    Fig. 3

    21 SNF Transportation cask considered in the study.

    • transportation cask submergence → O-ring breakage → seawater inflow → exposure of radioactive contents to seawater → leaching of contents → increase of radioactive material concentration → seawater outflow from transportation cask

    In this case, radioactive material can be released through diffusion, but the effect of convection is dominant. The driving force of the seawater flow into and out of the transportation cask is the buoyant force generated by the temperature difference. When the flow rate is not large, the concentration of radioactive material in the transportation cask reaches the solubility limit and the release rate is kept below the leaching rate of the radioactive material.

    JNFCWT-20-4-411_F4.gif
    Fig. 4

    Schematic of the barrier effect [7].

    2.2 Release Rate Calculation

    • - Buoyancy by natural convection

      F = Δ ρ g = ρ g β Δ T
      (1)

    • - Calculation of flow rate using energy conservation

      F δ = ρ 1 2 u m 2 ( 1 + λ f L d e )
      (2)

      u m = { 2 F δ ( 1 + λ f L d e ) ρ } 1 / 2
      (3)

      q = u m A
      (4)

    • - The nuclide release rate by natural convection

      R o = C q
      (5)

    • - The nuclide leaching rate from the contents into seawater

      R c = R p μ Q
      (6)

    • - Nuclide inventory decrease rate

      d Q d t = R c λ Q
      (7)

    • - Change of nuclides concentration

      d C d t = R C V R O V λ C Constraint: C 0 R c R o , C C s
      (8)

    CRIEPI’s barrier effect model is a computational model established based on the conservation of mass for individual nuclides assuming that flow occurs due to natural convection. The phenomenon in which the outflow rate of nuclides due to convection is smaller than the leaching rate of nuclides is called the barrier effect. Depending on the damage to the transportation cask, it is determined whether the barrier effect is applicable or not. As shown in Fig. 5, the transportation cask must maintain the integrity of the containment boundary against design-based accidents such as 9 m drop, fire, puncture, and water immersion [3-4]. When the containment boundary of the transportation cask is damaged and gap occurs, the size is expected to be very small and it is reasonable to consider the barrier effect. However, the CRIEPI’s model used simple formulas for the concentration of nuclides and the flow rate of seawater inside the transportation cask, rather than using precise computational analysis. Therefore, the following supplements to the CRIEPI’s barrier effect model are suggested.

    JNFCWT-20-4-411_F5.gif
    Fig. 5

    Hypothetical accident conditions for type-B transportation cask.

    Table 1

    Nuclear fuel storage ratio by nuclear power plant [10]

    JNFCWT-20-4-411_T1.gif
    • - Reflection of changes in seawater temperature due to changes in the ocean environment

    • - Reflecting the temperature change in the transportation cask considering the residual amount of nuclides

    • - Addition of a calculation module for input factors related to flow specifications according to the condition of the transportation cask

    • - Reflection of fuel damage ratio (FDR) in the calculation

    • - Calculation and reflection of the critical depth of seawater that can generate damages to the transportation cask and fuel cladding

    2.3 Modified Scenarios of Nuclide Release From a Submerged Cask

    In this section, a modified scenario of radio-nuclide release into seawater is proposed in Fig. 6. It is a regulatory requirement that a transportation cask should maintain its structural integrity under seawater of depth 200 m. Therefore, in the transportation cask whose containment boundary is not damaged, there is no release of contents until time passes the corrosion resistance life in the ocean environment at the depth of 200 m or less [4]. If the transportation cask is lost in the seawater of depth exceeding 200 m, it is necessary to obtain the critical depth at which the integrity of the containment boundary is compromised through the evaluation of the pressure limit of the transportation cask. If the transportation cask experienced a beyond-designbasis accident which causes a permanent damage to the containment function, the seawater can flow into the transportation cask and initiate the leaching of the nuclides which are exposed to the seawater. Then the concentration of nuclides in the seawater in the transportation cask cavity increases. When this seawater leaks to the outside, the nuclide is released into the ocean. If the release rate of nuclides through the containment boundary is less than the rate at which the nuclides are leached into the cavity seawater, the seawater in the cavity is saturated and the leaching rate is controlled by the flow rate of the seawater through the containment boundary.

    JNFCWT-20-4-411_F6.gif
    Fig. 6

    Radioactive material ocean release scenario.

    JNFCWT-20-4-411_F7.gif
    Fig. 7

    PLUS7 fuel assembly.

    3. Analysis of Internal Flow of SNF Transportation Cask

    3.1 SNF Model Selection

    The SNF, which is used as the basis for the development of the internal flow analysis model of the transportation cask, was selected. According to the paper [10] published by Cho et al., the PLUS7 occupies the biggest portion of 20.7% of the total SNF storage and it is reasonable to consider PLUS7 as the reference for the development.

    3.2 Development of a Simplified Model of Spent Fuel Assembly

    If the PLUS7 fuel assembly is modeled in detail and used for analysis, it is inefficient in terms of calculation time and cost. Therefore, a simplified model based on the porous media proposed by the U.S. DOE (Department of Energy) was used [11].

    The porous model is a model configured to conduct heat equally based on thermal conductivity as shown in Fig. 8. To develop a porous model, anisotropic thermal conductivity to each direction is required.

    JNFCWT-20-4-411_F8.gif
    Fig. 8

    Porous model [11].

    • - The formula for effective thermal conductivity in the x, y axis direction is as follows [11].

    : It is the heat diffusion equation for the steady temperature of a square media with heat generation. Assuming isotropic thermal conductivity it is as follows.

    T ( x , y ) = q ( a 2 x 2 ) 2 k 16 q a 2 k π 2 n = 0 ( 1 ) n cos [ ( 2 n + 1 ) π x 2 a ] cosh [ ( 2 n + 1 ) π y 2 a ] ( 2 n + 1 ) 3 cosh [ ( 2 n + 1 ) π 2 ]
    (9)

    : The maximum temperature occurring at the center of the cross-sectional area (x = y = 0) is as follows.

    T ( 0 , 0 ) = q ( a 2 ) 2 k 16 q a 2 k π 2 n = 0 ( 1 ) n ( 2 n + 1 ) 3 cosh [ ( 2 n + 1 ) π 2 ]
    (10)

    : Due to the strong divergence of cosh, the sum converges quickly and is as follows.

    T ( 0 , 0 ) = q ( a 2 ) k ( 0.2947 )
    (11)

    : Assuming the wall temperature is not zero, and substituting q''' = Q/4a2La gives as follows.

    k e = Q 4 L a ( T o T s ) ( 0.2947 )
    (12)

    : The thermal conductivity of PLUS7 is as follows.

    k e = 0.2501 4.3 × 10 4 T m + 2.38 × 10 6 T m 2 [ W m 1 K 1 ]
    (13)

    T m = T o + T s 2 [ K ]
    (14)

    • - An area-weighted average method was used for the formula for effective thermal conductivity in the z axis direction.

      k e , a x i a l = i = 1 n A i A k i
      (15)

    • - The finally obtained anisotropic effective thermal conductivity is summarized in Table 2.

    Table 2

    Anisotropic effective thermal conductivity

    JNFCWT-20-4-411_T2.gif

    Based on the obtained anisotropic thermal conductivity, SNF PLUS7 was simplified and used for CFD.

    3.3 CFD Model of Flow Inside the Spent Fuel Cask

    FLUENT 2021 R1, a general-purpose thermal fluid analysis code, was used for the analysis. Because the flow path created at the containment boundary due to a beyond design basis accident is expected to be very small and the flow rate through this gap is limited, the thermal fluidic analysis inside the transportation cask was performed assuming intact containment boundary. In this analysis model, the standard k-ε model commonly used for turbulence analysis is utilized. In addition to the Navier-Stokes equations of motion and continuity equations, related expression is [12]:

    • - Turbulent kinetic energy

      x i ( ρ k u i ) = x j [ ( μ + μ t σ k ) k x j ] + G k + G b ρ ε Y M + S k
      (16)

    • - Dissipation rate

      x i ( ρ u i ) = x j [ ( μ + μ t σ k ) x j ] + C 1 k ( G k + C 3 G b ) C 2 ρ 2 k + S
      (17)

    • - Coefficient of viscosity

      μ t = ρ C μ k 2 ε
      (18)

    In this study, as a preparatory stage for the development of the CFD model simulating the barrier effect, a three-dimensional 1/2 full model with symmetry was created by referring to the transportation cask capable of transporting 21 PWR (Pressurized Water Reactor) SNF as shown in Fig. 9.

    JNFCWT-20-4-411_F9.gif
    Fig. 9

    Shape (Left: model shape, Right: mesh shape /number of mesh: 3,923,189).

    The body of the referenced PWR transportation cask is made of carbon steel (SA-350 LF3), and the upper part of the transportation cask is sealed with a double lid. A neutron shielding material (NS-4-FR) is installed on the outer shell of the transportation cask. The material of the shock absorber installed on the upper and lower parts of the transportation cask was assumed to be Balsa wood. For the reference fuel considered in this transportation cask analysis, the data of decay heat according to the burnup and cooling period of PLUS7 presented in were referred as in Fig. 10 [13]. In this study, a burnup rate of 45 GWd/MTU and a cooling period of 10 years were assumed. The decay heat of 1 bundle of nuclear fuel assembly was 800 W, and since 21 assemblies were loaded, the total decay heat was set to 16.8 kW. The vertically upright posture of the transportation cask is assumed where the boundary layer is expected to develop the most. The natural convection situation was assumed by applying gravitational acceleration toward the bottom of the transportation cask. In addition, it was assumed that the inside of the transportation cask was filled with seawater due to the accident condition. The external condition is the temperature of the seawater outside the transportation cask (T = 15℃). 300 [W·m−2− K−1] was entered for the convective heat transfer coefficient to determine the temperature of the surface of the transportation cask [15]. The pressure inside and outside the container was assumed the same due to the breached containment boundary.

    JNFCWT-20-4-411_F10.gif
    Fig. 10

    Cooling Time-Decay Heat by PLUS7 burn rate [13].

    The material properties for the flow analysis inside the transportation cask are shown in Table 3~Table 8 [16-20].

    Table 3

    Thermal properties of SA-350 LF3

    JNFCWT-20-4-411_T3.gif
    Table 4

    Thermal properties of SA-240 Type-304

    JNFCWT-20-4-411_T4.gif
    Table 5

    Thermal properties of balsa wood

    JNFCWT-20-4-411_T5.gif
    Table 6

    Thermal properties of ceramics

    JNFCWT-20-4-411_T6.gif
    Table 7

    NS-4-FR Thermal Properties

    JNFCWT-20-4-411_T7.gif
    Table 8

    Thermal properties of seawater

    JNFCWT-20-4-411_T8.gif
    Table 9

    Surface Heat Coefficient

    JNFCWT-20-4-411_T9.gif

    4. Analysis of External Flow Around Transportation Cask

    The posture of the transportation cask submerged in the deep sea is random. Considering the shape of SNF transportation cask and the posture during the transportation, the horizontal position where the long side of cask is horizontal is most likely. In this study, the seawater flow and heat transfer around the submerged transportation cask are evaluated with CFD code, Fluent. Because the shape of the seabed is diverse, various postures can be made. However, as in the drop test of a transportation cask, vertical and horizontal postures are two limiting cases in terms of release rate. Therefore, the analyses are performed for the two postures, the horizontal and the vertical.

    Purpose of these analyses are as follows:

    • - Determine the analysis area for simulating the detailed seawater flow around the transportation cask submerged in the deep sea.

    • - Evaluate how the convective heat transfer coefficient between the transportation cask and the surrounding seawater changes depending on environmental factors (seawater temperature, seawater velocity, and cask posture).

    • - Creating a boundary layer in the vicinity of lid gap and identifying the detailed flow in this region.

    • - Identify severer posture in terms of release rate of contents.

    Finally, this analysis result is linked with the detailed heat flow evaluation model of the transportation cask in Section 3 and the assumptions made for the analyses are verified. Seawater flow into and out of the transportation cask will be evaluated in our future study and the results obtained in current work will be used as basic data for flow evaluation in microchannels created at the containment boundary.

    4.1 Simplification of Spent Nuclear Fuel Transportation Cask

    In this analysis, the transportation cask was modeled as a cylinder to which the surface temperature was given. The cuboidal seawater area around the transportation cask was selected as the area of analysis. In addition, the material properties for analysis are the same as those of the internal model [16-20].

    4.2 Setting of Analysis Area for Flow Evaluation Around Transportation Cask

    In general, when performing external flow analysis on a model with a height of H and a width of W, the height of the external analysis area should be at least 5 times the height of the model and 10 times the depth of the model [21]. After setting the size according to the above criteria, the value of the convective heat transfer coefficient between the transportation cask and seawater was calculated, and the size was reduced to the extent that this value does not change significantly.

    4.3 External Flow Model of the SNF Transportation Cask

    FLUENT 2021 R1, a general-purpose thermal fluid analysis code, was used for the analysis. For the flow model outside the transportation cask, the flow at the transportation cask surface is important. k-ω SST (shear stress transport) model with relatively accurate results in the wall boundary region was used. In addition to the Navier-Stokes equations of motion and continuity equations, related expression is [22-23]:

    • - Kinematic eddy viscosity

      v T = α 1 k m a x ( α 1 ω , S F 2 )
      (19)

    • - Turbulence kinetic energy

      k t + U j k x j = P k β * k ω + x j [ ( v + σ k v T ) k x j ]
      (20)

    • - Dissipation Rate

      ω t + U j ω x j = a S 2 β ω 2 + x j [ ( v + σ ω v T ) ω x j ] + 2 ( 1 F 1 ) σ ω 2 1 ω k x i ω x i
      (21)

    • - F1 mixing function

      JNFCWT-20-4-411_EQ22.gif
      (22)

    • - CD

      C D k ω = m a x ( 2 ρσ ω 2 1 ω k x i ω x i , 10 20 )
      (23)

    • - F2 mixing function

      F 2 = t a n h [ [ m a x ( 2 k β 2 ω y , 500 v y 2 ω ) ] 2 ]
      (24)

    • - Pk

      P k = m i n ( τ i j U i x j , 10 β * k ω )
      (25)

    This model is a mixture of the two models to apply the k-ω model for flows near the wall and the k-ε model for flows away from the wall. This model can be used without problems for multiple target problems. In general, it is known to provide good results near the wall boundary [24].

    In addition, a wall function [25] was used to determine whether all the flow on the wall is expressed in the boundary mesh layer to accurately simulate the flow on the wall as shown in Fig. 11. The wall function describes the empirical behavior on the wall as a dimensionless velocity and vertical distance from the wall and is expressed as follows.

    JNFCWT-20-4-411_F11.gif
    Fig. 11

    How to accurately determine the flow in a wall [21].

    The dimensionless normal distance on the surface is as follows [25].

    y + = y u τ v
    (26)

    • - The friction velocity is as follows

      u τ = τ w ρ
      (27)

    • - The dimensionless velocity is as follows

      U + = u u τ
      (28)

    Among them, the condition of y+ ≤ 5, which corresponds to a viscous sub-layer that reflects the effect on viscosity at the boundary layer wall well, was used as shown in Fig. 12. Mesh was created using an inflation layer on the boundary around the cylindrical surface. When checking the dimensionless distance from the boundary layer wall, it was confirmed that y+ = 0.648 and the condition of the viscous sub layer was satisfied.

    JNFCWT-20-4-411_F12.gif
    Fig. 12

    Wall function [25].

    In this model, a steady state in vertical and horizontal orientations was considered. For the condition of the outside of the transportation cask, the temperature of the external seawater T = 15℃ was assigned to the same as the internal model in Section 3. As shown in Fig. 13, the external seawater was assumed to be in a forced convection situation with a flow velocity of v = 0.1 m·s‒1.

    JNFCWT-20-4-411_F13.gif
    Fig. 13

    External flow model of transportation cask.

    JNFCWT-20-4-411_F14.gif
    Fig. 14

    Mesh shape (number of mesh: 13,707,670).

    5. Results

    5.1 Analysis Result of Internal Flow Model of SNF Transportation Cask

    The temperature distributions inside and on the surface of the transportation cask are depicted in Fig. 15. The maximum temperature of the internal seawater rises to 102℃ as shown in Fig. 15 and the temperature difference between the inside and outside of the transportation cask reaches 87℃. The average temperature of the transportation cask surface is evaluated as 20℃.

    JNFCWT-20-4-411_F15.gif
    Fig. 15

    Temperature distribution by natural convection (Left: inside, Right: outside).

    The flow in the transportation cask is shown in Fig. 16. The flow is generated by natural convection due to the decay heat of fuel assemblies. From the simulation results, it was confirmed that turbulence occurred in the middle of the transportation cask and the flow circulates to the top and bottom of the transportation cask. As a result, a flow velocity of up to 1.72 m·s‒1 occurred at the bottom of the transportation cask.

    JNFCWT-20-4-411_F16.gif
    Fig. 16

    Velocity distribution by natural convection (Left: pathlines, Right: vectors).

    5.2 Appropriate Size of Analysis Area for the External Flow Simulation

    The transportation cask has a large diameter and the flow around it has a big Reynolds number. This implies the development of turbulent flow paths. Therefore, the analysis area for the downstream was set large to account for the cylindrical wake well. If the analysis is performed in a flow field smaller than the above conditions, an error such as a backflow occurs at the outlet, and proper convergence may not be achieved. Therefore, for the flow field simulation, it is reasonable to set the height of analysis area to twice the height of the transportation cask, the depth to be eight times the height of the transportation cask, the upstream to twice the height of the transportation cask, and the downstream to twelve times. Therefore, the size of the flow field used in the analysis is Width: 28,000 mm, Depth: 16,000 mm, and Height: 10,000 mm.

    Table 10

    Heat transfer coefficient versus analysis area

    JNFCWT-20-4-411_T10.gif

    5.3 Analysis Result of External Flow Model of SNF Transportation Cask

    The results of the external flow simulation are shown in Figs. 1720. The surface heat transfer coefficient is calculated as h = 362.8 [W·m−2− K−1] in the vertical orientation, and h = 313.2 [W·m−2− K−1] in the horizontal orientation. The heat transfer coefficient is smaller when the transportation cask is in horizontal orientation and a bigger temperature difference between the inside and outside the transportation cask is expected. This will result in bigger buoyant force inside the transportation cask and eventually a large flow rate through the breached containment boundary of the transportation cask.

    JNFCWT-20-4-411_F17.gif
    Fig. 17

    Velocity pathlines when transportation cask is vertical.

    JNFCWT-20-4-411_F18.gif
    Fig. 18

    Velocity pathlines when transportation cask is horizontal.

    JNFCWT-20-4-411_F19.gif
    Fig. 19

    Temperature distribution when transportation cask is vertical.

    JNFCWT-20-4-411_F20.gif
    Fig. 20

    Temperature distribution when transportation cask is horizontal.

    In the internal flow model of the SNF transportation cask, the heat transfer coefficient was input as 300 [W·m−2− K−1] on the surface. When the decay heat of nuclear fuel is transferred to the surface of the transportation cask and thermal equilibrium is achieved with the temperature of the seawater outside the transportation cask with T = 15℃, the surface temperature is Ts = 20℃. The heat transfer coefficient calculated in this external flow model does not show significant difference with the input value 300 [W·m−2− K−1] for the internal flow model. Therefore, it is shown that the input parameters of the internal and external flow models are physically reasonable and can be used in the development of the detailed flow model around the gap created at the containment boundary of the transportation cask.

    It is expected that when the transportation cask is in horizontal position, bigger flow rate into and out of cask is expected than when in vertical position. In horizontal position, the bottom and top of the lid gap become the inlet and outlet of buoyant flow and the whole lid gap is used as flow passage. In vertical position, the inlet and outlet of flow share the lid gap and the effective flow area becomes smaller as shown in Fig. 21.

    JNFCWT-20-4-411_F21.gif
    Fig. 21

    Flow pattern at the lid gap of transportation cask.

    6. Discussion and Conclusions

    This study was carried out as part of the development of a CFD model simulating the barrier effect. In this study, a modified scenario of release of radio-nuclides from a submerged SNF transportation cask was presented. A detailed CFD model was built to simulate the fluid flow inside the transportation cask and simplified models were developed to simulate the fluid flow and heat transfer around the submerged transportation cask in the ocean.

    - A modified radioactive material release scenario was presented in consideration of the condition of the lost transportation cask and the depth of the sea where the cask is located. If the extent of damage is limited to the loss of O-ring function of the transportation cask, the barrier effect is significant, and a much smaller release rate than the leaching rate of nuclides is predicted.

    • - A detailed analysis model of the transportation cask was constructed as the preparatory step to perform the fluid analysis of the transportation cask with damage to the containment boundary.

    • - The external flow model of transportation cask was developed. The proper analysis area around the transportation cask was identified based on the flow characteristics and resulting heat transfer from the transportation cask.

    • - It is confirmed that the temperature difference between the inside and outside of the transportation cask is bigger in the case of the horizontal position because the heat transfer coefficient is lower. It will create a bigger buoyancy inside the transportation cask and a larger flow velocity through the damaged containment boundary. Therefore the horizontal posture is expected to be a conservative posture in terms of the release rate and this will be confirmed in our future study.

    The fluid flow through the gap in the containment boundary of the transportation cask should be evaluated by combining the fluid flow inside and outside the transportation cask. In our future research, the two models developed in this work will be combined to accurately simulate the fluid flow in the gap of micro-size and eventually calculate accurate release rate of nuclides from the transportation cask to the ocean.

    Acknowledgements

    This work was supported by the Nuclear Safety Research Program through the Korea Foundation Of Nuclear Safety (KoFONS) using the financial resource granted by the Nuclear Safety and Security Commission (NSSC) of the Republic of Korea (No. 2106042).

    Nomenclature

    • ρ : Density

    • g : Gravitational acceleration

    • k : Thermal conductivity

    • β : Coefficient of thermal expansion of seawater

    • T : Temperature

    • δ : Width of seawater channel

    • um : Average flow rate

    • λf : Friction loss coefficient of flow path (= 64 Re−1)

    • Re : Reynolds number (= umdeν−1)

    • ν : Kinematic viscosity coefficient of seawater (1.22×10−6 m2s−1)

    • L : Length of seawater channel

    • q : Flux

    • A : The projected area of the seawater channel

    • C : Nuclide concentration [Bq s−1]

    • Rp : Leaching rate [gm−2s−1]

    • μ : Specific surface area [m2g−1]

    • Q : Residual nuclide amount in container [Bq]

    • λ : Decay constant of nuclide

    • V : Volume m3

    • Cs : Solubility limit

    • ke : Transverse effective thermal conductivity [W·mK− 1]

    • ke,axial : Axial effective thermal conductivity [W·mK− 1]

    • To : SNF assembly center temperature [℃]

    • Ts : SNF assembly surface temperature [℃]

    • Tm : Assembly average temperature [℃]

    • La : SNF assembly length [m]

    • T : Seawater temperature

    • yp : Distance from the wall to the center point of the cell

    • Up : Velocity of the cell’s center point

    • y+1 : Dimensionless distance from the boundary layer wall

    • h : Convection heat transfer coefficient

    • μt : Represents eddy viscosity

    • Eij : Represents component of rate of deformation

    • υi : Represents velocity component in corresponding direction

    • Cμ : Experimental constant value (= 0.09)

    • σk : Experimental constant value (= 1.00)

    • σε : Experimental constant value (= 1.30)

    • C1ε : Experimental constant value (= 1.44)

    • C2ε : Experimental constant value (= 1.92)

    • β* : Experimental constant value (= 0.09)

    • α1 : Experimental constant value (= 0.31)

    Figures

    Tables

    References

    1. Public Data Portal Korea Hydro & Nuclear Power Co., Ltd. October 26 2022. “Status of Spent Nuclear Fuel Storage by Nuclear Power Plant.” Public Data Portal. Accessed Oct. 30 2022. Available from: https://www. data.go.kr/data/15060363/fileData.do.
    2. Korea Atomic Industry Forum. “Status of Nuclear Power Plants in Korea.” KAIF. Accessed Oct. 25 2022. Available from: https://www.kaif.or.kr/ko/ko/?c=189&s=188.
    3. International Atomic Energy Agency, Regulations for the Safe Transport of Radioactive Material, Specific Safety Requirements No. SSR-6 (Rev.1) (2018).
    4. Notice of Nuclear Safety Commission, Regulations on Packaging and Transport of Radioactive Materials, etc., No. 2011-49, NSC, Korea (2011).
    5. K.J. Connolly and R.B. Pope, A Historical Review of the Safe Transport of Spent Nuclear Fuel, FCRDNFST- 2016-000474, US Department of Energy (2016).
    6. J. Edwards, A. Hough, J.A.C. Marples, and T. Ohe, “Leaching of Unirradiated Mox Fuel in Sea Water”, Int. J. Radioac. Mater. Transp., 9(2), 147-156 (1998).
    7. D. Tsumune, T. Tsubono, K. Misumi, and Y. Yoshida, “Estimation of the Radiation Dose Equivalent for the Hypothetical Submergence of a Sea-Transport Package of Low-level Radioactive Waste”, J. Nucl. Sci. Technol., 54(6), 681-693 (2017).
    8. D. Tsumune, T. Saegusa, H. Suzuki, N. Watabe, H. Asano, K. Maruyama, and Y. Kinehara, “Dose Assessment for the Public Due to Packages Shipping Radioactive Materials Hypothetically Sunk on a Continental Shelf”, Int. J. Radioac. Mater. Transp., 11(4), 317-328 (2000).
    9. D. Tsumune, T. Saegusa, H. Suzuki, K. Maruyama, C. Ito, and N. Watabe, “Estimated Radiation Dose From a Mox Fuel Shipping Package That is Hypothetically Submerged Into Sea”, Int. J. Radioac. Mater. Transp., 11(3), 239-253 (2000).
    10. D.K. Cho, J. Kim, I.Y. Kim, and J.Y. Lee, “Investigation of PWR Spent Fuels for the Design of a Deep Geological Repository”, J. Nucl. Fuel Cycle Waste Technol., 17(3), 339-346 (2019).
    11. R.H. Bahney and T.L. Lotz. Spent Nuclear Fuel Effective Thermal Conductivity Report, US Department of Energy Technical Report, BBA000000-01717-5705- 00010 REV 00 (1996).
    12. B.E. Launder and D.B. Spalding, “The Numerical Computation of Turbulent Flows”, Comput. Methods Appl. Mech. Eng., 3(2), 269-289 (1974).
    13. Bechtel SAIC Company, Characteristics for the Representative Commercial Spent Fuel Assembly for Preclosure Normal Operations, 000-PSA-MGR0-00700- 000-00A (2007).
    14. O.M. Elmardi, Solution to Problems in Heat Transfer Transient Conduction or Unsteady Conduction, Anchor Academic Publishing, Hamburg (2017).
    15. International Atomic Energy Agency, Regulations for the Safe Transport of Radioactive Material, Safety Requirements No. TS-R-1 (2005).
    16. The American Society of Mechanical Engineers, ASME Boiler & Pressure Vessel Code Sec.Ⅱ Part D Properties, United States (2015).
    17. Holtec International. Safety Analysis Report the HISTAR 63 Package, Holtec Report, HI-2073777 (2009).
    18. A. Bejan, Convection Heat Transfer, 4th Ed., John Wiley & Sons Inc., New Jersey (2013).
    19. Korea Institute of Machinery & Materials. Development of Process Technology and Design of Al Composite Materials for Neutron Shielding, KIMM Report, UCN-5203.C (2019).
    20. H.S. Carslaw and J.C. Jaeger, Conduction of Heat in Solids, 2nd Ed., Oxford University Press, Oxford (1959).
    21. E.M. John, An Introduction to Ansys Fluent 2022, SDC Publication, Kansas (2022).
    22. F.R. Menter, “Zonal Two Equation k-ω Turbulence Models for Aerodynamic Flows”, Proc. of 23rd Fluid Dynamics, Plasmadynamics, and Lasers Conference, AIAA-93-2906, AIAA, Orlando, FL (1993).
    23. F.R. Menter, “Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications”, AIAA J., 32(8), 1598-1605 (1994).
    24. J. Paik and N.J. Lee, “Numerical Modeling of Free Surface Flow Over a Broad-Crested Rectangular Weir”, J. Korea Water Resour. Assoc., 48(4), 281-290 (2015).
    25. V.K. Th, Mechanische Ähnlichkeit und Turbulenz, Nachrichten von der Gesellschaft der Wissenschaften zu Göttingen, Mathematisch-Physikalische Klasse, 58-76 (1930).

    Editorial Office
    Contact Information

    - Tel: +82-42-861-5851, 866-4157
    - Fax: +82-42-861-5852
    - E-mail: krs@krs.or.kr

    SCImago Journal & Country Rank