Uncertainty quantification using SCALE 6.2 package and GPT techniques implemented in Serpent 2

Similar documents
SG39 Meeting May 16-17, Update on Continuous Energy Cross Section Adjustment. UC Berkeley / INL collaboration

Implementation of new adjoint-based methods for sensitivity analysis and uncertainty quantication in Serpent

Nuclear data sensitivity and uncertainty assessment of sodium voiding reactivity coefficients of an ASTRID-like Sodium Fast Reactor

Perturbation/sensitivity calculations with Serpent

Convergence Analysis and Criterion for Data Assimilation with Sensitivities from Monte Carlo Neutron Transport Codes

CASMO-5/5M Code and Library Status. J. Rhodes, K. Smith, D. Lee, Z. Xu, & N. Gheorghiu Arizona 2008

Neutronic analysis of SFR lattices: Serpent vs. HELIOS-2

Sensitivity and Uncertainty Analysis of the k eff and b eff for the ICSBEP and IRPhE Benchmarks

ANALYSIS OF THE COOLANT DENSITY REACTIVITY COEFFICIENT IN LFRs AND SFRs VIA MONTE CARLO PERTURBATION/SENSITIVITY

Potential Use of beta-eff and other Benchmarks for Adjustment

COMPARATIVE ANALYSIS OF WWER-440 REACTOR CORE WITH PARCS/HELIOS AND PARCS/SERPENT CODES

On the use of SERPENT code for few-group XS generation for Sodium Fast Reactors

TMS On-the-fly Temperature Treatment in Serpent

Sodium void coefficient map by Serpent

PIA: Progressive Incremental Adjustment

A.BIDAUD, I. KODELI, V.MASTRANGELO, E.SARTORI

A PERTURBATION ANALYSIS SCHEME IN WIMS USING TRANSPORT THEORY FLUX SOLUTIONS

USE OF LATTICE CODE DRAGON IN REACTOR CALUCLATIONS

New Approaches and Applications for Monte Carlo Perturbation Theory.

First ANDES annual meeting

A-priori and a-posteriori covariance data in nuclear cross section adjustments: issues and challenges

SENSITIVITY ANALYSIS OF ALLEGRO MOX CORE. Bratislava, Iľkovičova 3, Bratislava, Slovakia

On-the-fly Doppler Broadening in Serpent

Scope and Objectives. Codes and Relevance. Topics. Which is better? Project Supervisor(s)

IMPACT OF THE FISSION YIELD COVARIANCE DATA IN BURN-UP CALCULATIONS

New methods implemented in TRIPOLI-4. New methods implemented in TRIPOLI-4. J. Eduard Hoogenboom Delft University of Technology

CALCULATION OF TEMPERATURE REACTIVITY COEFFICIENTS IN KRITZ-2 CRITICAL EXPERIMENTS USING WIMS ABSTRACT

Sensitivity Computation with Monte Carlo Methods

Challenges in nuclear data evaluation of actinide nuclei

Homogenization Methods for Full Core Solution of the Pn Transport Equations with 3-D Cross Sections. Andrew Hall October 16, 2015

CONTROL ROD WORTH EVALUATION OF TRIGA MARK II REACTOR

In collaboration with NRG

Nuclear Data Uncertainty Quantification for Applications in Energy, Security, and Isotope Production

Cross Section Generation Strategy for High Conversion Light Water Reactors Bryan Herman and Eugene Shwageraus

WPEC Sub group 34 Coordinated evaluation of 239 Pu in the resonance region

PROPAGATION OF NUCLEAR DATA UNCERTAINTIES IN FUEL CYCLE USING MONTE-CARLO TECHNIQUE

Preliminary Uncertainty Analysis at ANL

Investigating Effects of Sensitivity Uncertainties

Demonstration of Full PWR Core Coupled Monte Carlo Neutron Transport and Thermal-Hydraulic Simulations Using Serpent 2/ SUBCHANFLOW

Excerpt from the Proceedings of the COMSOL Users Conference 2007 Grenoble

Consistent Code-to-Code Comparison of Pin-cell Depletion Benchmark Suite

PIA and REWIND: Two New Methodologies for Cross Section Adjustment. G. Palmiotti and M. Salvatores

Perspective on CIELO. Michael Dunn Nuclear Data & Criticality Safety Group Leader. NEMEA-7 / CIELO Workshop Geel, Belgium November 5-8, 2013

Resonance self-shielding methodology of new neutron transport code STREAM

CALCULATING UNCERTAINTY ON K-EFFECTIVE WITH MONK10

ABSTRACT 1 INTRODUCTION

Elastic scattering. Elastic scattering

Treatment of Implicit Effects with XSUSA.

DOPPLER COEFFICIENT OF REACTIVITY BENCHMARK CALCULATIONS FOR DIFFERENT ENRICHMENTS OF UO 2

The Effect of Burnup on Reactivity for VVER-1000 with MOXGD and UGD Fuel Assemblies Using MCNPX Code

Working Party on Pu-MOX fuel physics and innovative fuel cycles (WPPR)

On the Use of Serpent for SMR Modeling and Cross Section Generation

Assessment of the MCNP-ACAB code system for burnup credit analyses

Safety Analyses for Dynamical Events (SADE) SAFIR2018 Interim Seminar Ville Sahlberg

Investigation of Nuclear Data Accuracy for the Accelerator- Driven System with Minor Actinide Fuel

Monte Carlo Methods for Uncertainly Analysis Using the Bayesian R-Matrix Code SAMMY

Error Estimation for ADS Nuclear Properties by using Nuclear Data Covariances

Research Article Calculations for a BWR Lattice with Adjacent Gadolinium Pins Using the Monte Carlo Cell Code Serpent v.1.1.7

Serpent Monte Carlo Neutron Transport Code

2. The Steady State and the Diffusion Equation

Sensitivity Analysis of Gas-cooled Fast Reactor

Analysis of Neutron Thermal Scattering Data Uncertainties in PWRs

Nuclear Data for Innovative Fast Reactors: Impact of Uncertainties and New Requirements

VERIFICATION OF A REACTOR PHYSICS CALCULATION SCHEME FOR THE CROCUS REACTOR. Paul Scherrer Institut (PSI) CH-5232 Villigen-PSI 2

NEUTRON PHYSICAL ANALYSIS OF SIX ENERGETIC FAST REACTORS

Uncertainty propagation

Chapter 5: Applications Fission simulations

THE NEXT GENERATION WIMS LATTICE CODE : WIMS9

Needs for Nuclear Reactions on Actinides

FHR Neutronics Benchmarking White Paper. April, UC Berkeley Nuclear Engineering Department

Use of Monte Carlo and Deterministic Codes for Calculation of Plutonium Radial Distribution in a Fuel Cell

TENDL 2017: better cross sections, better covariances

Monte Carlo neutron transport and thermal-hydraulic simulations using Serpent 2/SUBCHANFLOW

The Use of Serpent 2 in Support of Modeling of the Transient Test Reactor at Idaho National Laboratory

Study of Predictor-corrector methods. for Monte Carlo Burnup Codes. Dan Kotlyar Dr. Eugene Shwageraus. Supervisor

Malcolm Bean AT THE MAY All Rights Reserved. Signature of Author: Malcolm Bean Department of Nuclear Science and Engineering

Calculation of uncertainties on DD, DT n/γ flux at potential irradiation positions (vertical ports) and KN2 U3 by TMC code (L11) Henrik Sjöstrand

Analysis of the TRIGA Reactor Benchmarks with TRIPOLI 4.4

CROSS SECTION WEIGHTING SPECTRUM FOR FAST REACTOR ANALYSIS

Uncertainties in activity calculations of different nuclides in reactor steels by neutron irradiation

Lecture 20 Reactor Theory-V

QUADRATIC DEPLETION MODEL FOR GADOLINIUM ISOTOPES IN CASMO-5

A Hybrid Deterministic / Stochastic Calculation Model for Transient Analysis

Development of Multigroup Cross Section Generation Code MC 2-3 for Fast Reactor Analysis

BURNUP CALCULATION CAPABILITY IN THE PSG2 / SERPENT MONTE CARLO REACTOR PHYSICS CODE

A Deterministic against Monte-Carlo Depletion Calculation Benchmark for JHR Core Configurations. A. Chambon, P. Vinai, C.

Implementation of the CLUTCH method in the MORET code. Alexis Jinaphanh

Complete activation data libraries for all incident particles, all energies and including covariance data

A PWR ASSEMBLY COMPUTATIONAL SCHEME BASED ON THE DRAGON V4 LATTICE CODE

VALIDATION OF VISWAM SQUARE LATTICE MODULE WITH MOX PIN CELL BENCHMARK

Lectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 1. Title: Neutron Life Cycle

A New MCNPX PTRAC Coincidence Capture File Capability: A Tool for Neutron Detector Design

High precision neutron inelastic cross section measurements

Computational and Experimental Benchmarking for Transient Fuel Testing: Neutronics Tasks

Symmetry in Monte Carlo. Dennis Mennerdahl OECD/NEA/NSC/WPNCS/AMCT EG, Paris, 18 September 2014

TRADE Preliminary comparison of codes for the activation calculation of the water cooling the target

Fuel cycle studies on minor actinide transmutation in Generation IV fast reactors

Present Status of JEFF-3.1 Qualification for LWR. Reactivity and Fuel Inventory Prediction

Covariance Generation using CONRAD and SAMMY Computer Codes

Solving the neutron slowing down equation

Transcription:

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 2 Giorgio BAIOCCO Alessandro PETRUZZI Marco CHERUBINI

Overview Introduction Codes and methods Test cases Results o UAM Exercise I-1 k eff sensitivity and uncertainty 1 group homogenized XS sensitivity and uncertainty Conclusions and future work 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 2/36

Introduction The main objective is to test Monte Carlo code SERPENT-2 in S/U analysis The program is tested on exercises I-1 of UAM benchmark o TMI fuel pin HZP Several cross sections (i.e. elastic & inelastic scattering, fission, capture etc.) of several isotopes (i.e. U, U, 16 O and 1 H) are perturbed 2 different cross section covariance libraries are used o 44 groups from SCALE 6.1 o 56 groups from SCALE 6.2 Uncertainty is assessed over different responses o k eff o 1 group s capt, s fiss, S capt and S fiss Results are compared with SCALE/TSUNAMI-1D 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 3/36

Codes and methods The effect of a perturbation of a parameter x on the response R is: dr/ dx S R, x / R x Uncertainty can be calculated using the sandwitch rule: R S Cov x S Var R, x R, x s R Var R 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 4/36

Codes and methods SERPENT-2 is a 3D continuous-energy Monte Carlo code for reactor physics calculation In this work a modified version of SERPENT-2 (provided by UC Berkeley), with implemented GPT is used to calculate sensitivity coefficients of different response functions (i.e. k eff, homogenized microscopic and macroscopic cross sections) related to the perturbation of input cross sections. Details about this program can be found in M. Aufiero et al. A collision historybased approach to sensitivity/perturbation calculations in the continuous energy Monte Carlo code SERPENT 2 different cross section covariance libraries are used for uncertainty quantification: o 44 groups from SCALE 6.1 o 56 groups from SCALE 6.2 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 5/36

Test cases Considered case o TMI pin at HZP 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 6/36

Test cases Isotopes considered in uncertainty quantification o U- o U- o O-16 o H-1 Perturbed cross sections o (n, elastic) o (n, n ) o (n, fission) o (n, g) o nu-bar o chi-bar 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 7/36

Test cases Uncertainty quantification k eff 1 group s cap U 1 group s fiss U TMI Pin 1 group s cap U 1 group s fiss U 1 group S cap Fuel 1 group S fiss Fuel 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 8/36

Results Uncertainty quantification k eff 1 group s cap U 1 group s fiss U TMI Pin 1 group s cap U 1 group s fiss U 1 group S cap Fuel 1 group S fiss Fuel 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 9/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff 9.0E-02 TMI Pin - k eff - U (n,fission) 8.0E-02 7.0E-02 6.0E-02 5.0E-02 4.0E-02 3.0E-02 2.0E-02 1.0E-02 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 SERPENT 44 groups SCALE 6.1-44 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 10/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff 9.0E-02 TMI Pin - k eff - U (n,fission) 8.0E-02 7.0E-02 6.0E-02 5.0E-02 4.0E-02 Integral Values Serpent = 0.25398 ± 2.36E-04 SCALE 6.1 = 0.25267 SCALE 6.2 = 0.25338 3.0E-02 2.0E-02 1.0E-02 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 SERPENT 44 groups SCALE 6.1-44 groups SCALE 6.2-56 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 11/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff -5.0E-03 TMI Pin - k eff - U (n,g) 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-1.0E-02-1.5E-02-2.0E-02-2.5E-02-3.0E-02-3.5E-02-4.0E-02-4.5E-02 SERPENT 44 groups SCALE 6.1-44 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 12/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff -5.0E-03 TMI Pin - k eff - U (n,g) 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-1.0E-02-1.5E-02-2.0E-02-2.5E-02-3.0E-02-3.5E-02 Integral Values Serpent = -0.15346 ± 8.90E-05 SCALE 6.1 = -0.15435 SCALE 6.2 = -0.15455-4.0E-02-4.5E-02 SERPENT 44 groups SCALE 6.1-44 groups SCALE 6.2-56 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 13/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff 3.5E-01 TMI Pin - k eff - U nubar 3.0E-01 2.5E-01 2.0E-01 1.5E-01 1.0E-01 5.0E-02 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 SERPENT 44 groups SCALE 6.1-44 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 14/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff 3.5E-01 TMI Pin - k eff - U nubar 3.0E-01 2.5E-01 2.0E-01 Integral Values Serpent = 0.93964 ± 5.64E-05 SCALE 6.1 = 0.94043 SCALE 6.2 = 0.94032 1.5E-01 1.0E-01 5.0E-02 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 SERPENT 44 groups SCALE 6.1-44 groups SCALE 6.2-56 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 15/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff 1.0E-03 TMI Pin - k eff - U (n,n') 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-1.0E-03-2.0E-03-3.0E-03-4.0E-03-5.0E-03-6.0E-03-7.0E-03 SERPENT 44 groups SCALE 6.1-44 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 16/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff 1.0E-03 TMI Pin - k eff - U (n,n') 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-1.0E-03-2.0E-03-3.0E-03-4.0E-03-5.0E-03-6.0E-03 Integral Values Serpent = -0.00581 ± 1.74E-04 SCALE 6.1 = -0.00588 SCALE 6.2 = -0.00593-7.0E-03 SERPENT 44 groups SCALE 6.1-44 groups SCALE 6.2-56 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 17/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff 3.0E-02 TMI Pin - k eff - U (n,fission) 2.5E-02 2.0E-02 1.5E-02 1.0E-02 5.0E-03 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-5.0E-03 SERPENT 44 groups SCALE 6.1-44 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 18/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff 3.0E-02 TMI Pin - k eff - U (n,fission) 2.5E-02 2.0E-02 1.5E-02 Integral Values Serpent = 0.02841 ± 6.42E-05 SCALE 6.1 = 0.02832 SCALE 6.2 = 0.02835 1.0E-02 5.0E-03 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-5.0E-03 SERPENT 44 groups SCALE 6.1-44 groups SCALE 6.2-56 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 19/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff TMI Pin - k eff - U (n,g) 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-1.0E-02-2.0E-02-3.0E-02-4.0E-02-5.0E-02-6.0E-02 SERPENT 44 groups SCALE 6.1-44 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 20/36

Sensitivity coefficient per unit lethargy Results TMI Pin - k eff -2.0E-02 TMI Pin - k eff - U (n,g) 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-4.0E-02-6.0E-02-8.0E-02-1.0E-01-1.2E-01 Integral Values Serpent = -0.19942 ± 1.08E-04 SCALE 6.1 = -0.21119 SCALE 6.2 = -0.20945-1.4E-01-1.6E-01-1.8E-01 SERPENT 44 groups SCALE 6.1-44 groups SCALE 6.2-56 groups 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 21/36

Results TMI Pin - k eff Uncertainty quantification TMI Pin - k eff SERPENT 44 groups SCALE 6.1 44 groups SCALE 6.2 56 groups Covariance Matrix Relative standard deviation in k eff (pcm) Relative standard deviation in k eff (pcm) Relative standard deviation in k eff (pcm) U nubar U nubar 264.63 264.78 341.22 U (n,g) U (n,g) 253.32 269.32 275.11 U (n,g) U (n,g) 208.99 210.53 196.52 U (n,n ) U (n,n ) 114.25 113.08 113.60 U c U c 93.73 92.72 152.32 U (n,fiss) U (n,fiss) 76.86 76.26 76.65 U nubar U nubar 70.47 69.66 70.41 1 H (n,el) 1 H (n,el) 23.97 24.90 16.15 SERPENT 44 groups SCALE 6.1 44 groups SCALE 6.2 56 groups k eff % Dk/k k eff % Dk/k k eff % Dk/k 1.43018 ± 4E-05 0.467% 1.42394 0.483% 1.42441 0.542% 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 22/36

Results Uncertainty quantification k eff 1 group s cap U 1 group s fiss U TMI Pin 1 group s cap U 1 group s fiss U 1 group S cap Fuel 1 group S fiss Fuel 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 23/36

Results Direct and indirect effects Sensitivity coefficients based on the Generalized Perturbation Theory expressions consist of 2 different effects: o Direct effects: due to the perturbation of cross sections appearing directly in the response function o Indirect effects: due to the perturbation of neutron flux (caused by the cross sections perturbation) These effects have been calculated separately with Serpent for the evaluation of the cross sections sensitivity coefficient 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 24/36

Sensitivity coefficient per unit lethargy Results TMI Pin s capt, U sensitivity to U (n,g) 6.0E-01 5.0E-01 44 groups 4.0E-01 3.0E-01 2.0E-01 1.0E-01 0.0E+00-1.0E-011.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-2.0E-01-3.0E-01-4.0E-01 Indirect term Direct term Total The indirect term is negative because an increase in the cross-section data enhances the local depression of the neutron flux, causing a decrease of the 1 group s The direct term is positive because an increase in the cross-section data directly increase the 1 group s 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 25/36

Sensitivity coefficient per unit lethargy Sensitivity coefficient per unit lethargy 3.00E-01 2.50E-01 2.00E-01 1.50E-01 1.00E-01 5.00E-02 6.0E-01 5.0E-01 4.0E-01 3.0E-01 2.0E-01 1.0E-01 0.0E+00-1.0E-01 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-2.0E-01-3.0E-01-4.0E-01 6.0E-01 5.0E-01 4.0E-01 3.0E-01 2.0E-01 1.0E-01 Results TMI Pin s capt, U sensitivity to U (n,g) 44 groups Integral Value Serpent = 0.67425 ± 1.62E-04 SCALE 6.1 = 0.69587 56 groups 0.00E+00 1.00E-03 0.0E+00 1.00E-02 1.00E-01 1.00E+00 1.00E+01 1.00E+02 1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07-5.00E-02 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 Integral Value SCALE 6.2 = 0.69558 56 groups Serpent - Indirect term Serpent - Direct term Serpent - Total SCALE 6.1 SCALE 6.2 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 26/36

Sensitivity coefficient per unit lethargy Sensitivity coefficient per unit lethargy Results TMI Pin s capt, U sensitivity to U (n,g) 2.5E-01 2.0E-01 44 groups 1.5E-01 1.0E-01 5.0E-02 Integral Value Serpent = 0.88167 ± 3.31E-04 SCALE 6.1 = 0.88332 3.00E-01 2.50E-01 2.00E-01 1.50E-01 1.00E-01 5.00E-02 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-5.0E-02 2.5E-01 2.0E-01 1.5E-01 1.0E-01 5.0E-02 56 groups 0.00E+00 1.00E-03 0.0E+00 1.00E-02 1.00E-01 1.00E+00 1.00E+01 1.00E+02 1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07-5.00E-02 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 Integral Value SCALE 6.2 = 0.88276 56 groups Serpent - Indirect term Serpent - Direct term Serpent - Total SCALE 6.1 SCALE 6.2 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 27/36

Sensitivity coefficient per unit lethargy Sensitivity coefficient per unit lethargy Results TMI Pin s fiss, U sensitivity to U (n,fiss) 3.5E-01 2.5E-01 1.5E-01 5.0E-02 44 groups Integral Value Serpent = 0.37806 ± 4.31E-04 SCALE 6.1 = 0.37480-5.0E-021.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-1.5E-01 3.00E-01 2.50E-01 2.00E-01-2.5E-01 3.5E-01 2.5E-01 56 groups 56 groups 1.50E-01 1.00E-01 1.5E-01 5.0E-02 Integral Value SCALE 6.2 = 0.37773 5.00E-02-5.0E-021.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 0.00E+00-1.5E-01 1.00E-03 1.00E-02 1.00E-01 1.00E+00 1.00E+01 1.00E+02 1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07-5.00E-02-2.5E-01 Serpent - Indirect term Serpent - Direct term Serpent - Total SCALE 6.1 SCALE 6.2 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 28/36

Sensitivity coefficient per unit lethargy Sensitivity coefficient per unit lethargy Results TMI Pin s fiss, U sensitivity to U (n,fiss) 1.2E+00 1.0E+00 44 groups 8.0E-01 6.0E-01 4.0E-01 2.0E-01 Integral Value Serpent = 0.97716 ± 5.83E-05 SCALE 6.1 = 0.97763 3.00E-01 2.50E-01 2.00E-01 1.50E-01 1.00E-01 5.00E-02 0.0E+00 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07-2.0E-01 1.2E+00 1.0E+00 8.0E-01 6.0E-01 4.0E-01 2.0E-01 56 groups Integral Value SCALE 6.2 = 0.97753 0.00E+00 1.00E-03 0.0E+00 1.00E-02 1.00E-01 1.00E+00 1.00E+01 1.00E+02 1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07-5.00E-02 1.E-03 1.E-02 1.E-01 1.E+00 1.E+01 1.E+02 1.E+03 1.E+04 1.E+05 1.E+06 1.E+07 56 groups Serpent - Indirect term Serpent - Direct term Serpent - Total SCALE 6.1 SCALE 6.2 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 29/36

Results TMI Pin - 1 group s Uncertainty quantification TMI Pin 1 group s s capt, U SERPENT 44 groups SCALE 6.1 44 groups SCALE 6.2 56 groups s fiss, U SERPENT 44 groups SCALE 6.1 44 groups SCALE 6.2 56 groups Covariance Matrix RSD (%) RSD (%) RSD (%) U (n,g) U (n,g) 1.284% 1.294% 1.171% Covariance Matrix RSD (%) RSD (%) RSD (%) U (n,n ) U (n,n ) 0.930% 0.918% 0.931% U (n,n ) U (n,n ) 0.962% 0.946% 0.960% U c U c 0.433% 0.430% 0.725% U c U c 0.447% 0.444% 0.750% 1 H (n,el) 1 H (n,el) 0.409% 0.407% 0.241% 1 H (n,el) 1 H (n,el) 0.413% 0.411% 0.215% U (n,g) U (n,g) 0.200% 0.211% 0.216% U (n,fiss) U (n,g) 0.340% 0.340% 0.340% U (n,g) U (n,g) 0.172% 0.172% 0.166% s capt, U SERPENT 44 groups SCALE 6.1 44 groups SCALE 6.2 56 groups s fiss, U SERPENT 44 groups SCALE 6.1 44 groups SCALE 6.2 56 groups Covariance Matrix RSD (%) RSD (%) RSD (%) U (n,n ) U (n,n ) 0.919% 0.902% 0.914% Covariance Matrix RSD (%) RSD (%) RSD (%) U (n,n ) U (n,n ) 3.273 % 3.222% 3.231% U (n,g) U (n,g) 0.844% 0.877% 0.901% U c U c 2.200 % 2.198% 3.623% U c U c 0.444% 0.437% 0.735% 1 H (n,el) 1 H (n,el) 0.362% 0.361% 0.149% U (n,g) U (n,g) 0.064% 0.065% 0.052% 1 H (n,el) 1 H (n,el) 0.518 % 0.513% 0.202% U (n,fiss) U (n,fiss) 0.509 % 0.510% 0.510% U c U c 0.166 % 0.166% 0.569% 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 30/36

Results TMI Pin - 1 group s Uncertainty quantification TMI Pin 1 group s SERPENT 44 groups s capt, U 8.263 ± 4E-05 Total RSD (%) SCALE 6.1 44 groups s capt, U Total RSD (%) SCALE 6.2 56 groups s capt, U Total RSD (%) 1.753 % 8.297E+00 1.765% 8.286E+00 1.758% SERPENT 44 groups s fiss, U 3.517E+01 ± 5E-05 Total RSD (%) SCALE 6.1 44 groups s fiss, U Total RSD (%) SCALE 6.2 56 groups s fiss, U Total RSD (%) 1.145 % 3.532E+01 1.140% 3.517E+01 1.252% Uncertainty in s fiss, U is big mainly due to the contribution of U (n,n ) and U c SERPENT 44 groups s capt, U 8.136E-01 ± 6E-05 Total RSD (%) SCALE 6.1 44 groups s capt, U Total RSD (%) SCALE 6.2 56 groups s capt, U Total RSD (%) 1.378 % 8.350E-01 1.385% 8.266E-01 1.500% SERPENT 44 groups s fiss, U 1.019E-01 ± 7E-05 Total RSD (%) SCALE 6.1 44 groups s fiss, U Total RSD (%) SCALE 6.2 56 groups s fiss, U Total RSD (%) 4.013 % 1.010E-01 3.963% 1.008E-01 4.913% 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 31/36

Results Uncertainty quantification k eff 1 group s cap U 1 group s fiss U TMI Pin 1 group s cap U 1 group s fiss U 1 group S cap Fuel 1 group S fiss Fuel 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 32/36

Results TMI Pin - 1 group S Uncertainty quantification TMI Pin 1 group S S capt, Fuel SERPENT 44 groups SCALE 6.1 44 groups SCALE 6.2 56 groups S fiss, Fuel SERPENT 44 groups SCALE 6.1 44 groups SCALE 6.2 56 groups Covariance Matrix RSD (%) RSD (%) RSD (%) U (n,n ) U (n,n ) 0.919% 0.917% 0.929% Covariance Matrix RSD (%) RSD (%) RSD (%) U (n,n ) U (n,n ) 0.707% 0.702% 0.713% U (n,g) U (n,g) 0.492% 0.516% 0.529% 1 H (n,el) 1 H (n,el) 0.360% 0.359% 0.218% U c U c 0.418% 0.440% 0.740% U c U c 0.293% 0.292% 0.502% U (n,g) U (n,g) 0.403% 0.401% 0.370% U (n,g) U (n,g) 0.185% 0.196% 0.200% 1 H (n,el) 1 H (n,el) 0.375% 0.377% 0.171% U (n,g) U (n,g) 0.160% 0.161% 0.155% SERPENT 44 groups S capt, Fuel 2.720E-02 ± 4E-05 Total RSD (%) SCALE 6.1 44 groups S capt, Fuel Total RSD (%) SCALE 6.2 56 groups S capt, Fuel Total RSD (%) 1.260 % 2.759E-02 1.278% 2.740E-02 1.380% SERPENT 44 groups S fiss, Fuel 4.184E-02 ± 4E-05 Total RSD (%) SCALE 6.1 44 groups S fiss, Fuel Total RSD (%) SCALE 6.2 56 groups S fiss, Fuel Total RSD (%) 0.893 % 4.199E-02 0.896% 4.182E-02 0.954% 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 33/36

Conclusions and future work SERPENT-2 GPT and SCALE 6.1 & 6.2 /TSUNAMI1D sensitivity calculations are in good agreement, regarding the TMI HZP pin analysis. In general uncertainty calculations with 44 & 56 groups covariance matrices are in good agreement (single contributions). Some discrepancies can be found in the contributions of U nubar, U & U c and 1 H (n,el), between 44 and 56 groups calculation. This is probably due to some update in the new 56 groups libraries of SCALE 6.2 (check Ds/s). U nubar and U (n,g) are the major contributors in k eff uncertainty. With the used options in SCALE, U (n,g) uncertainty may be slightly overestimated. U (n,n ) has a strong influence on the uncertainty of homogenized cross sections (both in pin and assembly). Also 1 H (n,el) plays an important role. 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 34/36

Include other sets o HFP o BWR and VVER cases Future work Conclusions and future work Extend the calculation to the fuel assembly and full core (exercises I-2 & I-3) Include other isotopes o Zr o Gd Use the new 56 and 252 groups covariance libraries of SCALE 6.2 with Serpent 2 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 35/36

Conclusions and future work THANK YOU FOR YOUR ATTENTION! 6th International Serpent User Group Meeting Politecnico di Milano, Milan, Italy - September 26 th - 30 th, 2016 36/36

2024 © sciencedocbox.com Privacy Policy | Terms of Service | Feedback