Vectorization and multitasking with a Monte Carlo code for neutron transport problems

2006 ◽  
pp. 234-244
Author(s):  
Yves Chauvet
Keyword(s):  
Monte Carlo ◽  
Neutron Transport ◽  
Monte Carlo Code ◽  
Transport Problems

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.

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 Simulating neutron transport in long beamlines at a spallation neutron source using Geant4

2020 ◽  
Vol 22 (2-3) ◽  
pp. 183-189
Author(s):  
Douglas D. DiJulio ◽  
Isak Svensson ◽  
Xiao Xiao Cai ◽  
Joakim Cederkall ◽  
Phillip M. Bentley
Keyword(s):  
Monte Carlo ◽  
Neutron Source ◽  
Variance Reduction ◽  
Neutron Transport ◽  
High Energy ◽  
Low Energy ◽  
Deep Penetration ◽  
Monte Carlo Code ◽  
Spallation Neutron Source ◽  
Unique Challenge

The transport of neutrons in long beamlines at spallation neutron sources presents a unique challenge for Monte-Carlo transport calculations. This is due to the need to accurately model the deep-penetration of high-energy neutrons through meters of thick dense shields close to the source and at the same time to model the transport of low- energy neutrons across distances up to around 150 m in length. Typically, such types of calculations may be carried out with MCNP-based codes or alternatively PHITS. However, in recent years there has been an increased interest in the suitability of Geant4 for such types of calculations. Therefore, we have implemented supermirror physics, a neutron chopper module and the duct-source variance reduction technique for low- energy neutron transport from the PHITS Monte-Carlo code into Geant4. In the current work, we present a series of benchmarks of these extensions with the PHITS software, which demonstrates the suitability of Geant4 for simulating long neutron beamlines at a spallation neutron source, such as the European Spallation Source, currently under construction in Lund, Sweden.

scholarly journals MULTI-PHYSICS TRANSIENT SIMULATIONS WITH TRIPOLI-4®

2021 ◽  
Vol 247 ◽  
pp. 07019 ◽  
Author(s):  
Margaux Faucher ◽  
Davide Mancusi ◽  
Andrea Zoia
Keyword(s):  
Monte Carlo ◽  
Neutron Transport ◽  
Coolant Temperature ◽  
Test Case ◽  
Monte Carlo Code ◽  
Thermal Hydraulics ◽  
Fuel Temperature ◽  
Control Rod ◽  
Transient Simulations ◽  
Neutron Power

In this work, we present the first dynamic calculations performed with the Monte Carlo neutron transport code TRIPOLI-4R with thermal-hydraulics feedback. For this purpose, the Monte Carlo code was extended for multi-physics capabilities and coupled to the thermal-hydraulics subchannel code SUBCHANFLOW. As a test case for the verification of transient simulation capabilities, a 3x3-assembly mini-core benchmark based on the TMI-1 reactor is considered with a pin-by-pin description. Two reactivity excursion scenarios initiated by control-rod movement are simulated starting from a critical state and compared to analogous simulations performed using the Serpent 2 Monte-Carlo code. The time evolution of the neutron power, fuel temperature, coolant temperature and coolant density are analysed to assess the multi-physics capabilities of TRIPOLI-4. The stabilizing e_ects of thermal-hydraulics on the neutron power appear to be well taken into account. The computational requirements for massively parallel calculations are also discussed.

Sign in / Sign up

Export Citation Format

Share Document