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

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

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

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

XV. Fission Product Poisoning

Lecture 28 Reactor Kinetics-IV

Reactivity Coefficients

Lecture 27 Reactor Kinetics-III

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

Reactor Operation with Feedback Effects

Lesson 14: Reactivity Variations and Control

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

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

Solving Bateman Equation for Xenon Transient Analysis Using Numerical Methods

Fundamentals of Nuclear Reactor Physics

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

Introduction to Reactivity and Reactor Control

Figure 22.1 Unflattened Flux Distribution

Reactor Operation Without Feedback Effects

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

Chem 481 Lecture Material 4/22/09

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

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

Physics Codes and Methods for CANDU Reactor

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

Delayed neutrons in nuclear fission chain reaction

17 Neutron Life Cycle

Reactor Kinetics and Operation

Chemical Engineering 412

NEUTRONIC CALCULATION SYSTEM FOR CANDU CORE BASED ON TRANSPORT METHODS

R.A. Chaplin Department of Chemical Engineering, University of New Brunswick, Canada

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

Comparative Safety Evaluation of Thorium Fuel to. Natural Uranium Fuel in a CANDU 6 Reactor

Coupled Neutronics Thermalhydraulics LC)CA Analysis

Nuclear Theory - Course 227 FISSION PRODUCT POISONING

Nuclear Fission. Q for 238 U + n 239 U is 4.??? MeV. E A for 239 U 6.6 MeV MeV neutrons are needed.

Numerical simulation of non-steady state neutron kinetics of the TRIGA Mark II reactor Vienna

Nuclear Theory - Course 227 EFFECT OF FUEL BURNUP

2. The Steady State and the Diffusion Equation

Year 11 Physics booklet Topic 1 Atomic structure and radioactivity Name:

LECTURE 4: CORE DESIGN ANALYSIS

Nuclear Theory - Course 227 REACTIVITY EFFECTS DUE TO TEMPERATURE CHANGES

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

PHYS-E0562 Ydinenergiatekniikan jatkokurssi Lecture 5 Burnup calculation

Reactivity Coefficients

CHAPTER NEUTRON DETECTORS

Nuclear Theory - Course 127 FISSION

Control of the fission chain reaction

Nuclear Theory - Course 227 NEUTRON FLUX DISTRI~UTION

Determination of research reactor fuel burnup

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

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

Power Changes in a Critical Reactor. The Critical Reactor

Subcritical Multiplication and Reactor Startup

FUEL PERFORMANCE EVALUATION THROUGH IODINE ACTIVITY MONITORING K.ANANTHARAMAN, RAJESH CHANDRA

NPP Simulators for Education Workshop - Passive PWR Models

THORIUM SELF-SUFFICIENT FUEL CYCLE OF CANDU POWER REACTOR

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

Name Date Class NUCLEAR RADIATION. alpha particle beta particle gamma ray

PhD Qualifying Exam Nuclear Engineering Program. Part 1 Core Courses

On the Use of Serpent for SMR Modeling and Cross Section Generation

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

General Introduction to Methods of Physics Analysis and Physics Computer Codes

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

TRANSMUTATION OF CESIUM-135 WITH FAST REACTORS

Reactor Physics: Basic Definitions and Perspectives. Table of Contents

Operational Reactor Safety

3. Detector Systems. This chapter describes the CANDU detector systems.

Chain Reactions. Table of Contents. List of Figures

Radioactivity. Radioactivity

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

Nuclear Fission. 1/v Fast neutrons. U thermal cross sections σ fission 584 b. σ scattering 9 b. σ radiative capture 97 b.

X. Neutron and Power Distribution

Binding Energy and Mass defect

Homework 06. Nuclear

Physics 3204 UNIT 3 Test Matter Energy Interface

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

P.M. THURSDAY, 16 January hour

Lecture 18 Neutron Kinetics Equations

A COMPARISON OF PEBBLE MIXING AND DEPLETION ALGORITHMS USED IN PEBBLE-BED REACTOR EQUILIBRIUM CYCLE SIMULATION

Neutronic Issues and Ways to Resolve Them. P.A. Fomichenko National Research Center Kurchatov Institute Yu.P. Sukharev JSC Afrikantov OKBM,

2017 Water Reactor Fuel Performance Meeting September 10 (Sun) ~ 14 (Thu), 2017 Ramada Plaza Jeju Jeju Island, Korea

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

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

Energy Dependence of Neutron Flux

Validation of Rattlesnake, BISON and RELAP-7 with TREAT

Information Nuclide = is an atomic species characterized by the specific constitution of its nucleus (protons and neutrons) Neutron

u d Fig. 6.1 (i) Identify the anti-proton from the table of particles shown in Fig [1]

Nuclear Chemistry. The Nucleus. Isotopes. Slide 1 / 43. Slide 2 / 43. Slide 3 / 43

Reactivity Effect of Fission and Absorption

SUB-CHAPTER D.1. SUMMARY DESCRIPTION

"Control Rod Calibration"

GASEOUS AND VOLATILE FISSION PRODUCT RELEASE FROM MOLTEN SALT NUCLEAR FUEL

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

Review A Z. a particle. proton. neutron. electron e -1. positron. e +1. Mass Number Atomic Number. Element Symbol

The discovery of nuclear reactions need not bring about the destruction of mankind any more than the discovery of matches - Albert Einstein

Introduction to Nuclear Reactor Kinetics

Research Article Calculations for a BWR Lattice with Adjacent Gadolinium Pins Using the Monte Carlo Cell Code Serpent v.1.1.7

Laplace Transform Problems

A Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis

Effect of Fission Fragments on the Properties of UO 2 Fuel of Pressurized Water Reactors

Transcription:

8. Effects of 135Xe The xenon isotope 135 Xe plays an important role in any power reactor. It has a very large absorption cross section for thermal neutrons and represents therefore a considerable load on the chain reaction. The 135 Xe concentration has an impact on the power distribution, and in turn is affected by the power distribution, by changes in power, and by movements of reactivity devices.

8.1 The Xe-I I Kinetics The 135 Xe/ 135 I kinetics are shown schematically in Figure 8.1. 135 Xe is produced to some degree directly in fission, but mostly as the result of the beta decay of its precursor 135 I (which has a half-life of 6.585 hours).

8.1 The Xe-I I Kinetics 135 Xe is destroyed in two ways: through its own radioactive decay ( 135 Xe has a half-life of 9.169 hours), and by absorption of neutrons to form 136 Xe, 135 I is a direct product of fission, but can also appear through the radioactive decay chain 135 Te to 135 Sb to 135 I. 135 Te and 135 Sb have half-lives which are very short (19.0 s and 1.71 s) compared to those of 135 I and of 135 Xe; it is sufficient to model the decay of 135 Te and 135 Sb as instantaneous, and add their fission yields to that of 135 I.

8.1 The Xe-I I Kinetics The 135 Xe/ 135 I kinetics in any particular fuel bundle can thus be represented by the following equations: di dt dx dt = γ Σ$ $ φ -λi i f F i = γ Σ$ $ φ + λi λ X σ$ X $ φ x f F i x X F (8.1) (8.2) where:

8.1 The Xe-I I Kinetics X = avge. conc. of 135 Xe in the bundle (atoms cm -3) I = avge. conc. of 135 I in the bundle (atoms cm -3) γ x = direct yield of 135 Xe per fission γ I = direct yield of 135 I in fission, including yields of 135 Te and 135 Sb λ x = decay constant of 135 Xe in s -1 λ I = decay constant of 135 I in s -1 $ φ = bundle-average flux in the fuel (n.cm -2 s -1 ) $Σ = macroscopic fission cross section of the fuel (cm -1 f ) and σ$ = microscopic 135 Xe cross section (cm 2 x )

8.1 The Xe-I I Kinetics $ $ Σ f F φ In the above equations, the term γ i gives the 135 I production rate, while λ i I gives the 135 I loss rate (and the production rate of 135 Xe by iodine decay). Similarly, the term γ x Σ$ $ f φ F gives the production rate of 135 Xe due to direct fission, while λ x X gives its decay rate. The term σ$ X $ x φf represents the destruction (burnout) rate of 135 Xe due to neutron capture. Because of the comparable magnitudes of the various terms, the 135 Xe concentration is very sensitive to changes in flux level.

8.1 The Xe-I I Kinetics The large absorption cross section of 135 Xe plays significant role in the overall neutron balance in the reactor, and thus directly affects the system reactivity, both in steady state and in transients. The 135 Xe/ 135 I kinetics also influences the spatial power distribution in the reactor.

8.2 Reactor Startup When a reactor is first started, or restarted after a long shutdown, the 135 Xe concentration will build up in all fuel bundles according to the equations derived above. In Figure 8.2, the variation of system reactivity as a function of time following startup is given for different final steady-state power levels. It can be seen that it takes ~ 40 hours for the 135 Xe concentration to fully reach equilibrium.

8.3 Steady-State State Xenon Load At steady state the time derivatives di/dt and dx/dt are zero. The above equations can then be solved to give the steady state concentrations of 135 I and 135 Xe (I ss and X ss ): I ss $ $ i = γ Σ φ f F λ i γ γ Σ$ $ φ ( + ) i x f X ss = λ + $ $ x x F σ φ F (8.3) (8.4)

8.3 Steady-State State Xenon Load It is obvious from these equations that, as a function of an increasing fuel flux, the steady-state concentration of 135 I increases indefinitely, while in contrast that of 135 Xe tends to an asymptotic value which will be denoted X ss, : ( ) $ X = γi + γ x Σ f ss, σ$ (8.5) x This asymptotic nature of the variation of X ss with is the reason why 135 Xe is termed a saturating fission product. (Other saturating fission products are 105 Rh, 149 Sm, 151 Sm, etc.) $ φ F

8.3 Steady-State State Xenon Load The limiting 135 Xe absorption rate at very high flux levels leads to a maximum reactivity of ~ -30 mk. In CANDU the equilibrium xenon load is approximately -28 mk. The flux level at full power in CANDU is such that the 135 Xe concentration is about 95% saturated, i.e., the average 135 Xe concentration is equal to about 95% of the value in an infinite flux. However, the steady-state 135 Xe concentration is not uniform in the core, but varies with flux according to Eq. (8.4). This is discussed in greater detail in Section 8.6 below.

8.4 Effect of Power Changes on Xenon Concentration Due to the presence of the term σ x X, the variation of the 135 Xe concentration with flux is non-linear. The 135 Xe reactivity following power (flux) changes will depend on: the starting power level, the time at that level, the new power level, and the time spent at the new power level. $ φ F

8.4 Effect of Power Changes on Xenon Concentration Generally speaking, when the power is reduced from a steady level, the 135 Xe concentration increases at first. This is due to the fact that 135 Xe is still being produced by the decay of 135 I, but its burnout rate (by neutron absorption) is decreased because of the reduced neutron flux (reduced power).

8.4 Effect of Power Changes on Xenon Concentration However, after a certain period (depending on the initial and final power and the rate of power reduction) the 135 I decay rate decreases sufficiently (due to the lower fission rate) that the rate of 135 Xe production drops below the rate of 135 Xe decay (and burnout). At this time, then, the 135 Xe concentration reaches a peak value and starts to decrease towards a new (lower) steady-state level.

8.4 Effect of Power Changes on Xenon Concentration Conversely, when the power is increased from a steady level, the 135 Xe concentration will first decrease, and then go through a minimum and start increasing again to a higher steady-state level.

8.4 Effect of Power Changes on Xenon Concentration Fig. 8.3 shows some typical reactivity variations due to 135 Xe following step reductions in power. Very similar variations, but in the opposite direction, ensue upon step increases in power. The quantitative effects will be different at different points in the core, due to the initial non-uniform distribution of 135 Xe. Thus, for an accurate assessment of xenon transients on the power distribution, a point-kinetics treatment is generally inadequate, and calculations in three dimensions will be required.

8.5 Xenon Transient Following a Shutdown Following a reactor shutdown, the burnout of 135 Xe stops, whereas the production by means of 135 I decay continues for several hours. The net result is that there is an initial increase in 135 Xe concentration and a decrease in core reactivity. If the reactor is required to be started up shortly after shutdown, extra positive reactivity must be supplied. The 135 Xe growth and decay following a shutdown for a typical CANDU is shown in Figure 8.4.

8.5 Xenon Transient Following a Shutdown It can be seen that, at about 10 hours after shutdown, the reactivity worth of 135 Xe increases to several times its equilibrium full-power value. At ~35-40 hours the 135 Xe has decayed back to its pre-shutdown level. If it were not possible to add positive reactivity during this period, every shutdown would necessarily last some 40 hours, when the reactor would again reach criticality.

8.5 Xenon Transient Following a Shutdown To achieve xenon override and permit power recovery following a shutdown (or reduction in reactor power), the CANDU-6 adjuster rods are withdrawn to provide positive reactivity. It is not possible to provide complete xenon override capability, this would require > 100 mk of positive reactivity. The CANDU-6 adjuster rods provide approximately 15 milli-k of reactivity, which is sufficient for about 30 minutes of xenon override following a shutdown.

8.6 Effects of Xenon on Power Distribution Xenon plays a role in the core power distribution. First, on the steady-state power distribution: Because the steady-state 135 Xe concentration depends on the flux (Eq. 8.2), high-power bundles will have a higher xenon load, and therefore a lower reactivity, than low-power bundles of the same irradiation. (cont d)

8.6 Effects of Xenon on Power Distribution The effect of xenon is therefore to flatten the steadystate power distribution: The reduction in the maximum bundle power due to the local 135 Xe concentration can be of the order of 5%. For maximum channel power, ~3%. Thus xenon helps to increase the margin between peak powers and licensed limits. The xenon effect should be taken into account when accurate calculations are desired.

8.6 Effects of Xenon on Power Distribution The xenon effect also plays a role after refuelling: Fresh bundles introduced into the reactor have no xenon (or other fission products) at first. Fig. 8.2 shows it takes a day or so for the xenon to build to near saturation. Thus, the power of these fresh bundles is high at first and subsequently decreases, by perhaps several per cent, to their equilibrium value.

8.6 Effects of Xenon on Power Distribution Also, bundles which are shifted along the refuelled channel experience a change in power. In these bundles, corresponding to the specific power change (decrease or increase), there is a xenon transient similar to those in Fig. 8.3 (or the opposite ). This effect is not limited to shifted bundles, but occurs in any bundle whose power is affected by the channel refuelling operation, e.g. bundles in neighbouring channels.

8.6 Effects of Xenon on Power Distribution Transient-xenon effects due to refuelling are perhaps not always taken into account in routine core-follow calculations, but need to be captured when it is desired to perform the most accurate simulations, and to properly monitor the margin between peak powers and licensed limits.

8.6 Effects of Xenon on Power Distribution The RFSP history-based local-parameter methodology can simulate the spatially non-uniform effects of saturatingfission-product kinetics for both Xe-I and other fission-product pairs (e.g. Sm-Pm) in both steady-state and transient situations, by recalculating the 135 Xe concentration individually in every bundle, and therefore the effect on the lattice properties and on the core power distribution.

Summary 135 Xe is an important poison, which makes itself felt at all times. Its spatial and time effects must be taken into account when accurate fuel-management calculations are required.