• Nuclear Engineering and Technology
• /
• 제49권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.

• 한국건축시공학회:학술대회논문집
• /
• 한국건축시공학회 2018년도 추계 학술논문 발표대회
• /
• pp.114-115
• /
• 2018

To develop an effective shielding material for spent fuel that emits fast neutrons is necessary. In this study, thermal neutron and fast neutron shielding performance of polyethylene coated boron carbide-incorporated cement paste was quantitatively analyzed by Monte Carlo N-Particle transport code (MCNP) simulations. As the results of the simulations, fast neutrons were effectively shielded through large quantity of hydrogen and boron elements in polyethylene and boron carbide.

• Nuclear Engineering and Technology
• /
• 제38권7호
• /
• pp.665-674
• /
• 2006

Three-dimensional neutron flux and fluence of Kori Unit 3 were evaluated using the synthesis technique described in Regulatory Guide 1.190 for all reactor geometry. For this purpose DORT neutron transport calculations from Cycle 1 to Cycle 15 were performed using BUGLE-96 cross-section library. The calculated flux and fluence were validated by comparing the calculated reaction rates to the measurement data from the dosimetry sensor set of the $5^{th}$ surveillance capsule withdrawn at the end of cycle 15 of Kori Unit 3. And then the best estimation of the neutron exposures for the reactor vessel beltline region was performed using the least square evaluation. These results can be used in the assessment of the state of embrittlement of Kori Unit 3 pressure vessel.

• Nuclear Engineering and Technology
• /
• 제56권3호
• /
• pp.916-921
• /
• 2024

The control rod worth sensitivity to the neutron detector position in the rod drop experiment is studied by the time-dependent Monte Carlo (TDMC) neutron transport calculations for AGN-201K educational reactor and the Kyoto University Critical Assembly. The TDMC simulations of the rod drop experiments are conducted by the Seoul National University Monte Carlo (MC) code, McCARD, yielding time-dependent neutron densities at detector positions. The detector-position-dependent results of the total control rod worth calculated by the extrapolation, the integral counting, and the inverse methods are compared with the numerical reference using the MC eigenvalue calculations and the experimental results. From these comparisons, it is observed that the total control rod worth can be estimated with a considerable difference depending on the detector position through the rod drop experiment. The proposed TDMC simulation of the rod drop experiment can be applied for searching a better detector position or quantifying a bias for the control rod worth measurement.

• Nuclear Engineering and Technology
• /
• 제47권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
• /
• 제54권9호
• /
• pp.3478-3487
• /
• 2022

In this study, we incorporate an anisotropic scattering scheme involving spherical harmonics into the method of characteristics (MOC). The neutron transport solution in a light water reactor can be significantly improved because of the impact of an anisotropic scattering source with the MOC flat source approximation. Several problems are selected to verify the proposed scheme and investigate its effects and accuracy. The MOC anisotropic scattering source is based on the expansion of spherical harmonics with Legendre polynomial functions. The angular flux, scattering source, and cross section are expanded in terms of the surface spherical harmonics. Later, the polynomial is expanded to achieve the odd and even parity of the source components. Ultimately, the MOC angular and scalar fluxes are calculated from a combination of two sources. This paper presents various numerical examples that represent the hot and cold conditions of a reactor core with boron concentration, burnable absorbers, and control rod materials, with and without a reflector or baffle. Moreover, a small critical core problem is considered which involves significant neutron leakage at room temperature. We demonstrate that an anisotropic scattering source significantly improves solution accuracy for the small core high-leakage problem, as well as for practical large core analyses.

• Nuclear Engineering and Technology
• /
• 제49권6호
• /
• pp.1211-1218
• /
• 2017

Monte Carlo neutron transport codes generally use the method of successive generations to converge the fission source distribution to-and then maintain it at-the fundamental mode. Recently, a phenomenon called "clustering" has been noted, which produces fission distributions that are very far from the fundamental mode. In this study, a mathematical model of clustering in Monte Carlo has been developed. The model draws on previous work for continuous-time birth-death processes, as well as methods from the field of population genetics.

• Nuclear Engineering and Technology
• /
• 제52권8호
• /
• pp.1784-1791
• /
• 2020

The 14 MeV neutron generator facility is being developed by the Institute for Plasma Research India to conduct the lab scale experiments related to Indian breeding blanket system for ITER and DEMO. It will also be utilized for material testing, shielding experiments and development of fusion diagnostics. Occupational radiation exposure control is necessary for the all kind of nuclear facilities to get the operational licensing from governing authorities and nuclear regulatory bodies. In the same way, the radiation exposure for the 14 MeV neutron generator facility at the occupational worker area and accessible zones for general workers should be under the permissible limit of AERB India. The generator is designed for the yield of 10¹² n/s. The shielding assessment has been made to estimate the radiation dose during the operational time of the neutron generator. The facility has many utilities and constraints like ventilation ducts, accessible doors, accessibility of neutron generator components and to conduct the experiments which make the shielding assessment challenging to provide proper safety for occupational workers and the general public. The neutron and gamma dose rates have been estimated using the MCNP radiation transport code and ENDF -VII nuclear data libraries. The ICRP-74 fluence to dose conversion coefficients has been used for the assessment. The annual radiation exposure has been assessed by considering 500 h per year operational time. The provision of local shield near to neutron generator has been also evaluated to reduce the annual radiation doses. The comprehensive results of radiation shielding capability of neutron generator building and local shield design have been presented in the paper along with detailed maps of radiation field.

• Nuclear Engineering and Technology
• /
• 제48권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
• /
• 제55권4호
• /
• pp.1280-1286
• /
• 2023

MgF₂ as a moderator material has been extensively used in the beam shaping assembly (BSA) that plays an important role in the boron neutron capture therapy (BNCT). Regarded as important for applications, the thermal neutron scattering data of MgF₂ were calculated, based on the phonon expansion model. The structural properties of MgF₂ were researched by the VASP code based on the ab-initio methods. The PHONOPY code was employed to calculate the phonon density of states. Furthermore, the NJOY code was used to calculate the thermal neutron scattering data of MgF₂. The calculated inelastic cross sections plus absorption cross sections are in agreement with the available experimental data. The neutron transport in the BSA has been simulated by using a hybrid Monte-Carlo-Deterministic code NECP-MCX. The results indicated that compared with the calculation of the free gas model, the thermal neutron flux and epithermal neutron flux at the BSA exit port calculated by using the thermal neutron scattering data of MgF₂ were reduced by 27.7% and 8.2%, respectively.