• Nuclear Engineering and Technology
• /
• v.54 no.5
• /
• pp.1769-1780
• /
• 2022

A new method of dose calculation algorithm, called GPU-accelerated Monte Carlo and collapsed cone Convolution (GMCC) was developed to improve the calculation speed of BNCT treatment planning system. The GPU-accelerated Monte Carlo routine in GMCC is used to simulate the neutron transport over whole energy range and the Collapsed Cone Convolution method is to calculate the gamma dose. Other dose components due to alpha particles and protons, are calculated using the calculated neutron flux and reaction data. The mathematical principle and the algorithm architecture are introduced. The accuracy and performance of the GMCC were verified by comparing with the FLUKA results. A water phantom and a head CT voxel model were simulated. The neutron flux and the absorbed dose obtained by the GMCC were consistent well with the FLUKA results. In the case of head CT voxel model, the mean absolute percentage error for the neutron flux and the absorbed dose were 3.98% and 3.91%, respectively. The calculation speed of the absorbed dose by the GMCC was 56 times faster than the FLUKA code. It was verified that the GMCC could be a good candidate tool instead of the Monte Carlo method in the BNCT dose calculations.

• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.499-507
• /
• 2020

Research reactors for radioisotope production, fuel and material testing and research activities are designed, constructed and operated based on the society's needs. In this study, neutronic and thermal hydraulic design of a high neutron flux research reactor core for radioisotope production is presented. Main parameters including core excess reactivity, reactivity variations, power and flux distribution during the cycle, axial and radial power peaking factors (PPF), Pu₂₃₉ production and minimum DNBR are calculated by nuclear deterministic codes. Core calculations performed by deterministic codes are validated with Monte Carlo code. Comparison of the neutronic parameters obtained from deterministic and Monte Carlo codes indicates good agreement. Finally, subchannel analysis performed for the hot channel to evaluate the maximum fuel and clad temperatures. The results show that the average thermal neutron flux at the beginning of cycle (BOC) is 1.0811 × 10¹⁴ n/㎠-s and at the end of cycle (EOC) is 1.229 × 10¹⁴ n/㎠-s. Total Plutonium (Pu₂₃₉) production at the EOC evaluated to be 0.9487 Kg with 83.64% grade when LEU (UO₂ with 3.7% enrichment) used as fuel. This designed reactor which uses LEU fuel and has high neutron flux and low plutonium production could be used for peaceful nuclear activities based on nuclear non-proliferation treaty concepts.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1165-1171
• /
• 2017

The traditional form of parallelism in Monte Carlo particle transport simulations, wherein each individual particle history is considered a unit of work, does not lend itself well to data-level parallelism. Event-based algorithms, which were originally used for simulations on vector processors, may offer a path toward better utilizing data-level parallelism in modern computer architectures. In this study, a simple model is developed for estimating the efficiency of the event-based particle transport algorithm under two sets of assumptions. Data collected from simulations of four reactor problems using OpenMC was then used in conjunction with the models to calculate the speedup due to vectorization as a function of the size of the particle bank and the vector width. When each event type is assumed to have constant execution time, the achievable speedup is directly related to the particle bank size. We observed that the bank size generally needs to be at least 20 times greater than vector size to achieve vector efficiency greater than 90%. When the execution times for events are allowed to vary, the vector speedup is also limited by differences in the execution time for events being carried out in a single event-iteration.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1999.10a
• /
• pp.49.1-49
• /
• 1999

• Nuclear Engineering and Technology
• /
• v.56 no.2
• /
• pp.575-600
• /
• 2024

During an earthquake on December 29th 2020, the Krško NPP automatically shutdown due to the trigger of the negative neutron flux rate signal on the power range nuclear instrumentation. From the time course of the detector signal, it can be concluded that the fluctuation in the detector signal may have been caused by the mechanical movement of the ex-core neutron detectors or the pressure vessel components rather than the actual change in reactor power. The objective of the analysis was to evaluate the sensitivity of the neutron flux at the ex-core detector position, if the detector is moved in the radial or axial direction. In addition, the effect of the core barrel movement and core inside the baffle movement in the radial direction were analysed. The analysis is complemented by the calculation of the thermal and total neutron flux gradient in radial, axial and azimuthal directions. The Monte Carlo particle transport code MCNP was used to study the changes in the response of the ex-core detector for the above-mentioned scenarios. Power and intermediate-range detectors were analysed separately, because they are designed differently, positioned at different locations, and have different response characteristics. It was found that the movement of the power range ex-core detector has a negligible effect on the value of the thermal neutron flux in the active part of the detector. However, the radial movement of the intermediate-range detector by 5 cm results in 7%-8% change in the thermal neutron flux in the active part of the intermediate-range detector. The analysis continued with an evaluation of the effects of moving the entire core barrel on the ex-core detector response. It was estimated that the 2 mm core barrel radial oscillation results in ~4% deviation in the power and intermediate-range detector signal. The movement of the reactor core inside baffle can contribute ~6% deviation in the ex-core neutron detector signal. The analysis showed that the mechanical movement of ex-core neutron detectors cannot explain the fluctuations in the ex-core detector signal. However, combined core barrel and reactor core inside baffle oscillations could be a probable reason for the observed fluctuations in the ex-core detector signal during an earthquake.

• Nuclear Engineering and Technology
• /
• v.55 no.5
• /
• pp.1593-1603
• /
• 2023

In deterministic and Monte Carlo transport codes, b-delayed emission is included using a group structure where all of the precursors are grouped together in 6 groups or families, but given the increase in computational power, nowadays there is no reason to keep this structure. Furthermore, there have been recent efforts to compile and evaluate all the available b-delayed neutron emission data and to measure new and improved data on individual precursors. In order to be able to perform a transient Monte Carlo simulation, data from individual precursors needs to be implemented in a transport code. This work is the first step towards the development of a tool to explore the effect of individual precursors in a fissile system. In concrete, individual precursor data is included by expanding the capabilities of the open source Monte Carlo code OpenMC. In the modified code - named Time Dependent OpenMC or OpenMC(TD)- time dependency related to β-delayed neutron emission was handled by using forced decay of precursors and combing of the particle population. The data for continuous energy neutron cross-sections was taken from JEFF-3.1.1 library. Regarding the data needed to include the individual precursors, cumulative yields were taken from JEFF-3.1.1 and delayed neutron emission probabilities and delayed neutron spectra were taken from ENDF-B/VIII.0. OpenMC(TD) was tested in a monoenergetic system, an energy dependent unmoderated system where the precursors were taken individually or in a group structure, and in a light-water moderated energy dependent system, using 6-groups, 50 and 40 individual precursors. Neutron flux as a function of time was obtained for each of the systems studied. These results show the potential of OpenMC(TD) as a tool to study the impact of individual precursor data on fissile systems, thus motivating further research to simulate more complex fissile systems.

• Journal of Energy Engineering
• /
• v.15 no.4 s.48
• /
• pp.250-255
• /
• 2006

In this study, the discrete ordinates method which has been widely used in the solution of neutron transport equation is applied to the solution of the time dependent radiative transfer equation. The self-adjoint form of the second order radiation intensity equation is used to enhance the stability of the solution, and a new multi-step linearization method is developed to avoid the nonlinearity in the material temperature equation. This new solution method is applied to the well known Marshak wave problem, and the numerical result is compared with that of the conventional Monte-Carlo method.

• Nuclear Engineering and Technology
• /
• v.48 no.3
• /
• pp.785-794
• /
• 2016

The Argonne National Laboratory of the United States and the Kharkov Institute of Physics and Technology of the Ukraine have been collaborating on the design, development and construction of a neutron source facility at Kharkov Institute of Physics and Technology utilizing an electron-accelerator-driven subcritical assembly. The electron beam power is 100 kW using 100-MeV electrons. The facility was designed to perform basic and applied nuclear research, produce medical isotopes, and train nuclear specialists. The biological shield of the accelerator building was designed to reduce the biological dose to less than 5.0e-03 mSv/h during operation. The main source of the biological dose for the accelerator building is the photons and neutrons generated from different interactions of leaked electrons from the electron gun and the accelerator sections with the surrounding components and materials. The Monte Carlo N-particle extended code (MCNPX) was used for the shielding calculations because of its capability to perform electron-, photon-, and neutron-coupled transport simulations. The photon dose was tallied using the MCNPX calculation, starting with the leaked electrons. However, it is difficult to accurately tally the neutron dose directly from the leaked electrons. The neutron yield per electron from the interactions with the surrounding components is very small, ~0.01 neutron for 100-MeV electron and even smaller for lower-energy electrons. This causes difficulties for the Monte Carlo analyses and consumes tremendous computation resources for tallying the neutron dose outside the shield boundary with an acceptable accuracy. To avoid these difficulties, the SOURCE and TALLYX user subroutines of MCNPX were utilized for this study. The generated neutrons were banked, together with all related parameters, for a subsequent MCNPX calculation to obtain the neutron dose. The weight windows variance reduction technique was also utilized for both neutron and photon dose calculations. Two shielding materials, heavy concrete and ordinary concrete, were considered for the shield design. The main goal is to maintain the total dose outside the shield boundary less than 5.0e-03 mSv/h during operation. The shield configuration and parameters of the accelerator building were determined and are presented in this paper.

• Nuclear Engineering and Technology
• /
• v.47 no.2
• /
• pp.165-175
• /
• 2015

The application of neutron noise analysis (NNA) to the ex-core neutron detector signal for monitoring the vibration characteristics of a reactor core support barrel (CSB) was investigated. Ex-core flux data were generated by using a nonanalog Monte Carlo neutron transport method in a simulated CSB model where the implicit capture and Russian roulette technique were utilized. First and third order beam and shell modes of CSB vibration were modeled based on parallel processing simulation. A NNA module was developed to analyze the ex-core flux data based on its time variation, normalized power spectral density, normalized cross-power spectral density, coherence, and phase differences. The data were then analyzed with a fuzzy logic module to determine the vibration characteristics. The ex-core neutron signal fluctuation was directly proportional to the CSB's vibration observed at 8Hz and15Hzin the beam mode vibration, and at 8Hz in the shell mode vibration. The coherence result between flux pairs was unity at the vibration peak frequencies. A distinct pattern of phase differences was observed for each of the vibration models. The developed fuzzy logic module demonstrated successful recognition of the vibration frequencies, modes, orders, directions, and phase differences within 0.4 ms for the beam and shell mode vibrations.

• Nuclear Engineering and Technology
• /
• v.51 no.8
• /
• pp.2018-2025
• /
• 2019

This study systematically compared two hybrid deterministic/Monte Carlo transport codes, ADVANTG/MCNP and MAVRIC, in solving a difficult shielding problem for a real-world spent fuel storage cask. Both hybrid codes were developed based on the consistent adjoint driven importance sampling (CADIS) methodology but with different implementations. The dose rate distributions on the cask surface were of primary interest and their predicted results were compared with each other and with a straightforward MCNP calculation as a baseline case. Forward-Weighted CADIS was applied for optimization toward uniform statistical uncertainties for all tallies on the cask surface. Both ADVANTG/MCNP and MAVRIC achieved substantial improvements in overall computational efficiencies, especially for gamma-ray transport. Compared with the continuous-energy ADVANTG/MCNP calculations, the coarse-group MAVRIC calculations underestimated the neutron dose rates on the cask's side surface by an approximate factor of two and slightly overestimated the dose rates on the cask's top and side surfaces for fuel gamma and hardware gamma sources because of the impact of multigroup approximation. The fine-group MAVRIC calculations improved to a certain extent and the addition of continuous-energy treatment to the Monte Carlo code in the latest MAVRIC sequence greatly reduced these discrepancies. For the two continuous-energy calculations of ADVANTG/MCNP and MAVRIC, a remaining difference of approximately 30% between the neutron dose rates on the cask's side surface resulted from inconsistent use of thermal scattering treatment of hydrogen in concrete.