scholarly journals Burnup calculation of the OECD VVER-1000 LEU benchmark assembly using MCNP6 and SRAC2006

2021 ◽  
Vol 8 (4) ◽  
pp. 10-19
Author(s):  
Tiep Nguyen Huu ◽  
Dung Nguyen Thi ◽  
Phu Tran Viet ◽  
Thanh Tran Vinh ◽  
Ha Pham Nhu Viet

The present work aims to perform burnup calculation of the OECD VVER-1000 LEU (lowenriched uranium) computational benchmark assembly using the Monte Carlo code MCNP6 and the deterministic code SRAC2006. The new depletion capability of MCNP6 was applied in the burnup calculation of the VVER-1000 LEU benchmark assembly. The OTF (on-the-fly) methodology of MCNP6, which involves high precision fitting of Doppler broadened cross sections over a wide temperature range, was utilized to handle temperature variation for heavy isotopes. The collision probability method based PIJ module of SRAC2006 was also used in this burnup calculation. The reactivity of the fuel assembly, the isotopic concentrations and the shielding effect due to the presenceof the gadolinium isotopes were determined with burnup using MCNP6 and SRAC2006 incomparison with the available published benchmark data. This study is therefore expected to reveal the capabilities of MCNP6 and SRAC2006 in burnup calculation of VVER-1000 fuel assemblies.

2016 ◽  
Vol 6 (3) ◽  
pp. 16-30
Author(s):  
Huy Hiep Nguyen ◽  
Huu Tiep Nguyen ◽  
Viet Phu Tran ◽  
Tuan Khai Nguyen

The paper aims to develop an MCNP5-ORIGEN2 coupling scheme for burnup calculation. Specifically, the Monte Carlo neutron transport code (MCNP5) and the nuclides depletion and decay calculation code (ORIGEN2) are combined by data processing and linking files written in the PERL programming language. The validity and applicability of the developed coupling scheme are tested through predicting the neutronic and isotopic behavior of the “VVER-1000 LEU Assembly Computational Benchmark”. The MCNP5-ORIGEN2 coupling results showed a good agreement with the k-inf benchmark values within 600 pcm during the entire burnup history. In addition, the differences of isotopes concentration at the end of the burnup (40 MWd/kgHM) when compared with benchmark values were reasonable and generally within 6.5%. The developed coupling scheme also considered the shielding effect due to gadolinium isotopes and simulated well the depletion of isotopes as a function of the radial position in gadolinium bearing fuel rods.


2021 ◽  
Vol 247 ◽  
pp. 06011
Author(s):  
A. Bernal ◽  
M. Pecchia ◽  
D. Rochman ◽  
A. Vasiliev ◽  
H. Ferroukhi

The main goal of this work is to perform pin-by-pin calculations of Swiss LWR fuel assemblies with neutron transport deterministic methods. At Paul Scherrer Institut (PSI), LWR calculations are performed with the core management system CMSYS, which is based on the Studsvik suite of codes. CMSYS includes models for all the Swiss reactors validated against a database of experimental information. Moreover, PSI has improved the pin power calculations by developing models of Swiss fuel assemblies for the Monte Carlo code MCNP, with the isotopic compositions obtained from the In-Core Fuel Management data of the Studsvik suite of codes, by using the SNF code. A step forward is to use a neutron code based on fast deterministic neutron transport methods. The method used in this work is based on a planar Method of Characteristics in which the axial coupling is solved by 1D SP3 method. The neutron code used is nTRACER. Thus, the methodology of this work develops nTRACER models of Swiss PWR fuel assemblies, in which the fuel of each pin and axial level is modelled with the isotopic composition obtained from SNF. This methodology was applied to 2D and 3D calculations of a Swiss PWR fuel assembly. However, this method has two main limitations. First, the cross sections libraries of nTRACER lack some of the isotopes obtained by SNF. Fortunately, this work proves that the missing isotopes do not have a strong effect on keff and the power distribution. Second, the 3D models require high computational memory resources, that is, more than 260 Gb. Thus, the nTRACER code was modified, so now it uses only 8 Gb, without any loss of accuracy. Finally, the keff and power results are compared with Monte Carlo calculations obtained by Serpent.


2020 ◽  
Vol 225 ◽  
pp. 03007
Author(s):  
Tanja Goričanec ◽  
Domen Kotnik ◽  
Žiga Štancar ◽  
Luka Snoj ◽  
Marjan Kromar

An approach for calculating ex-core detector response using Monte Carlo code MCNP was developed. As a first step towards ex-core detector response prediction a detailed MCNP model of the reactor core was made. A script called McCord was developed as a link between deterministic program package CORD-2 and Monte Carlo code MCNP. It automatically generates an MCNP input from the CORD-2 data. A detailed MCNP core model was used to calculate 3D power distributions inside the core. Calculated power distributions were verified by comparison to the CORD-2 calculations, which is currently used for core design calculation verification of the Krško nuclea power plant. For the hot zero power configuration, the deviations are within 3 % for majority of fuel assemblies and slightly higher for fuel assemblies located at the core periphery. The computational model was further verified by comparing the calculated control rod worth to the CORD-2 results. The deviations were within 50 pcm and considered acceptable. The research will in future be supplemented with the in-core and ex-core detector signal calculations and neutron transport outside the reactor core.


2013 ◽  
Vol 2013 ◽  
pp. 1-9 ◽  
Author(s):  
Mario Matijević ◽  
Dubravko Pevec ◽  
Krešimir Trontl

Revised guidelines with the support of computational benchmarks are needed for the regulation of the allowed neutron irradiation to reactor structures during power plant lifetime. Currently, US NRC Regulatory Guide 1.190 is the effective guideline for reactor dosimetry calculations. A well known international shielding database SINBAD contains large selection of models for benchmarking neutron transport methods. In this paper a PCA benchmark has been chosen from SINBAD for qualification of our methodology for pressure vessel neutron fluence calculations, as required by the Regulatory Guide 1.190. The SCALE6.0 code package, developed at Oak Ridge National Laboratory, was used for modeling of the PCA benchmark. The CSAS6 criticality sequence of the SCALE6.0 code package, which includes KENO-VI Monte Carlo code, as well as MAVRIC/Monaco hybrid shielding sequence, was utilized for calculation of equivalent fission fluxes. The shielding analysis was performed using multigroup shielding library v7_200n47g derived from general purpose ENDF/B-VII.0 library. As a source of response functions for reaction rate calculations with MAVRIC we used international reactor dosimetry libraries (IRDF-2002 and IRDF-90.v2) and appropriate cross-sections from transport library v7_200n47g. The comparison of calculational results and benchmark data showed a good agreement of the calculated and measured equivalent fission fluxes.


2021 ◽  
Vol 247 ◽  
pp. 06017
Author(s):  
Cheng Zhang ◽  
Liangzhi Cao ◽  
Yunzhao Li ◽  
Guowei Hua

In this paper, the modeling and simulation of the PWRs loaded with hexagonal fuel assemblies has been implemented with the NECP-Bamboo code. NECP-Bamboo, consisting of a 2D lattice code named Bamboo-Lattice and a 3D steady-state core code named Bamboo-Core, was primitively designed for the PWRs loaded with the rectangular fuel assemblies. As the capability extension for PWRs with hexagonal fuel assemblies, four aspects of improvement have been implemented in NECP-Bamboo. Firstly, the Constructive Solid Geometry (CSG) has been implemented in Bamboo-Lattice for the lattice modeling. Secondly, the explicit modeling of the reflector assembly has been applied to provide more reliable few-group constants, compared with the conventional 1D model for the reflector assembly. Thirdly, the assembly-homogenization capability has been extended to the hexagonal assembly. Fourthly, the diffusion solver in Bamboo-Core based on the Variational Nodal Method (VNM) has been extended to handle hexagonal geometry. With application of the capability-extended NECP-Bamboo, the modeling and simulations for the VVER-1000 benchmark loaded with MOX fuel has been implemented. It can be observed that the numerical results provided by NECP-Bamboo can agree well with corresponding results by the Monte-Carlo code.


2021 ◽  
Vol 9 ◽  
Author(s):  
Francesc Salvat ◽  
José Manuel Quesada

After a summary description of the theory of elastic collisions of nucleons with atoms, we present the calculation of a generic database of differential and integrated cross sections for the simulation of multiple elastic collisions of protons and neutrons with kinetic energies larger than 100 keV. The relativistic plane-wave Born approximation, with binding and Coulomb-deflection corrections, has been used to calculate a database of proton-impact ionization of K-shell and L-, M-, and N-subshells of neutral atoms These databases cover the whole energy range of interest for all the elements in the periodic system, from hydrogen to einsteinium (Z = 1–99); they are provided as part of the penh distribution package. The Monte Carlo code system penh for the simulation of coupled electron-photon-proton transport is extended to account for the effect of the transport of neutrons (released in proton-induced nuclear reactions) in calculations of dose distributions from proton beams. A simplified description of neutron transport, in which neutron-induced nuclear reactions are described as a fractionally absorbing process, is shown to give simulated depth-dose distributions in good agreement with those generated by the Geant4 code. The proton-impact ionization database, combined with the description of atomic relaxation data and electron transport in penelope, allows the simulation of proton-induced x-ray emission spectra from targets with complex geometries.


2020 ◽  
Vol 239 ◽  
pp. 14006
Author(s):  
Tim Ware ◽  
David Hanlon ◽  
Tara Hanlon ◽  
Richard Hiles ◽  
Malcolm Lingard ◽  
...  

Until recently, criticality safety assessment codes had a minimum temperature at which calculations can be performed. Where criticality assessment has been required for lower temperatures, indirect methods, including reasoned argument or extrapolation, have been required to assess reactivity changes associated with these temperatures. The ANSWERS Software Service MONK® version 10B Monte Carlo criticality code, is capable of performing criticality calculations at any temperature, within the temperature limits of the underlying nuclear data in the BINGO continuous energy library. The temperature range of the nuclear data has been extended below the traditional lower limit of 293.6 K to 193 K in a prototype BINGO library, primarily based on JEFF-3.1.2 data. The temperature range of the thermal bound scattering data of the key moderator materials was extended by reprocessing the NJOY LEAPR inputs used to produce bound data for JEFF-3.1.2 and ENDF/B-VIII.0. To give confidence in the low temperature nuclear data, a series of MONK and MCBEND calculations have been performed and results compared against external data sources. MCBEND is a Monte Carlo code for shielding and dosimetry and shares commonalities to its sister code MONK including the BINGO nuclear data library. Good agreement has been achieved between calculated and experimental cross sections for ice, k-effective results for low temperature criticality benchmarks and calculated and experimentally determined eigenvalues for thermal neutron diffusion in ice. To quantify the differences between ice and water bound scattering data a number of MONK criticality calculations were performed for nuclear fuel transport flask configurations. The results obtained demonstrate good agreement with extrapolation methods. There is a discernible difference in the use of ice and water data.


2021 ◽  
Vol 40 (3) ◽  
pp. 62-67
Author(s):  
H.A. Tanash ◽  
D.A.Solovyev undefined ◽  
V.G.Zimin undefined ◽  
A.A.Semenov undefined ◽  
N.V. Schukin ◽  
...  

2005 ◽  
Vol 48 (spe2) ◽  
pp. 191-199 ◽  
Author(s):  
Christophe Champion

When living cells are irradiated by charged particles, a wide variety of interactions occurs that leads to a deep modification of the biological material. To understand the fine structure of the microscopic distribution of the energy deposits, Monte Carlo event-by-event simulations are particularly suitable. However, the development of these track structure codes needs accurate interaction cross sections for all the electronic processes: ionization, excitation, Positronium formation (for incident positrons) and even elastic scattering. Under these conditions, we have recently developed a Monte Carlo code for electrons and positrons in water, this latter being commonly used to simulate the biological medium. All the processes are studied in detail via theoretical differential and total cross sections calculated by using partial wave methods. Comparisons with existing theoretical and experimental data show very good agreements. Moreover, this kind of detailed description allows one access to a useful microdosimetry, which can be coupled to a geometrical modelling of the target organ and then provide a detailed dose calculation at the nanometric scale.


Sign in / Sign up

Export Citation Format

Share Document