Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 4. Title: Control Rods and Sub-critical Systems

Similar documents
Chem 481 Lecture Material 4/22/09

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 5. Title: Reactor Kinetics and Reactor Operation

Reactor Kinetics and Operation

Chapter 7 & 8 Control Rods Fission Product Poisons. Ryan Schow

Shutdown Margin. Xenon-Free Xenon removes neutrons from the life-cycle. So, xenonfree is the most reactive condition.

Introduction to Reactivity and Reactor Control

20.1 Xenon Production Xe-135 is produced directly in only 0.3% of all U-235 fissions. The following example is typical:

Operational Reactor Safety

XV. Fission Product Poisoning

Lecture 28 Reactor Kinetics-IV

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

Lesson 14: Reactivity Variations and Control

Reactivity Coefficients

NUCLEAR SCIENCE ACAD BASIC CURRICULUM CHAPTER 5 NEUTRON LIFE CYCLE STUDENT TEXT REV 2. L th. L f U-235 FUEL MODERATOR START CYCLE HERE THERMAL NEUTRON

(1) The time t required for N generations to elapse is merely:

Lecture 27 Reactor Kinetics-III

Subcritical Multiplication and Reactor Startup

CRITICAL AND SUBCRITICAL EXPERIMENTS USING THE TRAINING NUCLEAR REACTOR OF THE BUDAPEST UNIVERSITY OF TECHNOLOGY AND ECONOMICS

X. Neutron and Power Distribution

Reactivity Balance & Reactor Control System

PHYSICS AND KINETICS OF TRIGA REACTOR. H. Böck and M. Villa AIAU 27307

Chemical Engineering 412

The moderator temperature coefficient MTC is defined as the change in reactivity per degree change in moderator temperature.

Fundamentals of Nuclear Reactor Physics

REACTOR PHYSICS FOR NON-NUCLEAR ENGINEERS

but mostly as the result of the beta decay of its precursor 135 I (which has a half-life of hours).

"Control Rod Calibration"

Xenon Effects. B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec.

Reactor Operation Without Feedback Effects

VI. Chain Reaction. Two basic requirements must be filled in order to produce power in a reactor:

3. State each of the four types of inelastic collisions, giving an example of each (zaa type example is acceptable)

N U C L : R E A C T O R O P E R A T I O N A N D R E G U L A T O R Y P O L I C Y, I

Hybrid Low-Power Research Reactor with Separable Core Concept

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

Lecture 20 Reactor Theory-V

Nuclear Theory - Course 227 REACTIVITY EFFECTS DUE TO TEMPERATURE CHANGES

Power Changes in a Critical Reactor. The Critical Reactor

MEASUREMENT AND COMPARISON OF CONTROL ROD WORTH OF BTRR USING INHOUR EQUATION AND PERIOD REACTIVITY CONVERSION TABLE

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

Lesson 8: Slowing Down Spectra, p, Fermi Age

CONTROL ROD WORTH EVALUATION OF TRIGA MARK II REACTOR

Application of the Dynamic Control Rod Reactivity Measurement Method to Korea Standard Nuclear Power Plants

Lamarsh, "Introduction to Nuclear Reactor Theory", Addison Wesley (1972), Ch. 12 & 13

A Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis

Estimation of Control Rods Worth for WWR-S Research Reactor Using WIMS-D4 and CITATION Codes

Verification of Predicted Powers and Shim Positions for Advanced Test Reactor Cycle 159A-1

2. The Steady State and the Diffusion Equation

Neutron reproduction. factor ε. k eff = Neutron Life Cycle. x η

Lesson 9: Multiplying Media (Reactors)

Lecture 18 Neutron Kinetics Equations

Reactivity Coefficients

Study of Control rod worth in the TMSR

Control Rod Calibration and Worth Calculation for Optimized Power Reactor 1000 (OPR-1000) Using Core Simulator OPR1000

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

Reactor Operation with Feedback Effects

A Dummy Core for V&V and Education & Training Purposes at TechnicAtome: In and Ex-Core Calculations

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 1. Title: Neutron Life Cycle

Solving Bateman Equation for Xenon Transient Analysis Using Numerical Methods

SUB-CHAPTER D.1. SUMMARY DESCRIPTION

17 Neutron Life Cycle

NEUTRONIC CALCULATION SYSTEM FOR CANDU CORE BASED ON TRANSPORT METHODS

Neutronic Calculations of Ghana Research Reactor-1 LEU Core

OECD/NEA Transient Benchmark Analysis with PARCS - THERMIX

A Method For the Burnup Analysis of Power Reactors in Equilibrium Operation Cycles

Control of the fission chain reaction

Fuel Element Burnup Determination in HEU - LEU Mixed TRIGA Research Reactor Core

PWR CONTROL ROD EJECTION ANALYSIS WITH THE MOC CODE DECART

Chain Reactions. Table of Contents. List of Figures

Calculation of Control Rod Worth of TRIGA Mark II Reactor Using Evaluated Nuclear Data Library JEFF-3.1.2

Simple Method to Predict Power Level and Core Flow Rate of Boiling Water Reactors by Using One-Point Core Model

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

REACTOR PHYSICS ASPECTS OF PLUTONIUM RECYCLING IN PWRs

VALIDATION OF BATAN'S STANDARD 3-D DIFFUSION CODE, BAT AN-3DIFF, ON THE FIRST CORE OF RSG GAS. Liem Peng Hong"

EXPERIMENTAL DETERMINATION OF NEUTRONIC PARAMETERS IN THE IPR-R1 TRIGA REACTOR CORE

Considerations for Measurements in Support of Thermal Scattering Data Evaluations. Ayman I. Hawari

USA HTR NEUTRONIC CHARACTERIZATION OF THE SAFARI-1 MATERIAL TESTING REACTOR

Radioactivity III: Measurement of Half Life.

Reactivity Effect of Fission and Absorption

MONTE CALRLO MODELLING OF VOID COEFFICIENT OF REACTIVITY EXPERIMENT

Reactors and Fuels. Allen G. Croff Oak Ridge National Laboratory (ret.) NNSA/DOE Nevada Support Facility 232 Energy Way Las Vegas, NV

The Effect of 99 Mo Production on the Neutronic Safety Parameters of Research Reactors

NEUTRON MODERATION. LIST three desirable characteristics of a moderator.

The Effect of Burnup on Reactivity for VVER-1000 with MOXGD and UGD Fuel Assemblies Using MCNPX Code

Reactivity Power and Temperature Coefficients Determination of the TRR

High-Order Finite Difference Nodal Method for Neutron Diffusion Equation

Excerpt from the Proceedings of the COMSOL Users Conference 2007 Grenoble

Chapter 2 Nuclear Reactor Calculations

Physics Codes and Methods for CANDU Reactor

Chemical Engineering 693R

Study of boron dilution phenomenon in the core and fuel assemblies of Bushehr VVER-1000 reactor in normal operating conditions

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

Applied Reactor Technology and Nuclear Power Safety, 4A1627; 4 cp. Course Description

Fundamentals of Nuclear Power. Original slides provided by Dr. Daniel Holland

Cold Critical Pre-Experiment Simulations of KRUSTy

LECTURE 4: CORE DESIGN ANALYSIS

Dynamic Analysis on the Safety Criteria of the Conceptual Core Design in MTR-type Research Reactor

General Introduction to Methods of Physics Analysis and Physics Computer Codes

Figure 22.1 Unflattened Flux Distribution

Metropolitan Community College COURSE OUTLINE FORM LAB: 3.0

INTRODUCTION TO NUCLEAR REACTORS AND NUCLEAR POWER GENERATION. Atsushi TAKEDA & Hisao EDA

Transcription:

Lectures on Nuclear Power Safety Lecture No 4 Title: Control Rods and Sub-critical Systems Department of Energy Technology KTH Spring 2005 Slide No 1

Outline of the Lecture Control Rods Selection of Control Rod Materials Types of Control Rods Control Rod Effectiveness Integral and Differential Control Rod Worth Rod Control Mechanisms Sub-critical systems Subcritical Multiplication Factor Effect of Reactivity Changes on Sub-critical Multiplication Use of 1/M Plots Slide No 2

Selection of Control Rod Material (1) Rods of neutron-absorbing material are installed in most reactors to provide precise, adjustable control of reactivity These rods are able to be moved into or out of the reactor core and typically contain elements such as silver, indium, cadmium, boron, or hafnium The material used for the control rods varies depending on reactor design Generally, the material selected should have a good absorption cross section for neutrons and have a long lifetime as an absorber (not burn out rapidly) Slide No 3

Selection of Control Rod Material (2) The ability of a control rod to absorb neutrons can be adjusted during manufacture A control rod that is referred to as a "black" absorber absorbs essentially all incident neutrons. A "grey" absorber absorbs only a part of them While it takes more grey rods than black rods for a given reactivity effect, the grey rods are often preferred because they cause smaller depressions in the neutron flux and power in the vicinity of the rod This leads to a flatter neutron flux profile and more even power distribution in the core Slide No 4

Selection of Control Rod Material (3) If grey rods are desired, the amount of material with a high absorption cross section that is loaded in the rod is limited Material with a very high absorption cross section may not be desired for use in a control rod, because it will burn out rapidly due to its high absorption cross section The same amount of reactivity worth can be achieved by manufacturing the control rod from material with a slightly lower cross section and by loading more of the material This also results in a rod that does not burn out as rapidly Slide No 5

Selection of Control Rod Material (4) Another factor in control rod material selection is that materials that resonantly absorb neutrons are often preferred to those that merely have high thermal neutron absorption cross sections Resonance neutron absorbers absorb neutrons in the epithermal energy range The path length traveled by the epithermal neutrons in a reactor is greater than the path length traveled by thermal neutrons Therefore, a resonance absorber absorbs neutrons that have their last collision farther (on the average) from the control rod than a thermal absorber This has the effect of making the area of influence around a resonance absorber larger than around a thermal absorber and is useful in maintaining a flatter flux profile Slide No 6

Types of Control Rods (1) There are several ways to classify the types of control rods One classification method is by the purpose of the control rods Three purposes of control rods are listed below Shim rods - used for coarse control and/or to remove reactivity in relatively large amounts Regulating rods - used for fine adjustments and to maintain desired power or temperature Safety rods - provide a means for very fast shutdown in the event of an unsafe condition. Addition of a large amount of negative reactivity by rapidly inserting the safety rods is referred to as a "scram" or "trip" Slide No 7

Types of Control Rods (2) Not all reactors have different control rods to serve the purposes mentioned above Depending upon the type of reactor and the controls necessary, it is possible to use dual-purpose or even triple-purpose rods For example, consider a set of control rods that can insert enough reactivity to be used as shim rods If the same rods can be operated at slow speeds, they will function as regulating rods Additionally, these same rods can be designed for rapid insertion, or scram These rods serve a triple function yet meet other specifications such as precise control, range of control, and efficiency Slide No 8

Types of Control Rods (3) Examples of BWR control rods Slide No 9

Types of Control Rods (5) Example of PWR control rods Slide No 10

Control Rod Effectiveness (1) The effectiveness of a control rod depends largely upon the value of the ratio of the neutron flux at the location of the rod to the average neutron flux in the reactor The control rod has maximum effect (inserts the most negative reactivity) if it is placed in the reactor where the flux is maximum If a reactor has only one control rod, the rod should be placed in the center of the reactor core The effect of such a rod on the flux is illustrated in Figure on the next slide Slide No 11

Control Rod Effectiveness (2) Slide No 12

Control Rod Effectiveness (3) If additional rods are added to this simple reactor, the most effective location is where the flux is maximum, that is, at point A Numerous control rods are required for a reactor that has a large amount of excess reactivity (that amount of reactivity in excess of that needed to be critical) The exact amount of reactivity that each control rod inserts depends upon the reactor design The change in reactivity caused by control rod motion is referred to as control rod worth Slide No 13

Control Rod Worth (1) The exact effect of control rods on reactivity can be determined experimentally For example, a control rod can be withdrawn in small increments, such as 1 cm, and the change in reactivity can be determined following each increment of withdrawal By plotting the resulting reactivity versus the rod position, a graph similar to Figure on the next slide is obtained The graph depicts integral control rod worth over the full range of withdrawal The integral control rod worth is the total reactivity worth of the rod at that particular degree of withdrawal and is usually defined to be the greatest when the rod is fully withdrawn Slide No 14

Control Rod Worth (2) Slide No 15

Control Rod Worth (3) The slope of the curve ( ρ/ x), and therefore the amount of reactivity inserted per unit of withdrawal, is greatest when the control rod is midway out of the core This occurs because the area of greatest neutron flux is near the center of the core; therefore, the amount of change in neutron absorption is greatest in this area If the slope of the curve for integral rod worth in Figure on the previous slide is taken, the result is a value for rate of change of control rod worth as a function of control rod position Slide No 16

Control Rod Worth (4) A plot of the slope of the integral rod worth curve, also called the differential control rod worth, is shown in Figure on the next slide At the bottom of the core, where there are few neutrons, rod movement has little effect so the change in rod worth per inch varies little As the rod approaches the center of the core its effect becomes greater, and the change in rod worth per unit length is greater At the center of the core the differential rod worth is greatest and varies little with rod motion From the center of the core to the top, the rod worth per unit length is basically the inverse of the rod worth per unit length from the center to the bottom Slide No 17

Control Rod Worth (5) Slide No 18

Control Rod Worth (6) Differential control rod worth is the reactivity change per unit movement of a rod and is normally expressed as ρ/inch, k/k per inch, or pcm/inch The integral rod worth at a given withdrawal is merely the summation of all the differential rod worths up to that point of withdrawal It is also the area under the differential rod worth curve at any given withdrawal position Slide No 19

Control Rod Worth (7) Consider a reactor with a central partially inserted control rod Applying a one-group diffusion approximation and using the perturbation theory it can be shown that the reactivity change ρ(x) depends on the insertion length as follows: ρ x H 1 2π 2πx H ( x) = ρ( H ) sin z a x H Where ρ(h) is the reactivity change due to full insertion of the control rod R Slide No 20

Control Rod Worth (8) When the control rod is placed at r = r 0 distance from the centerline of the reactor, the reactivity change for such a rod can be estimated as R ρ 2.4048r0 R 2 ( r ) = J () 0 0 0 ρ r 0 Where ρ(r 0 ) is the reactivity change for off-centerline rod and ρ(0) is the reactivity change for a rod inserted at the centerline Slide No 21

Rod Control Mechanisms (1) The control rod insertion rates on a scram are designed to be sufficient to protect the reactor against damage in all transients that are expected to occur during the life of the reactor During normal rod motion, the control rods must be able to move rapidly enough to compensate for the most rapid rate at which positive reactivity is expected to build within the reactor in order to provide positive control Slide No 22

Rod Control Mechanisms (2) The transient that is normally considered when setting this minimum rod speed is the burnout of maximum peak xenon while at full power Xenon burnout is usually the most rapid, non-accident transient expected The maximum rod speed is normally limited in order to reduce the severity of an accident involving the continuous withdrawal of control rods Slide No 23

Rod Control Mechanisms (3) Example of BWR control rod drive system System allows two independent ways of insertion of control rods: Normal operation with fine motion Hydraulic system for fast scram Slide No 24

Sub-critical Multiplication Factor (1) When a reactor is in a shutdown condition, neutrons are still present to interact with the fuel These source neutrons are produced by a variety of methods that were discussed earlier If neutrons and fissionable material are present in the reactor, fission will take place Therefore, a reactor will always be producing a small number of fissions even when it is shutdown Slide No 25

Sub-critical Multiplication Factor (2) Consider a reactor in which k eff is 0.6 If 100 neutrons are suddenly introduced into the reactor, these 100 neutrons that start the current generation will produce 60 neutrons (100 x 0.6) from fission to start the next generation The 60 neutrons that start the second generation will produce 36 neutrons (60 x 0.6) to start the third generation The number of neutrons produced by fission in subsequent generations due to the introduction of 100 source neutrons into the reactor is shown below Generation 1st 2nd 3rd 4th 5th 6th 7th 8th 9th 10th 11th 12th Neutrons 100 60 36 22 13 8 5 3 2 1 0 0 Slide No 26

Sub-critical Multiplication Factor (3) Because the reactor is subcritical, neutrons introduced in the reactor will have a decreasing effect on each subsequent generation The addition of source neutrons to the reactor containing fissionable material has the effect of maintaining a much higher stable neutron level due to the fissions occurring than the neutron level that would result from the source neutrons alone The effects of adding source neutrons at a rate of 100 neutrons per generation to a reactor with a k eff of 0.6 are shown in table on the next slide Slide No 27

Sub-critical Multiplication Factor (4) Generation 1st 2nd 3rd 4th 5th 6th 7th 8th 9th 10th 11th 12th 100 60 36 22 13 8 5 3 2 1 0 0 100 60 36 22 13 8 5 3 2 1 0 100 60 36 22 13 8 5 3 2 1 100 60 36 22 13 8 5 3 2 100 60 36 22 13 8 5 3 100 60 36 22 13 8 5 100 60 36 22 13 8 100 60 36 22 13 100 60 36 22 100 60 36 100 60 100 Total n 100 160 196 218 231 239 244 247 249 250 250... Slide No 28

Sub-critical Multiplication Factor (5) A neutron source strength of 100 neutrons per generation will result in 250 neutrons per generation being produced from a combination of sources and fission in a shutdown reactor with a k eff of 0.6 If the value of k eff were higher, the source neutrons would produce a greater number of fission neutrons and their effects would be felt for a larger number of subsequent generations after their addition to the reactor Slide No 29

Sub-critical Multiplication Factor (6) The effect of fissions in the fuel increasing the effective source strength of a reactor with a k eff of less than one is subcritical multiplication For a given value of k eff there exists a subcritical multiplication factor (M) that relates the source level to the steady-state neutron level of the core If the value of k eff is known, the amount that the neutron source strength will be multiplied (M) can easily be determined by Equation: M = 1/(1 k eff ) Slide No 30

Sub-critical Multiplication Factor (7) Example: calculate the subcritical multiplication factor for keff = 0.6 and keff = 0.986 Solution: M = 1/(1 k eff ) = 1/(1-0.6) = 2.5 and M = 1/(1 k eff ) = 1/(1-0.986) = 71.4 The example above illustrates that the subcritical multiplication factor will increase as positive reactivity is added to a shutdown reactor, increasing the value of k eff If the source strength of this reactor were 1000 neutrons/sec, the neutron level would increase from 2500 neutrons/second at a k eff of 0.6 to a neutron level of 71,400 neutrons/sec at a k eff of 0.986 Slide No 31

Effect of Reactivity Changes on Subcritical Multiplication (1) In a subcritical reactor, the neutron level is related to the source strength as: N = S x M N neutron level S neutron source strength M subcritical multiplication factor Or: N = S /(1 k eff ) Slide No 32

Effect of Reactivity Changes on Subcritical Multiplication (2) To this point it has been necessary to know the neutron source strength of the reactor in order to use the concept of subcritical multiplication In most reactors the actual strength of the neutron sources is difficult, if not impossible, to determine Even though the actual source strength may not be known, it is still possible to relate the change in reactivity to a change in neutron level Slide No 33

Effect of Reactivity Changes on Subcritical Multiplication (3) Consider a reactor at two different times when k eff is two different values, k 1 and k 2 The neutron level at each time can be determined based on the neutron source strength and the subcritical multiplication factor as follows N 1 = S / (1 k 1 ), N 2 = S / (1 k 2 ) Thus: N 1 /N 2 = (1 k 1 )/ (1 k 2 ) Slide No 34

Effect of Reactivity Changes on Subcritical Multiplication (4) The neutron level in a shutdown reactor is typically monitored using instruments that measure the neutron leakage out of the reactor The neutron leakage is proportional to the neutron level in the reactor Typical units for displaying the instrument reading are counts per second (cps) Because the instrument count rate (CR) is proportional to the neutron level, the above equation can be restated as shown below: CR 1 /CR 2 = (1 k 1 )/ (1 k 2 ) Slide No 35

Use of 1/M Plots (1) Because the subcritical multiplication factor is related to the value of k eff, it is possible to monitor the approach to criticality through the use of the subcritical multiplication factor As positive reactivity is added to a subcritical reactor, k eff will get nearer to one As k eff gets nearer to one, the subcritical multiplication factor (M) gets larger Slide No 36

Use of 1/M Plots (2) The closer the reactor is to criticality, the faster M will increase for equal step insertions of positive reactivity When the reactor becomes critical, M will be infinitely large For this reason, monitoring and plotting M during an approach to criticality is impractical because there is no value of M at which the reactor clearly becomes critical Slide No 37

Use of 1/M Plots (3) Instead of plotting M directly, its inverse (1/M) is plotted on a graph of 1/M versus rod height. 1/M = 1 - k eff As control rods are withdrawn and k eff approaches one and M approaches infinity, 1/M approaches zero Slide No 38

Use of 1/M Plots (4) For a critical reactor, 1/M is equal to zero A true 1/M plot requires knowledge of the neutron source strength Because the actual source strength is usually unknown, a reference count rate (CR) is substituted, and the calculation of the factor 1/M is through the use of the following Equation (where CR 0 reference count rate) 1/M = CR 0 /CR Slide No 39

Use of 1/M Plots (5) In practice, the reference count rate used is the count rate prior to the beginning of the reactivity change The startup procedures for many reactors include instructions to insert positive reactivity in incremental steps with delays between the reactivity insertions to allow time for subcritical multiplication to increase the steady-state neutron population to a new, higher level and allow more accurate plotting of 1/M The neutron population will typically reach its new steady-state value within 1-2 minutes, but the closer the reactor is to criticality, the longer the time will be to stabilize the neutron population Slide No 40

Use of 1/M Plots (6) Example: Given the following rod withdrawal data, construct a 1/M plot and estimate the rod position when criticality would occur. The initial count rate on the nuclear instrumentation prior to rod withdrawal is 50 cps. Rod Withdrawal (inches) Count Rate (cps) 2 55 4 67 6 86 8 120 10 192 12 500 Slide No 41

Use of 1/M Plots (7) Solution: step 1 Calculate 1/M for each of the rod positions The reference count rate is 50 cps at a rod position of zero Rod Withdrawal Count Rate CR 0 /CR (inches) (cps) 0 50 1 2 55 0.909 4 67 0.746 6 86 0.581 8 120 0.417 10 192 0.260 12 500 0.100 Slide No 42

Use of 1/M Plots (8) Solution: step 2 Plotting these values, as shown in Figure below, and extrapolating to a 1/M value of 0 reveals that the reactor will go critical at approximately 13 inches of rod withdrawal Slide No 43

Exercises Exercise 11: In a cylindrical reactor with hight H = 3.6 m and radius R = 2 m a centrally-placed control rod is fully inserted into the reactor. Calculate the reactivity increase if the rod is withdrawn 15 cm from the reactor. The control rod worth ρ(h) is equal to 2%. Exercise 12: In a cylindrical reactor with hight H = 3.6 m and radius R = 2 m a control rod is placed at r 0 = 1.2 m from the centerline. The rod is half-way withdrawn. Calculate the reactivity increase if the rod is withdrawn with additional 15 cm from the reactor. The centrallyplaced, identical control rod worth ρ(h) is equal to 2%. Slide No 44