Presentation for the CIELO Meeting of the NEA 9-11 May 2016 Paris, France. The Chinese work on 56 Fe

Similar documents
ORNL Nuclear Data Evaluation Accomplishments for FY 2013

The updated version of the Chinese Evaluated Nuclear Data Library (CENDL-3.1) and China nuclear data evaluation activities

CIELO Project. M. Herman 1) D. Brown 1), R. Capote 2), G. Nobre 1), A. Trkov 2) for the CIELO Collaboration 3)

The Updated Version of Chinese Evaluated Nuclear Data Library (CENDL-3.1)

Invited. ENDF/B-VII data testing with ICSBEP benchmarks. 1 Introduction. 2 Discussion

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

PROGRESS OF NUCLEAR DATA MEASUREMENT IN CHINA

Resonance Evaluations of 235 U for the CIELO Project

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

Sensitivity Computation with Monte Carlo Methods

NUCLEAR DATA VERIFICATION USING GALILEE AND TRIPOLI-4

Neutron Spectra Measurement and Calculations Using Data Libraries CIELO, JEFF-3.2 and ENDF/B-VII.1 in Spherical Iron Benchmark Assemblies

B. Morillon, L. Leal?, G. Noguere y, P. Romain, H. Duarte. April U, 238 U and 239 Pu JEFF-3.3T1 evaluations

ENDF8. M.B. Chadwick LANL. Mark Chadwick, LANL. WPEC-SG40 (CIELO), Paris, May17, Operated by Los Alamos National Security, LLC for NNSA

Potential Use of beta-eff and other Benchmarks for Adjustment

Review of nuclear data of major actinides and 56 Fe in JENDL-4.0

Benchmark of ENDF/B-VII.1 and JENDL-4.0 on Reflector Effects

VERIFICATION OFENDF/B-VII.0, ENDF/B-VII.1 AND JENDL-4.0 NUCLEAR DATA LIBRARIES FOR CRITICALITY CALCULATIONS USING NEA/NSC BENCHMARKS

Sensitivity and Uncertainty Analysis Methodologies for Fast Reactor Physics and Design at JAEA

Introduction to Nuclear Engineering

EVALUATION OF NEUTRON INDUCED REACTIONS ON 238 U NUCLEUS

A TEMPERATURE DEPENDENT ENDF/B-VI.8 ACE LIBRARY FOR UO2, THO2, ZIRC4, SS AISI-348, H2O, B4C AND AG-IN-CD

Presentation for the 28 th Meeting of the NEA WPEC 9-13 May 2016 Paris, France. GE Zhigang

JRPR. Current Status of ACE Format Libraries for MCNP at Nuclear Data Center of KAERI. Original Research. Introduction

Invited. The latest BROND-3 developments. 1 Introduction. 2 Cross section evaluations for actinides

Analysis of the TRIGA Reactor Benchmarks with TRIPOLI 4.4

TENDL-TMC for dpa and pka

Preliminary Uncertainty Analysis at ANL

1 Introduction. É É â. EPJ Web of Conferences 42, (2013) C Owned by the authors, published by EDP Sciences, 2013

Benchmark Tests of Gamma-Ray Production Data in JENDL-3 for Some Important Nuclides

TENDL-2011 processing and criticality benchmarking

Benchmark Test of JENDL High Energy File with MCNP

Nuclear Data Uncertainty Analysis in Criticality Safety. Oliver Buss, Axel Hoefer, Jens-Christian Neuber AREVA NP GmbH, PEPA-G (Offenbach, Germany)

Iron Neutron Data Benchmarking for 14 MeV Source

R&D in Nuclear Data for Reactor Physics Applications in CNL (CNL = Canadian Nuclear Laboratories) D. Roubtsov

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

International Evaluation

PIA: Progressive Incremental Adjustment

Experimental Activities in China

Needs for Nuclear Reactions on Actinides

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

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

High Energy Neutron Scattering Benchmark of Monte Carlo Computations

Michael Dunn Nuclear Data Group Leader Nuclear Science & Technology Division Medical Physics Working Group Meeting October 26, 2005

INTERCOMPARISON OF CALCULATIONS FOR GODIVA AND JEZEBEL

Testing of Nuclear Data Libraries for Fission Products

Vladimir Sobes 2, Luiz Leal 3, Andrej Trkov 4 and Matt Falk 5

Experimental Activities in China

Considerations for Measurements in Support of Thermal Scattering Data Evaluations. Ayman I. Hawari

Methods and Issues for the Combined Use of Integral Experiments and Covariance Data: Results of a NEA International Collaborative Study

Introduction to Nuclear Data

Benchmark Experiment for Fast Neutron Spectrum Potassium Worth Validation in Space Power Reactor Design

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

Fusion Neutronics, Nuclear Data, Design & Analyses - Overview of Recent FZK Activities -

MCNP6: Fission Multiplicity Model Usage in Criticality Calculations

Fundamentals in Nuclear Physics

Challenges in nuclear data evaluation of actinide nuclei

National Centre of Sciences, Energy and Nuclear Techniques (CNESTEN/CENM), POB 1382, Rabat, Morocco.

Recent Advances in the ENDF/B Library

User Perspective and Analysis of the Continuous- Energy Sensitivity Methods in SCALE 6.2 using TSUNAMI-3D

BENCHMARK CALCULATIONS FOR URANIUM 235

CIELO: Status of Cross Section Progress

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

International Evaluation

New Neutron-Induced Cross-Section Measurements for Weak s-process Studies

Improved nuclear data for material damage applications in LWR spectra

Principles of neutron TOF cross section measurements

Appendix C DESCRIPTION OF EVALUATION METHODS

Nuclear Data Section Department of Nuclear Sciences and Applications

Upcoming features in Serpent photon transport mode

Nuclear data uncertainty propagation using a Total Monte Carlo approach

FRENDY: A new nuclear data processing system being developed at JAEA

Theoretical Analysis of Neutron Double-Differential Cross Section of n + 19 F at 14.2 MeV

CALCULATING UNCERTAINTY ON K-EFFECTIVE WITH MONK10

Neutron-induced reactions on U and Th a new approach via AMS in collaboration with: KIT (Karlsruhe): F. Käppeler, I. Dillmann

Target accuracy of MA nuclear data and progress in validation by post irradiation experiments with the fast reactor JOYO

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

PROGRESS OF NUCLEAR DATA MEASUREMENT IN CHINA

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

Status of evaluated data for neutron induced reactions on 238 U in the resonance region

Integral Benchmark Experiments of the Japanese Evaluated Nuclear Data Library (JENDL)-3.3 for the Fusion Reactor Design

Benchmark Experiments of Accelerator Driven Systems (ADS) in Kyoto University Critical Assembly (KUCA)

Status and future plan of JENDL. Osamu Iwamoto Nuclear Data Center Japan Atomic Energy Agency

IRMM - Institute for Reference Materials and Measurements

Low Energy Neutron Verification in GEANT4: to 4.9.5

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

Nuclear data uncertainty quantification and data assimilation for a lead-cooled fast reactor

Error Estimation for ADS Nuclear Properties by using Nuclear Data Covariances

MUSE-4 BENCHMARK CALCULATIONS USING MCNP-4C AND DIFFERENT NUCLEAR DATA LIBRARIES

arxiv: v2 [hep-ex] 3 Jun 2018

TENDL 2017: better cross sections, better covariances

Present Status of the JENDL Project. Osamu IWAMOTO and Kenji YOKOYAMA Japan Atomic Energy Agency

Statistical Model Calculations for Neutron Radiative Capture Process

Neutron Capture Experiments with DANCE

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

INVESTIGATION OF THE NATURE OF WEAK NEUTRON CAPTURE RESONANCES OF U-238

Measurements of the Neutron Scattering Spectrum from 238 U and Comparison of the Results with a Calculation at the 36.

Present Status and Plans of JENDL FP Data Evaluation Project

Extension of the MCBEND Monte Carlo Code to Perform Adjoint Calculations using Point Energy Data

Coordinated evaluation of 239Pu in the resonance region

Transcription:

The Chinese work on 56 Fe Jing QIAN, Zhigang GE, Tingjin LIU, Hanlin LU, Xichao RUAN Guochang CHEN,Huanyu Huanyu ZHANG,Yangbo NIE China Nuclear Data Center(CNDC) China Institute of Atomic Energy(CIAE) P.O.Box 275-41,Beijing 102413, P.R.China E-Mail: Mail:qianj@ciae.ac.cn

Evaluation of all cross sections in fast range Integral Benchmarking Summary

Evaluation of all cross sections in fast range Experimental data evaluation Theoretical calculation by UNF series code Compilation of evaluated data(checking, CRECTJ6) New Resonance range evaluation (adopted d from CIELO-Iron-beta1) I were merged with the new fast-neutron evaluation ACE file processed by NJOY2012.50 (by M.Herman and Ramon) Comparison with other evaluated data (CIELO-Iron-beta1 B7.1 JEFF3.2 RSF JENDL4.0 CENDL-3.1)

Theoretical calculation result by code UNF UNF series Code APMN06 UNF2008 NDDX,TT DWUCK output A set of suitable OMP were given Calculate the contribution from direct part(dwba) Out put ENDF file DAE and DE were polted Level density parameters and pair correction: Level density parameters OMP V 0 Presentation 54.52155685 for the CIELO Meeting of Uthe 2 NEA -0.00040197 V 1-0.23792717 r r 1.09132423 V 2-0.00590990 r s 1.27747381 V 3-24.87143517 r v 1.49487329 V 4-0.00619497 α r 0.65029401 W 0 10.84276485 α s 0.50812095 W 1-0.25981072 α v 0.44999999 W 2-12.94733429 U SO 6.23620844 U 0-1.39829016 r so 1.08222973 U 1-0.19923522 a so 0.62750721 β param meters 1 so level Energy Jπ β 1 0.4868 2+ 0.260 2 2.0851 4+ 0.44 3 2.6576 2+ 0.06 5 2.9599 2+ 0.08 8 3.3698 2+ 0.12 11 3.4453 3+ 0.03 20 3.8564 3+ 0.03 24 4.1199 3+ 0.05 29 4.3948 3+ 0.01 30 4.4014 2+ 0.05 38 4.6583 3+ 0.12 (n,γ) (n,n ) (n,p) (n,α) (n, 3 He) (n,d) (n,t) (n,2n) (n,n α) (n,np) (n,3n) 7.650 5.152 6.395 5.680 6.076 6.334 6.337 5.261 5.416 7.157 5.577 Pair correction 0.100 0.25-0.75 0.40 0.35 0.15 0.00 0.50 0.40 0.20 0.50

MF MT Data Source 2 151 3 1 ENDF/B-VII.1 4 Cielo Iron-beta1 A resonance evaluation was performed at ORNL (Leal) from 1.0e-5 ev to 2 MeV using transmission data taken at the LINAC (Danon) and inelastic cross section measurements done GEEL (Plompen). Adopt GEEL2014 Exp Data(~6MeV)+Fit Nelson2004Exp Data(6~20MeV) 3,51 Adopt GEEL2014 Exp Data(~4.577MeV) + Theoretical calculated result(unf) (4.577~20MeV) 3,52~88 Theoretical calculated result(unf) 3,91 (3,4)-(3,51:88) 16 Fit Exp Data and ENDF/B-VII.1 103 ENDF/B-VII.1 102, 104,105, Theoretical calculated result(unf) 106,107 3 (3,3)=(3,102)+(3,4)+(3,16)+(3,103)+(3,104)+(3,105)+(3,106)+(3,107) 2 (3,2)=(3,1)-(3,3) 3) 4 2 Theoretical calculated result(unf) 5,6 Theoretical calculated result(unf)

(n,tot) (n,inl) (n,2n) (b) Section (n,p) (n,α) Cross (n,gamma)

The experimental data of r-ray production cross sections were corrected and converted to the inelastic cross section. Presentation for the CIELO Meeting of the NEA

Adopt GEEL experimental data below 6MeV and fit Nelson s data up to 20MeV

Exp.data given by GEEL fromthresholdto45mev to 4.5 Presentation for the CIELO Meeting of the NEA

Adopt GEEL Exp Data(~4.577MeV) + Theoretical calculated result(unf) (4.577~20MeV)

MF MT DataPresentation Source for the CIELO Meeting of the NEA 2 151 3 Cielo Iron-beta1 A resonance evaluation was performed at ORNL (Leal) from 1.0e-5 ev to 2 MeV using transmission data taken at the LINAC (Danon) and inelastic cross section measurements done GEEL (Plompen). 1 ENDF/B-VII.1 4 Adopt GEEL2014 Exp Data(~6MeV)+Fit Nelson2004Exp Data(6~20MeV) 3,51 Fit GEEL2014 Exp Data(~4.577MeV) + Theoretical calculated result(unf) (4.577~20MeV) 3,52~88 Theoretical calculated result(unf) 3,91 16 Fit Exp Data and ENDF/B-VII.1 103 ENDF/B-VII.1 102, 104,105, Theoretical calculated result(unf) 106,107 (3,4)-(3,51:88) (3,5 The continuous inelastic cross section was given by subtracting the all discrete levels from total inelastic cross section 3 (3,3)=(3,102)+(3,4)+(3,16)+(3,103)+(3,104)+(3,105)+(3,106)+(3,107) 2 (3,2)=(3,1)-(3,3) 4 2 Theoretical calculated result(unf) 5,6 Theoretical calculated result(unf)

56 Fe(n,p) The (n, p) cross section is a secondary standard cross section, there are no substantial difference for all evaluated libraries, and adopt evaluated data from ENDF/B-7.1 directly.

(n,non) (n,el) (n,2n) (n,p) (n,α)

MF MT Data Source 2 151 3 Cielo Iron-beta1 A resonance evaluation was performed at ORNL (Leal) from 1.0e-5 ev to 2 MeV using transmission data taken at the LINAC (Danon) and inelastic cross section measurements done GEEL (Plompen). 1 ENDF/B-VII.1 4 Fit GEEL2014 Exp Data(~6MeV)+Nelson2004Exp Data(6~20MeV) 3,51 Fit GEEL2014 Exp Data(~4.577MeV) + Theoretical calculated result(unf) (4.577~20MeV) 3,52~88 Theoretical calculated result(unf) 3,91 (3,4)-(3,51:88) 16 Fit Exp Data and ENDF/B-VII.1 103 ENDF/B-VII.1 102, 104,105, Theoretical calculated result(unf) 106,107 3 (3,3)=(3,102)+(3,4)+(3,16)+(3,103)+(3,104)+(3,105)+(3,106)+(3,107) 2 (3,2)=(3,1)-(3,3) 3) 4 2 Theoretical calculated result(unf) 5,6 Theoretical calculated result(unf)

Elastic scattering angular distribution

The relationship between each reactions Presentation for the CIELO Meeting of the NEA This work is reasonable and (n,tot) in good agreement with evaluated experimental data (n,el) (n,non) (n,inl) (n,2n) (n,p) (n,α)

MF MT Data Source 2 151 3 1 ENDF/B-VII.1 4 Cielo Iron-beta1 A resonance evaluation was performed at ORNL (Leal) from 1.0e-5 ev to 2 MeV using transmission data taken at the LINAC (Danon) and inelastic cross section measurements done GEEL (Plompen). Adopt GEEL2014 Exp Data(~6MeV)+ Fit Nelson2004Exp Data(6~20MeV) 3,51 Fit GEEL2014 Exp Data(~4.577MeV) + Theoretical calculated result(unf) (4.577~20MeV) 3,52~88 Theoretical calculated result(unf) 3,91 (3,4)-(3,51:88) 16 Fit Exp Data and ENDF/B-VII.1 103 ENDF/B-VII.1 102, 104,105, Theoretical calculated result(unf) 106,107 3 (3,3)=(3,102)+(3,4)+(3,16)+(3,103)+(3,104)+(3,105)+(3,106)+(3,107) 2 (3,2)=(3,1)-(3,3) 3) 4 2 Theoretical calculated result(unf) 5,6 Theoretical calculated result(unf)

0 Fe(n,xn) DE 14.1MeV Presentation for the CIELO Meeting of the NEA In good agreement with experimental data, especially in the energy region where the neutrons come from inelastic scattering of discrete levels. The results of (n, n ), (n, p), (n, 2n) and (n, α) cross sections are reasonable.

0 Fe(n,xn) DE 14.1MeV Presentation for the CIELO Meeting of the NEA This work is more close to the experimental p data in the energy region 4-7 MeV than others.

0 Fe(n,xn) DE 14.1MeV Presentation for the CIELO Meeting of the NEA C i t t ith RSF d JEFF b tt th Consistent with RSF and JEFF, better than CENDL-3.1

Evaluation of all cross sections in fast range Integral Benchmarking Summary

Integral Benchmarking Presentation for the CIELO Meeting of the NEA (By Huanyu Zhang, CNDC, China) The new evaluated Fe-56 of this work were tested t with shielding and criticality benchmarks by using MCNP5. Integral experiments Lawrence Livermore National libratory pulse sphere experiment Neutron leakage spectra measurement in CIAE 168 critical benchmarks was selected from ICSBEP (Fe-56 collision weight simulated by MCNP > 0.1). The evaluated Fe-56 CNDC (this work) CIELO(Cielo-beta1) Comparison of benchmarking results Other nuclides except Fe-56 in the benchmarks are from ENDF/B-VII.1.

Comparison of Shielding Benchmarks Result Presentation for the CIELO Meeting of the NEA 14.1MeV (thiswork) The new evaluated inelastic scattering cross sections improve the leakage neutron spectra above 7 MeV.

Comparison of Shielding Benchmarks Result (2) Presentation for the CIELO Meeting of the NEA 14.1MeV (thiswork) But lower than experiments, CIELO and ENDF/B-VII.1 in the range of 9-10MeV

Comparison of Shielding Benchmarks Result (3) 14.1MeV (thiswork)

Neutron leakage spectra measurement induced by 14 MeV neutrons (By Yangbo Nie, CIAE, China) Slab sample 10*10*5 cm D-T neutron source Neutron detector

The experimental data is a preliminary result which is consistent with evaluated data. Our new evaluation will be compared with data later.

Integral Benchmarking Presentation for the CIELO Meeting of the NEA (By Huanyu Zhang, CNDC, China) The new evaluated Fe-56 of this work were tested t with shielding and criticality benchmarks by using MCNP5. Integral experiments Lawrence Livermore National libratory pulse sphere experiment Neutron leakage spectra measurement in CIAE 168 critical benchmarks was selected from ICSBEP (Fe-56 collision weight simulated by MCNP > 0.1). The evaluated Fe-56 CNDC (this work) CIELO(Cielo-beta1) ) Comparison of benchmarking results Other nuclides except Fe-56 in the benchmarks are from ENDF/B-VII.1.

Comparison of Critical Benchmarks Result 2.5 2 1.5 1 ENDF/B-VII.1 JENDL-4.0 CENDL-3.1 CIELO CNDC Delta k 10 3 eff (1 pcm) 0.5 0 0.5 1 1.5 2 HMF HMI HMT HST ICI ICT IMF LCT LST Benchmark MMF PMF PST

Highly Enriched Uranium System HEU-MET-FAST C/E K-eff 1.008 1.006 1.004 1.002 1 0.998 0.996 0.994 0.992 0.99 exp ENDF/B VII.1 JENDL 4.0 40 CENDL 3.1 CIELO CNDC 0.988 1.E+03 1.E+04 1.E+05 1.E+06 EALF (ev) HMF021_01,HMF024_01,HMF034_03,HMF048_16,HMF048_17

C/E K-eff 1.02 1.015 1.01 1.005 1 exp ENDF/B VII.1 JENDL 4.0 CENDL 3.1 CIELO CNDC HEU-MET-THERM HMT013_2 HMT013_1 0.995 HMT015_1 0.06 0.065 0.07 0.075 0.08 0.085 EALF (ev) HMT013_01,HMT013_02,HMT015_01

C/E K-ef ff 1.03 1.025 1.02 1.015 1.01 1.005 1 0.995 0.99 0.985 exp ENDF/B VII.1 JENDL 4.0 CENDL 3.1 CIELO CNDC HEU-SOL-THERM HST016_3 0.98 1.E-02 1.E-01 1.E+00 1.E+01 EALF (ev) HST004_03~06, HST006_09, HST016_01, HST016_02, HST019_01, HST019_02, HST019_03, HST020_05, HST038_15, HST049_01, HST049_09

Intermediate Enriched Uranium System 1.015 IEU-COMP System 1.01 1.005 1 C/E K-eff 0.995 099 0.99 0.985 0.98 0.975 0.97 ICT05_4 exp ENDF/B VII.1 JENDL 4.0 CENDL 3.1 CIELO CNDC 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 EALF (ev) ICF001_01,ICI001_02,ICI001_03,ICT002_01,ICT002_03,ICT002_05

Low Enriched Uranium System 1.02 1.015 1.01 LEU-COMP-THERM exp ENDF/B VII.1 JENDL 4.0 CENDL 3.1 C/E K-eff 1.005 1 0.995 CIELO CNDC 099 0.99 0.985 0 0.2 0.4 0.6 0.8 1 1.2 1.4 EALF (ev) LCT018_01,LCT019_01~03,LCT020_01~07,LCT021_01,LCT022_01~06, LCT023_ 01~06,LCT024_ 01~02,LCT025_ 01~04,LCT026_ 01,LCT026_ 03, LCT026_05,LCT048_01

LEU-SOL-THERM Presentation for the CIELO Meeting of the NEA C/E 1.007 1.006 1.005 1.004 1.003 1002 1.002 1.001 1 exp ENDF/B VII.1 JENDL 4.0 CENDL 3.1 CIELO CNDC K-eff 0.999 0.998 0.03 0.035 0.04 0.045 0.05 0.055 EALF (ev) LST004_01~07,LST016_01~07,LST017_01~06,LST018_01~06, LST019_01~05,LST020_01~0404 For LST system, the k eff bias against EALF is related to Fe-56. The new evaluated RRPs cause the slope of the bias get a litter large.

ZPPR-21 Assembles 1.004 exp MIX-MET-FAST C/E K-eff 1.002 1 0.998 0.996 ENDF/B VII.1 JENDL 4.0 CENDL 3.1 CIELO CNDC 0.994 0.992 1.0E+05 1.5E+05 2.0E+05 2.5E+05 3.0E+05 3.5E+05 EALF (ev) MMF011_01 ~04 It is sensitive to 56 Fe(n,inl) and (n,el) cross sections. The calculated results get worse or better than B71

1.01 exp Plutonium system PU-ME-FAST (Fe reflect) C/E K-eff 1.005 1 0.995 0.99 ENDF/B VII.1 JENDL 4.0 CENDL 3.1 CIELO CNDC 0.985 0.98 6.0E+05 7.0E+05 8.0E+05 9.0E+05 1.0E+06 1.1E+06 1.2E+06 1.3E+06 1.4E+06 EALF (ev) PMF045_07,PMF045_06,PMF045_02,PMF015,PMF045_03,PMF028,PMF026,PMF032,PMF025

PU-SOL-THERM C/E K-eff 1.025 1.02 1.015 1.01 1.005 1 exp ENDF/B VII.1 JENDL 4.0 CENDL 3.1 CIELO CNDC 0995 0.995 0.99 0.0E+00 2.0E-02 4.0E-02 6.0E-02 8.0E-02 1.0E-01 1.2E-01 1.4E-01 1.6E-01 1.8E-01 EALF (ev) PST010_01~14,PST012_06~12,PST023_01~34

Evaluation of all cross sections in fast range Integral Benchmarking Summary

summary Evaluation of all cross sections from CNDC in fast range was finished. The result is in good agreement with evaluated experimental data. Energy spectra and double differential cross sections for neutron emission show the cross sections are reasonable. Integral Benchmarking shows the new evaluated inelastic scattering cross sections improve the leakage neutron spectra above 7 MeV. The results of CNDC show good agreements with those results from CIELO. But no significant ifi improvement was found in criticality benchmarks.