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

Size: px
Start display at page:

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

Transcription

1 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

2 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

3 ,, 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

4 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

5 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

6 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

7 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 , (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 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

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

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS http://www.paper.edu.c Iteratioal Joural of Bifurcatio ad Chaos, Vol. 1, No. 5 () 119 15 c World Scietific Publishig Compay AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC

More information

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

2. Neutronic calculations at uranium powered cylindrical reactor by using Bessel differential equation Trasworld Research Network 37/661 (), Fort P.O. Trivadrum-695 03 Kerala, Idia Nuclear Sciece ad Techology, 01: 15-4 ISBN: 978-81-7895-546-9 Editor: Turgay Korkut. Neutroic calculatios at uraium powered

More information

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

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

Appendix: The Laplace Transform

Appendix: The Laplace Transform Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio

More information

Kinetics of Complex Reactions

Kinetics of Complex Reactions Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet

More information

Information-based Feature Selection

Information-based Feature Selection Iformatio-based Feature Selectio Farza Faria, Abbas Kazeroui, Afshi Babveyh Email: {faria,abbask,afshib}@staford.edu 1 Itroductio Feature selectio is a topic of great iterest i applicatios dealig with

More information

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

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would

More information

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

ChE 471 Lecture 10 Fall 2005 SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS SAFE OPERATION OF TUBULAR (PFR) ADIABATIC REACTORS I a exothermic reactio the temperature will cotiue to rise as oe moves alog a plug flow reactor util all of the limitig reactat is exhausted. Schematically

More information

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

Mechatronics. Time Response & Frequency Response 2 nd -Order Dynamic System 2-Pole, Low-Pass, Active Filter Time Respose & Frequecy Respose d -Order Dyamic System -Pole, Low-Pass, Active Filter R 4 R 7 C 5 e i R 1 C R 3 - + R 6 - + e out Assigmet: Perform a Complete Dyamic System Ivestigatio of the Two-Pole,

More information

ADVANCED SOFTWARE ENGINEERING

ADVANCED SOFTWARE ENGINEERING ADVANCED SOFTWARE ENGINEERING COMP 3705 Exercise Usage-based Testig ad Reliability Versio 1.0-040406 Departmet of Computer Ssciece Sada Narayaappa, Aeliese Adrews Versio 1.1-050405 Departmet of Commuicatio

More information

Infinite Sequences and Series

Infinite Sequences and Series Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet

More information

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

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series

More information

Numerical Methods in Fourier Series Applications

Numerical Methods in Fourier Series Applications Numerical Methods i Fourier Series Applicatios Recall that the basic relatios i usig the Trigoometric Fourier Series represetatio were give by f ( x) a o ( a x cos b x si ) () where the Fourier coefficiets

More information

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

The Method of Least Squares. To understand least squares fitting of data. The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve

More information

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

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows

More information

EECE 301 Signals & Systems

EECE 301 Signals & Systems EECE 301 Sigals & Systems Prof. Mark Fowler Note Set #8 D-T Covolutio: The Tool for Fidig the Zero-State Respose Readig Assigmet: Sectio 2.1-2.2 of Kame ad Heck 1/14 Course Flow Diagram The arrows here

More information

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

Addition: Property Name Property Description Examples. a+b = b+a. a+(b+c) = (a+b)+c Notes for March 31 Fields: A field is a set of umbers with two (biary) operatios (usually called additio [+] ad multiplicatio [ ]) such that the followig properties hold: Additio: Name Descriptio Commutativity

More information

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.

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. .3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(

More information

Some New Iterative Methods for Solving Nonlinear Equations

Some New Iterative Methods for Solving Nonlinear Equations World Applied Scieces Joural 0 (6): 870-874, 01 ISSN 1818-495 IDOSI Publicatios, 01 DOI: 10.589/idosi.wasj.01.0.06.830 Some New Iterative Methods for Solvig Noliear Equatios Muhammad Aslam Noor, Khalida

More information

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 +

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 + 62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

COMM 602: Digital Signal Processing

COMM 602: Digital Signal Processing COMM 60: Digital Sigal Processig Lecture 4 -Properties of LTIS Usig Z-Trasform -Iverse Z-Trasform Properties of LTIS Usig Z-Trasform Properties of LTIS Usig Z-Trasform -ve +ve Properties of LTIS Usig Z-Trasform

More information

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE UPB Sci Bull, Series A, Vol 79, Iss, 207 ISSN 22-7027 NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE Gabriel Bercu We itroduce two ew sequeces of Euler-Mascheroi type which have fast covergece

More information

Time-Domain Representations of LTI Systems

Time-Domain Representations of LTI Systems 2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable

More information

Ma 530 Introduction to Power Series

Ma 530 Introduction to Power Series Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power

More information

Chapter 6 Infinite Series

Chapter 6 Infinite Series Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat

More information

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig

More information

Chapter 2: Numerical Methods

Chapter 2: Numerical Methods Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,

More information

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

Classification of problem & problem solving strategies. classification of time complexities (linear, logarithmic etc) Classificatio of problem & problem solvig strategies classificatio of time complexities (liear, arithmic etc) Problem subdivisio Divide ad Coquer strategy. Asymptotic otatios, lower boud ad upper boud:

More information

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

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET Ray Optics Theory ad Mode Theory Dr. Mohammad Faisal Dept. of, BUT Optical Fiber WG For light to be trasmitted through fiber core, i.e., for total iteral reflectio i medium, > Ray Theory Trasmissio Ray

More information

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

1. Linearization of a nonlinear system given in the form of a system of ordinary differential equations . Liearizatio of a oliear system give i the form of a system of ordiary differetial equatios We ow show how to determie a liear model which approximates the behavior of a time-ivariat oliear system i a

More information

Recurrence Relations

Recurrence Relations Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The

More information

ENGI Series Page 6-01

ENGI Series Page 6-01 ENGI 3425 6 Series Page 6-01 6. Series Cotets: 6.01 Sequeces; geeral term, limits, covergece 6.02 Series; summatio otatio, covergece, divergece test 6.03 Stadard Series; telescopig series, geometric series,

More information

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

1 6 = 1 6 = + Factorials and Euler s Gamma function Royal Holloway Uiversity of Lodo Departmet of Physics Factorials ad Euler s Gamma fuctio Itroductio The is a self-cotaied part of the course dealig, essetially, with the factorial fuctio ad its geeralizatio

More information

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

Higher-order iterative methods by using Householder's method for solving certain nonlinear equations Math Sci Lett, No, 7- ( 7 Mathematical Sciece Letters A Iteratioal Joural http://dxdoiorg/785/msl/5 Higher-order iterative methods by usig Householder's method for solvig certai oliear equatios Waseem

More information

ABSTRACT 1. INTRODUCTION

ABSTRACT 1. INTRODUCTION Iwamura, Y., T. Itoh, ad M. SakaNo. Nuclear Products ad Their Time Depedece Iduced by Cotiuous Diffusio of Deuterium Through Multi-layer Palladium Cotaiig Low Work Fuctio Material. i 8th Iteratioal Coferece

More information

CHAPTER 10 INFINITE SEQUENCES AND SERIES

CHAPTER 10 INFINITE SEQUENCES AND SERIES CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece

More information

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

Temperature Distribution in a Cylindrical Core and Coil Assembly with Heat Generation Temperature Distributio i a Cylidrical Core ad Coil Assembly with Heat Geeratio P. S. Yeh, Ph.D. Abstract A residetial distributio trasformer cosists of a cylidrical steel tak ad a core-ad-coil assembly

More information

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS

A NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet

More information

Orthogonal Gaussian Filters for Signal Processing

Orthogonal Gaussian Filters for Signal Processing Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios

More information

Notes on iteration and Newton s method. Iteration

Notes on iteration and Newton s method. Iteration Notes o iteratio ad Newto s method Iteratio Iteratio meas doig somethig over ad over. I our cotet, a iteratio is a sequece of umbers, vectors, fuctios, etc. geerated by a iteratio rule of the type 1 f

More information

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

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0, 4. Series Solutios of Differetial Equatios:Special Fuctios 4.. Illustrative examples.. 5. Obtai the geeral solutio of each of the followig differetial equatios i terms of Maclauri series: d y (a dx = xy,

More information

Some properties of Boubaker polynomials and applications

Some properties of Boubaker polynomials and applications Some properties of Boubaker polyomials ad applicatios Gradimir V. Milovaović ad Duša Joksimović Citatio: AIP Cof. Proc. 179, 1050 (2012); doi: 10.1063/1.756326 View olie: http://dx.doi.org/10.1063/1.756326

More information

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

CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS C.PRAX ad H.SADAT Laboratoire d'etudes Thermiques,URA CNRS 403 40, Aveue du Recteur Pieau 86022 Poitiers Cedex,

More information

1 Approximating Integrals using Taylor Polynomials

1 Approximating Integrals using Taylor Polynomials Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................

More information

8. Applications To Linear Differential Equations

8. Applications To Linear Differential Equations 8. Applicatios To Liear Differetial Equatios 8.. Itroductio 8.. Review Of Results Cocerig Liear Differetial Equatios Of First Ad Secod Orders 8.3. Eercises 8.4. Liear Differetial Equatios Of Order N 8.5.

More information

Design and Analysis of Algorithms

Design and Analysis of Algorithms Desig ad Aalysis of Algorithms Probabilistic aalysis ad Radomized algorithms Referece: CLRS Chapter 5 Topics: Hirig problem Idicatio radom variables Radomized algorithms Huo Hogwei 1 The hirig problem

More information

Chapter 10: Power Series

Chapter 10: Power Series Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because

More information

Chapter 5 Vibrational Motion

Chapter 5 Vibrational Motion Fall 4 Chapter 5 Vibratioal Motio... 65 Potetial Eergy Surfaces, Rotatios ad Vibratios... 65 Harmoic Oscillator... 67 Geeral Solutio for H.O.: Operator Techique... 68 Vibratioal Selectio Rules... 7 Polyatomic

More information

NUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK

NUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK NUMERICAL METHODS COURSEWORK INFORMAL NOTES ON NUMERICAL INTEGRATION COURSEWORK For this piece of coursework studets must use the methods for umerical itegratio they meet i the Numerical Methods module

More information

6.003 Homework #3 Solutions

6.003 Homework #3 Solutions 6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the

More information

Physics Oct Reading

Physics Oct Reading Physics 301 21-Oct-2002 17-1 Readig Fiish K&K chapter 7 ad start o chapter 8. Also, I m passig out several Physics Today articles. The first is by Graham P. Collis, August, 1995, vol. 48, o. 8, p. 17,

More information

Massachusetts Institute of Technology

Massachusetts Institute of Technology 6.0/6.3: Probabilistic Systems Aalysis (Fall 00) Problem Set 8: Solutios. (a) We cosider a Markov chai with states 0,,, 3,, 5, where state i idicates that there are i shoes available at the frot door i

More information

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

A numerical Technique Finite Volume Method for Solving Diffusion 2D Problem The Iteratioal Joural Of Egieerig d Sciece (IJES) Volume 4 Issue 10 Pages PP -35-41 2015 ISSN (e): 2319 1813 ISSN (p): 2319 1805 umerical Techique Fiite Volume Method for Solvig Diffusio 2D Problem 1 Mohammed

More information

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

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001. Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;

More information

Complex Analysis Spring 2001 Homework I Solution

Complex Analysis Spring 2001 Homework I Solution Complex Aalysis Sprig 2001 Homework I Solutio 1. Coway, Chapter 1, sectio 3, problem 3. Describe the set of poits satisfyig the equatio z a z + a = 2c, where c > 0 ad a R. To begi, we see from the triagle

More information

Castiel, Supernatural, Season 6, Episode 18

Castiel, Supernatural, Season 6, Episode 18 13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio

More information

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

Proof of Fermat s Last Theorem by Algebra Identities and Linear Algebra Proof of Fermat s Last Theorem by Algebra Idetities ad Liear Algebra Javad Babaee Ragai Youg Researchers ad Elite Club, Qaemshahr Brach, Islamic Azad Uiversity, Qaemshahr, Ira Departmet of Civil Egieerig,

More information

SINGLE-CHANNEL QUEUING PROBLEMS APPROACH

SINGLE-CHANNEL QUEUING PROBLEMS APPROACH SINGLE-CHANNEL QUEUING ROBLEMS AROACH Abdurrzzag TAMTAM, Doctoral Degree rogramme () Dept. of Telecommuicatios, FEEC, BUT E-mail: xtamta@stud.feec.vutbr.cz Supervised by: Dr. Karol Molár ABSTRACT The paper

More information

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

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan Mathematical ad Computatioal Applicatios, Vol. 9, No. 3, pp. 30-40, 04 DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS Muhammad Aslam Noor, Khalida Iayat Noor ad Asif Waheed

More information

MATH 129 FINAL EXAM REVIEW PACKET (Revised Spring 2008)

MATH 129 FINAL EXAM REVIEW PACKET (Revised Spring 2008) MATH 9 FINAL EXAM REVIEW PACKET (Revised Sprig 8) The followig questios ca be used as a review for Math 9. These questios are ot actual samples of questios that will appear o the fial exam, but they will

More information

Chapter 2 Feedback Control Theory Continued

Chapter 2 Feedback Control Theory Continued Chapter Feedback Cotrol Theor Cotiued. Itroductio I the previous chapter, the respose characteristic of simple first ad secod order trasfer fuctios were studied. It was show that first order trasfer fuctio,

More information

Signals & Systems Chapter3

Signals & Systems Chapter3 Sigals & Systems Chapter3 1.2 Discrete-Time (D-T) Sigals Electroic systems do most of the processig of a sigal usig a computer. A computer ca t directly process a C-T sigal but istead eeds a stream of

More information

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

Taylor polynomial solution of difference equation with constant coefficients via time scales calculus TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu

More information

Chapter 6 Sampling Distributions

Chapter 6 Sampling Distributions Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to

More information

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

Chapter 6 Overview: Sequences and Numerical Series. For the purposes of AP, this topic is broken into four basic subtopics: Chapter 6 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals (which is what most studets

More information

TEACHER CERTIFICATION STUDY GUIDE

TEACHER CERTIFICATION STUDY GUIDE COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra

More information

Analytic Theory of Probabilities

Analytic Theory of Probabilities Aalytic Theory of Probabilities PS Laplace Book II Chapter II, 4 pp 94 03 4 A lottery beig composed of umbered tickets of which r exit at each drawig, oe requires the probability that after i drawigs all

More information

Sequences of Definite Integrals, Factorials and Double Factorials

Sequences of Definite Integrals, Factorials and Double Factorials 47 6 Joural of Iteger Sequeces, Vol. 8 (5), Article 5.4.6 Sequeces of Defiite Itegrals, Factorials ad Double Factorials Thierry Daa-Picard Departmet of Applied Mathematics Jerusalem College of Techology

More information

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

Thermal characterization of an insulating material through a tri-layer transient method Advaced Sprig School «Thermal Measuremets & Iverse techiques», Domaie de Fraço, Biarritz, March -6 05 T6 Thermal characterizatio of a isulatig material through a tri-layer trasiet method V. Félix, Y. Jaot,

More information

SECTION 2 Electrostatics

SECTION 2 Electrostatics SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial

More information

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

Monte Carlo Optimization to Solve a Two-Dimensional Inverse Heat Conduction Problem Australia Joural of Basic Applied Scieces, 5(): 097-05, 0 ISSN 99-878 Mote Carlo Optimizatio to Solve a Two-Dimesioal Iverse Heat Coductio Problem M Ebrahimi Departmet of Mathematics, Karaj Brach, Islamic

More information

Chapter 7: Numerical Series

Chapter 7: Numerical Series Chapter 7: Numerical Series Chapter 7 Overview: Sequeces ad Numerical Series I most texts, the topic of sequeces ad series appears, at first, to be a side topic. There are almost o derivatives or itegrals

More information

ON POINTWISE BINOMIAL APPROXIMATION

ON POINTWISE BINOMIAL APPROXIMATION Iteratioal Joural of Pure ad Applied Mathematics Volume 71 No. 1 2011, 57-66 ON POINTWISE BINOMIAL APPROXIMATION BY w-functions K. Teerapabolar 1, P. Wogkasem 2 Departmet of Mathematics Faculty of Sciece

More information

MOMENT-METHOD ESTIMATION BASED ON CENSORED SAMPLE

MOMENT-METHOD ESTIMATION BASED ON CENSORED SAMPLE Vol. 8 o. Joural of Systems Sciece ad Complexity Apr., 5 MOMET-METHOD ESTIMATIO BASED O CESORED SAMPLE I Zhogxi Departmet of Mathematics, East Chia Uiversity of Sciece ad Techology, Shaghai 37, Chia. Email:

More information

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

Discrete Orthogonal Moment Features Using Chebyshev Polynomials Discrete Orthogoal Momet Features Usig Chebyshev Polyomials R. Mukuda, 1 S.H.Og ad P.A. Lee 3 1 Faculty of Iformatio Sciece ad Techology, Multimedia Uiversity 75450 Malacca, Malaysia. Istitute of Mathematical

More information

: Transforms and Partial Differential Equations

: Transforms and Partial Differential Equations Trasforms ad Partial Differetial Equatios 018 SUBJECT NAME : Trasforms ad Partial Differetial Equatios SUBJECT CODE : MA 6351 MATERIAL NAME : Part A questios REGULATION : R013 WEBSITE : wwwharigaeshcom

More information

PAPER : IIT-JAM 2010

PAPER : IIT-JAM 2010 MATHEMATICS-MA (CODE A) Q.-Q.5: Oly oe optio is correct for each questio. Each questio carries (+6) marks for correct aswer ad ( ) marks for icorrect aswer.. Which of the followig coditios does NOT esure

More information

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

Helium-like Ions Below. the N=2-9 Hydrogenic Thresholds Adv. Studies Theor. Phys., Vol. 6, 01, o. 7, 1345-1354 The 1,3 S e Itershell Autoioizig States of Helium-like Ios Below the N=-9 Hydrogeic Thresholds M. Dieg *, A. Wagué Laboratory Atoms Lasers, Departmet

More information

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

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy Liear Differetial Equatios of Higher Order Basic Theory: Iitial-Value Problems d y d y dy Solve: a( ) + a ( )... a ( ) a0( ) y g( ) + + + = d d d ( ) Subject to: y( 0) = y0, y ( 0) = y,..., y ( 0) = y

More information

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

CS 270 Algorithms. Oliver Kullmann. Growth of Functions. Divide-and- Conquer Min-Max- Problem. Tutorial. Reading from CLRS for week 2 Geeral remarks Week 2 1 Divide ad First we cosider a importat tool for the aalysis of algorithms: Big-Oh. The we itroduce a importat algorithmic paradigm:. We coclude by presetig ad aalysig two examples.

More information

THE KALMAN FILTER RAUL ROJAS

THE KALMAN FILTER RAUL ROJAS THE KALMAN FILTER RAUL ROJAS Abstract. This paper provides a getle itroductio to the Kalma filter, a umerical method that ca be used for sesor fusio or for calculatio of trajectories. First, we cosider

More information

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

*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully X0/70 NATIONAL QUALIFICATIONS 006 MONDAY, MAY.00 PM.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.

More information

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

THE APPEARANCE OF FIBONACCI AND LUCAS NUMBERS IN THE SIMULATION OF ELECTRICAL POWER LINES SUPPLIED BY TWO SIDES THE APPEARANCE OF FIBONACCI AND LUCAS NUMBERS IN THE SIMULATION OF ELECTRICAL POWER LINES SUPPLIED BY TWO SIDES Giuseppe Ferri Dipartimeto di Ihgegeria Elettrica-FacoM di Igegeria, Uiversita di L'Aquila

More information

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

MID-YEAR EXAMINATION 2018 H2 MATHEMATICS 9758/01. Paper 1 JUNE 2018 MID-YEAR EXAMINATION 08 H MATHEMATICS 9758/0 Paper JUNE 08 Additioal Materials: Writig Paper, MF6 Duratio: hours DO NOT OPEN THIS BOOKLET UNTIL YOU ARE TOLD TO DO SO READ THESE INSTRUCTIONS FIRST Write

More information

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

The optimal online control of the instantaneous power and the multiphase source s current BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES, Vol. 65, No. 6, 27 DOI:.55/bpasts-27-9 The optimal olie cotrol of the istataeous power ad the multiphase source s curret M. SIWCZYŃSKI ad

More information

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

Exact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method Exact Solutios for a Class of Noliear Sigular Two-Poit Boudary Value Problems: The Decompositio Method Abd Elhalim Ebaid Departmet of Mathematics, Faculty of Sciece, Tabuk Uiversity, P O Box 741, Tabuki

More information

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

Numerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem Leoardo Joural of Scieces ISSN 1583-0233 Issue 9, July-December 2006 p. 25-32 Numerical Simulatio of Thermomechaical Problems i Applied Mechaics: Applicatio to Solidificatio Problem Vicet Obiajulu OGWUAGWU

More information

Unit 4: Polynomial and Rational Functions

Unit 4: Polynomial and Rational Functions 48 Uit 4: Polyomial ad Ratioal Fuctios Polyomial Fuctios A polyomial fuctio y px ( ) is a fuctio of the form p( x) ax + a x + a x +... + ax + ax+ a 1 1 1 0 where a, a 1,..., a, a1, a0are real costats ad

More information

μ are complex parameters. Other

μ are complex parameters. Other A New Numerical Itegrator for the Solutio of Iitial Value Problems i Ordiary Differetial Equatios. J. Suday * ad M.R. Odekule Departmet of Mathematical Scieces, Adamawa State Uiversity, Mubi, Nigeria.

More information

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

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at

More information

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

Exact scattering and bound states solutions for novel hyperbolic potentials with inverse square singularity Exact scatterig ad boud states solutios for ovel hyperbolic potetials with iverse square sigularity A. D. Alhaidari Saudi Ceter for Theoretical Physics, P. O. Box 37, Jeddah 38, Saudi Arabia Abstract:

More information

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

( ) ( ), (S3) ( ). (S4) Ultrasesitivity i phosphorylatio-dephosphorylatio cycles with little substrate: Supportig Iformatio Bruo M.C. Martis, eter S. Swai 1. Derivatio of the equatios associated with the mai model From the differetial

More information

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS

POWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity

More information

The Riemann Zeta Function

The Riemann Zeta Function Physics 6A Witer 6 The Riema Zeta Fuctio I this ote, I will sketch some of the mai properties of the Riema zeta fuctio, ζ(x). For x >, we defie ζ(x) =, x >. () x = For x, this sum diverges. However, we

More information

SCORE. Exam 2. MA 114 Exam 2 Fall 2017

SCORE. Exam 2. MA 114 Exam 2 Fall 2017 Exam Name: Sectio ad/or TA: Do ot remove this aswer page you will retur the whole exam. You will be allowed two hours to complete this test. No books or otes may be used. You may use a graphig calculator

More information

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

A New Method to Order Functions by Asymptotic Growth Rates Charlie Obimbo Dept. of Computing and Information Science University of Guelph A New Method to Order Fuctios by Asymptotic Growth Rates Charlie Obimbo Dept. of Computig ad Iformatio Sciece Uiversity of Guelph ABSTRACT A ew method is described to determie the complexity classes of

More information

CS322: Network Analysis. Problem Set 2 - Fall 2009

CS322: Network Analysis. Problem Set 2 - Fall 2009 Due October 9 009 i class CS3: Network Aalysis Problem Set - Fall 009 If you have ay questios regardig the problems set, sed a email to the course assistats: simlac@staford.edu ad peleato@staford.edu.

More information

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

Frank Laboratory of Neutron Physics, JINR, Dubna, Russia. Institute of Heavy Ion Physics, Peking University, Beijing, P.R. China. 1. SYSTEMATICS OF (,α) CROSS SECTIONS FOR 4-6 MeV NEUTRONS G.Khuukhekhuu 1, Yu.M.Gledeov, Guohui Zhag 3, M.V.Sedysheva, J.Mukhsaikha 1, M.Odsure 1 1) Nuclear Research Ceter, Natioal Uiversity of Mogolia,

More information

Chapter 4 : Laplace Transform

Chapter 4 : Laplace Transform 4. Itroductio Laplace trasform is a alterative to solve the differetial equatio by the complex frequecy domai ( s = σ + jω), istead of the usual time domai. The DE ca be easily trasformed ito a algebraic

More information