HWR Moderator Sub-cooling Requirements to Demonstrate Back-up Capabilities of Moderator During Accidents

Similar documents
In-Vessel Retention Analysis for Pressurised Heavy Water Reactors (PHWR) under Severe Core Damage Accident (SCDA)

CFX SIMULATION OF A HORIZONTAL HEATER RODS TEST

In-Vessel Retention Analysis for PHWR Calandria under Severe Core Damage Accident Condition using ASTEC

EFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION

Reactivity Coefficients

Development of 3D Space Time Kinetics Model for Coupled Neutron Kinetics and Thermal hydraulics

4th International CANDU In-service Inspection Workshop and NDT in Canada 2012 Conference, 2012 June 18-21, Toronto, Ontario

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

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

VVER-1000 Reflooding Scenario Simulation with MELCOR Code in Comparison with MELCOR Simulation

ENGINEERING OF NUCLEAR REACTORS

Comparison of 2D and 3D FEM Models of Eddy Current Pressure Tube to Calandria Tube Gap Measurement

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

Chapter 8 - Technology of Accident Analysis

1. INTRODUCTION 2. EAEA EXISTING CAPABILITIES AND FACILITIES

Institute of Atomic Energy POLATOM OTWOCK-SWIERK POLAND. Irradiations of HEU targets in MARIA RR for Mo-99 production. G.

Chapter 10: Boiling and Condensation 1. Based on lecture by Yoav Peles, Mech. Aero. Nuc. Eng., RPI.

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

MECHANICAL VIBRATIONS OF CANDU FEEDER PIPES

Thermal-Hydraulic Design

CANDU Safety #3 - Nuclear Safety Characteristics Dr. V.G. Snell Director Safety & Licensing

SUB-CHAPTER D.1. SUMMARY DESCRIPTION

Tick the box next to those resources for which the Sun is also the source of energy.

ATLAS Facility Description Report

Comparison of Silicon Carbide and Zircaloy4 Cladding during LBLOCA

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

Safety Assessment on the Storage of Irradiated Graphite Waste Produced from the Decommissioning of KRR-2

Nuclear Research Facilities in Russia for innovative nuclear development N. Arkhangelskiy, ROSATOM IAEA Consultancy Meeting June 2013

Utilization of the Irradiation Holes in HANARO

NUCLEAR SAFETY AND RELIABILITY WEEK 8

Improvement of Critical Heat Flux Performance by Wire Spacer

SINGLE-PHASE AND TWO-PHASE NATURAL CONVECTION IN THE McMASTER NUCLEAR REACTOR

Characterisation of Deposits on Membrane Walls of Steam Generators by Heat Flux Density Measurement

Data analysis with uncertainty evaluation for the Ignalina NPP RBMK 1500 gas-gap closure evaluation

Available online at ScienceDirect. Energy Procedia 71 (2015 ) 52 61

Coupled Neutronics Thermalhydraulics LC)CA Analysis

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

INNER EDDY CURRENT TRANSDUCER WITH ROTATING MAGNETIC FIELD APPLICATION TO EXAMINATION OF FUEL CHANNELS OF PHWR NUCLEAR POWER PLANT

Topic 2: Heat Affects Matter in Different Ways

Heat Transfer with Phase Change

Application of System Codes to Void Fraction Prediction in Heated Vertical Subchannels

Structural Integrity and NDE Reliability II

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

Hybrid Low-Power Research Reactor with Separable Core Concept

Advanced Heavy Water Reactor. Amit Thakur Reactor Physics Design Division Bhabha Atomic Research Centre, INDIA

Prediction of Critical Heat Flux (CHF) for Vertical Round Tubes with Uniform Heat Flux in Medium Pressure Regime

Qualification of the system code AC² (submodule ATHLET) for the safety assessment of passive residual heat removal systems

If there is convective heat transfer from outer surface to fluid maintained at T W.

40P (2 x 60 x 60) = 2.5 x 10 6 (4200)(5) P = 1.82 x 10 5 W

M. R. Gartia Reactor Engineering Division Bhabha Atomic Research Centre, Trombay Mumbai P. K. Vijayan.

USE OF CFD TO PREDICT CRITICAL HEAT FLUX IN ROD BUNDLES

Analysis of Flow Reversal under Two-Phase. Natural Circulation in CANDU9 during Small Loss. of Coolant Accident with Loss of Class IV Power

DOE NNSA B&W Y-12, LLC Argonne National Lab University of Missouri INR Pitesti. IAEA Consultancy Meeting Vienna, August 24-27, 2010

PREDICTION OF THE RESPONSE OF THE CANADIAN SUPER CRITICAL WATER REACTOR TO POTENTIAL LOSS OF FORCED FLOW SCENARIOS

Loads on RPV Internals in a PWR due to Loss-of-Coolant Accident considering Fluid-Structure Interaction

Uncertainty Analysis on Containment Failure Frequency for a Japanese PWR Plant

Chapter 7. Design of Steam Generator

heat transfer process where a liquid undergoes a phase change into a vapor (gas)

Lecture 3: Light and Temperature

Module 4 : Deflection of Structures Lecture 4 : Strain Energy Method

VERIFICATION AND VALIDATION OF ONE DIMENSIONAL MODELS USED IN SUBCOOLED FLOW BOILING ANALYSIS

LOFT Experiment LP-02-6 Analysis by RELAP5/ MOD2 Code with Improved Minimum Film Boiling Temperature

CFD Simulation of Sodium Boiling in Heated Pipe using RPI Model

Coupling of thermal-mechanics and thermalhydraulics codes for the hot channel analysis of RIA events First steps in AEKI toward multiphysics

HEAT TRANSFER AT SUPERCRITICAL PRESSURES (SURVEY) 1

Investigation of CTF void fraction prediction by ENTEK BM experiment data

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

XA04C1727. Abstract. 1 Introduction

* Defining Temperature * Temperature is proportional to the kinetic energy of atoms and molecules. * Temperature * Internal energy

PHYSICAL MECHANISM OF NATURAL CONVECTION

CFD SIMULATIONS OF THE SPENT FUEL POOL IN THE LOSS OF COOLANT ACCIDENT

Effect Analysis of Volume Fraction of Nanofluid Al2O3-Water on Natural Convection Heat Transfer Coefficient in Small Modular Reactor

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

Effect of flow velocity on the process of air-steam condensation in a vertical tube condenser

SENSITIVITY ANALYSIS FOR ULOF OF PGSFR

Graduate Preliminary Examination, Thursday, January 6, Part I 2

The Research of Heat Transfer Area for 55/19 Steam Generator

EXPERIMENTAL INVESTIGATION OF CRITICAL HEAT FLUX FOR ZERO AND NATURAL CIRCULATION FLOW OF WATER IN THREE RODS BUNDLE NEAR ATMOSPHERIC PRESSURE

Demonstrate understanding of aspects of heat

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

ENGINEERING OF NUCLEAR REACTORS. Fall December 17, 2002 OPEN BOOK FINAL EXAM 3 HOURS

«CALCULATION OF ISOTOPE BURN-UP AND CHANGE IN EFFICIENCY OF ABSORBING ELEMENTS OF WWER-1000 CONTROL AND PROTECTION SYSTEM DURING BURN-UP».

4.1. Physics Module Form 4 Chapter 4 - Heat GCKL UNDERSTANDING THERMAL EQUILIBRIUM. What is thermal equilibrium?

A PWR HOT-ROD MODEL: RELAP5/MOD3.2.2Y AS A SUBCHANNEL CODE I.C. KIRSTEN (1), G.R. KIMBER (2), R. PAGE (3), J.R. JONES (1) ABSTRACT

A Fundamental Study on High Heat Flux Cooling using Subcooled Flow Boiling with Microbubble Emission

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

DELAY OF WETTING PROPAGATION DURING JET IMPINGEMENT QUENCHING FOR A HIGH TEMPERATURE SURFACE

Utilization of Egyptian Research Reactor and modes of collaboration

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

The three-level scaling approach with application to the Purdue University Multi-Dimensional Integral Test Assembly (PUMA)

Scaling Analysis as a part of Verification and Validation of Computational Fluid Dynamics and Thermal-Hydraulics software in Nuclear Industry

Profile SFR-64 BFS-2. RUSSIA

Some Results in studying NPP accidents using PCTRAN-2 LOOP and WWER-1000 simulators

The Pennsylvania State University. The Graduate School. College of Engineering

Lesson 14: Reactivity Variations and Control

Assessment of an isolation condenser performance in case of a LOHS using the RMPS+ Methodology

QUALIFICATION OF A CFD CODE FOR REACTOR APPLICATIONS

Sebastian Buchholz, Daniel von der Cron, Andreas Schaffrath. System codes improvements for modelling passive safety systems and their validation

English text only NUCLEAR ENERGY AGENCY COMMITTEE ON THE SAFETY OF NUCLEAR INSTALLATIONS

Work by Friction. A box slides 10 m across a surface. A frictional force of 20 N is acting on the box.

Transcription:

HWR Moderator Sub-cooling Requirements to Demonstrate Back-up Capabilities of Moderator During Accidents HEMANT KALRA NPCIL IAEA International Collaboration Standard Problem (ICSP) 1 st Workshop November, 19-21,2012 Ottawa CANADA

Working Team S. Hajela, Chief Participant Hemant Kalra V.V. Reddy

CONTENTS INTRODUCTION SALIENT FEATURE OF IPHWR ACCIDENT SCENARIO : LOCA & ECCS FAILURE CASE STUDY SUMMARY ICSP INPUT REQUIREMENTS

INTRODUCTION Nuclear Power Corporation of India Limited is a Public Sector Enterprise Under the administrative control of the Department of Atomic Energy (DAE), Government of India. NPCIL is responsible for design, construction, commissioning and operation of nuclear power reactors. NPCIL is presently operating 20 nuclear power reactors with an installed capacity of 4780 MWe. The reactor fleet comprises of 18 Pressurised Heavy Water Reactors (PHWRs) 7 units are under construction including four units of 700 MWe PHWRs

TYPICAL IPHWR 540MW(e) 7

Salient Features of IPHWR with respect to severe accident 392 Pressure tubes Horizontal orientation On power refueling Two independent diverse fast acting shutdown systems 8

Salient Features of IPHWR with respect to severe accident Heavy water coolant Fuel bundles lie inside the pressure tube Calandria tube separates the pressure tube from the cold moderator A gas gap between the pressure tube and the calandria tube provides thermal insulation Short bundles (0.5 meters) 9

LOCA AND ECCS FAILURE During accidents such as LOCA and ECCS Failure, the Pressure tube will Balloon/Sag or strain into contact with the calandria tube. The initial contact between hot pressure tube and cold calandria tube would result in a sudden increase in the heat flux to the moderator. The magnitude of peak heat flux combined with the temperature of moderator surrounding the calandria tube determines the boiling regime on the calandria tube surface and thus the rate of heat removal.

LOCA AND ECCS FAILURE If nucleate boiling is predicted to occur, the stored energy of the pressure tube is quickly removed and the fuel radiates to a low temperature pressure tube. If film boiling is the predicted mode of heat transfer, the stored energy in the pressure tube is not removed quickly. Heat generated in the fuel radiates to higher temperature pressure tube until the calandria tube surface rewets and the situation returns to nucleate boiling condition. NPCIL In-House developed Computer Code CONTACT simulates the Thermo-Mechanical Channel Behavior under such accident conditions.

LOCA AND ECCS FAILURE VOIDING IN CHANNEL FUEL TEMPERATURE INCREASES HEAT TRANSFER FUEL ELEMENT FUEL to PT to CT CT to MODERATOR BY CONDUCTION BY RADIATION BY CONVECTION HIGH TEMPERATURE METAL WATER REACTION D2 GENERATION CREEP BALLOON/SAG PT CT CONTACT HEAT FLUX MODERATOR ACTS AS A HEAT SINK

LOCA AND ECCS FAILURE FUEL CHANNEL MODEL UPPER SECTOR TOP BOTTOM LOWER SECTOR

LOCA AND ECCS FAILURE

HEAT TRANSFER MODELS Conduction Heat Transfer in Fuel Elements Radiation Heat Transfer Between Fuel Pressure Tube (PT) and PT-CT Gas Conduction between PT-CT Contact Conductance Between PT-CT Convective Heat Transfer between CT and Moderator BALLOONING & SAGGING MODEL Thermal Creep models for Sagging & Ballooning

16 HEAT TRANSFER MODELS S.No. Heat transfer Mode Boundaries 1. Natural convention TW <Tsat 2. Pool boiling TSat<TW<T 3. Transition boiling TCHF<TW<Tmin 4. Film boiling TW>Tmin LOCA AND ECCS FAILURE t T c q T k r r r T k r r p δ δ ρ θ δ θ δ δ δ δ δ δ δ = + + " 1 1 2 Conduction in Fuel Elements with Heat Flux at Surface as Boundary Condition Heat Transfer Between Calandria Tube to Moderator

Pressure Tube Deformation Ballooning Incase of LOCA where pressure tube sees high temperature and pressure especially in critical break LOCA Sagging In case of LOCA+ECC failure where the system pressure is low after blowdown 17

Modeling Concept Pressure tube Sagging Calandria Tube Roll Joint Garter spring Pressure Tube Initial Segment Deflection Axial Deformation Strained segment Note: There is a local deformation of pressure tube around each Garter spring 18

LOCA AND ECCS FAILURE METAL WATER REACTION Zr +2D 2 O ZrO 2 + 2D 2 +ENERGY Kinetics of oxidation dφ φ = dt K φ 2 Where K φ = A exp(-q/rt) where φ is mass of Zr, oxidized per unit area, [mg/cm2] Kφ is parabolic rate constant [ mg/cm2]2 S-1 The kinetic data of zircaloy steam reaction is considered for three corresponding temperature ranges 19

Effect of Moderator Temperature The type of deformation and rate of heat up determine the temperature at which the pressure tube contacts the calandria tube. The parameter which determines the boiling regime on the calandria tube when contact occurs is the moderator subcooling. This is because CHF increases with increase in degree of subcooling. The degree of sub cooling affects the time period during which dry out exists. Period of dry-out is one of the detrimental factor for PT-CT failure.

Results Sample case; IPHWR-540MW(e) LOCA + ECC failure + Moderator as a heat sink Using Computer code CONTACT Objective To prove moderator is an effective heat sink Fuel temperature is maintained below melting temperature Channel integrity is maintained so as to prevent loss of cooling/moderator inventory NOTE; During LOCA only one loop i.e. half the core is affected 21

LOCA AND ECCS FAILURE Fuel & PT Temperature In Upper Sector Temperature(c) 1800.00 1600.00 1400.00 1200.00 1000.00 800.00 600.00 CONTACT:SEVERE ACCIDENT ANALYSIS FOR TAPP3&4 LONG-TERM SYSTEM RESPONSE FOR 100% BREAK IN ROH WITH ECCS UNAVAILABLE Max.Rated Bundle of 5500Kw Channel Legend Title Center Bottom Center Top Second Bottom Second Top Third Bottom, Third Top Fourth Bottom Fourth Top Pr.TubeTop Cal. Tube Top 400.00 200.00 0.00 0.00 600.00 1200.00 1800.00 2400.00 3000.00 3600.00 Time(sec.) Figure F.21 a:sheath,pr. Tube and Cal. Tube Temperature Upward Sector HK,RSA 22

LOCA AND ECCS FAILURE Fuel & PT Temperature In Lower Sector Temperature(c) 1800.00 1600.00 1400.00 1200.00 1000.00 800.00 600.00 CONTACT:SEVERE ACCIDENT ANALYSIS FOR TAPP3&4 LONG-TERM SYSTEM RESPONSE FOR 100% BREAK IN ROH WITH ECCS UNAVAILABLE Max.Rated Bundle of 5500Kw Channel Legend Title Center Bottom Center Top Second Bottom Second Top Third Bottom, Third Top Fourth Bottom Fourth Top Pr.Tube Bottom Cal. Tube Bottom 400.00 200.00 0.00 0.00 600.00 1200.00 1800.00 2400.00 3000.00 3600.00 Time(sec.) Figure F.21 b:sheath, Pr.Tube and Cal. Tube Temperature Downward Sector HK,RSA 23

Dry Out Period

Summary Acceptance criteria for LOCA and ECCS failure LOCA AND ECCS FAILURE Channel integrity is to be maintained Channel integrity is maintained if the calandria tube remains intact after pressure tube contacts. This criterion is satisfied if the calandria tube outer surface does not go into prolonged film boiling. Analysis shows that dry-out occurs at few local points of calandria tube for a short period. 25

Requirements to be looked at LOCA AND ECCS FAILURE Relation between Channel integrity & dry out and its time period Note; An experiment is planned at CMRI, India using gas burner as heaters Gap conductance between PT-CT and its time dependency CHF on calandria tube and its effect on dry out period Note; Experiments are planned at IIT, Roorkee to create a data base of CHF on various modified surfaces Slumping of fuel and its effect on conduction 26

ICSP INPUT Requirements 1. Geometric Details 2. Loading Conditions 3. Heat Load 4. Thermal & Physical Properties of Materials 5. Creep parameters 6. Contact Conductance 7. Heat Transfer parameters