Study of Predictor-corrector methods. for Monte Carlo Burnup Codes. Dan Kotlyar Dr. Eugene Shwageraus. Supervisor
|
|
- Silvester Hunt
- 6 years ago
- Views:
Transcription
1 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
2 Introduction High Conversion LWRs Eigen value D (no T-H) Time, days 1
3 Introduction High Conversion LWRs Eigen value D (T-H) 2D (no T-H) Time, days 1
4 Introduction Advisor s response Power share in seed What? T [i] T[i+1] What? Why? Programming? Physical (model)? Numerical (coupling)? Height, cm Cause How? Solution! Is this problem common for all MC based codes? 2
5 Outline Simplify the model Generic PWR assembly in 3D No T-H feedback, but non-uniform coolant density fixed in time Study efficient schemes for coupling MC with multi-physics feedbacks Identify origins of the observed problem Focus of this work: Examine MC-burnup coupling schemes (& t size) Was this problem identified earlier? Recent burnup schemes sensitivity studies Description of the test case (3D) & Results Conclusions & future work 3
6 Introduction Examples of Integrated MC-burnup codes: MOCUP (Moore et al., 1995) MCODE (Xu et al., 2002) MONTEBURNS (Trellue, 2003) SERPENT (Leppänen, 2007) BGCore (Fridman et al., 2008) MCNPX (Hendricks et al., 2008) What about Integration of MC-burnup-TH BGCore SERPENT Issues related to Integrated MC codes Accuracy Computation requirements 4
7 Objectives Investigate stability and accuracy of current depletion calculations Monte-Carlo/Burnup coupling scheme Depletion timestep size Commonly used methods Explicit Euler predictor method Euler predictor-corrector method Propose an extension to the modified Euler predictor-corrector Accuracy as the predictor-corrector method Computation requirements as the predictor method 5
8 Possible multi-physics coupling schemes Integration scheme 1 Beginning of Step analysis (predictor) Inner loop: TH-neutronic Total exe time (Depletion/TH) n T mcnp Outer loop: Depletion Integration scheme 2 Inner loop: Depletion Outer loop: TH-neutronuc Important issues Combination of #1 and #2 Calculation time Source convergence Distribution of errors Variable convergence tolerances Neutronics MC (r, E) Depletion module New concentration for next step Temp & density Distribution Thermal feedback Power & BU Distribution (x n) (x 1) 6
9 Recent sensitivity studies Yamamoto et al., 2008 Projected predictor corrector method Linear correlation between the number density and the microscopic RR Tested for the Gd-bearing 2D fuel assembly Carpenter et al., 2010 (Bettis Atomic Power Laboratory) Modified Log Linear correlation (Yamamoto) Tested for the Gd-bearing 2D fuel assembly Isotalo et al., 2011 Higher order methods Use more (previous) BU points Tested on PWR pin cell and seed/blanket assembly Saadi et al., 2012 Burnup sensitivity analysis Tested on 1D UO2 PWR unit cell Dufek et al., 2009 Numerical stability of MC-burnup codes Instability demonstrated on infinitely reflected 3D unit cell 7
10 Burnup coupling scheme methodology Formulation of the burnup problem Explicit Euler Predictor Euler Predictor-Corrector Extended Predictor-Corrector
11 Burnup calculations n n+1 n+2 EOC Eigen-value transport problem Burnup equation Matrix Exponential solution L N r, t FN r, t r,, E, t0 dn r, t dt N 0 1 k r, E, tede Nr, t r N r t t n n n1 1 exp n 8
12 Burnup calculations n n+1 n+2 EOC Eigen-value equation Burnup equation L N r, t FN r, t r,, E, t0 dn r, t dt 0 1 k r, E, tede Nr, t Mat. Exp. solution N r N r t t n n n1 1 exp n Objective: coupled space-energy-time dependent solution: Coupling scheme: Independent neutronics and depletion solvers 8
13 Euler predictor-corrector methods 1. Explicit Euler Predictor method Reaction rate (RR) calculation (neutronics - ) at BOT Depletion with BOT RR s N n+1 N n const RR n N n+1 n n+1 const RR n+1 Time/BU 2. Modified Euler Predictor-Corrector method Reaction rate calculation ( ) with predicted N n+1 values Depletion with EOT (n+1) reaction rates to obtain corrected N n+1 values Finally, the predicted and corrected ND are averaged Used as initial values for the following step 9
14 Why is the P-C method is not enough? t! Power Share Normalized commulative energy (50 days) Quote: Predictor-corrector methods are numerically explicit The stable size may not be much larger than of that of the predictor Numerical methods for engineers and scientists, Joe D. Hoffman 0 50 d t E i t 0,50 P t t Predictor Corrector 50 days d 10 d 25 d 50 d Height, cm Height, cm 10
15 Suggested E-P-C method 3. Extended Predictor-Corrector method (E-P-C) Assume (power, RR, σ, N) are known for the previous [t n-1,t n ] interval Deplete with previous average RR N n+0.5 ( =N n exp[-rr n-1 t 0.5 ] ) Update transport solution (RR calculation ) at MOT (t n+0.5 ) Deplete with MOT RR s N n+1 Reaction rate calculation ( ) with predicted N n+1 values Depletion with EOT (n+1) reaction rates to obtain corrected N n+1 values Finally, the predicted and corrected ND are averaged N n-1 N n N n+0.5 const RR n+0.5 N n+1 n-1 75d n 100d n d const RR n+1 n+1 150d Time/BU 11
16 Different coupling scheme results
17 3D assembly test case Radial +20,000 layout Cylindrical pins With guide tubes Axial layout / Coolant density profile 21.5 cm Coolant Density, gr/cm 3 12
18 BGCore vs. SERPENT comparison Concentration of Pu239, #/b cm Eigen value Parameter k-eff Maximum difference, SERPENT vs. BGCore ~70 pcm Xe 135 ~0.5 % ~0.4 % BGCore SERPENT U 235 ~0.4 % Pu E E E-05 BGCore SERPENT Time, days 7.0E E E E E-05 Method : predictor-corrector Timestep : 5 days 2.0E E E Time, days 13
19 Burnup coupling scheme method & T Case no Method Predictor Predictor-Corrector a Predictor-Corrector Predictor-Corrector Description Extended-Predictor-Corrector a. Chosen to be the reference case T, days Designated P-C-0 (25d) P-C-1 ( 5d) P-C-1 (25d) P-C-1 (50d) E-P-C (50d) 14
20 Difference in eigen value, pcm Burnup coupling scheme method & T Case no Method Predictor Predictor-Corrector a Predictor-Corrector Predictor-Corrector Description Extended-Predictor-Corrector a. Chosen to be the reference case T, days Designated P-C-0 (25d) P-C-1 ( 5d) P-C-1 (25d) P-C-1 (50d) E-P-C (50d) P-C-0 (25d) P-C-1 (25d) P-C-1 (50d) Time, days 14
21 Concentration of Xe-135 Difference in eigen value, pcm Burnup coupling scheme method & T Case no. Method Description T, days Designated P-C-0 (25d) P-C-1 (25d) P-C-1 (50d) 1 Predictor 25 P-C-0 (25d) Predictor-Corrector a 5 P-C-1 ( 5d) Predictor-Corrector 25 P-C-1 (25d) Predictor-Corrector 50 P-C-1 (50d) 0 5 Extended-Predictor-Corrector 50 E-P-C (50d) Time, days a. Chosen to be the reference case 1.2E E E E E E E+00 P-C-0 (25d) P-C-1 (25d) P-C-1 (50d) P-C-1 (5d) Time, days 14
22 Flux distribution, n/s cm2 Concentration of Xe-135 Difference in eigen value, pcm Burnup coupling scheme method & T Case no. Method Description T, days Designated P-C-0 (25d) P-C-1 (25d) P-C-1 (50d) 1 Predictor 25 P-C-0 (25d) Predictor-Corrector a 5 P-C-1 ( 5d) Predictor-Corrector 25 P-C-1 (25d) Predictor-Corrector 50 P-C-1 (50d) 0 5 Extended-Predictor-Corrector 50 E-P-C (50d) Time, days a. Chosen to be the reference case 4.0E E+14 T(i) T(i+1) 1.2E E E E E E E E E Height, cm 8.0E E E E E+00 P-C-0 (25d) P-C-1 (25d) P-C-1 (50d) P-C-1 (5d) Time, days 14
23 Eigen value The effect of non-uniform coolant density Eigen value Realistic coolant density along the z-axis Averaged coolant density along the z-axis P-C-0 (25d) P-C-1 (25d) P-C-1 (50d) P-C-0 P-C Time, days Time, days 15
24 k-eff Comparison between P-C and E-P-C Difference [pcm] to P-C-1-25d Concentration of U-235 Relative difference to P-C-1-25d Maximum difference (compared to P-C-1) in different core parameters Parameter k-eff Maximum difference, P-C-1 (50d) vs. P-C-1 (25 d) ~500 pcm Maximum difference, E-P-C (50d) vs. P-C-1 (25 d) ~20 pcm Xe 135 ~2.0 % ~0.5 % U 235 ~1.9 % ~0.1 % Pu 239 ~2.5 % ~0.2 % P-C-1 (25d) P-C-1 (50d) E-P-C (50d) DIFF (P-C-1) DIFF (E-P-C) E E E E % 0.5% 0.0% -0.5% E % Time, days 6.0E E E E-04 P-C-1 (25d) P-C-1 (50d) E-P-C (50d) DIFF (P-C-1) DIFF (E-P-C) Time, days -1.5% -2.0% -2.5% -3.0% 16
25 Conclusions & future work
26 Summary Sensitivity studies of different coupling schemes were performed Explicit Euler predictor method may be unstable for some 3D problems Predictor-corrector method with Δt ~50d does not resolve the problem EPC method allows using larger timesteps Maintains sufficient accuracy Future work: How/When is the Thermal-hydraulic feedback should be applied? 17
27 Thank You for your attention!
FULL CORE POWER AND ISOTOPIC OSCILLATIONS WITH VARIOUS DEPLETION SCHEMES
FULL CORE POWER AND ISOTOPIC OSCILLATIONS WITH VARIOUS DEPLETION SCHEMES A N D R E W J O H N S O N C O M P U TAT I O N A L R E A C T O R E N G I N E E R I N G OUTLINE Governing depletion equations Summary
More informationPreventing xenon oscillations in Monte Carlo burnup calculations by forcing equilibrium
Preventing xenon oscillations in Monte Carlo burnup calculations by forcing equilibrium Aarno Isotaloa), Jaakko Leppänenb), Jan Dufekcc) a) Aalto University, Finland b) VTT Technical Research Centrte of
More informationImproved time integration methods for burnup calculations with Monte Carlo neutronics
Improved time integration methods for burnup calculations with Monte Carlo neutronics Aarno Isotalo 13.4.2010 Burnup calculations Solving time development of reactor core parameters Nuclide inventory,
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 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 informationBURNUP CALCULATION CAPABILITY IN THE PSG2 / SERPENT MONTE CARLO REACTOR PHYSICS CODE
International Conference on Mathematics, Computational Methods & Reactor Physics (M&C 29) Saratoga Springs, New York, May 3-7, 29, on CD-ROM, American Nuclear Society, LaGrange Park, IL (29) BURNUP CALCULATION
More informationConsistent Code-to-Code Comparison of Pin-cell Depletion Benchmark Suite
Consistent Code-to-Code Comparison of Pin-cell Depletion Benchmark Suite September 27, 26 Jinsu Park, Deokjung Lee * COmputational Reactor Physics & Experiment lab Contents VERA depletion benchmark suite
More informationOn-the-fly Doppler Broadening in Serpent
On-the-fly Doppler Broadening in Serpent 1st International Serpent User Group Meeting 16.9.2011, Dresden Tuomas Viitanen VTT Technical Research Centre of Finland Outline Fuel temperatures in neutronics
More informationTMS On-the-fly Temperature Treatment in Serpent
TMS On-the-fly Temperature Treatment in Serpent Tuomas Viitanen & Jaakko Leppänen Serpent User Group Meeting, Cambridge, UK September 17 19, 2014 Effects of thermal motion on neutron transport On reaction
More informationSerpent Monte Carlo Neutron Transport Code
Serpent Monte Carlo Neutron Transport Code NEA Expert Group on Advanced Monte Carlo Techniques, Meeting September 17 2012 Jaakko Leppänen / Tuomas Viitanen VTT Technical Research Centre of Finland Outline
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 informationUSE OF LATTICE CODE DRAGON IN REACTOR CALUCLATIONS
USE OF LATTICE CODE DRAGON IN REACTOR CALUCLATIONS ABSTRACT Dušan Ćalić ZEL-EN razvojni center Hočevarjev trg 1 Slovenia-SI8270, Krško, Slovenia dusan.calic@zel-en.si Andrej Trkov, Marjan Kromar J. Stefan
More informationUse of Monte Carlo and Deterministic Codes for Calculation of Plutonium Radial Distribution in a Fuel Cell
Use of Monte Carlo and Deterministic Codes for Calculation of Plutonium Radial Distribution in a Fuel Cell Dušan Ćalić, Marjan Kromar, Andrej Trkov Jožef Stefan Institute Jamova 39, SI-1000 Ljubljana,
More informationAssessment of the MCNP-ACAB code system for burnup credit analyses
Assessment of the MCNP-ACAB code system for burnup credit analyses N. García-Herranz, O. Cabellos, J. Sanz UPM - UNED International Workshop on Advances in Applications of Burnup Credit for Spent Fuel
More informationMOx Benchmark Calculations by Deterministic and Monte Carlo Codes
MOx Benchmark Calculations by Deterministic and Monte Carlo Codes G.Kotev, M. Pecchia C. Parisi, F. D Auria San Piero a Grado Nuclear Research Group (GRNSPG), University of Pisa via Diotisalvi 2, 56122
More informationIncineration of Plutonium in PWR Using Hydride Fuel
Incineration of Plutonium in PWR Using Hydride Fuel Francesco Ganda and Ehud Greenspan University of California, Berkeley ARWIF-2005 Oak-Ridge, TN February 16-18, 2005 Pu transmutation overview Many approaches
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 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 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 informationCriticality analysis of ALLEGRO Fuel Assemblies Configurations
Criticality analysis of ALLEGRO Fuel Assemblies Configurations Radoslav ZAJAC Vladimír CHRAPČIAK 13-16 October 2015 5th International Serpent User Group Meeting at Knoxville, Tennessee ALLEGRO Core - Fuel
More informationSolving Burnup Equations in Serpent:
Solving Burnup Equations in Serpent: Matrix Exponential Method CRAM Maria Pusa September 15th, 2011 Outline Burnup equations and matrix exponential solution Characteristics of burnup matrices Established
More informationDETERMINATION OF THE EQUILIBRIUM COMPOSITION OF CORES WITH CONTINUOUS FUEL FEED AND REMOVAL USING MOCUP
Supercomputing in Nuclear Applications (M&C + SNA 2007) Monterey, California, April 15-19, 2007, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) DETERMINATION OF THE EQUILIBRIUM COMPOSITION
More informationHigher-order Chebyshev Rational Approximation Method (CRAM) and application to Burnup Calculations
Higher-order Chebyshev Rational Approximation Method (CRAM) and application to Burnup Calculations Maria Pusa September 18th 2014 1 Outline Burnup equations and matrix exponential solution Characteristics
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 informationResearch Article Calculations for a BWR Lattice with Adjacent Gadolinium Pins Using the Monte Carlo Cell Code Serpent v.1.1.7
Science and Technology of Nuclear Installations Volume 2, Article ID 65946, 4 pages doi:.55/2/65946 Research Article Calculations for a BWR Lattice with Adjacent Gadolinium Pins Using the Monte Carlo Cell
More informationQUADRATIC DEPLETION MODEL FOR GADOLINIUM ISOTOPES IN CASMO-5
Advances in Nuclear Fuel Management IV (ANFM 2009) Hilton Head Island, South Carolina, USA, April 12-15, 2009, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2009) QUADRATIC DEPLETION MODEL FOR
More informationNew Capabilities for the Chebyshev Rational Approximation method (CRAM)
New Capabilities for the Chebyshev Rational Approximation method (CRAM) A. Isotaloa,b W. Wieselquista M. Pusac aoak Ridge National Laboratory PO Box 2008, Oak Ridge, TN 37831-6172, USA baalto University
More informationSome remarks on XS preparation with SERPENT
Some remarks on XS preparation wit SERPENT E. Fridman Text optional: Institutsname Prof. Dr. Hans Mustermann www.fzd.de Mitlied der Leibniz-Gemeinscaft Outline Leakae-corrected omoenized XS XS preparation
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 informationCharacterization of waste by R2S methodology: SEACAB system. Candan Töre 25/11/2017, RADKOR2017, ANKARA
Characterization of waste by R2S methodology: SEACAB system Candan Töre 25/11/2017, RADKOR2017, ANKARA SEA Ingeniería y Análisis de Blindajes Avda. de Atenas, 75, 106-107 28230 LAS ROZAS (Madrid) Tel:
More informationCASMO-5/5M Code and Library Status. J. Rhodes, K. Smith, D. Lee, Z. Xu, & N. Gheorghiu Arizona 2008
CASMO-5/5M Code and Library Status J. Rhodes, K. Smith, D. Lee, Z. Xu, & N. Gheorghiu Arizona 2008 CASMO Methodolgy Evolution CASMO-3 Homo. transmission probability/external Gd depletion CASMO-4 up to
More informationTOWARDS A COUPLED SIMULATION OF THERMAL HYDRAULICS AND NEUTRONICS IN A SIMPLIFIED PWR WITH A 3x3 PIN ASSEMBLY
TOWARDS A COUPLED SIMULATION OF THERMAL HYDRAULICS AND NEUTRONICS IN A SIMPLIFIED PWR WITH A 3x3 PIN ASSEMBLY Anni Schulze, Hans-Josef Allelein Institute for Reactor Safety and Reactor Technology, RWTH
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 informationTemperature treatment capabilites in Serpent 2(.1.24)
Temperature treatment capabilites in Serpent 2(.1.24) Tuomas Viitanen Serpent User Group Meeting Knoxville, TN, USA OCt. 13 16, 2015 Introduction Temperature affects the neutron transport through: - Reaction
More informationA Deterministic against Monte-Carlo Depletion Calculation Benchmark for JHR Core Configurations. A. Chambon, P. Vinai, C.
A Deterministic against Monte-Carlo Depletion Calculation Benchmark for JHR Core Configurations A. Chambon, P. Vinai, C. Demazière Chalmers University of Technology, Department of Physics, SE-412 96 Gothenburg,
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 informationClick to edit Master title style
Automated calculation sequence for group constant generation in Serpent 4th International Serpent UGM, Cambridge, UK, Sept. 17-19, 014 Jaakko Leppänen VTT Technical Research Center of Finland Click to
More information2017 Water Reactor Fuel Performance Meeting September 10 (Sun) ~ 14 (Thu), 2017 Ramada Plaza Jeju Jeju Island, Korea
Neutronic evaluation of thorium-uranium fuel in heavy water research reactor HADI SHAMORADIFAR 1,*, BEHZAD TEIMURI 2, PARVIZ PARVARESH 1, SAEED MOHAMMADI 1 1 Department of Nuclear physics, Payame Noor
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 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 informationPHYS-E0562 Ydinenergiatekniikan jatkokurssi Lecture 5 Burnup calculation
PHYS-E0562 Ydinenergiatekniikan jatkokurssi Lecture 5 Burnup calculation Jaakko Leppänen (Lecturer), Ville Valtavirta (Assistant) Department of Applied Physics Aalto University, School of Science Jaakko.Leppanen@aalto.fi
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 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 informationAdvanced Methods Development for Equilibrium Cycle Calculations of the RBWR. Andrew Hall 11/7/2013
Advanced Methods Development for Equilibrium Cycle Calculations of the RBWR Andrew Hall 11/7/2013 Outline RBWR Motivation and Desin Why use Serpent Cross Sections? Modelin the RBWR Axial Discontinuity
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 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 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 informationEvaluation of RAPID for a UNF cask benchmark problem
Evaluation of RAPID for a UNF cask benchmark problem Valerio Mascolino, Alireza Haghighat, and Nathan Roskoff Virginia Tech Nuclear Science and Engineering Laboratory Nuclear Engineering Program, Mechanical
More informationCASMO-5 Development and Applications. Abstract
Organized and hosted by the Canadian Nuclear Society. Vancouver, BC, Canada. 2006 September 10-14 CASMO-5 Development and Applications Joel Rhodes *1, Kord Smith 1, and Deokjung Lee 1 1 Studsvik Scandpower
More informationIMPACT OF THE FISSION YIELD COVARIANCE DATA IN BURN-UP CALCULATIONS
IMPACT OF THE FISSION YIELD COVARIANCE DATA IN BRN-P CALCLATIONS O. Cabellos, D. Piedra, Carlos J. Diez Department of Nuclear Engineering, niversidad Politécnica de Madrid, Spain E-mail: oscar.cabellos@upm.es
More informationTechnical workshop : Dynamic nuclear fuel cycle
Technical workshop : Dynamic nuclear fuel cycle Reactor description in CLASS Baptiste LENIAU* Institut d Astrophysique de Paris 6-8 July, 2016 Introduction Summary Summary The CLASS package : a brief overview
More informationDissolution, Reactor, and Environmental Behavior of ZrO2-MgO Inert Fuel Matrix: Neutronic Evaluation of ZrO2-MgO Inert Fuels
Fuels Campaign (RP) ransmutation Research Program Projects 1-2006 Dissolution, Reactor, and Environmental Behavior of ZrO2-MgO Inert Fuel Matrix: Neutronic Evaluation of ZrO2-MgO Inert Fuels E. Fridman
More informationThe Effect of Burnup on Reactivity for VVER-1000 with MOXGD and UGD Fuel Assemblies Using MCNPX Code
Journal of Nuclear and Particle Physics 2016, 6(3): 61-71 DOI: 10.5923/j.jnpp.20160603.03 The Effect of Burnup on Reactivity for VVER-1000 with MOXGD and UGD Fuel Assemblies Using MCNPX Code Heba K. Louis
More informationSerpent Neutronic Calculations of Nuclear Thermal Propulsion Reactors
Serpent Neutronic Calculations of Nuclear Thermal Propulsion Reactors Paolo Venneri 1,2, Yonghee Kim 1 Reactor Core Design and Transmutation, KAIST 1 Advanced Systems, Ultra Safe Nuclear Corp 2. Presented
More informationDevelopment of Multigroup Cross Section Generation Code MC 2-3 for Fast Reactor Analysis
Development o Multigroup Cross Section Generation Code MC 2-3 or Fast Reactor Analysis International Conerence on Fast Reactors and Related Fuel Cycles December 7-11, 2009 Kyoto, Japan Changho Lee and
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 informationInvestigation of Nuclear Data Accuracy for the Accelerator- Driven System with Minor Actinide Fuel
Investigation of Nuclear Data Accuracy for the Accelerator- Driven System with Minor Actinide Fuel Kenji Nishihara, Takanori Sugawara, Hiroki Iwamoto JAEA, Japan Francisco Alvarez Velarde CIEMAT, Spain
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 informationNuclear data sensitivity and uncertainty assessment of sodium voiding reactivity coefficients of an ASTRID-like Sodium Fast Reactor
Nuclear data sensitivity and uncertainty assessment of sodium voiding reactivity coefficients of an ASTRID-like Sodium Fast Reactor García-Herranz Nuria 1,*, Panadero Anne-Laurène 2, Martinez Ana 1, Pelloni
More informationPROPAGATION OF NUCLEAR DATA UNCERTAINTIES IN FUEL CYCLE USING MONTE-CARLO TECHNIQUE
PROPAGATION OF NUCLEAR DATA UNCERTAINTIES IN FUEL CYCLE CALCULATIONS USING MONTE-CARLO TECHNIQUE C.J. Díez (1), O. Cabellos (1), J.S. Martínez (1) (1) Universidad Politécnica de Madrid (UPM) International
More informationIntroduction. Problem Summary. Attila Code-to-Code Comparison Subcritical Spent Nuclear Fuel Canister with Primary Neutron and Gamma Sources
Introduction An example spent nuclear fuel canister calculation was performed with Attila, with results compared to the MCNPX Monte Carlo code. The canister contained representative subcritical neutron
More informationTHOMAS P. GERRITY III. June Signature of Author: Departmeff-of Nuclear Scieince and Engineering May 18, 2012
MCODE-3: TIME-DEPENDENT DEPLETION ISOTOPICS WITH MCNP-5 AND SCALE-6.1 by THOMAS P. GERRITY III a R1 ES Submitted to the Department of Nuclear Science and Engineering in Partial Fulfillment of the Requirements
More informationMalcolm Bean AT THE MAY All Rights Reserved. Signature of Author: Malcolm Bean Department of Nuclear Science and Engineering
COMPUTATIONAL NEUTRONICS ANALYSIS OF TRIGA REACTORS DURING POWER PULSING ARCHIIVE By Malcolm Bean SUBMITTED TO THE DEPARTMENT OF NUCLEAR SCIENCE AND ENGINEERING IN PARTIAL FULFILLMENT OF THE REQUIREMENT
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 informationAEGIS: AN ADVANCED LATTICE PHYSICS CODE FOR LIGHT WATER REACTOR ANALYSES
AEGIS: AN ADVANCED LATTICE PHYSICS CODE FOR LIGHT WATER REACTOR ANALYSES AKIO YAMAMOTO *, 1, TOMOHIRO ENDO 1, MASATO TABUCHI 2, NAOKI SUGIMURA 2, TADASHI USHIO 2, MASAAKI MORI 2, MASAHIRO TATSUMI 3 and
More informationImproved PWR Simulations by Monte-Carlo Uncertainty Analysis and Bayesian Inference
Improved PWR Simulations by Monte-Carlo Uncertainty Analysis and Bayesian Inference E. Castro, O. Buss, A. Hoefer PEPA1-G: Radiology & Criticality, AREVA GmbH, Germany Universidad Politécnica de Madrid
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 informationKord Smith Art DiGiovine Dan Hagrman Scott Palmtag. Studsvik Scandpower. CASMO User s Group May 2003
Kord Smith Art DiGiovine Dan Hagrman Scott Palmtag CASMO User s Group May 2003 TFU-related data is required input for: -CASMO-4 - SIMULATE-3 - SIMULATE-3K and SIMULATE-3R (implicit in XIMAGE and GARDEL)
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 informationNew methods implemented in TRIPOLI-4. New methods implemented in TRIPOLI-4. J. Eduard Hoogenboom Delft University of Technology
New methods implemented in TRIPOLI-4 New methods implemented in TRIPOLI-4 J. Eduard Hoogenboom Delft University of Technology on behalf of Cheikh Diop (WP1.1 leader) and all other contributors to WP1.1
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 informationNonlinear Iterative Solution of the Neutron Transport Equation
Nonlinear Iterative Solution of the Neutron Transport Equation Emiliano Masiello Commissariat à l Energie Atomique de Saclay /DANS//SERMA/LTSD emiliano.masiello@cea.fr 1/37 Outline - motivations and framework
More informationDevelopment of 3D Space Time Kinetics Model for Coupled Neutron Kinetics and Thermal hydraulics
Development of 3D Space Time Kinetics Model for Coupled Neutron Kinetics and Thermal hydraulics WORKSHOP ON ADVANCED CODE SUITE FOR DESIGN, SAFETY ANALYSIS AND OPERATION OF HEAVY WATER REACTORS October
More informationTHE NEXT GENERATION WIMS LATTICE CODE : WIMS9
THE NEXT GENERATION WIMS LATTICE CODE : WIMS9 T D Newton and J L Hutton Serco Assurance Winfrith Technology Centre Dorchester Dorset DT2 8ZE United Kingdom tim.newton@sercoassurance.com ABSTRACT The WIMS8
More informationThe Pennsylvania State University. The Graduate School HIGH ACCURACY MODELING FOR ADVANCED NUCLEAR REACTOR
The Pennsylvania State University The Graduate School Department of Mechanical and Nuclear Engineering HIGH ACCURACY MODELING FOR ADVANCED NUCLEAR REACTOR CORE DESIGNS USING MONTE CARLO BASED COUPLED CALCULATIONS
More informationA Hybrid Stochastic Deterministic Approach for Full Core Neutronics Seyed Rida Housseiny Milany, Guy Marleau
A Hybrid Stochastic Deterministic Approach for Full Core Neutronics Seyed Rida Housseiny Milany, Guy Marleau Institute of Nuclear Engineering, Ecole Polytechnique de Montreal, C.P. 6079 succ Centre-Ville,
More informationSafety Analyses for Dynamical Events (SADE) SAFIR2018 Interim Seminar Ville Sahlberg
VTT TECHNICAL RESEARCH CENTRE OF FINLAND LTD Safety Analyses for Dynamical Events (SADE) SAFIR2018 Interim Seminar Ville Sahlberg SADE Project SADE produces more reliable answers to safety requirements
More informationThe Pennsylvania State University. The Graduate School. Department of Mechanical and Nuclear Engineering
The Pennsylvania State University The Graduate School Department of Mechanical and Nuclear Engineering THE DEVELOPMENT OF A THERMAL HYDRAULIC FEEDBACK MECHANISM WITH A QUASI-FIXED POINT ITERATION SCHEME
More informationNEUTRONIC ANALYSIS STUDIES OF THE SPALLATION TARGET WINDOW FOR A GAS COOLED ADS CONCEPT.
NEUTRONIC ANALYSIS STUDIES OF THE SPALLATION TARGET WINDOW FOR A GAS COOLED ADS CONCEPT. A. Abánades, A. Blanco, A. Burgos, S. Cuesta, P.T. León, J. M. Martínez-Val, M. Perlado Universidad Politecnica
More informationA PWR ASSEMBLY COMPUTATIONAL SCHEME BASED ON THE DRAGON V4 LATTICE CODE
Supercomputing in Nuclear Applications (M&C + SNA 2007) Monterey, California, April 15-19, 2007, on CD-ROM, American Nuclear Society, LaGrange Park, IL (2007) A PWR ASSEMBLY COMPUTATIONAL SCHEME BASED
More informationStatus of MORET5 source convergence improvements and benchmark proposal for Monte Carlo depletion calculations
Status of MORET5 source convergence improvements and benchmark proposal for Monte Carlo depletion calculations Y. Richet ; W. Haeck ; J. Miss Criticality analysis department Study, Research, Codes Development
More informationImpact of the Hypothetical RCCA Rodlet Separation on the Nuclear Parameters of the NPP Krško core
International Conference Nuclear Energy for New Europe 2005 Bled, Slovenia, September 5-8, 2005 Impact of the Hypothetical RCCA Rodlet Separation on the Nuclear Parameters of the NPP Krško core ABSTRACT
More informationFirst ANDES annual meeting
First ANDES Annual meeting 3-5 May 011 CIEMAT, Madrid, Spain 1 / 0 *C.J. Díez e-mail: cj.diez@upm.es carlosjavier@denim.upm.es UNCERTAINTY METHODS IN ACTIVATION AND INVENTORY CALCULATIONS Carlos J. Díez*,
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 informationGB5 A LINKING CODE BETWEEN MCNP5 AND ORIGEN2.1 DEN/UFMG VERSION: INCLUDING RADIOACTIVITY HAZARD.
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 informationComparison of Silicon Carbide and Zircaloy4 Cladding during LBLOCA
Comparison of Silicon Carbide and Zircaloy4 Cladding during LBLOCA Prepared By: Kwangwon Ahn Prepared For: 22.314 Prepared On: December 7 th, 2006 Comparison of Silicon Carbide and Zircaloy4 Cladding during
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 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 informationOn-The-Fly Neutron Doppler Broadening for MCNP"
LA-UR-12-00700" 2012-03-26! On-The-Fly Neutron Doppler Broadening for MCNP" Forrest Brown 1, William Martin 2, " Gokhan Yesilyurt 3, Scott Wilderman 2" 1 Monte Carlo Methods (XCP-3), LANL" 2 University
More informationBenchmark Calculation of KRITZ-2 by DRAGON/PARCS. M. Choi, H. Choi, R. Hon
Benchmark Calculation of KRITZ-2 by DRAGON/PARCS M. Choi, H. Choi, R. Hon General Atomics: 3550 General Atomics Court, San Diego, CA 92121, USA, Hangbok.Choi@ga.com Abstract - Benchmark calculations have
More informationAnalysis of the Neutronic Characteristics of GFR-2400 Fast Reactor Using MCNPX Transport Code
Amr Ibrahim, et al. Arab J. Nucl. Sci. Appl, Vol 51, 1, 177-188 The Egyptian Arab Journal of Nuclear Sciences and Applications (2018) Society of Nuclear Vol 51, 1, (177-188) 2018 Sciences and Applications
More informationCzech Technical University in Prague. Faculty of Nuclear Sciences and Physical Engineering. Department of Nuclear Reactors
Czech Technical University in Prague Faculty of Nuclear Sciences and Physical Engineering Department of Nuclear Reactors Feasibility study on high-conversion Th-U233 fuel cycle for current generation of
More informationFHR Neutronics Benchmarking White Paper. April, UC Berkeley Nuclear Engineering Department
FHR Neutronics Benchmarking White Paper April, 2016 UC Berkeley Nuclear Engineering Department Contents Contents... 1 List of Figures... 3 List of Tables... 4 1 Introduction... 5 1.1 White paper outline...
More informationAdvanced Heavy Water Reactor. Amit Thakur Reactor Physics Design Division Bhabha Atomic Research Centre, INDIA
Advanced Heavy Water Reactor Amit Thakur Reactor Physics Design Division Bhabha Atomic Research Centre, INDIA Design objectives of AHWR The Advanced Heavy Water Reactor (AHWR) is a unique reactor designed
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 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 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 informationSystems Analysis of the Nuclear Fuel Cycle CASMO-4 1. CASMO-4
1. CASMO-4 1.1 INTRODUCTION 1.1.1 General The CASMO-4 code is a multi-group two-dimensional transport code developed by Studsvik, which is entirely written in FORTRAN 77. It is used for burnup calculations
More informationSteady-State and Transient Neutronic and Thermal-hydraulic Analysis of ETDR using the FAST code system
Steady-State and Transient Neutronic and Thermal-hydraulic Analysis of ETDR using the FAST code system Sandro Pelloni, Evaldas Bubelis and Paul Coddington Laboratory for Reactor Physics and Systems Behaviour,
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 information