Excerpt from the Proceedings of the COMSOL Users Conference 2007 Grenoble
|
|
- Lawrence Everett Conley
- 6 years ago
- Views:
Transcription
1 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 of Nuclear Engineering *Corresponding author: via Ponzio 34/3 ilano ilan (Italy), vito.memoli@mail.polimi.it Abstract: The moderator temperature coefficient (briefly TC) plays an important role within the power thermal reactor dynamics. In order to run a reactor safely, a negative moderator coefficient temperature is necessary to reach stability during changes in temperature that can be caused by a step insertion of reactivity. Thus, TC calculation is a key point in the reactor design process. The aim of the present work is to implement a simple model for the TC estimation in a PWR. The model is based on the well-known twogroup diffusion theory, which allows to perform criticality calculations with good approximation. The two-group diffusion PDE system has been defined and solved in COSOL ulti-physics which provides a finite-elements solution for any kind of geometry. The reactor has been modeled in -D as a cylindrical homogeneous core along with a radial reflector. The calculations of TC have been carried out performing several static criticality calculations for different temperature variations and different reflector thicknesses. The change temperature effect on group constants in the core zone has been evaluated by means of simple relations obtained starting from the definitions of resonance escape probability and absorption macroscopic cross section. Two core sizes have been analyzed: a large core, H=370cm and W=340cm, and a small core with H=W=00cm. Keywords: moderator temperature coefficient, PWR, two-group diffusion theory. Introduction A nuclear reactor should be designed with a negative TC. This requirement makes sure that, for example, positive excursion of power are quenched by the negative feedback. In particular in a PWR a negative moderator temperature coefficient has an important effect on the control of the reactor following changes in demand for power from the turbine. In a thermal reactor, an increase in the moderator temperature affects the multiplication factor in two ways: (a) the temperature at which thermal cross-sections are computed is changed, and (b) the physical density of the moderator changes because of thermal expansion. The main purpose of this work consists in implementing a simple model for the evaluation of the moderator temperature coefficient in a PWR, which accounts for the effect (b), ordinarily the more significant one. The model [,], which is discussed in the next section, is based on the two-group diffusion theory and the TC can be obtained by running several static criticality calculations. From the calculated multiplication factors the TC can be calculated as: δρ δt δt k k () where k and k are the multiplication factors for two different temperature (T, and T, ) and δt = T, -T,. In this paper we have focused our attention on the influence of the reflector thickness and the magnitude of the moderator temperature change on the moderator temperature coefficient estimation.. odeling The reactor model which has been developed consists of a homogeneous cylindrical multiplying core surrounded by a reflector. The effective multiplication factor is determined by solving the eigenvalue equations derived from the two-group diffusion theory. In order to reduce the computing time we have replaced the original 3-D geometry with a -D one. The loss of neutrons in the axial direction has been taken into account by means of the axial buckling B z which can be calculated as B z π = H 0 ()
2 Excerpt from the Proceedings of the COSOL Users Conference 007 Grenoble where H 0 is the height of the cylinder plus the reflector saving which has been taken as one diffusion length in the reflector medium [3,4,5]. That means: H 0 =H + δ (3) with δ L (4) b where L b is the neutron diffusion length in the blanket (reflector). The diffusion length has been calculated as the square root of ratio between the diffusion coefficient and the absorption cross section in the reflector. According to the previous hypothesis the twogroup diffusion equations, in the core region, become: D B z D x, y + a, + = ( ν f, + ν f, ) k eff D + = B z D x, y a, where and are respectively the fast and the thermal neutron fluxes. D i is the diffusion coefficient and a,i is the absorption cross section (where index i= stands for fast group and i= for thermal groupp), is the removal cross section (i.e. the probability per centimeter that a fast neutron is slowed down to the thermal group). The fudge factor /k eff, the inverse of the multiplication factor, makes the equations time independent. The two group diffusion equations have been processed in the PDE module of COSOL ulti-physics where the direct method UFPACK has been set as eigenvalue solver. In order to validate the numerical code the numerical solution has been compared with the well-known analytical one in the case of D reflected core[5]. In Fig. the calculated thermal and fast fluxes are shown. Figure. thermal and fast normalized neutron fluxes evaluated with Comsol for a -D reflected core (core diameter 40 cm, reflector thickness 0cm) As can be seen the agreement with the analytical solution is good. The fast flux shows a peak in the core region because of the presence of fuel which can be regarded as a source of fast neutrons due to fission of U-35. The thermal flux increases in nearby reflector region where the moderation is the main effect. Within the two-group diffusion theory, the effects of temperature change on nuclear quantities can be simply calculated from the definition of resonance escape probability, removal cross section from the fast group to the thermal one and the absorption cross section. oreover, the assumption of homogeneous core requires the estimation of the effect on the fission cross section. In fact for a homogeneous core, the density variation of the moderator due to a temperature change should have an impact on the fuel which is dispersed in the liquid moderator. In general the resonance escape probability, p, can be expressed by the following relation: p ρ I F ξρσ exp. (5) = s where ρ F and ρ are respectively the fuel and moderator densities, I is the resonance integral[3,4,5], ξ is the lethargy (average log energy loss) and σ s is the microscopic scattering cross section.
3 Excerpt from the Proceedings of the COSOL Users Conference 007 Grenoble The latter relation, which refers to the one energy group theory, reduces to the two-group resonance escape probability definition given by p= + a, (6) in the hypothesis of low neutron absorption. By calculating the variation of the resonance escape probability with respect to the temperature variation we obtain the following relation: δp δρ = log = log βδt (7) p p ρ p where β is the thermal expansion coefficient of the moderator. In order to evaluate the removal cross section variation with respect to temperature change, δt Μ, the two-group resonance escape probability definition can be rewritten as: p = a, (8) p Thus, making the derivative of the previous one and using (7) the removal cross section variation is given by: log( / p) δ β T = δ (9) p The temperature variation effect on the thermal absorption cross section is simply given by : δ a, = β a, δt. (0) Parameter Group Group a(cm - ) D(cm) (cm - ) νf(cm - ) Table. Two-Group diffusion constants in the core region Parameter Group Group a(cm - ) D(cm) (cm - ) Table. Two-Group diffusion constants in the reflector region. The above data have been used to evaluate the multiplication factor k. The new multiplication factor k is obtained by updating the cross sections, a,, f,i to the new ones i.e. + δ, a, +δ a,, f,i +δ f,i. Two core sizes have been taken into account. A large core with H=370cm and W=340cm and small core with H=00cm and W=00cm. For the reactor geometry an unstructured fine mesh has been used as it is shown in Fig. which corresponds to the large core with a reflector thickness of 0cm. In this case the number of the cells is equal to 95. In the same way we can evaluate the variation of the fission cross section in the thermal region: δ f, = β f, δt. () 3. Results In Table and are shown the typical nuclear constants for a PWR[3]. The thermal expansion coefficient [6] of the moderator, β = / C, corresponds to the usual PWR operating conditions i.e. pressure of 5Pa and water temperature of 30 C which corresponds to the water mean temperature between core inlet and core exit. Figure. esh of the core geometry. This choice can be justified by Fig.3 which shows the k eff estimation as a function of the element number. The k eff remains practically
4 Excerpt from the Proceedings of the COSOL Users Conference 007 Grenoble constant for meshes with an number of elelemnts greater or equal to oderator Temperature Coefficient as a function of Reflector Thickness ultiplication factor versus number of mesh elements TC(pcm/ C) Keff Number of elements Figure 3. k eff estimation as a function of the mesh size. In Fig. 4 the fast and thermal neutron fluxes as a function of radial distance is shown for the large core with a reflector thickness of 0cm Reflector Thickness (cm) Figure 5. oderator temperature coefficient versus reflector thickness for PWR big core for δt =4 C. In the case of small core the effect of the reflector thickness on the TC is much more significant. The TC %-variation between and reflected (0cm) reactor is about 7.7%. oderator Temperature Coefficient vs Reflector Thickness TC (pcm/ C) Reflector Thickness(cm) Figure 4. Fast and thermal neutron flux as a function of the radial distance in large core. The effect of the reflector thickness is shown in Figs. 5 and 6 for large and small core respectively. A moderator temperature change of 4 C has been considered as reference value. In the case of large core, the reflector thickness has a small influence on the moderator temperature coefficient. The variation of the TC between and reflected (0cm) reactor is about 0.4 %. Figure 6. oderator temperature coefficient versus reflector thickness for PWR small core for δt =4 C. In both cases the curve reaches a plateau due to the fact that after a certain thickness value (about 0cm) the reflector behaves practically as an infinite reflector [3,4,5]. The TC as function of temperature change is shown in Figs. 7 and 8 for and reflected reactor in the two core size cases.
5 Excerpt from the Proceedings of the COSOL Users Conference 007 Grenoble oderator Temperature Coefficient as a function of temperature chenge 4. Conclusions -55 TC(pcm/ C) δt ( C) reflected 0cm Figure 7. oderator temperature coefficient versus temperature change for large core. TC(pcm/ C) oderator Temperature Coefficient versus Temperature Change Reflected 0cm The influence of both the reflector thickness and temperature change magnitude on the moderator temperature coefficient in a PWR was investigated for two different size. To carry out the criticality calculations a -D homogeneous reactor has been modeled on the basis of the twogroup diffusion theory. Simple relations have been used to quantify the temperature change effect on the nuclear constants of absorption cross section, removal cross section and fission cross section. On the basis of the results it can be asserted that the reflector thickness has a negligible effect on the moderator temperature coefficient estimation for large core. Instead the small core introduces an uncertainty on the determination of the TC which can reach a value of about 8%. Within the range of 3 C to 5 C the TC estimation does not change significantly δt ( C) Figure 8. oderator temperature coefficient versus temperature change for small core. The percentage variation of the TC as a function of the temperature change is shown in Fig. 9. The results refer to the large core only, The percentage variation is practically negligible in a range 3 C δt 5 C which is the range of investigation recommended by the nuclear standard (ANS, 997). %-variation %-variation of the TC versus temperature change δt ( C) reflected 0cm 5. References. ANS,997. Calculation and measurement of the moderator temperature coefficient of reactivity for water moderated reactor power reactors. American Nuclear Society, American National ANSI/ANS C. Housidias et al., Calculation of the moderator temperature coefficient of reactivity for water moderated reactors, Annals of Nuclear Energy, 8, (00) 3. J.J. Duderstadt, L.J. Hamilton, 976. Nuclear Reactor Analysis. John Wiley & Sons, New York. 4. J. L. EE, 964. Two Group Reactor Theory, Gordon and Breach Science Publishers, Inc. 5. J.R. Lamarsh, 966. Introduction to nuclear Reactor Theory. Addison-Wesley Publishing Company Inc, Reading, assachusetts. 6. NIST, Standard Reference Database 3. Physical and Chemical Properties Division. Figure 9. Percentage Variation of TC as a function of moderator temperature change for and reflected reactor.
The 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 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 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 informationX. Assembling the Pieces
X. Assembling the Pieces 179 Introduction Our goal all along has been to gain an understanding of nuclear reactors. As we ve noted many times, this requires knowledge of how neutrons are produced and lost.
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 informationLamarsh, "Introduction to Nuclear Reactor Theory", Addison Wesley (1972), Ch. 12 & 13
NEEP 428 REACTOR PULSING Page 1 April 2003 References: Lamarsh, "Introduction to Nuclear Reactor Theory", Addison Wesley (1972), Ch. 12 & 13 Notation: Lamarsh (2), "Introduction to Nuclear Engineering",(1975),
More informationReactivity Power and Temperature Coefficients Determination of the TRR
Reactivity and Temperature Coefficients Determination of the TRR ABSTRACT Ahmad Lashkari Nuclear Science and Technology Research Institute (NSTRI), Atomic Energy Organization of Iran Tehran 14399-51113,
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 informationEnergy Dependence of Neutron Flux
Energy Dependence of Neutron Flux B. Rouben McMaster University Course EP 4D03/6D03 Nuclear Reactor Analysis (Reactor Physics) 2015 Sept.-Dec. 2015 September 1 Contents We start the discussion of the energy
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 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 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 informationA Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis
A Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis A. Aures 1,2, A. Pautz 2, K. Velkov 1, W. Zwermann 1 1 Gesellschaft für Anlagen- und Reaktorsicherheit (GRS) ggmbh Boltzmannstraße
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 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 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 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 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 informationCOMPARATIVE ANALYSIS OF WWER-440 REACTOR CORE WITH PARCS/HELIOS AND PARCS/SERPENT CODES
COMPARATIVE ANALYSIS OF WWER-440 REACTOR CORE WITH PARCS/HELIOS AND PARCS/SERPENT CODES S. Bznuni, A. Amirjanyan, N. Baghdasaryan Nuclear and Radiation Safety Center Yerevan, Armenia Email: s.bznuni@nrsc.am
More informationElastic scattering. Elastic scattering
Elastic scattering Now we have worked out how much energy is lost when a neutron is scattered through an angle, θ We would like to know how much energy, on average, is lost per collision In order to do
More informationPressurized Water Reactor (PWR)
Pressurized Water Reactor (PWR) Models: Involve neutron transport, thermal-hydraulics, chemistry Inherently multi-scale, multi-physics CRUD Measurements: Consist of low resolution images at limited number
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 informationNuclear Reactor Physics I Final Exam Solutions
.11 Nuclear Reactor Physics I Final Exam Solutions Author: Lulu Li Professor: Kord Smith May 5, 01 Prof. Smith wants to stress a couple of concepts that get people confused: Square cylinder means a cylinder
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 informationThe Use of Serpent 2 in Support of Modeling of the Transient Test Reactor at Idaho National Laboratory
The Use of Serpent 2 in Support of Modeling of the Transient Test Reactor at Idaho National Laboratory Sixth International Serpent User s Group Meeting Politecnico di Milano, Milan, Italy 26-29 September,
More informationA Dummy Core for V&V and Education & Training Purposes at TechnicAtome: In and Ex-Core Calculations
A Dummy Core for V&V and Education & Training Purposes at TechnicAtome: In and Ex-Core Calculations S. Nicolas, A. Noguès, L. Manifacier, L. Chabert TechnicAtome, CS 50497, 13593 Aix-en-Provence Cedex
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 informationUncertainty quantification using SCALE 6.2 package and GPT techniques implemented in Serpent 2
6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy September 26 th -30 th, 2016 Uncertainty quantification using SCALE 6.2 package and GPT techniques implemented in Serpent
More informationMODERATING RATIO PARAMETER EVALUATION FOR DIFFERENT MATERIALS BY MEANS OF MONTE CARLO CALCULATIONS AND REACTIVITY DIRECT MEASUREMENTS
MODERING RIO PRMEER EVLUION FOR DIFFEREN MERIL BY MEN OF MONE CRLO CLCULION ND RECIVIY DIREC MEUREMEN. Borio, M. Cagnazzo, F. Marchetti, P. Pappalardo,. alvini. Laboratorio Energia Nucleare pplicata (L.E.N..)
More informationHomogenization Methods for Full Core Solution of the Pn Transport Equations with 3-D Cross Sections. Andrew Hall October 16, 2015
Homogenization Methods for Full Core Solution of the Pn Transport Equations with 3-D Cross Sections Andrew Hall October 16, 2015 Outline Resource-Renewable Boiling Water Reactor (RBWR) Current Neutron
More informationNEUTRONIC CALCULATION SYSTEM FOR CANDU CORE BASED ON TRANSPORT METHODS
Romanian Reports in Physics, Vol. 63, No. 4, P. 948 960, 2011 NEUTRONIC CALCULATION SYSTEM FOR CANDU CORE BASED ON TRANSPORT METHODS V. BALACEANU 1, M. PAVELESCU 2 1 Institute for Nuclear Research, PO
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 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 information240 ETSEIB School of Industrial Engineering of Barcelona
Name of the subject: Reactor Physics and Thermal hydraulics Code: 240NU013 ECTS Credits: 7.5 Unit responsible: 240 ETSEIB School of Industrial Engineering of Barcelona Department: 721, Physics and Nuclear
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 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 informationApplicability of Point Reactor Kinetic Equations for Subcriticality Monitoring
Proceedings of Korean Nuclear Society Spring Meeting Cheji, Korea, May 21 Applicability of Point Reactor Kinetic Equations for Subcriticality Monitoring Yoichiro Shimazu, Masashi Tsuji Hokkaido University
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 1. Title: Neutron Life Cycle
Lectures on Nuclear Power Safety Lecture No 1 Title: Neutron Life Cycle Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture Infinite Multiplication Factor, k Four Factor Formula
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 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 informationRe-Evaluation of SEFOR Doppler Experiments and Analyses with JNC and ERANOS systems
PHYSOR 2004 -The Physics of Fuel Cycles and Advanced Nuclear Systems: Global evelopments Chicago, Illinois, April 25-29, 2004, on C-ROM, American Nuclear Society, Lagrange Park, IL. (2004) Re-Evaluation
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 informationNeutron reproduction. factor ε. k eff = Neutron Life Cycle. x η
Neutron reproduction factor k eff = 1.000 What is: Migration length? Critical size? How does the geometry affect the reproduction factor? x 0.9 Thermal utilization factor f x 0.9 Resonance escape probability
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 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 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 Enineerin Department, Faculty of Technoloy and Vocational Education Indonesia University of Education
More informationA Method For the Burnup Analysis of Power Reactors in Equilibrium Operation Cycles
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
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 informationDOPPLER COEFFICIENT OF REACTIVITY BENCHMARK CALCULATIONS FOR DIFFERENT ENRICHMENTS OF UO 2
Supercomputing in Nuclear Applications (M&C + SNA 2007) Monterey, California, April 15-19, 2007, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) DOPPLER COEFFICIENT OF REACTIVITY BENCHMARK
More informationSolving the neutron slowing down equation
Solving the neutron slowing down equation Bertrand Mercier, Jinghan Peng To cite this version: Bertrand Mercier, Jinghan Peng. Solving the neutron slowing down equation. 2014. HAL Id: hal-01081772
More informationInstitution : Majmaa'h University Academic Department : Physics Programme :
Institution : Majmaa'h University Academic Department : Physics Programme : Bachelor of Science in Physics Course : Neutron Physics and Reactors (PHYS485) Course Coordinator : Elassaad Jemii Programme
More informationNUCLEAR SCIENCE ACAD BASIC CURRICULUM CHAPTER 5 NEUTRON LIFE CYCLE STUDENT TEXT REV 2. L th. L f U-235 FUEL MODERATOR START CYCLE HERE THERMAL NEUTRON
ACAD BASIC CURRICULUM NUCLEAR SCIENCE CHAPTER 5 NEUTRON LIFE CYCLE 346 RESONANCE LOSSES p 038 THERMAL NEUTRON 2 THERMAL NEUTRON LEAKAGE 52 THERMAL ABSORBED BY NON-FUEL ATOMS L th 07 THERMAL f 965 THERMAL
More informationANALYSIS OF THE OECD PEACH BOTTOM TURBINE TRIP 2 TRANSIENT BENCHMARK WITH THE COUPLED NEUTRONIC AND THERMAL-HYDRAULICS CODE TRAC-M/PARCS
ANALYSIS OF THE OECD PEACH BOTTOM TURBINE TRIP 2 TRANSIENT BENCHMARK WITH THE COUPLED NEUTRONIC AND THERMAL-HYDRAULICS CODE TRAC-M/PARCS Deokjung Lee and Thomas J. Downar School of Nuclear Engineering
More informationSteady State Analysis of Small Molten Salt Reactor (Effect of Fuel Salt Flow on Reactor Characteristics)
610 Steady State Analysis of Small Molten Salt Reactor (Effect of Fuel Salt Flow on Reactor Characteristics) Takahisa YAMAMOTO,KoshiMITACHI and Takashi SUZUKI The Molten Salt Reactor (MSR) is a thermal
More informationReactivity Balance & Reactor Control System
Reactivity Balance & Reactor Control System K.S. Rajan Professor, School of Chemical & Biotechnology SASTRA University Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 6 Table of Contents 1 MULTIPLICATION
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 informationPWR CONTROL ROD EJECTION ANALYSIS WITH THE MOC CODE DECART
PWR CONTROL ROD EJECTION ANALYSIS WITH THE MOC CODE DECART Mathieu Hursin and Thomas Downar University of California Berkeley, USA mhursin@nuc.berkeley.edu,downar@nuc.berkeley.edu ABSTRACT During the past
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 informationPotential Use of beta-eff and other Benchmarks for Adjustment
Potential Use of beta-eff and other Benchmarks for Adjustment Ivo Kodeli SG39 Meeting, NEA, May 19, 2015 Analysed benchmarks from ICSBEP and IRPhE SNEAK-7A & -7B: MOX fuel reflected by metallic depleted
More informationNuclear Fission. 1/v Fast neutrons. U thermal cross sections σ fission 584 b. σ scattering 9 b. σ radiative capture 97 b.
Nuclear Fission 1/v Fast neutrons should be moderated. 235 U thermal cross sections σ fission 584 b. σ scattering 9 b. σ radiative capture 97 b. Fission Barriers 1 Nuclear Fission Q for 235 U + n 236 U
More informationDemonstration of Full PWR Core Coupled Monte Carlo Neutron Transport and Thermal-Hydraulic Simulations Using Serpent 2/ SUBCHANFLOW
Demonstration of Full PWR Core Coupled Monte Carlo Neutron Transport and Thermal-Hydraulic Simulations Using Serpent 2/ SUBCHANFLOW M. Daeubler Institute for Neutron Physics and Reactor Technology (INR)
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 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 informationTABLE OF CONTENTS CHAPTER TITLE PAGE
v TABLE OF CONTENTS CHAPTER TITLE PAGE TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES LIST OF SYMBOLS LIST OF APPENDICES v viii ix xii xiv CHAPTER 1 INTRODUCTION 1.1 Introduction 1 1.2 Literature Review
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 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 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 informationLesson 6: Diffusion Theory (cf. Transport), Applications
Lesson 6: Diffusion Theory (cf. Transport), Applications Transport Equation Diffusion Theory as Special Case Multi-zone Problems (Passive Media) Self-shielding Effects Diffusion Kernels Typical Values
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 informationApplication of the next generation of the OSCAR code system to the ETRR-2 multi-cycle depletion benchmark
Application of the next generation of the OSCAR code system to the ETRR-2 multi-cycle depletion benchmark M. Mashau 1, S.A. Groenewald 1, F.A. van Heerden 1 1) The South African Nuclear Energy Corporation
More informationA PERTURBATION ANALYSIS SCHEME IN WIMS USING TRANSPORT THEORY FLUX SOLUTIONS
A PERTURBATION ANALYSIS SCHEME IN WIMS USING TRANSPORT THEORY FLUX SOLUTIONS J G Hosking, T D Newton, B A Lindley, P J Smith and R P Hiles Amec Foster Wheeler Dorchester, Dorset, UK glynn.hosking@amecfw.com
More informationTHREE-DIMENSIONAL INTEGRAL NEUTRON TRANSPORT CELL CALCULATIONS FOR THE DETERMINATION OF MEAN CELL CROSS SECTIONS
THREE-DIMENSIONAL INTEGRAL NEUTRON TRANSPORT CELL CALCULATIONS FOR THE DETERMINATION OF MEAN CELL CROSS SECTIONS Carsten Beckert 1. Introduction To calculate the neutron transport in a reactor, it is often
More informationA Taylor series solution of the reactor point kinetics equations
1 A Taylor series solution of the reactor point kinetics equations David McMahon and Adam Pierson Department of Nuclear Safety Analysis,Sandia National Laboratories, Albuquerque, NM 87185-1141 E-mail address:
More informationStudy of Predictor-corrector methods. for Monte Carlo Burnup Codes. Dan Kotlyar Dr. Eugene Shwageraus. Supervisor
Serpent International Users Group Meeting Madrid, Spain, September 19-21, 2012 Study of Predictor-corrector methods for Monte Carlo Burnup Codes By Supervisor Dan Kotlyar Dr. Eugene Shwageraus Introduction
More informationComparison of the Monte Carlo Adjoint-Weighted and Differential Operator Perturbation Methods
Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vol., pp.836-841 (011) ARTICLE Comparison of the Monte Carlo Adjoint-Weighted and Differential Operator Perturbation Methods Brian C. KIEDROWSKI * and Forrest
More informationSodium void coefficient map by Serpent
Wir schaffen Wissen heute für morgen Paul Scheer Institut 4th Annual Serpent Users Group Meeting (Cambridge Sept 17-19, 2014): Sandro Pelloni Sodium void coefficient map by Serpent PSI, 8. September 2014
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 information"Control Rod Calibration"
TECHNICAL UNIVERSITY DRESDEN Institute of Power Engineering Training Reactor Reactor Training Course Experiment "Control Rod Calibration" Instruction for Experiment Control Rod Calibration Content: 1...
More informationUSA HTR NEUTRONIC CHARACTERIZATION OF THE SAFARI-1 MATERIAL TESTING REACTOR
Proceedings of HTR2008 4 th International Topical Meeting on High Temperature Reactors September 28-October 1, 2008, Washington, D.C, USA HTR2008-58155 NEUTRONIC CHARACTERIZATION OF THE SAFARI-1 MATERIAL
More informationPhD Qualifying Exam Nuclear Engineering Program. Part 1 Core Courses
PhD Qualifying Exam Nuclear Engineering Program Part 1 Core Courses 9:00 am 12:00 noon, November 19, 2016 (1) Nuclear Reactor Analysis During the startup of a one-region, homogeneous slab reactor of size
More informationControl Rod Calibration and Worth Calculation for Optimized Power Reactor 1000 (OPR-1000) Using Core Simulator OPR1000
World Journal of Nuclear Science and Technology, 017, 7, 15-3 http://www.scirp.org/journal/wjnst ISSN Online: 161-6809 ISSN Print: 161-6795 Control Rod Calibration and Worth Calculation for Optimized Power
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 informationSymmetry in Monte Carlo. Dennis Mennerdahl OECD/NEA/NSC/WPNCS/AMCT EG, Paris, 18 September 2014
Symmetry in Monte Carlo Dennis Mennerdahl OECD/NEA/NSC/WPNCS/AMCT EG, Paris, 18 September 2014 OVERVIEW Identical events - Full model results contain everything and more Symmetry to improve convergence?
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 informationA.BIDAUD, I. KODELI, V.MASTRANGELO, E.SARTORI
SENSITIVITY TO NUCLEAR DATA AND UNCERTAINTY ANALYSIS: THE EXPERIENCE OF VENUS2 OECD/NEA BENCHMARKS. A.BIDAUD, I. KODELI, V.MASTRANGELO, E.SARTORI IPN Orsay CNAM PARIS OECD/NEA Data Bank, Issy les moulineaux
More informationNeutronic analysis of SFR lattices: Serpent vs. HELIOS-2
Neutronic analysis of SFR lattices: Serpent vs. HELIOS-2 E. Fridman 1, R. Rachamin 1, C. Wemple 2 1 Helmholtz Zentrum Dresden Rossendorf 2 Studsvik Scandpower Inc. Text optional: Institutsname Prof. Dr.
More informationEfficient Energy Conversion of the 14MeV Neutrons in DT Inertial Confinement Fusion. By F. Winterberg University of Nevada, Reno
Efficient Energy Conversion of the 14MeV Neutrons in DT Inertial Confinement Fusion By F. Winterberg University of Nevada, Reno Abstract In DT fusion 80% of the energy released goes into 14MeV neutrons,
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 informationAn alternative experimental approach for subcritical configurations of the IPEN/MB-01 nuclear reactor
Journal of Physics: Conference Series PAPER OPEN ACCESS An alternative experimental approach for subcritical configurations of the IPEN/MB-01 nuclear reactor To cite this article: E Gonnelli et al 2015
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 informationLecture 18 Neutron Kinetics Equations
24.505 Lecture 18 Neutron Kinetics Equations Prof. Dean Wang For a nuclear reactor to operate at a constant power level, the rate of neutron production via fission reactions should be exactly balanced
More informationPassive reactivity control with 10 B burnable poison in the U-battery
Passive reactivity control with B burnable poison in the U-battery Simon de Die April 2, 27 Foreword This thesis represents my graduation thesis for the Master in Applied Physics at Delft 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 informationExtension of the MCBEND Monte Carlo Code to Perform Adjoint Calculations using Point Energy Data
Extension of the MCBEND Monte Carlo Code to Perform Adjoint Calculations using Point Energy Data Malcolm Grimstone Abstract In radiation transport calculations there are many situations where the adjoint
More informationOECD/NEA Transient Benchmark Analysis with PARCS - THERMIX
OECD/NEA Transient Benchmark Analysis with PARCS - THERMIX Volkan Seker Thomas J. Downar OECD/NEA PBMR Workshop Paris, France June 16, 2005 Introduction Motivation of the benchmark Code-to-code comparisons.
More informationComparison of PWR burnup calculations with SCALE 5.0/TRITON other burnup codes and experimental results. Abstract
Comparison of PWR burnup calculations with SCALE 5.0/TRITON other burnup codes and experimental results Ph.Oberle, C.H.M.Broeders, R.Dagan Forschungszentrum Karlsruhe, Institut for Reactor Safety Hermann-von-Helmholtz-Platz-1,
More informationModeling of the Multi-SERTTA Experiment with MAMMOTH
INL/MIS-17-43729 INL/MIS-16-40269 Approved for for public release; distribution is is unlimited. Modeling of the Multi-SERTTA Experiment with MAMMOTH Javier Ortensi, Ph.D. P.E. R&D Scientist Nuclear Science
More informationComputational and Experimental Benchmarking for Transient Fuel Testing: Neutronics Tasks
Computational and Experimental Benchmarking for Transient Fuel Testing: Neutronics Tasks T. Downar W. Martin University of Michigan C. Lee Argonne National Laboratory November 19, 2015 Objective of Neutronics
More informationNucE 497A RAMP Class #1
1 COURSE OBJECTIVES This course is designed as an intensive course providing an introduction to nuclear engineering (NucE) for graduate students with non-nuce background and to returning students. After
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 information