A Method For the Burnup Analysis of Power Reactors in Equilibrium Operation Cycles
|
|
- Logan Todd
- 5 years ago
- Views:
Transcription
1 Journal of NUCLEAR SCIENCE and TECHNOLOGY, 3[5], p.184~188 (May 1966). A Method For the Burnup Analysis of Power Reactors in Equilibrium Operation Cycles Shoichiro NAKAMURA* Received February 7, 1966 This paper describes a method based on the consistent thermal neutron flux concept for burnup analysis of power reactors in equilibrium operation cycle. Radial flux distributions are discussed, as well as the variation of isotopic densities in the fuel, and considerations are given on control rod distribution for maintaining a consistent thermal flux distribution. The present method can be applied to most power reactors like BWR and PWR in which a multiple-batch refueling scheme is adopted. I. INTRODUCTION The main objective of burnup calculation is to determine the characteristics of the core that vary with operation time. Burnup calculation consists of the main parts: Flux distribution calculation with a diffusion code and isotopic density calculation with burnup equation. From the burnup calculations we obtain the variations of flux distribution, the isotopic densities, the excess reactivity, the available operation time and the total energy output produced by a unit volume of fuel. Since the eigenvalue (keff) of the reactor should be maintained at unity, excess reactivity is invariably cancelled by the control rods during full power operation. The control rods having as they do a strong influence on flux distribution, flux calculations in burnup analysis must be performed with consideration given to the effect of the control rods. When the operation cycles of the core are stabilized, refueling is performed periodically. We define gone operation cycle h as the period of operation between one refueling and the next. If an identical process of physical property variation is repeated in each cycle, calculations for only one cycle should suffice. Excess reactivity being large at the beginning of the cycle, a large number of control rods are inserted, and it is easy to modify the flux distribution. The flux distribution decreases its flexibility with operation time, and finally, it becomes quite rigid at the end of the cycle when all the control rods are fully withdrawn. The flux distribution at the end of the cycle is uniquely determined by the history of the flux distribution during the cycle. It has been shown since our previous study(1) that this fact is compatible with burnup calculations based on consistent thermal flux distribution. By gconsistent flux his meant in this instance the regid shape of the thermal flux distribution maintained throughout the cycle. Using the two-group theory, this paper discusses the method based on the consistent flux distribution of radial coordinates. The control rod patterns required for maintaining the consistent thermal flux are also discussed. The present concept not only simplifies the burnup calculations as a whole, but also has another merit in that the consistent thermal flux is the flattest among all the possible flux distributions during the cycle, as it can easily be proved. II. THEORY In the one dimensional two-group diffusion approximation, the flux distributions on radial coordinates satisfy the equation * Nuclear Reactor Design Section, Hitachi Works, Hitachi Ltd., Hitachi-shi, Ibaraki-ken. 14
2 Vol. 3, No. 5 (May 1966) 185 (1) into a number of zones bounded by concentric circles, each numbered (represented by the symbol k) sequentially from center to periphery. The spent fuel is discharged and replaced by fresh fuel in some of these where (2) In the above equations, Sa is the absorption cross section, SIS is the absorption cross section of Xe and Sm. Assuming that the control rods can be approximated by a homogeneous poison, the absorption is represented by Sc. We consider that the variations of the absorption and fission cross sections are functions of exposure, so that zones. Shuffling of partly irradiated fuel is also undertaken. It is assumed that the shuffled fuel is completely mixed and therefore uniform in in each zone. The refueling schemes are represented by the matrix designated as gshuffling matrix h, the element Gkk Œbeing unity when the fuel in the k-th zone has been previously irradiated in the k Œ-th zone, and otherwise zero. If we assume that, the thermal flux where t is the thermal neutron flux time of exposure of the fuel, N the isotopic density, and SMa the absorption cross section of materials other than fuel. The dependence of the fast neutron cross sections on exposure is neglected. The relation of the isotopic densities to is obtained(2) from the burnup equations (3) (4) (5) (6) t distribution is maintained constant throughout, in conformity with the consistent flux concept, the exposure of the fuel at the end of a cycle is given by (12) where g is the position in the k-th zone, T the duration of one cycle, and <P2>k Œthe average thermal flux in the k Œ-th zone. With usual power reactors, the reactivity of the core at the end of a cycle is a decreasing function of the operation time, so that it may be assumed that the reactor is just critical at the full power with Xe and Sm at equilibrium, and with all the control rods fully withdrawn. Therefore, we can obtain the consistent thermal flux P02 by solving (13) (7) (8) (9) (10) (11) where P1 is the non-leakage factor for the neutrons slowing down from fission energy to resonance. All the refueling schemes considered here are multiple-batch. The core is divided which can be done by letting keff=1.0 and Sc =0 in Eq.(1). Here we consider that the neutron flux is normalized so as to satisfy where Pt, is the thermal output in MW. (14) If the neutron flux distribution is constant, the total power varies with time depending upon the variation of the fission cross section. Constant power can nevertheless be maintained if desired by changing the neutron flux level and without changing the shape of the flux distribution. 15
3 and 186 J. Nucl. Sci. Technol. In order to obtain the solution of Eq. (13), it is convenient to let Sc=0 in Eq. (1). Two iteration schemes are used in order to obtain the solution of Eq. (13): the first being an ginner iteration hto determine eff for a given value of T, and the k second an gouter iteration hto obtain the value of T that gives keff=1 for Sc=0 in Eq. (1). In the inner iteration, a normalized flux 2 is assumed in obtaining tt from Eq. (12), P upon which the exposure-dependent constants can be calculated. Solving Eq. (1) for this tt, a new q2 and eigenvalue are obtained. Then the new P2 is used for the next value of tt. This procedure is repeated until tt2 and keff converges for the given T. For some inputs, however, oscillation of the eigenvalue is experienced during the iterations performed by the above procedure. In order to remove this oscillation, we replace P2(g-1) (g being the iteration number) by the formula (15) for calculating tt in Eq. (12). It has been found that this modification not only removes the oscillation but also reduces the computing time for all cases to approximately one third of that required when P(g-1)2 is directly used. In the first and second steps of the outer iterations, certain values are assumed for T1 and T2 for each of which the corresponding and k2, - are calculated keff-i. with e k1 Eq. (1). After the second iteration, letting e be the iteration number, Te is extrapolated by where (17) (18) (19) In the above equations, t is the time measured from the beginning of the cycle. Assuming that koo is a decreasing function of exposure, it can be proved that the consistent flux has the lowest HSF (Hot Spot Factor). If the core is operated with a flux with HSF lower than that of the consistent flux, koo at flux peaking becomes higher than the case of consistent flux at the end of the cycle, so that the peaking level becomes higher than the consistent flux. Conversely, in order to obtain a lower peaking at the end of the cycle, a higher flux peaking must be maintained at the beginning of the cycle. III. NUMERICAL ILLUSTRATIONS For the purpose of illustrating the application of the present method, calculations were performed for a large power reactor with the characteristics shown in Table 1. Table 1 Characteristics of the Core (16) and the corresponding ke is calculated. This procedure is repeated until ke converges to unity. Once the solution of Eq. (13) is obtained, it can be maintained throughout the cycle by varying the control rod distribution. The time dependence of the control rod during the cycle can be obtained by solving of Eq. (1), letting P2=P02 and keff=1.0: Sc Two types of refueling schemes were considered: (1) 5-batch out-in The core is divided into five concentric circular or annular zones of equal area, During each refueling, the fuel in the first zone is discharged, the fuel in the second to fifth zones are shifted to the next inner zones, and fresh fuel is 16
4 Vol. 3 No. 5 (May 1966) 187 introduced in the fifth zone. The shuffling matrix for this scheme is (2) 5-batch scatter The core is divided into ten concentric zones of equal area for instance. Each of the zones are further divided into five subzones of equal area. During each refueling, fuel in the highest exposure subzone in each zone is discharged, the fuel in the remaining subzones shifted inwards, and the outermost subzone loaded with new fuel, as shown in Fig. 1. Fig. 2 The Consistent Thermal Flux and Control Rod Distributions for the 5-batch Out-in Scheme Fig. 1 Approximation for the 5-batch Scatter Scheme This refueling scheme is none other than an out-in refueling scheme performed in each zone, and can be expressed in terms of a diagonal matrix where This approximation for the scatter scheme is valid when the widths of the subzones are so small that the variation of thermal flux within the five adjacent subzones is negligible. Figures 2 and 3 show, for the two types of refueling schemes being considered, the consistent thermal flux distribution obtained from Eq. (13), as well as the control rod distribution required to maintain these flux distributions. Despite the appearance of Fig. 3 The Consistent Thermal Flux and Control Rod Distributions for the 5-batch Scatter Scheme (The shape of the small flux peakings will be changed if other calculation model is used.) negative values of Sc in Fig. 3, the calculated can still be used for control Sc rod programming, based on the following considerations: (1) Local variations in Sc hardly influence the results of the burnup calculation, so that Sc can be averaged in each of the small regions. When Sc is averaged for each zone, only positive values are obtained. (2) The actual control rods used in practice are heterogeneous, whereas in this paper they have been treated as homogeneous. The control rod pattern is determined in such manner as to satisfy the condition that the locally averaged values of the 17
5 188 J. Nucl. Sci Technol effective absorption cross section of the actual control rods are equal to the locally averaged values of Sc. ACKNOWLEDGEMENT The author wishes to thank his colleagues Mr. T. Kawai, Mr. M. Yokomi and Mr. T. Morita for their stimulating and valuable discussions. The numerical calculations were performed by Mr. H. Kobayashi. REFERENCES (1) NAKAMURA, S.: Effect of Fuel; Shuffling on Burnup of BWR, Preprint of 1961 National Nuclear Congress, Tokyo, Japan, A25, (1961). (2) BENEDICT, M., PIGFORD, T.H.: gnuclear Chemical Engineering h, 91 (1957), McGraw-Hill Book Co. 18
Simple Method to Predict Power Level and Core Flow Rate of Boiling Water Reactors by Using One-Point Core Model
Journal of NUCLEAR SCIENCE and TECHNOLOGY, 19[7], pp. 513~520 (July 1982). 513 Simple Method to Predict Power Level and Core Flow Rate of Boiling Water Reactors by Using One-Point Core Model Kikuo UMEGAKI,
More informationLesson 14: Reactivity Variations and Control
Lesson 14: Reactivity Variations and Control Reactivity Variations External, Internal Short-term Variations Reactivity Feedbacks Reactivity Coefficients and Safety Medium-term Variations Xe 135 Poisoning
More informationHigh-Order Finite Difference Nodal Method for Neutron Diffusion Equation
Journal of NUCLEAR SCIENCE and TECHNOLOGY, 28[4], pp. 285~292 (April 1991) 285 High-Order Finite Difference Nodal Method for Neutron Diffusion Equation Kazuo AZEKURA and Kunitoshi KURIHARA Energy Research
More informationLesson 9: Multiplying Media (Reactors)
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,
More informationbut mostly as the result of the beta decay of its precursor 135 I (which has a half-life of hours).
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
More information2. The Steady State and the Diffusion Equation
2. The Steady State and the Diffusion Equation The Neutron Field Basic field quantity in reactor physics is the neutron angular flux density distribution: Φ( r r, E, r Ω,t) = v(e)n( r r, E, r Ω,t) -- distribution
More informationAvailable online at ScienceDirect. Energy Procedia 71 (2015 )
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 71 (2015 ) 97 105 The Fourth International Symposium on Innovative Nuclear Energy Systems, INES-4 High-Safety Fast Reactor Core Concepts
More informationLecture 20 Reactor Theory-V
Objectives In this lecture you will learn the following We will discuss the criticality condition and then introduce the concept of k eff.. We then will introduce the four factor formula and two group
More informationReactivity Effect of Fission and Absorption
Journal of NUCLEAR SCIENCE and TECHNOLOGY, 1, No. 7, p.246~254 (October, 1964). J. Nucl. Sci. Tech. Space Dependent Reactivity Effect of Fission and Absorption Tohru HAGA* and Iwao KOBAYASHI* Received
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 4. Title: Control Rods and Sub-critical Systems
Lectures on Nuclear Power Safety Lecture No 4 Title: Control Rods and Sub-critical Systems Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture Control Rods Selection of Control
More informationTRANSMUTATION OF CESIUM-135 WITH FAST REACTORS
TRANSMUTATION OF CESIUM-3 WITH FAST REACTORS Shigeo Ohki and Naoyuki Takaki O-arai Engineering Center Japan Nuclear Cycle Development Institute (JNC) 42, Narita-cho, O-arai-machi, Higashi-Ibaraki-gun,
More informationREACTOR PHYSICS ASPECTS OF PLUTONIUM RECYCLING IN PWRs
REACTOR PHYSICS ASPECTS OF PLUTONIUM RECYCLING IN s Present address: J.L. Kloosterman Interfaculty Reactor Institute Delft University of Technology Mekelweg 15, NL-2629 JB Delft, the Netherlands Fax: ++31
More informationLesson 8: Slowing Down Spectra, p, Fermi Age
Lesson 8: Slowing Down Spectra, p, Fermi Age Slowing Down Spectra in Infinite Homogeneous Media Resonance Escape Probability ( p ) Resonance Integral ( I, I eff ) p, for a Reactor Lattice Semi-empirical
More informationJOYO MK-III Performance Test at Low Power and Its Analysis
PHYSOR 200 -The Physics of Fuel Cycles and Advanced Nuclear Systems: Global Developments Chicago, Illinois, April 25-29, 200, on CD-ROM, American Nuclear Society, Lagrange Park, IL. (200) JOYO MK-III Performance
More informationO-arai Engineering Center Power Reactor and Nuclear Fuel Development Corporation 4002 Narita, O-arai-machi, Ibaraki-ken JAPAN ABSTRACT
Characteristics of TRU Transmutation in an LMFBR M. Yamaoka, M. Ishikawa, and T. Wakabayashi O-arai Engineering Center Power Reactor and Nuclear Fuel Development Corporation 4002 Narita, O-arai-machi,
More informationChapter 7 & 8 Control Rods Fission Product Poisons. Ryan Schow
Chapter 7 & 8 Control Rods Fission Product Poisons Ryan Schow Ch. 7 OBJECTIVES 1. Define rod shadow and describe its causes and effects. 2. Sketch typical differential and integral rod worth curves and
More informationStudy of Burnup Reactivity and Isotopic Inventories in REBUS Program
Study of Burnup Reactivity and Isotopic Inventories in REBUS Program T. Yamamoto 1, Y. Ando 1, K. Sakurada 2, Y. Hayashi 2, and K. Azekura 3 1 Japan Nuclear Energy Safety Organization (JNES) 2 Toshiba
More informationStudy of Control rod worth in the TMSR
Nuclear Science and Techniques 24 (2013) 010601 Study of Control rod worth in the TMSR ZHOU Xuemei * LIU Guimin 1 Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China
More informationXenon Effects. B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec.
enon Effects B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec. 2015 September 1 Contents We study the importance of e-135 in the operation of
More informationFigure 22.1 Unflattened Flux Distribution
22 Neutron Flux Control If nothing were done to flatten the flux in our reactors, it would look something like Figure 22.1. The flux would be a maximum in the efuel of the reactor (where neutrons are moving
More informationChemical Engineering 412
Chemical Engineering 412 Introductory Nuclear Engineering Lecture 18 Nuclear Reactor Theory IV Reactivity Insertions 1 Spiritual Thought 2 Mosiah 2:33 33 For behold, there is a wo pronounced upon him who
More informationReactor Operation with Feedback Effects
22.05 Reactor Physics - Part Twenty-Nine Reactor Operation with Feedback Effects 1. Reference Material: See pp. 368 372 in Light Water Reactor Control Systems, in Wiley Encyclopedia of Electrical and Electronics
More informationXV. Fission Product Poisoning
XV. Fission Product Poisoning XV.1. Xe 135 Buil-Up As we already know, temperature changes bring short-term effects. That is to say, once a power change is produced it is rapidly manifested as a change
More informationOn the Use of Serpent for SMR Modeling and Cross Section Generation
On the Use of Serpent for SMR Modeling and Cross Section Generation Yousef Alzaben, Victor. H. Sánchez-Espinoza, Robert Stieglitz INSTITUTE for NEUTRON PHYSICS and REACTOR TECHNOLOGY (INR) KIT The Research
More informationX. Neutron and Power Distribution
X. Neutron and Power Distribution X.1. Distribution of the Neutron Flux in the Reactor In order for the power generated by the fission reactions to be maintained at a constant level, the fission rate must
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 5. Title: Reactor Kinetics and Reactor Operation
Lectures on Nuclear Power Safety Lecture No 5 Title: Reactor Kinetics and Reactor Operation Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture (1) Reactor Kinetics Reactor
More informationNeutronic Issues and Ways to Resolve Them. P.A. Fomichenko National Research Center Kurchatov Institute Yu.P. Sukharev JSC Afrikantov OKBM,
GT-MHR Project High-Temperature Reactor Neutronic Issues and Ways to Resolve Them P.A. Fomichenko National Research Center Kurchatov Institute Yu.P. Sukharev JSC Afrikantov OKBM, GT-MHR PROJECT MISSION
More informationChem 481 Lecture Material 4/22/09
Chem 481 Lecture Material 4/22/09 Nuclear Reactors Poisons The neutron population in an operating reactor is controlled by the use of poisons in the form of control rods. A poison is any substance that
More informationTarget accuracy of MA nuclear data and progress in validation by post irradiation experiments with the fast reactor JOYO
Target accuracy of MA nuclear data and progress in validation by post irradiation experiments with the fast reactor JOYO Shigeo OHKI, Kenji YOKOYAMA, Kazuyuki NUMATA *, and Tomoyuki JIN * Oarai Engineering
More informationCross Section Generation Strategy for High Conversion Light Water Reactors Bryan Herman and Eugene Shwageraus
Cross Section Generation Strategy for High Conversion Light Water Reactors Bryan Herman and Eugene Shwageraus 1 Department of Nuclear Science and Engineering Massachusetts Institute of Technology 77 Massachusetts
More informationIdaho National Laboratory Reactor Analysis Applications of the Serpent Lattice Physics Code
Idaho National Laboratory Reactor Analysis Applications of the Serpent Lattice Physics Code By Frederick N. Gleicher II, Javier Ortensi, Benjamin Baker, and Mark DeHart Outline Intra-Pin Power and Flux
More informationLecture 27 Reactor Kinetics-III
Objectives In this lecture you will learn the following In this lecture we will understand some general concepts on control. We will learn about reactivity coefficients and their general nature. Finally,
More informationTHORIUM SELF-SUFFICIENT FUEL CYCLE OF CANDU POWER REACTOR
International Conference Nuclear Energy for New Europe 2005 Bled, Slovenia, September 5-8, 2005 ABSTRACT THORIUM SELF-SUFFICIENT FUEL CYCLE OF CANDU POWER REACTOR Boris Bergelson, Alexander Gerasimov Institute
More informationAnalytic Solution Technique for Solving One-Group Diffusion Equations for Core Simulations
Journal of NUCLEAR SCIENCE and TECHNOLOGY, 20[8] pp. 627~644 (August 1983). 627 Analytic Solution Technique for Solving One-Group Diffusion Equations for Core Simulations Masafumi ITAGAKI Japan Nuclear
More informationError Estimation for ADS Nuclear Properties by using Nuclear Data Covariances
Error Estimation for ADS Nuclear Properties by using Nuclear Data Covariances Kasufumi TSUJIMOTO Center for Proton Accelerator Facilities, Japan Atomic Energy Research Institute Tokai-mura, Naka-gun, Ibaraki-ken
More informationReactivity Coefficients
Reactivity Coefficients B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec. 2015 September 1 Reactivity Changes In studying kinetics, we have seen
More informationReactivity Coefficients
Revision 1 December 2014 Reactivity Coefficients Student Guide GENERAL DISTRIBUTION GENERAL DISTRIBUTION: Copyright 2014 by the National Academy for Nuclear Training. Not for sale or for commercial use.
More informationPIA and REWIND: Two New Methodologies for Cross Section Adjustment. G. Palmiotti and M. Salvatores
PIA and REWIND: Two New Methodologies for Cross Section Adjustment G. Palmiotti and M. Salvatores Nuclear Systems Design and Analysis Division, Idaho National Laboratory, P.O. Box 1625, MS 3860, Idaho
More informationA COMPARISON OF PEBBLE MIXING AND DEPLETION ALGORITHMS USED IN PEBBLE-BED REACTOR EQUILIBRIUM CYCLE SIMULATION
International Conference on Mathematics, Computational Methods & Reactor Physics (M&C 2009) Saratoga Springs, New York, May 3-7, 2009, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2009) A COMPARISON
More informationNEUTRONIC ANALYSIS OF HE-EFIT EFIT ADS - SOME RESULTS -
NEUTRONIC ANALYSIS OF HE-EFIT EFIT ADS - SOME RESULTS - Alan Takibayev & Danas Ridikas CEA Saclay / DSM / IRFU Atelier GEDEPEON 'Accelerator Driven System' Aix-en-Provence 15-10-2008 HE-EFIT MAIN CHARACTERISTICS
More informationNuclear Theory - Course 127 EFFECTS OF FUEL BURNUP
Nuclear Theory - Course 127 EFFECTS OF FUEL BURNUP The effect of fuel burnup wa~ considered, to some extent, in a previous lesson. During fuel burnup, U-235 is used up and plutonium is produced and later
More informationOn the use of SERPENT code for few-group XS generation for Sodium Fast Reactors
On the use of SERPENT code for few-group XS generation for Sodium Fast Reactors Raquel Ochoa Nuclear Engineering Department UPM CONTENTS: 1. Introduction 2. Comparison with ERANOS 3. Parameters required
More informationULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor
ULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor Ade afar Abdullah Electrical Engineering Department, Faculty of Technology and Vocational Education Indonesia University of
More information3. State each of the four types of inelastic collisions, giving an example of each (zaa type example is acceptable)
Nuclear Theory - Course 227 OBJECTIVES to: At the conclusion of this course the trainee will be able 227.00-1 Nuclear Structure 1. Explain and use the ZXA notation. 2. Explain the concept of binding energy.
More informationSUB-CHAPTER D.1. SUMMARY DESCRIPTION
PAGE : 1 / 12 CHAPTER D. REACTOR AND CORE SUB-CHAPTER D.1. SUMMARY DESCRIPTION Chapter D describes the nuclear, hydraulic and thermal characteristics of the reactor, the proposals made at the present stage
More informationDESIGN OF B 4 C BURNABLE PARTICLES MIXED IN LEU FUEL FOR HTRS
DESIGN OF B 4 C BURNABLE PARTICLES MIXED IN LEU FUEL FOR HTRS V. Berthou, J.L. Kloosterman, H. Van Dam, T.H.J.J. Van der Hagen. Delft University of Technology Interfaculty Reactor Institute Mekelweg 5,
More informationSolving Bateman Equation for Xenon Transient Analysis Using Numerical Methods
Solving Bateman Equation for Xenon Transient Analysis Using Numerical Methods Zechuan Ding Illume Research, 405 Xintianshiji Business Center, 5 Shixia Road, Shenzhen, China Abstract. After a nuclear reactor
More informationSENSITIVITY ANALYSIS OF ALLEGRO MOX CORE. Bratislava, Iľkovičova 3, Bratislava, Slovakia
SENSITIVITY ANALYSIS OF ALLEGRO MOX CORE Jakub Lüley 1, Ján Haščík 1, Vladimír Slugeň 1, Vladimír Nečas 1 1 Institute of Nuclear and Physical Engineering, Slovak University of Technology in Bratislava,
More information17 Neutron Life Cycle
17 Neutron Life Cycle A typical neutron, from birth as a prompt fission neutron to absorption in the fuel, survives for about 0.001 s (the neutron lifetime) in a CANDU. During this short lifetime, it travels
More informationMONTE CALRLO MODELLING OF VOID COEFFICIENT OF REACTIVITY EXPERIMENT
MONTE CALRLO MODELLING OF VOID COEFFICIENT OF REACTIVITY EXPERIMENT R. KHAN, M. VILLA, H. BÖCK Vienna University of Technology Atominstitute Stadionallee 2, A-1020, Vienna, Austria ABSTRACT The Atominstitute
More informationThe moderator temperature coefficient MTC is defined as the change in reactivity per degree change in moderator temperature.
Moderator Temperature Coefficient MTC 1 Moderator Temperature Coefficient The moderator temperature coefficient MTC is defined as the change in reactivity per degree change in moderator temperature. α
More informationHTR Spherical Super Lattice Model For Equilibrium Fuel Cycle Analysis. Gray S. Chang. September 12-15, 2005
INEEL/CON-05-02655, Revision 3 PREPRINT HTR Spherical Super Lattice Model For Equilibrium Fuel Cycle Analysis Gray S. Chang September 12-15, 2005 Mathematics And Computation, Supercomputing, Reactor Physics
More informationNuclear Theory - Course 227 EFFECT OF FUEL BURNUP
Nuclear Theory - Course 227 EFFECT OF FUEL BURNUP The changes in the composition of the fuel as it is depleted give rise to a number of effects which may be described under the following headings: 1) Long
More informationStudy on SiC Components to Improve the Neutron Economy in HTGR
Study on SiC Components to Improve the Neutron Economy in HTGR Piyatida TRINURUK and Assoc.Prof.Dr. Toru OBARA Department of Nuclear Engineering Research Laboratory for Nuclear Reactors Tokyo Institute
More informationReactor Kinetics and Operation
Reactor Kinetics and Operation Course No: N03-002 Credit: 3 PDH Gilbert Gedeon, P.E. Continuing Education and Development, Inc. 9 Greyridge Farm Court Stony Point, NY 0980 P: (877) 322-5800 F: (877) 322-4774
More informationNeutronic Calculations of Ghana Research Reactor-1 LEU Core
Neutronic Calculations of Ghana Research Reactor-1 LEU Core Manowogbor VC*, Odoi HC and Abrefah RG Department of Nuclear Engineering, School of Nuclear Allied Sciences, University of Ghana Commentary Received
More informationNEUTRON MODERATION. LIST three desirable characteristics of a moderator.
Reactor Theory (eutron Characteristics) DOE-HDBK-1019/1-93 EUTRO MODERATIO EUTRO MODERATIO In thermal reactors, the neutrons that cause fission are at a much lower energy than the energy level at which
More informationShutdown Margin. Xenon-Free Xenon removes neutrons from the life-cycle. So, xenonfree is the most reactive condition.
22.05 Reactor Physics - Part Thirty-One Shutdown Margin 1. Shutdown Margin: Shutdown margin (abbreviated here as SDM) is defined as the amount of reactivity by which a reactor is subcritical from a given
More informationIntroduction to Reactivity and Reactor Control
Introduction to Reactivity and Reactor Control Larry Foulke Adjunct Professor Director of Nuclear Education Outreach University of Pittsburgh IAEA Workshop on Desktop Simulation October 2011 Learning Objectives
More informationSPACE-DEPENDENT DYNAMICS OF PWR. T. Suzudo Japan Atomic Energy Research Institute, JAERI Tokai-Mura, Naka-Gun Japan
SPACE-DEPENDENT DYNAMICS OF PWR T. Suzudo Japan Atomic Energy Research Institute, JAERI Tokai-Mura, Naka-Gun 319-11 Japan E. Türkcan and J.P. Verhoef Netherlands Energy Research Foundation P.O. Box 1,
More informationEVALUATION OF PWR AND BWR CALCULATIONAL BENCHMARKS FROM NUREG/CR-6115 USING THE TRANSFX NUCLEAR ANALYSIS SOFTWARE
ANS MC2015 - Joint International Conference on Mathematics and Computation (M&C), Supercomputing in Nuclear Applications (SNA) and the Monte Carlo (MC) Method Nashville, Tennessee April 19 23, 2015, on
More informationChallenges in Prismatic HTR Reactor Physics
Challenges in Prismatic HTR Reactor Physics Javier Ortensi R&D Scientist - Idaho National Laboratory www.inl.gov Advanced Reactor Concepts Workshop, PHYSOR 2012 April 15, 2012 Outline HTR reactor physics
More informationAdaptation of Pb-Bi Cooled, Metal Fuel Subcritical Reactor for Use with a Tokamak Fusion Neutron Source
Adaptation of Pb-Bi Cooled, Metal Fuel Subcritical Reactor for Use with a Tokamak Fusion Neutron Source E. Hoffman, W. Stacey, G. Kessler, D. Ulevich, J. Mandrekas, A. Mauer, C. Kirby, D. Stopp, J. Noble
More informationCritical Experiment Analyses by CHAPLET-3D Code in Two- and Three-Dimensional Core Models
Journal of NUCLEAR SCIENCE and TECHNOLOGY, Vol. 42, No. 1, p. 101 108 (January 2005) TECHNICAL REPORT Critical Experiment Analyses by CHAPLET-3D Code in Two- and Three-Dimensional Core Models Shinya KOSAKA
More informationEstimation of Control Rods Worth for WWR-S Research Reactor Using WIMS-D4 and CITATION Codes
Estimation of Control Rods Worth for WWR-S Research Reactor Using WIMS-D4 and CITATION Codes M. S. El-Nagdy 1, M. S. El-Koliel 2, D. H. Daher 1,2 )1( Department of Physics, Faculty of Science, Halwan University,
More informationFuel Element Burnup Determination in HEU - LEU Mixed TRIGA Research Reactor Core
Fuel Element Burnup Determination in HEU - LEU Mixed TRIGA Research Reactor Core Tomaž Žagar, Matjaž Ravnik Institute "Jožef Stefan", Jamova 39, Ljubljana, Slovenia Tomaz.Zagar@ijs.si Abstract This paper
More information(1) SCK CEN, Boeretang 200, B-2400 Mol, Belgium (2) Belgonucléaire, Av. Arianelaan 4, B-1200 Brussels, Belgium
The REBUS Experimental Programme for Burn-up Credit Peter Baeten (1)*, Pierre D'hondt (1), Leo Sannen (1), Daniel Marloye (2), Benoit Lance (2), Alfred Renard (2), Jacques Basselier (2) (1) SCK CEN, Boeretang
More informationVIII. Neutron Moderation and the Six Factors
Introduction VIII. Neutron Moderation and the Six Factors 130 We continue our quest to calculate the multiplication factor (keff) and the neutron distribution (in position and energy) in nuclear reactors.
More informationExcerpt from the Proceedings of the COMSOL Users Conference 2007 Grenoble
Excerpt from the Proceedings of the COSOL Users Conference 007 Grenoble Evaluation of the moderator temperature coefficient of reactivity in a PWR V. emoli *,, A. Cammi Politecnico di ilano, Department
More informationSpectral History Correction of Microscopic Cross Sections for the PBR Using the Slowing Down Balance. Abstract
Organized and hosted by the Canadian Nuclear Society. Vancouver, BC, Canada. 2006 September 10-14 Spectral History Correction of Microscopic Cross Sections for the PBR Using the Slowing Down Balance Nathanael
More informationFuel cycle studies on minor actinide transmutation in Generation IV fast reactors
Fuel cycle studies on minor actinide transmutation in Generation IV fast reactors M. Halász, M. Szieberth, S. Fehér Budapest University of Technology and Economics, Institute of Nuclear Techniques Contents
More informationOperational Reactor Safety
Operational Reactor Safety 22.091/22.903 Professor Andrew C. Kadak Professor of the Practice Lecture 3 Reactor Kinetics and Control Page 1 Topics to Be Covered Time Dependent Diffusion Equation Prompt
More informationPHYSICS AND KINETICS OF TRIGA REACTOR. H. Böck and M. Villa AIAU 27307
PHYSICS AND KINETICS OF TRIGA REACTOR H. Böck and M. Villa AIAU 27307 *prepared for NTEC Overview This training module is written as an introduction to reactor physics for reactor operators. It assumes
More informationCRITICAL AND SUBCRITICAL EXPERIMENTS USING THE TRAINING NUCLEAR REACTOR OF THE BUDAPEST UNIVERSITY OF TECHNOLOGY AND ECONOMICS
CRITICAL AND SUBCRITICAL EXPERIMENTS USING THE TRAINING NUCLEAR REACTOR OF THE BUDAPEST UNIVERSITY OF TECHNOLOGY AND ECONOMICS É. M. Zsolnay Department of Nuclear Techniques, Budapest University of Technology
More informationPower Changes in a Critical Reactor. The Critical Reactor
Chapter 8 Power Changes in a Critical Reactor n For very small reactivity increases n For small reactivity increases n For large reactivity increases/decreases The Critical Reactor < k = hfpel f L t =
More informationNuclear Fission. Q for 238 U + n 239 U is 4.??? MeV. E A for 239 U 6.6 MeV MeV neutrons are needed.
Q for 235 U + n 236 U is 6.54478 MeV. Table 13.11 in Krane: Activation energy E A for 236 U 6.2 MeV (Liquid drop + shell) 235 U can be fissioned with zero-energy neutrons. Q for 238 U + n 239 U is 4.???
More informationSpatially Dependent Self-Shielding Method with Temperature Distribution for the Two-Dimensional
Journal of Nuclear Science and Technology ISSN: 0022-3131 (Print) 1881-1248 (Online) Journal homepage: http://www.tandfonline.com/loi/tnst20 Spatially Dependent Self-Shielding Method with Temperature Distribution
More informationReactor Operation Without Feedback Effects
22.05 Reactor Physics - Part Twenty-Six Reactor Operation Without Feedback Effects 1. Reference Material: See pp. 363-368 of the article, Light Water Reactor Control Systems, in Wiley Encyclopedia of Electrical
More informationDelayed neutrons in nuclear fission chain reaction
Delayed neutrons in nuclear fission chain reaction 1 Critical state Temporal flow Loss by leakage Loss by Absorption If the number of neutrons (the number of fission reactions) is practically constant
More informationChapter 6 Development of the Method to Assay Barely Measurable Elements in Spent Nuclear Fuel and Application to BWR 9 9 Fuel
Chapter 6 Development of the Method to Assay Barely Measurable Elements in Spent Nuclear Fuel and Application to BWR 9 9 Fuel Kenya Suyama, Gunzo Uchiyama, Hiroyuki Fukaya, Miki Umeda, Toru Yamamoto, and
More informationCALCULATION OF TEMPERATURE REACTIVITY COEFFICIENTS IN KRITZ-2 CRITICAL EXPERIMENTS USING WIMS ABSTRACT
CALCULATION OF TEMPERATURE REACTIVITY COEFFICIENTS IN KRITZ-2 CRITICAL EXPERIMENTS USING WIMS D J Powney AEA Technology, Nuclear Science, Winfrith Technology Centre, Dorchester, Dorset DT2 8DH United Kingdom
More informationVerification of Core Monitoring System with Gamma Thermometer
GENES4/ANP2003, Sep. 15-19, 2003, Kyoto, JAPAN Paper 1006 Verification of Core Monitoring System with Gamma Thermometer Hisashi SHIRAGA 1*, Hiromi MARUYAMA 1, Atsushi FUSHIMI 2, Yoshiji KARINO 3 and Hiroyuki
More informationMA/LLFP Transmutation Experiment Options in the Future Monju Core
MA/LLFP Transmutation Experiment Options in the Future Monju Core Akihiro KITANO 1, Hiroshi NISHI 1*, Junichi ISHIBASHI 1 and Mitsuaki YAMAOKA 2 1 International Cooperation and Technology Development Center,
More informationVI. Chain Reaction. Two basic requirements must be filled in order to produce power in a reactor:
VI. Chain Reaction VI.1. Basic of Chain Reaction Two basic requirements must be filled in order to produce power in a reactor: The fission rate should be high. This rate must be continuously maintained.
More informationDevelopment of depletion models for radionuclide inventory, decay heat and source term estimation in discharged fuel
Development of depletion models for radionuclide inventory, decay heat and source term estimation in discharged fuel S. Caruso, A. Shama, M. M. Gutierrez National Cooperative for the Disposal of Radioactive
More informationTesting the EPRI Reactivity Depletion Decrement Uncertainty Methods
Testing the EPRI Reactivity Depletion Decrement Uncertainty Methods by Elliot M. Sykora B.S. Physics, Massachusetts Institute of Technology (0) Submitted to the Department of Nuclear Science and Engineering
More informationLecture 28 Reactor Kinetics-IV
Objectives In this lecture you will learn the following In this lecture we will understand the transient build up of Xenon. This can lead to dead time in reactors. Xenon also induces power oscillations
More informationDETAILED CORE DESIGN AND FLOW COOLANT CONDITIONS FOR NEUTRON FLUX MAXIMIZATION IN RESEARCH REACTORS
Proceedings of ICONE14 14th International Conference on Nuclear Engineering June 17-21, 2006, Miami, FL, USA ICONE14-89547 DETAILED CORE DESIGN AND FLOW COOLANT CONDITIONS FOR NEUTRON FLUX MAXIMIZATION
More informationNumerical analysis on element creation by nuclear transmutation of fission products
NUCLEAR SCIENCE AND TECHNIQUES 26, S10311 (2015) Numerical analysis on element creation by nuclear transmutation of fission products Atsunori Terashima 1, and Masaki Ozawa 2 1 Department of Nuclear Engineering,
More informationNuclear Data for Emergency Preparedness of Nuclear Power Plants Evaluation of Radioactivity Inventory in PWR using JENDL 3.3
Nuclear Data for Emergency Preparedness of Nuclear Power Plants Evaluation of Radioactivity Inventory in PWR using JENDL 3.3 Yoshitaka Yoshida, Itsuro Kimura Institute of Nuclear Technology, Institute
More informationREACTOR PHYSICS FOR NON-NUCLEAR ENGINEERS
International Conference on Mathematics and Computational Methods Applied to Nuclear Science and Engineering (M&C 2011) Rio de Janeiro, RJ, Brazil, May 8-12, 2011, on CD-ROM, Latin American Section (LAS)
More informationInvariability of neutron flux in internal and external irradiation sites in MNSRs for both HEU and LEU cores
Indian Journal of Pure & Applied Physics Vol. 49, February 2011, pp. 83-90 Invariability of neutron flux in internal and external irradiation sites in MNSRs for both HEU and LEU cores M Albarhoum Department
More informationNeutronic Analysis of Moroccan TRIGA MARK-II Research Reactor using the DRAGON.v5 and TRIVAC.v5 codes
Physics AUC, vol. 27, 41-49 (2017) PHYSICS AUC Neutronic Analysis of Moroccan TRIGA MARK-II Research Reactor using the DRAGON.v5 and TRIVAC.v5 codes DARIF Abdelaziz, CHETAINE Abdelouahed, KABACH Ouadie,
More informationR.A. Chaplin Department of Chemical Engineering, University of New Brunswick, Canada
REACTVTY CHANGES R.A. Chaplin Department of Chemical Engineering, University of New Brunswick, Canada Keywords: Reactivity Coefficient, Fuel Burnup, Xenon, Samarium, Temperature Contents 1. ntroduction
More informationFundamentals of Nuclear Reactor Physics
Fundamentals of Nuclear Reactor Physics E. E. Lewis Professor of Mechanical Engineering McCormick School of Engineering and Applied Science Northwestern University AMSTERDAM BOSTON HEIDELBERG LONDON NEW
More informationRecycling Spent Nuclear Fuel Option for Nuclear Sustainability and more proliferation resistance In FBR
Recycling Spent Nuclear Fuel Option for Nuclear Sustainability and more proliferation resistance In FBR SIDIK PERMANA a, DWI IRWANTO a, MITSUTOSHI SUZUKI b, MASAKI SAITO c, ZAKI SUUD a a Nuclear Physics
More informationABSTRACT 1 INTRODUCTION
A NODAL SP 3 APPROACH FOR REACTORS WITH HEXAGONAL FUEL ASSEMBLIES S. Duerigen, U. Grundmann, S. Mittag, B. Merk, S. Kliem Forschungszentrum Dresden-Rossendorf e.v. Institute of Safety Research P.O. Box
More informationModule 6. Neutron Flux Distribution 6.1 MODULE OVERVIEW '!' MODULE OBJECTIVES OVERALL FLUX SHAPE FLUX FLATTENING 6
Module 6 Neutron Flux Distribution 6.1 MODULE OVERVIEW '!' 2 6.2 MODULE OBJECTIVES 2 6.3 OVERALL FLUX SHAPE 3 6.4 FLUX FLATTENING 6 6.4.1 The Reflector 7 6.4.2 Bi-
More informationRequests on Nuclear Data in the Backend Field through PIE Analysis
Requests on Nuclear Data in the Backend Field through PIE Analysis Yoshihira Ando 1), Yasushi Ohkawachi 2) 1) TOSHIBA Corporation Power System & Services Company Power & Industrial Systems Research & Development
More informationThe Effect of 99 Mo Production on the Neutronic Safety Parameters of Research Reactors
The Effect of 99 Mo Production on the Neutronic Safety Parameters of Research Reactors Riham M. Refeat and Heba K. Louis Safety Engineering Department, Nuclear and Radiological Regulation Authority (NRRA),
More information