scholarly journals METHODOLOGY FOR HIGH-FIDELITY DETERMINISTIC MODELLING OF SWISS LWR FUEL ASSEMBLIES

2021 ◽  
Vol 247 ◽  
pp. 06011
Author(s):  
A. Bernal ◽  
M. Pecchia ◽  
D. Rochman ◽  
A. Vasiliev ◽  
H. Ferroukhi
Keyword(s):  
Monte Carlo ◽  
Cross Sections ◽  
Power Distribution ◽  
Neutron Transport ◽  
3D Models ◽  
Experimental Information ◽  
Monte Carlo Code ◽  
Deterministic Modelling ◽  
Fuel Assemblies ◽  
Paul Scherrer

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.

scholarly journals Predicting Ex-core Detector Response in a PWR with Monte Carlo Neutron Transport Methods

2020 ◽  
Vol 225 ◽  
pp. 03007
Author(s):  
Tanja Goričanec ◽  
Domen Kotnik ◽  
Žiga Štancar ◽  
Luka Snoj ◽  
Marjan Kromar
Keyword(s):  
Monte Carlo ◽  
Neutron Transport ◽  
Reactor Core ◽  
Response Prediction ◽  
Detector Response ◽  
Design Calculation ◽  
Monte Carlo Code ◽  
Control Rod ◽  
The Core ◽  
Fuel Assemblies

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.

scholarly journals Collisions of Nucleons with Atoms: Calculated Cross Sections and Monte Carlo Simulation

2021 ◽  
Vol 9 ◽  
Author(s):  
Francesc Salvat ◽  
José Manuel Quesada
Keyword(s):  
Monte Carlo ◽  
Cross Sections ◽  
Nuclear Reactions ◽  
Impact Ionization ◽  
Neutron Transport ◽  
Emission Spectra ◽  
Monte Carlo Code ◽  
Dose Distributions ◽  
Elastic Collisions ◽  
Proton Impact

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.

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
Keyword(s):  
Monte Carlo ◽  
Cross Sections ◽  
Wide Temperature Range ◽  
Collision Probability ◽  
Shielding Effect ◽  
Probability Method ◽  
Monte Carlo Code ◽  
Benchmark Data ◽  
Fuel Assemblies ◽  
Burnup Calculation

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.

scholarly journals Processing of the multigroup cross-sections for MCNP calculations

2019 ◽  
Vol 9 (2) ◽  
pp. 17-24
Author(s):  
Jakub Lüley ◽  
Branislav Vrban ◽  
Štefan Čerba ◽  
Filip Osuský ◽  
Vladimír Nečas
Keyword(s):  
Monte Carlo ◽  
Cross Section ◽  
Cross Sections ◽  
Power Distribution ◽  
Neutron Transport ◽  
Resonance Energy ◽  
Radiation Shielding ◽  
Processing Scheme ◽  
Section Data ◽  
Model Complex

Stochastic Monte Carlo (MC) neutron transport codes are widely used in various reactorphysics applications, traditionally related to criticality safety analyses, radiation shielding and validation of deterministic transport codes. The main advantage of Monte Carlo codes lies in their ability to model complex and detail geometries without the need of simplifications. Currently, one of the most accurate and developed stochastic MC code for particle transport simulation is MCNP. To achieve the best real world approximations, continuous-energy (CE) cross-section (XS) libraries are often used. These CE libraries consider the rapid changes of XS in the resonance energy range; however, computing-intensive simulations must be performed to utilize this feature. To broaden ourcomputation abilities for industrial application and partially to allow the comparison withdeterministic codes, the CE cross section library of the MCNP code is replaced by the multigroup (MG) cross-section data. This paper is devoted to the cross-section processing scheme involving modified versions of TRANSX and CRSRD codes. Following this approach, the same data may be used in deterministic and stochastic codes. Moreover, using formerly developed and upgraded crosssection processing scheme, new MG libraries may be tailored to the user specific applications. For demonstration of the proposed cross-section processing scheme, the VVER-440 benchmark devoted to fuel assembly and pip-by-pin power distribution was selected. The obtained results are compared with continues energy MCNP calculation and multigroup KENO-VI calculation.

A Monte Carlo study on burnup treatment in sodium-cooled reactor with Th fuel

Kerntechnik ◽  
2021 ◽  
Vol 86 (4) ◽  
pp. 302-311
Author(s):  
M. E. Korkmaz ◽  
N. K. Arslan
Keyword(s):  
Monte Carlo ◽  
Cross Sections ◽  
Power Distribution ◽  
Reactor Core ◽  
Monte Carlo Study ◽  
Minor Actinides ◽  
Monte Carlo Code ◽  
Effective Multiplication Factor ◽  
Conceptual Core ◽  
Effective Multiplication

Abstract Sodium Cooled Reactors is one of the Generation-IV plants selected to manage the long-lived minor actinides and to transmute the long-life radioactive elements. This study presents the comparison between two-designed SFR cores with 600 and 800 MWth total heating power. We have analyzed a conceptual core design and nuclear characteristic of SFR. Monte Carlo depletion calculations have been performed to investigate essential characteristics of the SFR core. The core calculations were performed by using the Serpent Monte Carlo code for determining the burnup behavior of the SFR, the power distribution and the effective multiplication factor. The neutronic and burn-up calculations were done by means of Serpent-2 Code with the ENDF-7 cross-sections library. Sodium Cooled Fast Reactor core was taken as the reference core for Th-232 burnup calculations. The results showed that SFR is an important option to deplete the minor actinides as well as for transmutation from Th-232 to U-233.

scholarly journals KYLIN-V2.0 CODE CALCULATION ABILITY VERIFICATION BASED ON VERA BENCHMAR

2021 ◽  
Vol 247 ◽  
pp. 10020
Author(s):  
Dongyong Wang ◽  
Yingrui Yu ◽  
Xingjie Peng ◽  
Chenlin Wang ◽  
Kun Liu ◽  
...  
Keyword(s):  
Fuel Assembly ◽  
Power Distribution ◽  
Nuclear Power ◽  
Neutron Transport ◽  
Monte Carlo Code ◽  
Neutron Transport Equation ◽  
Control Rod ◽  
Fuel Assemblies ◽  
Fuel Pin ◽  
Calculation Ability

Virtual Environmental for Reactor Analysis (VERA) benchmark was released by the Consortium for Advanced Simulation of Light water reactors (CASL) project in 2012. VERA benchmark includes more than ten problems at different levels, from 2D fuel pin case to 2D fuel assembly case to 3D core refuelling case, in addition, reference results and experimental measured data of some problems were provided by CASL. Fuel assemblies in VERA benchmark are various, including control rod assemblies, Pyrex assembly, IFBA assembly, WABA assembly and gadolinium poison assembly, and so on. In this paper, various fuel assembly models in the VERA benchmark have been built by using KYIIN-V2.0 code to verify its calculation ability from 2D fuel pin case to 2D fuel assembly case to 2D 3x3 fuel assembly case, and making a comparative analysis on the reference results in VERA benchmark, as well as the calculation results of the Monte Carlo code RMC. KYLIN-V2.0 is an advanced neutron transport lattice code developed by Nuclear Power Institute of China (NPIC). The subgroup resonance calculation method is used in KYIIN-V2.0 to obtain effective resonance selfshielding cross section, method of modular characteristics (MOC) is adopted to solve the neutron transport equation, and CRAM method and PPC method is adopted to solve the depletion equation. The numerical results show that KYLIN-V2.0 code has the reliable capability of direct heterogeneous calculation of 2D fuel assembly, and the effective multiplication factor, assembly power distribution, rod power distribution and control rod reactivity worths of various fuel assemblies that are calculated by KYLIN-V2.0 are in better agreement with the reference.

3D Calculations of PWR Vessels Neutron Fluence With EFLUVE 3D Code

2013 ◽  
Author(s):  
Cécile-Aline Gosmain ◽  
Sylvain Rollet ◽  
Damien Schmitt
Keyword(s):  
Monte Carlo ◽  
Fast Neutron ◽  
Pressure Vessel ◽  
Reaction Rate ◽  
Neutron Fluence ◽  
Neutron Transport ◽  
Reaction Rates ◽  
Surveillance Program ◽  
Monte Carlo Code ◽  
Neutron Sources

In the framework of surveillance program dosimetry, the main parameter in the determination of the fracture toughness and the integrity of the reactor pressure vessel (RPV) is the fast neutron fluence on pressure vessel. Its calculated value is extrapolated using neutron transport codes from measured reaction rate value on dosimeters located on the core barrel. EDF R&D has developed a new 3D tool called EFLUVE3D based on the adjoint flux theory. This tool is able to reproduce on a given configuration the neutron flux, fast neutron fluence and reaction rate or dpa results of an exact Monte Carlo calculation with nearly the same accuracy. These EFLUVE3D calculations does the Source*Importance product which allows the calculation of the flux, the neutronic fluence (flux over 1MeV integrated on time) received at any point of the interface between the skin and the pressure vessel but also at the capsules of the pressurized water reactor vessels surveillance program and the dpa and reaction rates at different axial positions and different azimuthal positions of the vessel as well as at the surveillance capsules. Moreover, these calculations can be carried out monthly for each of the 58 reactors of the French current fleet in challenging time (less than 10mn for the total fluence and reaction rates calculations considering 14 different neutron sources of a classical power plant unit compared to more than 2 days for a classic Monte Carlo flux calculation at a given neutron source). The code needs as input: - for each reaction rate, the geometric importance matrix produced for a 3D pin by pin mesh on the basis of Green’s functions calculated by the Monte Carlo code TRIPOLI; - the neutron sources calculated on assemblies data (enrichment, position, fission fraction as a function of evolution), pin by pin power and irradiation. These last terms are based on local in-core activities measurements extrapolated to the whole core by use of the EDF core calculation scheme and a pin by pin power reconstruction methodology. This paper presents the fundamental principles of the code and its validation comparing its results to the direct Monte Carlo TRIPOLI results. Theses comparisons show a discrepancy of less than 0,5% between the two codes equivalent to the order of magnitude of the stochastic convergence of Monte Carlo results.

Validating the Serpent-Ants Calculation Chain Using BEAVRS Fresh Core HZP Data

2021 ◽  
Author(s):  
Ville Valtavirta ◽  
Antti Rintala ◽  
Unna Lauranto
Keyword(s):  
Monte Carlo ◽  
Best Practices ◽  
Power Distribution ◽  
Measured Data ◽  
Monte Carlo Code ◽  
Generation Process ◽  
Isothermal Temperature ◽  
Control Rod ◽  
Temperature Coefficients ◽  
Group Constant

Abstract The Serpent Monte Carlo code and the Serpent-Ants two step calculation chain are used to model the hot zero power physics tests described in the BEAVRS benchmark. The predicted critical boron concentrations, control rod group worths and isothermal temperature coefficients are compared between Serpent and Serpent-Ants as well as against the experimental measurements. Furthermore, radial power distributions in the unrodded and rodded core configurations are compared between Serpent and Serpent-Ants. In addition to providing results using a best practices calculation chain, the effects of several simplifications or omissions in the group constant generation process on the results are estimated. Both the direct and two-step neutronics solutions provide results close to the measured values. Comparison between the measured data and the direct Serpent Monte Carlo solution yields RMS differences of 12.1 mg/kg, 25.1 × 10-5 and 0.67 × 10-5 K-1 for boron, control rod worths and temperature coefficients respectively. The two-step Serpent-Ants solution reaches a similar level of accuracy with RMS differences of 17.4 mg/kg, 23.6 × 10-5 and 0.29 × 10-5 K-1. The match in the radial power distribution between Serpent and Serpent-Ants was very good with the RMS and maximum for pin power errors being 1.31 % and 4.99 % respectively in the unrodded core and 1.67 %(RMS) and 8.39 % (MAX) in the rodded core.

scholarly journals Comparative Neutronics Analysis of DIMPLE S06 Criticality Benchmark with Contemporary Reactor Core Analysis Computer Code Systems

2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
Wonkyeong Kim ◽  
Jinsu Park ◽  
Tomasz Kozlowski ◽  
Hyun Chul Lee ◽  
Deokjung Lee
Keyword(s):  
Comparative Analysis ◽  
Cross Sections ◽  
Power Distribution ◽  
Reactor Core ◽  
Computer Code ◽  
Anisotropic Scattering ◽  
Monte Carlo Code ◽  
Core Analysis ◽  
Modeling Approach ◽  
Significant Difference

A high-leakage core has been known to be a challenging problem not only for a two-step homogenization approach but also for a direct heterogeneous approach. In this paper the DIMPLE S06 core, which is a small high-leakage core, has been analyzed by a direct heterogeneous modeling approach and by a two-step homogenization modeling approach, using contemporary code systems developed for reactor core analysis. The focus of this work is a comprehensive comparative analysis of the conventional approaches and codes with a small core design, DIMPLE S06 critical experiment. The calculation procedure for the two approaches is explicitly presented in this paper. Comprehensive comparative analysis is performed by neutronics parameters: multiplication factor and assembly power distribution. Comparison of two-group homogenized cross sections from each lattice physics codes shows that the generated transport cross section has significant difference according to the transport approximation to treat anisotropic scattering effect. The necessity of the ADF to correct the discontinuity at the assembly interfaces is clearly presented by the flux distributions and the result of two-step approach. Finally, the two approaches show consistent results for all codes, while the comparison with the reference generated by MCNP shows significant error except for another Monte Carlo code, SERPENT2.

scholarly journals Modeling of the ORNL PCA Benchmark Using SCALE6.0 Hybrid Deterministic-Stochastic Methodology

2013 ◽  
Vol 2013 ◽  
pp. 1-9 ◽  
Author(s):  
Mario Matijević ◽  
Dubravko Pevec ◽  
Krešimir Trontl
Keyword(s):  
Cross Sections ◽  
Neutron Transport ◽  
National Laboratory ◽  
General Purpose ◽  
Monte Carlo Code ◽  
Oak Ridge ◽  
Benchmark Data ◽  
Selection Of ◽  
Good Agreement ◽  
Reactor Dosimetry

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.

Sign in / Sign up

Export Citation Format

Share Document