Nuclear Engineering and Design

Similar documents
Hideki Kamide 1. A. Ono 1, N. Kimura 1, J. Endoh 2 and O. Watanabe 2

CFD STUDIES IN THE PREDICTION OF THERMAL STRIPING IN AN LMFBR

Application of Grid Convergence Index Estimation for Uncertainty Quantification in V&V of Multidimensional Thermal-Hydraulic Simulation

NUMERICAL APPROACH OF SELF-WASTAGE PHENOMENA IN STEAM GENERATOR OF SODIUM-COOLED FAST REACTOR

Department of Engineering and System Science, National Tsing Hua University,

NATURAL CONVECTION HEAT TRANSFER CHARACTERISTICS OF KUR FUEL ASSEMBLY DURING LOSS OF COOLANT ACCIDENT

ICONE18- ICONE NUMERICAL STUDY ON PASSIVE CONTROL OF THERMAL STRIPING PHENOMENON USING LORENTZ FORCE IN FAST REACTOR

Analysis and interpretation of the LIVE-L6 experiment

EFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION

Numerical Simulation of the MYRRHA reactor: development of the appropriate flow solver Dr. Lilla Koloszár, Philippe Planquart

HEAT TRANSFER CAPABILITY OF A THERMOSYPHON HEAT TRANSPORT DEVICE WITH EXPERIMENTAL AND CFD STUDIES

Comparison between Numerical and Experimental for UVP Measurement in Double Bent Pipe with Out-of-Plane Angle

PREDICTION OF MASS FLOW RATE AND PRESSURE DROP IN THE COOLANT CHANNEL OF THE TRIGA 2000 REACTOR CORE

Coolant Flow and Heat Transfer in PBMR Core With CFD

3D CFD and FEM Evaluations of RPV Stress Intensity Factor during PTS Loading

Developments and Applications of TRACE/CFD Model of. Maanshan PWR Pressure Vessel

Instability Analysis in Peach Bottom NPP Using a Whole Core Thermalhydraulic-Neutronic Model with RELAP5/PARCS v2.7

A 2-D Test Problem for CFD Modeling Heat Transfer in Spent Fuel Transfer Cask Neutron Shields. G Zigh and J Solis U.S. Nuclear Regulatory Commission

QUALIFICATION OF A CFD CODE FOR REACTOR APPLICATIONS

Onset of Flow Instability in a Rectangular Channel Under Transversely Uniform and Non-uniform Heating

A Study on Hydraulic Resistance of Porous Media Approach for CANDU-6 Moderator Analysis

Analysis of the Cooling Design in Electrical Transformer

Theoretical and Experimental Studies on Transient Heat Transfer for Forced Convection Flow of Helium Gas over a Horizontal Cylinder

Journal of NUCLEAR SCIENCE and TECHNOLOGY, Vol. 41, No. 7, p (July 2004)

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 6

CFD Simulation of Sodium Boiling in Heated Pipe using RPI Model

THERMAL HYDRAULIC REACTOR CORE CALCULATIONS BASED ON COUPLING THE CFD CODE ANSYS CFX WITH THE 3D NEUTRON KINETIC CORE MODEL DYN3D

Visualization of Distribution of Shear Stress due to Water Vortex Flow with SSLCC

Safety Analysis of Loss of Flow Transients in a Typical Research Reactor by RELAP5/MOD3.3

Uncertainty Quantification of EBR-II Loss of Heat Sink Simulations with SAS4A/SASSYS-1 and DAKOTA

CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel

EXPERIMENTAL AND NUMERICAL STUDIES OF A SPIRAL PLATE HEAT EXCHANGER

Numerical Simulation on Flow and Heat Transfer in Oscillating Heat Pipes

Numerical Investigation of Sodium-Water Reaction Phenomenon in a Tube Bundle Configuration

Mechanical Engineering Journal

THERMAL HYDRAULIC ANALYSIS IN REACTOR VESSEL INTERNALS CONSIDERING IRRADIATION HEAT

A COUPLED CFD-FEM ANALYSIS ON THE SAFETY INJECTION PIPING SUBJECTED TO THERMAL STRATIFICATION

Analysis of Heat Transfer Enhancement in Spiral Plate Heat Exchanger

NTHAS7: The Seventh Korea-Japan Symposium on Nuclear Thermal Hydraulics and Safety Chuncheon, Korea, November 14-17, 2010.

CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE

ANALYSIS OF METASTABLE REGIMES IN A PARALLEL CHANNEL SINGLE PHASE NATURAL CIRCULATION SYSTEM WITH RELAP5/ MOD 3.2

The Dynamical Loading of the WWER440/V213 Reactor Pressure Vessel Internals during LOCA Accident in Hot and Cold Leg of the Primary Circuit

Computation of turbulent natural convection with buoyancy corrected second moment closure models

SIMULATION OF FLOW IN A RADIAL FLOW FIXED BED REACTOR (RFBR)

STUDY ON BUBBLE BEHAVIOR OF INERT GASES AT ENTRANCE NOZZLE IN SODIUM-COOLED FAST REACTOR

Available online at ScienceDirect. Procedia Engineering 157 (2016 )

Experimental and Numerical Investigation of Two- Phase Flow through Enlarging Singularity

CONSTRUCTION OF HIERACHICAL MODEL BASED ON FACTOR ANALYSIS FOR WASTAGE RATE PREDICTION IN SODIUM-WATER REACTION

Unsteady RANS and LES Analyses of Hooper s Hydraulics Experiment in a Tight Lattice Bare Rod-bundle

Research Article Analysis of NEA-NSC PWR Uncontrolled Control Rod Withdrawal at Zero Power Benchmark Cases with NODAL3 Code

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

Numerical simulation of high pressure gas quenching of H13 steel

Computational fluid dynamics study of flow depth in an open Venturi channel for Newtonian fluid

Simplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis

VHTR Thermal Fluids: Issues and Phenomena

CFX SIMULATION OF A HORIZONTAL HEATER RODS TEST

Numerical Simulation Of Fluid Flow And Heat Transfer Of Supercritical Pressure

Computational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor

THE EFFECT OF TWO PHASE (AIR-WATER) FLOW CHARACTERISTICS ON MOMENTUM FLUX DUE TO FLOW TURNING ELEMENTS AT ATMOSPHERIC CONDITIONS

Development and Applicability Evaluation of Frequency Response Function of Structures to Fluctuations of Thermal Stratification

Transport phenomenon in a jet type mold cooling pipe

Measuring of velocity and temperature field in a model of reactor vessel downcomer and cold-legs inlets

Vertical Mantle Heat Exchangers for Solar Water Heaters

ULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor

ANALYSIS OF THE OECD MSLB BENCHMARK WITH THE COUPLED NEUTRONIC AND THERMAL-HYDRAULICS CODE RELAP5/PARCS

Studies on flow through and around a porous permeable sphere: II. Heat Transfer

THE USE OF PB-BI EUTECTIC AS THE COOLANT OF AN ACCELERATOR DRIVEN SYSTEM. Joint research Centre of the European Commission Ispra, Italy.

Flow analysis in centrifugal compressor vaneless diffusers

Mechanical Engineering Journal

Keywords: Spiral plate heat exchanger, Heat transfer, Nusselt number

Increase Productivity Using CFD Analysis

Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions

FLOW CHARACTERIZATION WITHIN A SPHERE-PACKED BED USING PIV MEASUREMENT

JOYO MK-III Performance Test at Low Power and Its Analysis

APPLICATION OF THE COUPLED THREE DIMENSIONAL THERMAL- HYDRAULICS AND NEUTRON KINETICS MODELS TO PWR STEAM LINE BREAK ANALYSIS

HIGH TEMPERATURE THERMAL HYDRAULICS MODELING

Available online at ScienceDirect. Energy Procedia 71 (2015 )

Study on the flow and thermal characteristics of a heat storage system

SUB-CHAPTER D.1. SUMMARY DESCRIPTION

Tutorial for the supercritical pressure pipe with STAR-CCM+

Status and Future Challenges of CFD for Liquid Metal Cooled Reactors

Sensitivity Analyses of the Peach Bottom Turbine Trip 2 Experiment

The influence of a magnetic field on turbulent heat transfer of a high Prandtl number fluid

ANALYSIS OF THE OECD PEACH BOTTOM TURBINE TRIP 2 TRANSIENT BENCHMARK WITH THE COUPLED NEUTRONIC AND THERMAL-HYDRAULICS CODE TRAC-M/PARCS

International Journal of Research in Advent Technology, Vol.6, No.11, November 2018 E-ISSN: Available online at

CFD ANALYSIS OF A THERMAL MIXING IN A SUBCOOLED WATER POOL UNDER A HIGH STEAM MASS FLUX

CFD-Modeling of Turbulent Flows in a 3x3 Rod Bundle and Comparison to Experiments

Application of Computational Fluid Dynamics to the Flow Mixing and Heat Transfer in Rod Bundle

Steady-State and Transient Neutronic and Thermal-hydraulic Analysis of ETDR using the FAST code system

Investigation of Flow Profile in Open Channels using CFD

EasyChair Preprint. Numerical Simulation of Fluid Flow and Heat Transfer of the Supercritical Water in Different Fuel Rod Channels

Coupled CFD-STH analysis of liquid metal flows: STAR-CCM+ - RELAP5

Condition monitoring of steam generator by estimating the overall heat transfer coefficient

Thermal Hydraulic System Codes Performance in Simulating Buoyancy Flow Mixing Experiment in ROCOM Test Facility

CONDENSATION HEAT TRANSFER COEFFICIENT CORRELATION BASED ON SLIP RATIO MODEL IN A HORIZONTAL HEAT EXCHANGER

Design optimization of first wall and breeder unit module size for the Indian HCCB blanket module

Fluid structure interaction dynamic analysis of a mixed-flow waterjet pump

ARIES-AT Blanket and Divertor Design (The Final Stretch)

Experiences of TRAC-P code at INS/NUPEC

Heat Transfer Experiment and Simulation of the Verification Facility for High Power Rotating Tritium Target System

Transcription:

Nuclear Engineering and Design 258 (2013) 226 234 Contents lists available at SciVerse ScienceDirect Nuclear Engineering and Design j ourna l ho me pag e: www.elsevier.com/locate/nucengdes Numerical study on thermal stratification phenomena in upper plenum of LMFBR MONJU Makoto Shibahara, Takashi Takata, Akira Yamaguchi Graduate School of Engineering, Osaka University, 2-1 Yamada-Oka, Suita, Osaka, Japan h i g h l i g h t s Three-dimensional analysis of thermal stratification in the upper plenum of MONJU is conducted using a commercial CFD code, FLUENT. A comparison between the numerical simulation and the trip test of MONJU is made. This study shows that the interface of thermal stratification is influenced by the flow pattern in the upper plenum of MONJU. a r t i c l e i n f o Article history: Received 29 February 2012 Received in revised form 28 December 2012 Accepted 1 February 2013 a b s t r a c t The three-dimensional analysis of thermal stratification in the upper plenum of MONJU is conducted using the commercial CFD code, FLUENT ver.12.1. Since the temperature gradient near the thermal stratification interface would cause thermal stress in the reactor components, it is important to understand the characteristics of thermal stratification for evaluating the structural integrity in the upper plenum. As the result of numerical analysis, it is understood that the interface of thermal stratification is influenced by the flow pattern in the upper plenum. After the jet from the core outlet impinges on the upper core structure, the hot fluid flows obliquely upward to the inner barrel under the 40% electric output and flow condition. On the other hand, the jet from the core outlet flows to the lower part of the upper plenum, and then cold fluid flows through the s after the turbine trip. Hence, the flow structure has changed from the initial condition as the flow rate and temperature of the core outlet decrease due to the turbine trip. It is considered that the flow path of the s has been maintained for a long duration since the thermal stratification changes the distribution of buoyancy forces. 2013 Elsevier B.V. All rights reserved. 1. Introduction The knowledge of thermal stratification phenomena may be important for safety assessment in the upper plenum of some liquid metal fast breeder reactors (LMFBRs). The thermal stratification interface in the upper plenum of MONJU was observed due to the flow coasting down after the turbine trip (Doi et al., 1996). Since the temperature gradient near the thermal stratification interface could cause thermal stress in the reactor components, it is important to understand the characteristics of thermal stratification to evaluate the impact on the structural integrity in the upper plenum. In order to predict the temperature distribution after the turbine trip, many numerical analyses have been reported under various conditions. Ieda et al. (1990) conducted the analysis of the thermal stratification for water and sodium with 1/10 scaled models of the upper plenum using a multi-dimensional thermal hydraulic code, AQUA. And then, Doi and Muramatsu (1996) performed the Corresponding author. Tel.: +81 6 6879 7891. E-mail address: takata t@see.eng.osaka-u.ac.jp (T. Takata). numerical analysis of the thermal stratification for the turbine trip test was conducted using AQUA. Since it was complicated to represent the reactor components such as the s of the inner barrel and the upper core structure (UCS) accurately, the empirical correlations for the pressure drop were given in the numerical analysis. Kimura et al. also carried out the numerical analysis using AQUA on the thermal stratification of plant dynamics test loop integrated with direct heat exchanger (PLANDTL-DHX) (Kamide et al., 1998; Nishimura et al., 2000) which has an upper plenum of 2 m diameter (Kimura et al., 1999). They evaluated the thermal stratification interface and clarified the effect of Richardson number and Reynolds number. Moreover, Ohno et al. (2007) conducted the numerical analysis for thermal stratification in the 1/10 scale model using FLUENT (ANSYS, 2012) and STAR-CD (STAR, 2012) to validate the commercial codes. They clarified the influence of mesh arrangement, discretized scheme, and various turbulence models on the thermal stratification. Although the axial temperature distribution and the thermal stratification interface were evaluated in the 1/10 scale model, the full scale model including the reactor components was not simulated. As a benchmark problem of the International Atomic Energy Agency (IAEA), Yoshikawa and 0029-5493/$ see front matter 2013 Elsevier B.V. All rights reserved. http://dx.doi.org/10.1016/j.nucengdes.2013.02.007

M. Shibahara et al. / Nuclear Engineering and Design 258 (2013) 226 234 227 Fig. 1. Schematic of the upper plenum. Minami (2008) summarized detailed engineering data of the turbine trip test under the 40% rated power operation of MONJU. Honda et al. (2010) performed the numerical analysis of the upper plenum of MONJU under the 40% rated power operational condition using the commercial code, FrontFlow/Red (Tubokura et al., 2009), which was developed at the University of Tokyo under the national project Revolutionary Simulation Software (RSS). As a parametric study on the thermal stratification in the upper plenum of MONJU, Sakamoto et al. (2010) carried out the numerical analysis using FLUENT code in order to investigate the influence of mesh arrangement on the thermal stratification. It was clarified that thermal stratification interface was affected by the mesh arrangement of the s. In this research, we focused on the flow dynamic structure in the upper plenum of MONJU after the turbine trip. The three-dimensional analysis of thermal stratification in the upper plenum of MONJU is conducted using the commercial CFD code, FLUENT ver. 12.1.2 (ANSYS, 2012). The purposes of this study are: (1) to clarify the flow pattern of cold sodium which flows to the upper plenum after the turbine trip and (2) to compare the axial temperature distribution and the temperature gradient with those of experimental data (Doi et al., 1996). Fig. 2. UCS overall view. flow guide tubes and honeycomb structure in Fig. 2 is assumed as the porous media body, which is applied with the porosity of 0.5. Fig. 3 shows the top view of the core subassemblies of MONJU (Yoshikawa and Minami, 2008). The entire core channels are divided into eighteen regions. Table 1 shows the flow channel ID at the core outlet. The core subassemblies of MONJU consist of inner diver core, outer diver core, neutron source, blanket, neutron 2. Outline of turbine trip test of MONJU In the turbine trip test of MONJU, the 40% rated power operational condition was set as an initial condition and the turbine trip was conducted manually. The steady-state condition means the 40% electric output and flow condition. When the turbine trip test started, the circulation pump stopped with a flow coasting down characteristic and a pony motor was activated. This operation after the turbine trip was defined as the transient condition. As the result of the turbine trip test, the flow rate of the primary coolant decreased with the flow coasting down characteristic and was kept to be approximately 10% of the rated operation under the transient condition. The temperature distribution in the axial direction of the upper plenum of MONJU was measured with thermocouple (T/C) plug installed in the upper plenum. The upper plenum consists of upper instrument structure (UIS), T/C plug, core barrel, three hot leg pipes (H/L), fuel handling machine (FHM) and inner barrel with the s as shown in Fig. 1. 3. Numerical method 3.1. Analytical model Fig. 1 shows the analytical model of the upper plenum. In the numerical analysis, the upper core structure (UCS) such as fingers, Fig. 3. Top view of the core subassemblies.

228 M. Shibahara et al. / Nuclear Engineering and Design 258 (2013) 226 234 Table 1 Flow channel ID at the core outlet. Core region Flow channel ID Inner core 1 2 3 4 5 Outer core 6 7 8 Blanket 9 10 11 Neutron shielding 12 13 FC Fig. 4. Cross-sectional view of the mesh arrangement. shielding, and control rods such as coarse control rods, fine control rods and backup control rods. Fig. 4 shows the cross-sectional view of the mesh arrangement with C-loop H/L for the analytical model. In the numerical analysis, the mesh number of approximately 5.11 million was determined based on the previous study (Sakamoto et al., 2010). To clarify the flow distribution in the s and overflow of the inner barrel, the fine mesh is arranged in the s and top of the inner barrel as shown in Fig. 5. Since it was complicated to model the flow holes on the inner barrel accurately, the pressure drop correlations were applied in the previous analysis (Doi and Muramatsu, 1996) as mentioned in Section 1. Since the thermal stratified layer would affect on the coefficient of pressure drop, the s have been represented by as built-geometry in this numerical analysis. 3.2. Boundary condition and numerical method The flow rate and the sodium temperature for the inlet boundary condition are consistent with the IAEA benchmark condition (Yoshikawa and Minami, 2008). Fig. 6 shows the typical flow rate of channels at the inlet boundary condition. As mentioned in Section 2, the flow rate from the core outlet decreases with time and kept approximately 10% of the rated operation. The turbulence intensity, I, is estimated to be 5% at the inlet boundary conditions, and Flow rate [kg/s] Nutron source N Coarse CR rods C Fine CR rods F Backup CR rods B 400 350 300 250 200 150 100 50 0 Channel 1 (Inner core) Channel 6 (Outer core) Channel 9 (Blanket) Channel 12 (Neutron Shielding) 0 100 200 300 400 500 600 Time [s] Fig. 6. Flow rate of channels after the turbine trip. the turbulent kinematic energy, k, is calculated as k = 3/2(U avg I) 2, where, U avg is the mean outlet velocity for each channel. And also, Fig. 7 shows the typical sodium outlet temperature of each channel (Yoshikawa and Minami, 2008). As shown in Fig. 7, the sodium Fig. 5. Mesh arrangement for (a) and (b) top of inner barrel.

M. Shibahara et al. / Nuclear Engineering and Design 258 (2013) 226 234 229 Temperature [ o C] 550 500 450 Channel 1 (Inner core) Channel 6 (Outer core) Channel 9 (Blanket) Channel 12 (Neutron Shielding) of gravity is considered in the numerical analysis. The basic conservation equations are discretized with finite volume method (FVM) in FLUENT code (ANSYS, 2012) The QUICK scheme is applied for convective terms so as to reduce numerical diffusion (Maekawa, 1990), and the time marching term is calculated by the second order Euler implicit method. The time step is 2 s based on a previous study (Sakamoto et al., 2010). The analysis time is ranged from 0 s to 600 s. Also, the SIMPLE method and standard k ε model are applied in the numerical analysis. 400 4. Result and discussion 350 0 100 200 300 400 500 600 Time [s] Fig. 7. Sodium temperature of channels after the turbine trip. temperature decreases with time after the trip. The outlets of three H/L pipes are bounded by satisfying a constant pressure condition. The walls of reactor components are adiabatic in this numerical analysis. The density, isobaric specific heat, and kinematic viscosity of sodium are expressed by quadratic or cubic polynomial of sodium temperature (Sodium Technology Education Committee, 2005). The thermal conductivity of sodium is defined with a straight-line approximation polynomial, which depends on sodium temperature. Since the buoyancy force would affect stratified flow, the effect 4.1. Steady-state analysis Fig. 8 shows the velocity (a) and temperature (b) profiles under the steady-state condition. It can be seen from Fig. 8(a), the jet from the core outlet ascends through the porous media zone which is assumed as the UCS, and it impinges on UIS. Then, the hot fluid flows obliquely upward to the inner barrel. After impinging on the inner barrel, the flow is branched off into two directions; upward flow and downward flow. The upward flow along the inner barrel overflows the top of the inner barrel, and then it descends through the annular gap between the inner barrel and the reactor vessel wall. In the numerical analysis, the ratio of flow rate of s is evaluated to be 10% compared with the overflow rate of the inner barrel. Hence, the oblique flow comes to govern under the steady-state condition. On the other hand, as shown in Fig. 8(b), the temperature profile is not uniform at the lower side of the upper plenum. Fig. 8. Velocity and temperature profiles under the steady state condition. Fig. 9. Velocity profiles for the depth of (a) 5300 mm, (c) 4400 mm, and (d) 1000 mm from the liquid surface.

230 M. Shibahara et al. / Nuclear Engineering and Design 258 (2013) 226 234 Fig. 10. Temperature profiles for the depth of (a) 5300 mm, (c) 4400 mm, and (d) 1000 mm from the liquid surface. This is because the coolant temperature of blanket fuel subassembly is lower than that of the inner core fuel subassembly, so that the lower temperature fluid descends along the core barrel and accumulates between the support plate and the lower s. Fig. 9 shows the velocity profiles at the depth of (a) 5300 mm, (b) 4400 mm, and (c) 1000 mm from the liquid surface. These correspond to the vertical cross-sectional view of figure in Fig. 8(a). The circumference direction of velocity profile is almost uniform at each depth from the liquid surface even though the reactor components such as the FHM and the T/C plug are installed, while there are circumferential variations in the annular flow regions due to the outlet pipes as shown in Fig. 7(b). Fig. 10 shows the circumference direction of the temperature profiles at the depth of (a) 5300 mm, (b) 4400 mm and (c) 1000 mm from the liquid surface. The circumference direction of temperature profiles is also uniform at each depth from the liquid surface. Hence, the asymmetric structure of the upper plenum has relatively little effect on the flow structure and temperature distribution under the steady-state condition. 4.2. Transient analysis Fig. 11 shows the flow pattern during transient operation ranged from 60 s to 600 s. The flow pattern has changed from the Fig. 11. Flow pattern during the transient operation ranged from 60 s to 600 s.

M. Shibahara et al. / Nuclear Engineering and Design 258 (2013) 226 234 231 Fig. 12. Temperature profiles during the transient operation ranged from 60 s to 600 s. steady-state condition due to the flow coasting down. As shown in Fig. 11(a), the jet is spread to the upper plenum horizontally after it flows between UIS and the core outlet. Then, the jet flows to the support plate as shown in Fig. 11(b). It is considered that the inertial force of jet becomes smaller than the plugging effect, which means that the hot sodium acts as the plug due to the buoyancy force. After the coolant flows on the support plate, it flows to the lower s. For the transient operation ranged from 240 s to 600 s as shown in Fig. 11(c) (f), the flow pattern has not changed since the plugging effect is greater than the inertial force. Fig. 12 shows the vertical cross-sectional view of temperature profiles during the transient operation ranged from 60 s to 600 s. It can be seen from Fig. 12(a) that the cold sodium begins to penetrate the upper plenum since the core outlet temperature has reduced rapidly due to scram. The temperature variations result in density changes that affect the flow conditions. As shown in Fig. 12(b), the thermal stratification interface is formed near the upper flow holes as cold sodium fills with lower part of upper plenum. It seems that the thermal stratification is characterized by the flow distribution. Once the thermal stratification is established, it would last for a long time. Moreover, the axial temperature gradient near the interface of thermal stratification may be of significant concern as mentioned in Section 1. Fig. 12(c) (f) shows that the interface of thermal stratification ascends due to bottom-up flow from the core Fig. 13. Temperature profiles in horizontal cross section for the elapsed time of (a) 60 s, (b) 120 s, and (c)240 s at 5300 mm.

232 M. Shibahara et al. / Nuclear Engineering and Design 258 (2013) 226 234 100 90 80 70 Annular gap Upper s Lower s Flow ratio [%] 60 50 40 30 20 10 0 0 100 200 300 400 500 600 Time [s] Fig. 14. Temperature distribution under the transient condition. Fig. 15. Flow ratio of annular gap and s in the transient analysis. thermal stratification interface, the non-dimensional temperature T * is defined as following equation (Ohno et al., 2007): outlet since the flow rate through the s is limited. In other words, the position of thermal stratification interface moves slowly under the transient condition. As the result of numerical analysis of flow dynamic structure and thermal stratification in the upper plenum of MONJU, it is understood that the interface of thermal stratification is influenced by the flow pattern after the turbine trip. Fig. 13 shows the temperature profiles in horizontal crosssection at the depth of 5300 mm from the liquid surface, where the lower s exist, for the elapsed time of (a) 60 s, (b) 120 s, and (c) 240 s. These correspond to the vertical cross-sectional view of figures in Fig. 12(a) (c), respectively. As shown in Fig. 13(a), a uniform temperature distribution exists at 60 s. When the time is longer than 60 s, the temperature difference appears between the annular gap and the inside of the inner barrel since the coolant flows in the upper plenum as shown in Fig. 13(b). There are temperature variations inside the inner barrel at 120 s, however, since the T/C plug is far from the modalities of the temperature change, it seems that they have a little influence on the long-time temperature profiles at the T/C plug. As the thermal stratification interface rises with time, the coolant inside the inner barrel flows through the s as shown in Fig. 13(c). And then, after that the coolant in the inner barrel and coolant in the annular gap mix and flow toward H/L pipes. In short, the temperature distribution changes along with the radial direction, in contrast, it is uniform on the circumference direction under the transient condition. Fig. 14 shows the comparison of the axial temperature distributions between this analysis and the experimental data at the T/C plug (see Fig. 1). The symbols show the experimental data for the turbine trip under the 40% rated power operational condition while the solid lines show the numerical result. The numerical result shows good agreement with the experimental data except for the depth lower than 6000 mm under the initial condition of the turbine trip (0 s). After the turbine trip, the calculated temperature accords with the experimental data at the s, however, it is lower than that of the experiment above the upper s after 240 s. Therefore, we investigated the flow ratio of s to the annulus gap as shown in Fig. 15. As a result, it is found that the fraction of the flow rate through the s increases after the trip, while the flow rate of the annular gap decreases with time. However, since the flow ratio of annular gap remains about 18% in Fig. 15, the fluctuation of annular flow has a possibility to delay the raise of the thermal stratification layer. Moreover, it remains possible that the flow rate in s is lower in the numerical analysis than in the experiment. Fig. 16 shows the position of the thermal stratification interface after the turbine trip. In order to evaluate the position of the T = T T c (1) T h T c where T c is minimum temperature, and T h is maximum temperature of liquid surface at the trip for the T/C plug. When the value of T * become 0.5 both in the numerical analysis and the experiment, the position will be defined as the interface of thermal stratification. As shown in Fig. 16, it is understood that the interface of the thermal stratification raises as the coolant flows into the upper plenum. The trend of numerical values is consistent with the experimental data at the early time ranged from the turbine trip (0 s) to 240 s, while both the numerical result and the experiment curve after 240 s. The slope of numerical result is higher that of experiment, however, the numerical result shows similar curvature to the experiment. As shown in Fig. 11, there is no obvious change of the flow pattern after 240 s in the numerical result, while the interface of the thermal stratification would be not liner during 240 600 s. Since the flow rate through the s increase with time as shown in Fig. 15, the thermal stratification layer would be affected by the bypass flow through the s. Fig. 17 shows the comparison of the axial temperature gradients between the numerical result and experimental data for the elapsed time of (a) 120 s, (b) 240 s, and (c) 600 s. Each temperature gradient (dt/dz) is determined from the T/C temperature distribution as shown in Fig. 14. As mentioned in Section 1, the evaluation of temperature gradient is important for determining the structural integrity of the reactor component such as the inner barrel. In the experiment, the temperature gradient is obtained from the [mm] Depth from liquid surface 0-1000 -2000-3000 -4000-5000 -6000 Numerical result Experimental data -7000 0 50 100 150 200 250 300 350 400 450 500 550 600 650 Time [s] Fig. 16. Thermal stratification interface after the turbine trip. Liquid surface Upper Lower

M. Shibahara et al. / Nuclear Engineering and Design 258 (2013) 226 234 233 (a) 0.00 Time 120s Liquid surface 0.00 Time 240s Liquid surface Depth from liquid surface [m] -1.00-2.00-3.00-4.00-5.00 Experimental data Numerical result Top of inner barrel Upper Lower Depth from liquid surface [m] -1.00-2.00-3.00-4.00-5.00 Experimental data Numerical result Top of inner barrel Upper Lower -6.00-200 -150-100 -50 0 50 100 150 200 Temperature gradient [ ºC/m] -6.00-200 -150-100 -50 0 50 100 150 200 Temperature gradient [ºC/m] (c) 0.00 Time 600s Liquid surface Depth from liquid surface [m] -1.00-2.00-3.00-4.00-5.00 Experimental data Numerical result Top of inner barrel Upper Lower -6.00-200 -150-100 -50 0 50 100 150 200 Temperature Gradient [ºC/m] Fig. 17. Temperature gradients for the elapsed time of (a) 120 s, (b) 240 s, and (c) 600 s. readings from thermocouples and the arrangement. The symbols show the experimental data and numerical result, at each elevation defined by the arrangement of the thermocouple in the experiment. As shown in Fig. 17(a), the maximum temperature gradient in the numerical result is about 80 C lower than that of the experimental data. After 120 s, it appears that the maximum temperature gradient decreases due to the mixing effect on the thermal stratification. The numerical result is approximately consistent with the experimental data. Moreover, the maximum temperature gradient of the numerical simulation shows good agreement with that of experiment. In Fig. 17(c), the elevation of the maximum temperature gradient at 600 s is higher than that of experiment, although the maximum temperature gradient approximately accords with that of the experiment. As mentioned in Fig. 14, the thermal stratification interface is overestimated after 240 s, therefore, the elevation of the maximum temperature gradient becomes higher than that of experiment. And also, the spike of the numerical result was caused by a discrepancy in the temperature distribution at the top of the inner barrel as shown in Fig. 14. Moreover, since the interface of thermal stratification has risen at 600 s in the numerical result, the temperature change becomes lower than that of the experiment at the upper s. 5. Conclusions The three-dimensional analysis of thermal stratification in the upper plenum of MONJU is conducted using the commercial CFD code, FLUENT ver. 12.1.2 to analyze the flow dynamic structure on the thermal stratification phenomena. As a result, thermal stratification interface is influenced by the flow pattern in the upper plenum of MONJU. The flow structure of the upper plenum has changed from obliquely flow beyond the top of inner barrel to downward flow through the s after the turbine trip. The numerical result is consistent with that of the experimental data until 240 s, while the thermal stratification interface after 240 s ascends faster than that of experimental data. The rise of the thermal stratification layer would be affected by the flow ratio of the s to the annular gap. Therefore, it may be important to clarify a fluctuation of the flow ratio by the related boundary conditions of the s and annular flow. Acknowledgements Present study is the result of 2010 R&D core program for practical realization of fast breeder reactor by using MONJU entrusted to University of Fukui by the Ministry of Education, Culture, Sports, Science and Technology of Japan (MEXT). References ANSYS FLUENT, ANSYS, Inc., http://www.ansys.com/products/simulation+ Technology/Fluid+Dynamics/. Fluid + Dynamics + Products/ANSYS + Fluent (accessed 20-12-2012). Doi, Y., Muramatsu, T., 1996. Investigation of analytical methods in thermal stratification analysis evaluation of flow rates through s for normal and scram conditions of 40% power operation with AQUA code. Power reactor and Nuclear Fuel Development Corporation, PNC TN9410 97-083 (in Japanese). Doi, Y., Muramatsu, T., Kunoki, K., 1996. Thermal stratification tests in MONJU upper plenum (I) temperature distribution under normal and scram conditions with 40% power operation. Power Reactor and Nuclear Fuel Development Corporation, PNC TN9410 96-117 (in Japanese). Honda, K., Ohira, H., Sotsu, M., Yhoshikawa, S., 2010. Thermal-hydraulic analysis of MONJU upper plenum under 40% rated power operational condition. In: Proc. Int. Topical Meeting on Nuclear Thermal-Hydraulics, Operation and Safety (NUTHOS-8). Ieda, Y., Maekawa, I., Muramatsu, T., Nakanishi, S., 1990. Experimental and analytical studies of the thermal stratification phenomenon in the outlet plenum of fast breeder reactor. Nucl. Eng. Des. 120, 403 414.

234 M. Shibahara et al. / Nuclear Engineering and Design 258 (2013) 226 234 Kamide, H., Hayashi, K., Toda, S., 1998. An experimental study of inter-subassembly heat transfer during natural circulation decay heat removal in fast breeder reactors. Nucl. Eng. Des. 183, 97 106. Kimura, N., Nishimura, M., Hayashi, K., Kamide, H., 1999. An investigation of thermal stratification phenomena in fast reactors sodium experiments and analyses. In: Proc. Int. Conf. on Nucl. Eng. (ICONE7), ICONE7-7107. Maekawa, I., 1990. Numerical diffusion in single-phase multi-dimensional thermalhydraulic analysis. Nucl. Eng. Des. 120, 323 339. Nishimura, M., Kamide, H., Hayashi, K., Momoi, K., 2000. Transient experiments on fast reactor core thermal-hydraulics and its numerical analysis: intersubassembly heat transfer and inter-wrapper flow under natural circulation conditions. Nucl. Eng. Des. 200, 157 175. Ohno, S., Ohki, H., Sugahara, A., Ohshima, H., 2007. Validation of a computational simulation method for evaluating thermal stratification in the reactor vessel upper plenum of fast reactors. J. Nucl. Sci. Technol. 48 (2), 205 214. Sakamoto, T., Shibahara, M., Takata, T., Yamaguchi, A., 2010. Numerical study of three dimensional thermal hydraulics effect on thermal stratification phenomena in upper plenum of MONJU. In: Proc. Korea-Japan Symposium on Nuclear Thermal Hydraulics and Safety (NTHAS7), N7P0032. Sodium Technology Education Committee ed., 2005. Sodium Technology Handbook, Japan Nuclear Cycle Development Institute. JNC-TN9410 2005-011 (in Japanese). STAR-CD, CD-adapco, Inc., http://www.cd-adapco.com/products/star cd/ (accessed 20-12-2012). Tubokura, M., Kobayashi, K., Nakashima, T., Nouzawa, T., Nakamura, T., Zhang, H., Onishi, K., Oshima, N., 2009. Computational visualization of unsteady flow around vehicles using high performance computing. Comput. Fluids 38, 981 990. Yoshikawa, S., Minami, M., 2008. Data description for coordinated research project on benchmark analyses of sodium natural convection in the upper plenum of the MONJU reactor vessel under supervisory of technical working group on fast reactors. International Atomic Energy Agency, JAEA-Data/Code 2008-024, Japan Atomic energy Agency.