Natural Circulation Level Optimization and the Effect during ULOF Accident in the SPINNOR Reactor Ade Gafar Abdullah 1,2,*, Zaki Su ud 2, Rizal Kurniadi 2, Neny Kurniasih 2, Yanti Yulianti 2,3 1 Electrical Engineering Dept. Universitas Pendidikan Indonesia (UPI), Jl. Dr. Setiabudhi 229 Bandung, Indonesia 2 Physics Dept. Institut Teknologi Bandung (ITB), Jl. Ganesha 10 Bandung, Indonesia 3 Physics Dept. Universitas Lampung (UNILA), Jl. Sumantri Brojonegoro 1 Bandar Lampung, Indonesia *e-mail :ade_gaffar@students.itb.ac.id Abstract. Natural circulation level optimization and the effect during loss of flow accident in the 250 MWt MOX fuelled small Pb-Bi Cooled non-refuelling nuclear reactors (SPINNOR) have been performed. The simulation begin with steady state calculation of neutron flux, power distribution and temperature distribution accros the core, hot and cool pool, and also steam generator. Then the accident started by the loss of pumping power. The sequence is then the transient flow calculation across the core, the core temperature distribution, hot pool and cool pool fluid dynamic calculation and temperature change simulation, and the transient flow and temperature calculation across the steam generator. Then the reactivity feedback calculation is conducted, followed by kinetic calculation, and then the sequence repeated. The results show that the SPINNOR reactor has inherent safety capability against this accident. Keywords: natural circulation, ULOF accident, SPINNOR reactor PACS: 28.20.Gd INTRODUCTION Safety aspect is the most important part of designing a nuclear reactor. After the Chernobyl Accident and Tree Mile Island II, there was a change of paradigm in the design and reactor safety. Fourthgeneration nuclear reactors must have inherent safety criteria. The SPINNOR (Small Power Reactor, Indonesia, No On-Site Refuelling) are concepts of small leadbismuth cooled nuclear power reactors with fast neutron spectrum that could be operated for more than 15 years without on-site refuelling. They are based on the concept of a long-life core reactor developed in Indonesia since early 1990 in collaboration with the Research Laboratory for Nuclear Reactors of The Tokyo Institute of Technology [2]. The reactor has been used in relatively isolated areas, small island and operated to the end of its life without refueling or fuel shuffling. Some important characteristics required are easy operation, easy maintenance, transportability, inherent/passive safety and nuclear proliferation resistance [2]. During ULOF accident, natural circulation level, coolant-fuel temperature difference, and reactivity feedback play important role to reach inherent safety capabilities. The paper discuss the optimization natural circulation and the effect during ULOF accident in the SPINNOR reactor. COMPUTATIONAL METHOD The model and computational aspect of accident analysis performed in this research are mainly based on the method described in [8].The general calculations are neutronic, thermal hydraulic, and feedback calculations. Calculation of distribution fluxs neutron If the reactor is evaluated in a steady state and the material core region is assumed homogen, the diffusion equation can is writen as [6] :
x, Dgg ( r) rg g ( r) k G ( r) ( r) 1 1 1 1 1 fg g 1 sg g g g 1 g 1 G g eff (1) The neutronic calculation yields neutron flux distribution and multiplication factor (K eff ). In adiabatic model, the amplitude function p( is satisfied by the following equations [8]: dp( ( p( dt 6 i1 ic i (2) dci ( i p( ici ( (3) dt Thermal Hydraulic Calculation Calculation of thermal hydraulic covers calculation on distribution of temperature in all parts of reactor reactor core and steam generator ( SG) that is fuel temperature, coolant temperature and temperature cladding and gap, so pressure in reactor and all coolant circulation aspects in the reactor. CENTER OF COOLANT CHANNEL SPATIAL MESH FOR THERMAL HYDRAULIC CALCULATION FIGURE 1. The model for thermo-hydraulic calculation Figure 1. shows thermal hydraulic model in this calculations. Reactor is divided into concentric ring, where the cross-flow between two adjacent rings is assumed zero. For coolant, calculation uses mass and energy conservation equation: T T ''' c p wc p Q (4) t z Pressure drop calculation is solved by applying momentum conservation equation. G P G fg g t z z 2D 2 CENTER OF PELLET FUEL-CLADDING GAP (COOLANT) FUEL PELLET WIRE SPACER COOLANT e 2 CLADDING (5) Hot pool and Cool Pool During this analysis, average temperature in both hot pool and cool pool are used. For hot pool, the following equation used : T hp h hp A hp hp G psg C php G core C pcore T [( h A G C T old hp hp hp hp psg php hp G tc T ] /[( h A core pcore core hp hp hp G t G C ] psg core php h ( h A ( G G ) i 1 j hp hp hp hp psg core / ( A ) hp hp = Temperature of hot pool = Height of hot pool = Area of hot pool = Mass density of hot pool = Total mass flow rate in primary SG = Specific heat of hot pool = Total mass flow rate of core = Specific heat of core (6) (7) Similarly the temperature cool pool is solved using the following equations: old T [( h A G C T T cp h cp A cp cp cp cp cp cp psg pcp cp G tc T ] /[( h A core pcore core cp cp cp G t G C ] psg core pcp h ( h A ( G G ) i 1 j cp cp cp cp psg core / ( A ) cp cp = Temperature of cool pool = Height of cool pool = Area of cool pool = Mass density of cool pool Feedback Calculations (8) (9) ( in Eq. (3) is the summation of external reactivity and feedback reactivity. Feedback reactivity included Doppler, axial expansion, radial expansion, and void reactivity. Feedback reactivity is the effect of the core temperature T [6]. Can be express as follow the temperature coefficient of reactivity: (11) T Where is a negative value that is inserted in calculation since increasing temperature will cause decreasing. T is average temperature in side the core (fuel and coolan.
General Calculation Algorithm steam Reactor Vessel Subcooled water Hot pool Pump shielding reflector core Steam Generator (SG) inlet FIGURE 3. Overview of the SPINNOR reactor FIGURE 2. Calculational flow diagram For the first step, the steady state neutronic and the thermal-hydraulic calculation are performed, then the accident condition are taken.the ULOF accident is started by failure of the primary pumping system. Temperature in the core will be changed depend on how much flow rate has been inserted. The change core temperature will change average hot pool temperature and steam generator temperature. In this research, reactivity external of reactor is assumed constant, then heat will be transferred to water circulation in steam generator. RESULT AND DISCUSSION A schematic view of the system being analyzed is shown in Fig.3. The overall system includes a reactor core, hot pool, cool pool, steam generator and pump. The intermediate heat exchanger (IHX) is eliminated, and heat from primary coolant system is transferred directly to the steam-water loop through the steam generator. The coolant flows trough the core, removes heat generated in the core, and then flows up to the hot pool. From the hot pool, coolant flows into the steam generator, transferring the heat into the steam-water loop, and then goes down to the cool pool. From the cool pool the coolant is pumped back to the core. TABLE 1. Main reactor design parameter No Parameter Value 1 Reactor Power 250 MWth 2 Fuel UO 2 - PUO 2 3 4 5 6 7 8 Shielding Material Coolant Pin/Pitch diameter Pin diameter Cladding Thickness Total Flowrate B4C Pb-Bi 1.2 cm 1.0 cm 0.8 cm 8000 kg/s Table 1. shows the main reactor design parameter. In this section, calculation results including power reactor, reactivity feedback, max of fuel and coolant temperature and total flowrate, will be discussed. Fig.4 shows the core flow and steam generator primary side flow pattern change during ULOF accident. As reactors loose pumping power at t=0, the core flow rate decreases significantly and goes to the natural circulation level. Primary SG flow rate Core flow rate FIGURE 4. Flow change during ULOF accident
Natural circulation can drive the coolant circulation, and the flow rate remain about 20-30% of normal value at the latter part of the accident. There is a small oscillation of the total coolant flow rate in the primary side of the steam generator. which are mainly caused by the fact that flow rate in the core and primary side SG in the absence of primary pumping is mainly controlled by level and temperature of coolant in the hot and cool pool which in turn depend on the flowrate of the core and primary SG. Primary SG flow rate Core flow rate FIGURE 5. Primary SG flow rate change during ULOF accident with up chimney level Fig.5 shows the pattern of change core flow rate and primary SG flow rate after optimizing the pin pitch, pin diameter and height the chimney. This optimization contributes to increase natural circulation level. After optimization, the natural circulation is about 45% of the initial flow rate. FIGURE 7. Reactivity component change during ULOF accident Fig. 7 shows the reactivity component change during ULOF accident.the most important feedback in this accident is core radial expansion. At the beginning stage of the accident, the Doppler component is negative but then move to positive value near the asymptotic condition. At the asymptotic condition, a new equilibrium condition occurs at reduced power level (40% of nominal power) balancing the natural circulation level. Doppler reactivity feedback, core radial expansion reactivity feedback, coolant density feedback and axial expansion feedback, all give negative feedback contribution at the beginning of accident sequence to force the power down. Fuel temperature Cladding temperature Coolant temperature FIGURE 6. Hot spot temperature during ULOF accident Fig.6 show temperature of fuel, cladding and coolant. The maximum of fuel temperature, cladding temperature and coolant temperature, respectively o C, 603,32 o C, and 564,22 o C. Fuel, 791,61 cladding and coolant temperature are still below the maximum temperature limits. FIGURE 8. Power change during ULOF accident Fig. 8 shows the reactor power change pattern during accident. The power increase at the beginning of the accident but then decrease and goes towards a new stable condition which is about 45% of initial power. CONCLUSION The conclusions are summarized as follows, 1. In general, the reactor can survive the ULOF accident inherently.
2. The relatively high natural circulation component contributes important factor to keep the reactor survive ULOF accident. 3. To exchange from the natural circulation level, optimization can be done by adjusting the ratio between pin pitch and pin diameter and height of the chimney. 4. The result of analysis shows that the maximum temperatures of coolant, fuel and cladding are still below each temperature limit, and margin to coolant and fuel are large. ACKNOWLEDGMENTS The Authors would like to acknowledge Dr. Eng Zaki Su ud, Dr. Rizal Kurniadi, Dr. Neny Kurniasih and Ms. Yanti Yulianti comments and suggestions to this paper. REFERENCES 1. A. G. Abdullah, Z. Su ud, and Y. Yulianti, ULOF Accident Analysis for 300 MWth Pb-Bi Cooled MOX Fuelled SPINNOR Reactor: Proc. of the Intern. Conf. on Advances in Nuclear Science and Engineering 2009.] 2. Z. Su ud, Advanced SPINNORs Concept and The Prospect of Their Deployment in Remote Area: International Conference on Advances in Nuclear Science and Engineering, Bandung, Indonesia, 2007, pp. 199-207. 3. Z. Su ud and H. Sekimoto, Safety Aspect of Long- Life Small Safe Power Reactors: Ann. Nucl. Energy Vol. 22, No 11,1995, pp. 711-722. 4. Z. Su ud, Comparative Study on Safety Performance of Nitride Fueled Lead Bismuth Cooled Fast Reactor With Various Power Levels: Progress in Nuclear Energy, Vol. 32, No. 3/4, 1998, pp. 571-577. 5. Z. Su ud and H. Sekimoto, Accident Analysis of Lead-Bismuth Cooled Small Safe Long-Life Fast Reactor Using Metallic or Nitride Fuel: Nuclear Engineering and Design 162, 1996, pp. 205-222. 6. Y. Yulianti and Z. Su ud, Development of Three Dimensional Accident Analysis Code for Pb-Bi Cooled Tank type Fast Reactors: Proc. of the Intern. Conf. on Advances in Nuclear Science and Engineering 2007. 7. Waltar A. E. and Reynolds A. B. Fast Breeder Reactor: Pergamon Press, 1981 8. Duderstadt J.J. and Hamilton L. J. Nuclear Reactor Analysis: Joh Wiley and Sons,1076