REAL-TIME MONITORING OF POWER AND NEUTRON CAPTURE CROSS SECTION IN NUCLEAR RESEARCH REACTOR

Similar documents
AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

2. Neutronic calculations at uranium powered cylindrical reactor by using Bessel differential equation

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Appendix: The Laplace Transform

Kinetics of Complex Reactions

Information-based Feature Selection

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

ChE 471 Lecture 10 Fall 2005 SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter

ADVANCED SOFTWARE ENGINEERING

Infinite Sequences and Series

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Numerical Methods in Fourier Series Applications

The Method of Least Squares. To understand least squares fitting of data.

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

EECE 301 Signals & Systems

Addition: Property Name Property Description Examples. a+b = b+a. a+(b+c) = (a+b)+c

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

Some New Iterative Methods for Solving Nonlinear Equations

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

6.3 Testing Series With Positive Terms

COMM 602: Digital Signal Processing

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

Time-Domain Representations of LTI Systems

Ma 530 Introduction to Power Series

Chapter 6 Infinite Series

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Chapter 2: Numerical Methods

Classification of problem & problem solving strategies. classification of time complexities (linear, logarithmic etc)

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations

Recurrence Relations

ENGI Series Page 6-01

1 6 = 1 6 = + Factorials and Euler s Gamma function

Higher-order iterative methods by using Householder's method for solving certain nonlinear equations

ABSTRACT 1. INTRODUCTION

CHAPTER 10 INFINITE SEQUENCES AND SERIES

Temperature Distribution in a Cylindrical Core and Coil Assembly with Heat Generation

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

Orthogonal Gaussian Filters for Signal Processing

Notes on iteration and Newton s method. Iteration

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,

Some properties of Boubaker polynomials and applications

CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS

1 Approximating Integrals using Taylor Polynomials

8. Applications To Linear Differential Equations

Design and Analysis of Algorithms

Chapter 10: Power Series

Chapter 5 Vibrational Motion

NUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK

6.003 Homework #3 Solutions

Physics Oct Reading

Massachusetts Institute of Technology

A numerical Technique Finite Volume Method for Solving Diffusion 2D Problem

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

Complex Analysis Spring 2001 Homework I Solution

Castiel, Supernatural, Season 6, Episode 18

Proof of Fermat s Last Theorem by Algebra Identities and Linear Algebra

SINGLE-CHANNEL QUEUING PROBLEMS APPROACH

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

MATH 129 FINAL EXAM REVIEW PACKET (Revised Spring 2008)

Chapter 2 Feedback Control Theory Continued

Signals & Systems Chapter3

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus

Chapter 6 Sampling Distributions

Chapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics:

TEACHER CERTIFICATION STUDY GUIDE

Analytic Theory of Probabilities

Sequences of Definite Integrals, Factorials and Double Factorials

Thermal characterization of an insulating material through a tri-layer transient method

SECTION 2 Electrostatics

Monte Carlo Optimization to Solve a Two-Dimensional Inverse Heat Conduction Problem

Chapter 7: Numerical Series

ON POINTWISE BINOMIAL APPROXIMATION

MOMENT-METHOD ESTIMATION BASED ON CENSORED SAMPLE

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

: Transforms and Partial Differential Equations

PAPER : IIT-JAM 2010

Helium-like Ions Below. the N=2-9 Hydrogenic Thresholds

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy

CS 270 Algorithms. Oliver Kullmann. Growth of Functions. Divide-and- Conquer Min-Max- Problem. Tutorial. Reading from CLRS for week 2

THE KALMAN FILTER RAUL ROJAS

*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully

THE APPEARANCE OF FIBONACCI AND LUCAS NUMBERS IN THE SIMULATION OF ELECTRICAL POWER LINES SUPPLIED BY TWO SIDES

MID-YEAR EXAMINATION 2018 H2 MATHEMATICS 9758/01. Paper 1 JUNE 2018

The optimal online control of the instantaneous power and the multiphase source s current

Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method

Numerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem

Unit 4: Polynomial and Rational Functions

μ are complex parameters. Other

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

Exact scattering and bound states solutions for novel hyperbolic potentials with inverse square singularity

( ) ( ), (S3) ( ). (S4)

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

The Riemann Zeta Function

SCORE. Exam 2. MA 114 Exam 2 Fall 2017

A New Method to Order Functions by Asymptotic Growth Rates Charlie Obimbo Dept. of Computing and Information Science University of Guelph

CS322: Network Analysis. Problem Set 2 - Fall 2009

Frank Laboratory of Neutron Physics, JINR, Dubna, Russia. Institute of Heavy Ion Physics, Peking University, Beijing, P.R. China. 1.

Chapter 4 : Laplace Transform

Transcription:

D5 REAL-TIME MONITORING OF POWER AND NEUTRON CAPTURE CROSS SECTION IN NUCLEAR RESEARCH REACTOR D.A.P. PALMA Brazilia Nuclear Eergy Commissio (CNEN), Rio de Jaeiro, Brazil A.Z. MESQUITA, R.M.G.P. SOUZA Nuclear Techology Developmet Cetre (CDTN), Belo Horizote, Brazil A.S. MARTINEZ Nuclear Eergy Program, PEN COPPE/UFRJ, Rio de Jaeiro, Brazil Abstract Nuclear reactor operators have to moitor the behaviour of differet uclear ad desig parameters that vary i time so to esure the safety operatio of the reactor. I recet years several operatig parameters for the IPR-R1 TRIGA research reactor were moitored ad idicated i real-time by the data acquisitio system developed for the reactor, with all the data beig stored i a hard drive i the data acquisitio computer, to build i this way a database. The goal of this work is to isert i the set of parameters already collected the eutro capture cross sectios for the fuel, from the uclear power obtaied i real-time. The experimetal data was obtaied by usig a fuel elemet istrumeted with temperature sesors, located i the core of the IPR-R1 TRIGA research reactor at the Nuclear Techology Developmet Cetre (CDTN). This iformatio is useful for the cotiuous moitorig of the eutro capture reactio rate. For that, a ew aalytical formulatio is used for the Doppler broadeig fuctio proposed by Palma ad Martiez, which is free from special fuctios i its fuctioal form ad with easy computig implemetatio. The results obtaied were satisfactory from the stadpoit of accuracy i compariso with the umerical referece model ad idicate that it is possible to carry out real-time moitorig of the eutro capture cross sectio i the fuel. 1. INTRODUCTION The movemet of thermal agitatio of the uclei is adequately represeted i the microscopic cross sectios of eutro-ucleus iteractio through the Doppler broadeig fuctio ad of iterferece term. This fuctio is calculated umerically i moder systems for the calculatio of the macrogroup costats, eeded to determie the power distributio of a uclear reactor ad is writte by: beig: dy 1 y 4 (1) x, exp x y E E x () Iteratioal Atomic Eergy Agecy 1

4 E kt / A 1/ D.A.P. Palma et al., (3) where E is the eergy of the icidet eutro, E the eergy of the resoace peak, the total width of the resoace, T the temperature of the medium ad k the Boltzma costat. Cosiderig a medium thermal equilibrium at temperature T where the target uclei are i movemet ad their velocities give by the Maxwell-Boltzma distributio ad the expressios for effective absorptio cross sectios are obtaied from the Briet-Wiger formalism [1], E E, T x,. E (4) Due the itegral form i the equatio (1), the equatio (4) is ot adequate to a real time moitorig of the absorptio cross sectio.. THE IPR-R1 TRIGA REACTOR The IPR-R1 TRIGA (Traiig, Research, Isotopes, ad Geeral Atomic), showed i Fig. 1 is a pool type uclear research reactor, with a ope water surface ad the core has a cylidrical cofiguratio []. The maximum core power is 5 kw, cooled by light water ad with graphite reflectors. The fuel is a alloy of zircoium hydride ad uraium eriched at % i 35 U. The reactor core has 63 cylidrical fuel elemets, 58 alumium- clad fuel elemets, ad 5 stailess steel-clad fuel elemets with 8.5 wt % uraium. FIG. 1. IPR-R1 TRIGA reactor pool ad core. I the ext sectio it will be preseted a ew aalytical approximatio for the Doppler x, that will allow your calculatio fast eough to make a real time broadeig fuctio moitorig of the ET, from the data of temperature ad power comig from istrumeted fuel. 3. AN ANALYTICAL APPROXIMATION FOR THE DOPPLER BROADENING x, FUNCTION The Doppler broadeig fuctio x, equatio [3], ca be defied as the solutio of differetial

,, D.A.P. Palma et al. 4 d x d x x x x,, dx dx 4 4 (5) subjected to the iitial coditios, x, x exp 1 erf 4 x, x. x (6) A aalytical expressio for the Doppler broadeig fuctio ca be obtaied by solvig the equatio (5), usig a cojuctio betwee the Frobeius method to obtai the homogeeous solutio ad the method of parameter variatio to obtai a particular solutio. The Frobeius method cosists basically of seekig solutios i the form of: h s x, c x, with c, (7) where s is a o-determied parameter which provides the method with flexibility. The replacemet of equatio (7) ad its derived elemets i equatio (11), oe obtais the followig solutio for the homogeeous equatio: h x x x x, exp k1cos ksi, 4 (8) where k ad k are costats to be determied after the impositio of the iitial coditios expressed by equatio (6). Havig obtaied the liearly idepedet solutios of the homogeeous equatio associated to Eq. (11), x x 1 x, exp cos 4 x x x, exp si, 4 (9) (1) oe ca apply the method of parameter variatio i obtaiig a particular solutio for Eq. (5), seekig solutios as: where,,, (11) p x u1 x 1 x u x x u1 x ad u x are obtaied from the coditios expressed by equatio (6). I doig so oe has the followig aalytical expressio for the Doppler broadeig fuctio: x x x, cos exp x 11 Re x, ta Im x,, 4 (1) 3

ix where x, erf. The equatio (1) provides a simple ad precise expressio for the Doppler broadeig fuctio but is possible to separate the error fuctio with complex argumet that exists i the i its real ad imagiary parts accordig to the fuctioal form of fuctio Frob x, expasio proposed by Salzer [3]. Accordig to this approximatio, the error fuctio with complex argumet ca be expaded i ifiite series: where, Re. Im., erf a bi F F i (13) a 1 cos( ab) max exp exp / 4 FRe. erf a f a, b, a 1 4a a si( ab) max exp exp / 4 FIm. g a, b, a 1 4a ad the auxiliary fuctios (14) (15) f a, b a acosh b cos ab sih b si ab, (16) g a, b acosh b si ab sih b cos ab. (17) Deotig a ad b ax ad imagiary parts: it is possible to separate fuctio x, ito its real max 1 x exp / 4 Re x, erf exp cos 1 f x,, Ap. 4 1 (18) max 1 x exp / 4 Im x, exp si g x,, Ap. 4 1 where the auxiliary fuctios f x, ad g, (19) x are re-writte by: x x x x f x, cosh cos sih si, x x x x g x, cosh si sih cos. () (1) 4

Whe replacig equatios (18) ad (19) i the expressio of the Doppler broadeig fuctio (x,ξ), equatio (1), it is possible to write the followig expressio for the Doppler broadeig fuctio: 1 x x x, exp x 1 cos 1 Re x, ta Im x, Ap.. () Ap. 4 x,, give by equatio (), is a valid expressio for all the values of Expressio for parameters x ad. I the preset paper N max =1 was used. 4. RESULTS 4.1 Thermal power moitorig Experimetal studies have bee performed i the IPR-R1 reactor by to fid out the o-lie moitorig of the core thermal power [4] ad [5]. Oe method uses the temperature differece betwee a istrumeted fuel elemet with thermocouples, ad the pool water below the reactor core. Aother method cosists of the steady-state eergy balace of the primary ad secodary reactor coolig loops. A further method is by calorimetric procedure whereby the reactor power is moitored as a fuctio of the temperature-rise rate ad the system heat capacity costat. Figure shows oe of the video-scree displays of the digital moitorig system computer that cosolidates iformatio for the reactor power status i real time. This scree moitors the power measured by the eutroic chaels ad by the ew thermal chaels. 4.. Capture cross sectio FIG. Power moitorig o the scree of the data acquisitio system. The relatio betwee the temperature i the uclear fuel as a fuctio of the power durig the operatio of the IPR-R1 TRIGA reactor ca be visualized i Fig. 3. 5

FIG. 3. Temperature i the uclear fuel as a fuctio of the power durig the operatio of the IPR-R1 TRIGA This iformatio of the temperature is useful for the cotiuous moitorig of the reactio rate i eutro capture. For that, a ew aalytical formulatio is used for the Doppler broadeig fuctio proposed by Palma ad Martiez to calculate the absorptio cross sectio that ca be visualized i the figure 4. FIG. 4. Capture cross sectio calculated from the equatio (4), i fuctio of the power measured ad the eergy of the icidet eutro. 5. CONCLUSIONS Two importat criteria for operatioal parameters measuremet i uclear reactors are redudacy ad diversity. Other criteria such as accuracy, reliability ad speed i respose are also of major cocer. The greater the umber of chaels for measurig power, the greater is the reliability ad safety of reactor operatios. The aim of this research is to develop ew methodologies for o-lie moitorig of uclear reactor power ad isert i the set of parameters already collected the eutro 6

capture cross sectios for the fuel, from the power obtaied olie. The system was implemeted, with success, i the IPR-R1 TRIGA research reactor at the Nuclear Techology Developmet Cetre (CDTN). ACKNOWLEDGEMENTS These experimets are part of a research project supported by the Brazilia Coucil for Scietific ad Techological Developmet (CNPq) ad the Research Support Foudatio of the State of Mias Gerais (FAPEMIG). REFERENCES [1] DUDERSTADT, J.J., HAMILTON, L. J. Nuclear Reactor Aalysis, Ed. Joh Wiley ad Sos, Ic. Nova York, 65 p. (1976). [] MESQUITA, A.Z. Experimetal Heat Trasfer Aalysis of the IPR-R1 TRIGA Reactor. Iteratioal Atomic Eergy Agecy (IAEA-CN-156). IAEA Publicatio - Research Reactors: Safe Maagemet ad Effective Utilizatio, v. IAEA, pp. 1-1, 8. x, [3] PALMA, D. A. P., Martiez, A.S. A Faster Procedure for the Calculatio of the Aals of Nuclear Eergy, Vol. 36, pp. 1469-1471, (9). [4] MESQUITA, A.Z.; REZENDE, H.C. Thermal Methods for O-Lie Power Moitorig of the IPR-R1 TRIGA Reactor. Progress i Nuclear Eergy, v. 5, p. 68-7. 1. doi:1.116/j.pucee.9.5.6, (9). [5] MESQUITA, A.Z.; REZENDE, H.C.; SOUZA, R.M.G.P. Thermal Power Calibratios of the IPR-R1 TRIGA Reactor by the Calorimetric ad the Heat Balace Methods. Progress i Nuclear Eergy (i press), doi:1.116/j.pucee, (11). 7