Lesson 9: Multiplying Media (Reactors)

Similar documents
Lesson 6: Diffusion Theory (cf. Transport), Applications

Lesson 8: Slowing Down Spectra, p, Fermi Age

Lecture 20 Reactor Theory-V

X. Assembling the Pieces

X. Neutron and Power Distribution

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

2. The Steady State and the Diffusion Equation

VIII. Neutron Moderation and the Six Factors

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

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

Reactivity Balance & Reactor Control System

Serco Assurance. Resonance Theory and Transport Theory in WIMSD J L Hutton

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

Fundamentals of Nuclear Reactor Physics

Nuclear Reactor Physics I Final Exam Solutions

Excerpt from the Proceedings of the COMSOL Users Conference 2007 Grenoble

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

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

Control of the fission chain reaction

Lesson 14: Reactivity Variations and Control

17 Neutron Life Cycle

Chain Reactions. Table of Contents. List of Figures

POWER DENSITY DISTRIBUTION BY GAMMA SCANNING OF FUEL RODS MEASUREMENT TECHNIQUE IN RA-8 CRITICAL FACILITY

Reactivity Effect of Fission and Absorption

DESIGN OF B 4 C BURNABLE PARTICLES MIXED IN LEU FUEL FOR HTRS

Gauss s Law. 3.1 Quiz. Conference 3. Physics 102 Conference 3. Physics 102 General Physics II. Monday, February 10th, Problem 3.

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

The Henryk Niewodniczański INSTITUTE OF NUCLEAR PHYSICS Polish Academy of Sciences ul. Radzikowskiego 152, Kraków, Poland.

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

Delayed neutrons in nuclear fission chain reaction

Power Series. Part 1. J. Gonzalez-Zugasti, University of Massachusetts - Lowell

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

Lecture 27 Reactor Kinetics-III

Subcritical Multiplication and Reactor Startup

Source Neutron Amplification and "Criticality Safety Margin" of an Energy Amplifier

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

MNFFY 221/SIF 4082 Energy and Environmental Physics. 2002

III. Chain Reactions and Criticality

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

REACTOR PHYSICS ASPECTS OF PLUTONIUM RECYCLING IN PWRs

"Control Rod Calibration"

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

ENEN/Reactor Theory/ Laboratory Session 1 DETERMINATION OF BASIC STATIC REACTOR PARAMETERS IN THE GRAPHITE PILE AT THE VENUS FACILITY

On the use of SERPENT code for few-group XS generation for Sodium Fast Reactors

Potential & Potential Energy

Experimental and computational determination of Alpha modes of the Delphi subcritical assembly

Applicability of Point Reactor Kinetic Equations for Subcriticality Monitoring

NEUTRONIC CALCULATION SYSTEM FOR CANDU CORE BASED ON TRANSPORT METHODS

Breeding K.S. Rajan Professor, School of Chemical & Biotechnology SASTRA University

CHAPTER 4 Reactor Statics. Table of Contents

Boundary value problems

Lecture 13: Electromagnetic Theory Professor D. K. Ghosh, Physics Department, I.I.T., Bombay. Poisson s and Laplace s Equations

PHYS-E0562 Ydinenergiatekniikan jatkokurssi Lecture 4 Diffusion theory

High performance computing for neutron diffusion and transport equations

MONTE CALRLO MODELLING OF VOID COEFFICIENT OF REACTIVITY EXPERIMENT

Nuclear Theory - Course 227

CALCULATION OF TEMPERATURE REACTIVITY COEFFICIENTS IN KRITZ-2 CRITICAL EXPERIMENTS USING WIMS ABSTRACT

New methods implemented in TRIPOLI-4. New methods implemented in TRIPOLI-4. J. Eduard Hoogenboom Delft University of Technology

THREE-DIMENSIONAL INTEGRAL NEUTRON TRANSPORT CELL CALCULATIONS FOR THE DETERMINATION OF MEAN CELL CROSS SECTIONS

Nuclear Physics (13 th lecture)

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

Nuclear Theory - Course 127 EFFECTS OF FUEL BURNUP

in U. S. Department of Energy Nuclear Criticality Technology and Safety Project 1995 Annual Meeting SanDiego, CA

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

Study of the End Flux Peaking for the CANDU Fuel Bundle Types by Transport Methods

Low-Grade Nuclear Materials as Possible Threats to the Nonproliferation Regime. (Report under CRDF Project RX0-1333)

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

A MONTE CARLO METHOD FOR CALCULATING INITIATION PROBABILITY

Chem 481 Lecture Material 4/22/09

Student-level numerical simulation of conditions inside an exploding fission-bomb core

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

Characterization of a Portable Neutron Coincidence Counter Angela Thornton and W. Charlton Texas A&M University College Station, TX

Efficient Energy Conversion of the 14MeV Neutrons in DT Inertial Confinement Fusion. By F. Winterberg University of Nevada, Reno

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

AP Physics C. Gauss s Law. Free Response Problems

The collision probability method in 1D part 2

Effect of WIMSD4 libraries on Bushehr VVER-1000 Core Fuel Burn-up

1 Introduction. Green s function notes 2018

Critical Experiment Analyses by CHAPLET-3D Code in Two- and Three-Dimensional Core Models

Solving Bateman Equation for Xenon Transient Analysis Using Numerical Methods

PhD Qualifying Exam Nuclear Engineering Program. Part 1 Core Courses

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

Lecture 18 Neutron Kinetics Equations

A Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis

MURI teleconference 28 May Optical Antimatter. John Pendry and Sebastien Guenneau Imperial College London. 24 May 2004 page 1

Nuclear Theory - Course 127 NEUTRON BALANCE AND THE FOUR FACTOR FORMULA

Lecture 10. Central potential

Energy Dependence of Neutron Flux

Postprint.

Form #221 Page 1 of 7

Theoretical Task 3 (T-3) : Solutions 1 of 9

REFLECTOR FEEDBACK COEFFICIENT: NEUTRONIC CONSIDERATIONS IN MTR-TYPE RESEARCH REACTORS ABSTRACT

Fundamental Stellar Parameters. Radiative Transfer. Stellar Atmospheres

Module 6. Neutron Flux Distribution 6.1 MODULE OVERVIEW '!' MODULE OBJECTIVES OVERALL FLUX SHAPE FLUX FLATTENING 6

ELECTRO MAGNETIC FIELDS

Nuclear Theory - Course 227 NEUTRON FLUX DISTRI~UTION

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

THORIUM SELF-SUFFICIENT FUEL CYCLE OF CANDU POWER REACTOR

Chemical Reaction between Solids x 2 = Kt

Reactivity Coefficients

Transcription:

Lesson 9: Multiplying Media (Reactors) Laboratory for Reactor Physics and Systems Behaviour Multiplication Factors Reactor Equation for a Bare, Homogeneous Reactor Geometrical, Material Buckling Spherical, Slab Reactors Cylindrical Reactor Absolute Value of Flux Maximum- to-average Flux Ratio Comments on Criticality Condition Multiplying Media.. 1

Media Containing Fuel Neutron propagation considered in passive media till now Characterised by Σ a, Σ t (or Σ s ), External neutron source(s) For a medium containing fissile material, one has a source per unit volume: arising from fissions Here, Q f is not a part of the data provided, but rather is dependent on the flux, i.e. on the sought solution itself The neutronics description yields one or several homogeneous equations Stationary solution exists only if a certain criticality condition is satisfied by the geometrical and material characteristics of the system Multiplying Media.. 2

Media Containing Fuel (contd.) For the critical state: productions = absorptions + leakage Condition ~ independent of the flux (all three reaction rates vary in about the same proportions) We will first consider a reactor which is homogeneous and bare Isolated zone, without a reflector or a blanket (i.e. just the core ) Treatment will first be one-group Monoenergetic neutrons (e.g. thermal with appropriately averaged cross-sections) Later will be considered Multizone reactors (e.g. with a reflector, i.e. an outer zone of pure moderator) Multigroup theory (basis used, in most practical cases, for numerical calculations) Certain aspects of the heterogeneity of the reactor lattice (heterogeneous unit-cell: fuel/moderator/ ) Multiplying Media.. 3

Multiplication Factors For the infinite medium, in general Multiplying Media.. 4

Multiplication Factors (contd.) η c fuel multiplication factor f utilisation factor (for the homogeneous mixture) Laboratory for Reactor Physics and Systems Behaviour For a heterogeneous unit-cell, one needs to consider that the flux is depressed within the fuel region Multiplying Media.. 5

System of Finite Size For the infinite medium: k = η c. f This corresponds to k = η c. f. ε. p, with ε = p = 1 All the neutrons have the same energy There are neither resonance absorptions, nor fast fissions For the finite system, one has the effective multiplication factor Multiplying Media.. 6

System of Finite Size (contd.) Laboratory for Reactor Physics and Systems Behaviour Three cases to be considered: k eff = 1 critical reactor (self-sustaining chain reaction, constant flux) k eff > 1 supercritical reactor (divergent reaction, increasing flux) k eff < 1 subcritical reactor (convergent reaction, decreasing flux) k eff most important single parameter for the functioning of a reactor Important to note: k eff < k < η c < ν Multiplying Media.. 7

Equation for the Critical Reactor Laboratory for Reactor Physics and Systems Behaviour One uses the 1-group diffusion equation for the stationary case (critical reactor) Thereby: Thus, with One-group Reactor Equation Eqn. homogeneous Can be solved for a given geometry, Condition to be satisfied can be identified Simplest systems homogeneous spherical reactor, slab reactor Material Buckling (depends only material properties) Multiplying Media.. 8

Spherical Reactor One has (positive sign, cf. passive medium with point source), i.e. For the finite system: 1. Φ at ρ = 0 C = 0 2. Φ = 0 at ρ = R+d = R+ 0.71λ t, i.e. sin {B m (R+d)} = 0 Multiplying Media.. 9

Spherical Reactor (contd.) From the condition sin {B m (R+d)} = 0 : For the critical reactor, only i = 1 is valid Smallest eigenvalue Fundamental Mode Other solutions: higher harmonics Exist only in a subcritical system E.g. near the external source B m = B i = iπ R + d for i = 1, 2, Critical condition for the spherical reactor is thus: Multiplying Media.. 10

Spherical Reactor (contd.2) The flux distribution is: For a given medium (specific values of criticality condition Critical radius: ), R is determined by the Critical mass: Conversely, if the size (R) is fixed, the material properties need to be identified which yield the appropriate (material buckling) E.g. adjust the enrichment, i.e. k Multiplying Media.. 11

Comments For a bare homogeneous reactor, the criticality condition demands an eigenvalue search for the reactor equation: Eigenvalues need to go to zero at the extrapolated outer surface The square of the smallest eigenvalue: B 2 (geometrical buckling) Criticality condition: Flux distribution: proportional to the eigenfunction corresponding to B 2 Absolute level of the flux not yet known A constant ( A, in the example considered) remains undetermined, depends on the power (determined, in turn, by technological constraints ) allows us to determine the criticality condition and the spatial distribution of the flux, but does not fix the latter s absolute value Multiplying Media.. 12

Slab Reactor System infinite in the X, Y directions Height H Flux, function only of z Reactor Equation: General solution: Flux, symmetric relative to plane z = 0 ( dφ/dz = 0 at z = 0 ) C = 0 Multiplying Media.. 13

Slab Reactor (contd.) Thus, Condition Φ = 0 at yields as eigenvalues: with Square of the smallest eigenvalue ( i = 0 ) (geometrical buckling) Flux distribution: Multiplying Media.. 14

Cylindical Reactor For a cylindrical reactor of height H and radius a, Critical Reactor Equation: With the asumption, Multiplying Media.. 15

Cylindical Reactor (contd.) Multiplying Media.. 16

Cylindical Reactor (contd.2) Comments: The criticality condition implies: (i) If H is so small that, even a does not suffice System remains subcritical (ii) Similarly, if a is too small, Even with H, Multiplying Media.. 17

Cylindical Reactor (contd.3) (iii) There are minimal values for H, a (iv) Many combinations of H, a for the critical state - Shaded area corresponds to supercritical states (e.g. power reactor at start-of-cycle ) Multiplying Media.. 18

Absolute Flux value For a spherical reactor, if one neglects d extrapolated radius Constant A is determined by the reactor power watts joules/fission Multiplying Media.. 19

Absolute Flux value (contd.) For a slab reactor with Per cm 2 of the slab, watts/cm 2 For a cylindrical reactor, one can show: Multiplying Media.. 20

( Φ max / Φ ) Ratio The flux distribution determines the power distribution Reactor homogeneous constant Σ f The maximum-to-average ratio, same for flux, power For bare, homogeneous reactors, flux always maximum at centre, varies strongly (goin to zero at extrapolated surface) In practice, one has a reflector and/ or several different zones in the reactor, which render the flux distribution flatter Multiplying Media.. 21

( Φ max / Φ ) Ratio (Examples) Laboratory for Reactor Physics and Systems Behaviour Spherical Reactor Slab Reactor Multiplying Media.. 22

Comments on Criticality Condition Laboratory for Reactor Physics and Systems Behaviour (a) Rewriting above equation: (1) Previously, one had : (2) From (1) (criticality condition): (3) Comparing (2), (3): Multiplying Media.. 23

Comments on Criticality Condition (contd.) (b) In general, for given values of k, B 2 : (system either subctritical or supercritical) One may nevertheless use the formalism of a critical system (stationary flux) - We consider a fictitious medium in which the number of n s produced per fission For this medium, yields (criticality condition) For the actual system, one needs to search for the modification which leads to. - Change in B 2 or in k. For a bare homogeneous reactor, the search is direct For a multizone system (heterogeneous layout), one needs to adopt an iterative approach - The dimensions (or the material characteristics) are varied until. Multiplying Media.. 24

Comments on Criticality Condition (contd.) (c) We have Laboratory for Reactor Physics and Systems Behaviour Comparing this with: Leakage term DB 2 is like a supplementary Σ a One may consider a reactor of finite dimensions (with a geometrical buckling of B 2 ), as though it were an infinite medium with a poisoning of (However, reactor equation still needs to be solved for obtaining B 2 ) Multiplying Media.. 25

Summary, Lesson 9 Multiplying media, multiplication factors 1-group, diffusion equation (Reactor Eq.) Material and geometrical buckling Bare homogeneous reactors (sphere, slab, cylinder, etc.) Absolute flux and reactor power Maximum- to-average flux ratio Comments on criticality condition Non-leakage probability, criticality search, leakage as absorptions (DB 2 ) Multiplying Media.. 26