• Nuclear Engineering and Technology
• /
• v.49 no.5
• /
• pp.920-927
• /
• 2017

With ever-advancing computer technology, the Monte Carlo (MC) neutron transport calculation is expanding its application area to nuclear reactor transient analysis. Dynamic MC (DMC) neutron tracking for transient analysis requires efficient algorithms for delayed neutron generation, neutron population control, and initial condition modeling. In this paper, a new MC steady-state simulation method based on time-dependent MC neutron tracking is proposed for steady-state initial condition modeling; during this process, prompt neutron sources and delayed neutron precursors for the DMC transient simulation can easily be sampled. The DMC method, including the proposed time-dependent DMC steady-state simulation method, has been implemented in McCARD and applied for two-dimensional core kinetics problems in the time-dependent neutron transport benchmark C5G7-TD. The McCARD DMC calculation results show good agreement with results of a deterministic transport analysis code, nTRACER.

• Nuclear Engineering and Technology
• /
• v.54 no.9
• /
• pp.3181-3204
• /
• 2022

The variational nodal method for solving the neutron transport equation has evolved over 40 years. Based on a functional form of the Boltzmann neutron transport equation, the method now comprises a complete set of variants that can be employed for different problems. This paper presents an extensive review of the development of the variational nodal method. The emphasis is on summarizing the whole theoretical system rather than validating the methodologies. The paper covers the variational nodal formulation of the Boltzmann neutron transport equation, the Ritz procedure for various application purposes, the derivation of boundary conditions, the extension for adjoint and perturbation calculations, and treatments for anisotropic scattering sources. Acceleration approaches for constructing response matrices and solving the resulting system of algebraic equations are also presented.

• Nuclear Engineering and Technology
• /
• v.52 no.12
• /
• pp.2852-2859
• /
• 2020

This paper deals with the simulation-based design optimization and experimental validation of the characteristics of an in-core silver Self-Powered Neutron Detector (SPND). Optimized dimensions of the SPND are determined by combining Monte Carlo simulations and analytical methods. As a first step, the Monte Carlo transport code MCNPX is used to follow the trajectory and fate of the neutrons emitted from an external source. This simulation is able to seamlessly integrate various phenomena, including neutron slowing-down and shielding effects. Then, the expected number of beta particles and their energy spectrum following a neutron capture reaction in the silver emitter are fetched from the TENDEL database using the JANIS software interface and integrated with the data from the first step to yield the origin and spectrum of the source electrons. Eventually, the MCNPX transport code is used for the Monte Carlo calculation of the ballistic current of beta particles in the various regions of the SPND. Then, the output current and the maximum insulator thickness to avoid breakdown are determined. The optimum design of the SPND is then manufactured and experimental tests are conducted. The calculated design parameters of this detector have been found in good agreement with the obtained experimental results.

• Nuclear Engineering and Technology
• /
• v.53 no.1
• /
• pp.30-43
• /
• 2021

In the 2D/1D method, a global adjoint CMFD based on the generalized equivalence theory is built to synthesize the 2D radial MOC adjoint and 1D axial NEM adjoint calculation and also to accelerate the iteration convergence of 3D whole-core adjoint transport calculation. Even more important, an advanced yet accurate two-level (TL) CMFD acceleration technique is proposed, in which an equivalent one-group adjoint CMFD is established to accelerate the multi-group adjoint CMFD and then to accelerate the 3D whole-core adjoint transport calculation efficiently. Based on these method, a new code is developed to perform 3D adjoint neutron flux calculation. Then a set of VERA and C5G7 benchmark problems are chosen to verify the capability of the 3D adjoint calculations and the effectiveness of TL CMFD acceleration. The numerical results demonstrate that acceptable accuracy of 2D/1D adjoint calculations and superior acceleration of TL CMFD are achievable.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3506-3513
• /
• 2023

In this paper, we demonstrate the neutron attenuation of dry cask shielding materials using the S670e helium-4 detector manufactured by Arktis Radiation Ltd. In particular, two materials expected to be applied to the TN-32 dry cask manufactured by ORANO Korea and KORAD-21 by the Korea Radioactive Waste Agency (KORAD) were utilized. The measured neutron attenuation was compared with our Monte Carlo N-Particle Transport simulation results, and the difference is given as the root mean square (RMS). For the fast neutron case, a rapid decline in neutron counts was observed as a function of increasing material thickness, exhibiting an exponential relationship. The discrepancy between the experimentally acquired data and simulation results for the fast neutron was maintained within a 2.3% RMS. In contrast, the observed thermal neutron count demonstrated an initial rise, attained a maximum value, and exhibited an exponential decline as a function of increasing thickness. In particular, the discrepancy between the measured and simulated peak locations for thermal neutrons displayed an RMS deviation of approximately 17.3-22.4%. Finally, the results suggest that a minimum thickness of 5 cm for Li-6 is necessary to achieve a sufficiently significant cross-section, effectively capturing incoming thermal neutrons within the dry cask.

• Geophysics and Geophysical Exploration
• /
• v.15 no.4
• /
• pp.219-226
• /
• 2012

This study provides an analysis on results of neutron logging for various borehole environments through numerical simulation based on a Monte Carlo N-Particle (MCNP) code developed and maintained by Los Alamos National Laboratory. MCNP is suitable for the simulation of neutron logging since the algorithm can simulate transport of nuclear particles in three-dimensional geometry. Rather than simulating a specific tool of a particular service company between many commercial neutron tools, we have constructed a generic thermal neutron tool characterizing commercial tools. This study makes calibration chart of the neutron logging tool for materials (e.g., limestone, sandstone and dolomite) with various porosities. Further, we provides correction charts for the generic neutron logging tool to analyze responses of the tool under various borehole conditions by considering brine-filled borehole fluid and void water, and presence of borehole fluid.

• Nuclear Engineering and Technology
• /
• v.55 no.2
• /
• pp.756-768
• /
• 2023

Neutron transport calculations are extremely challenging due to the high computational cost of large and complex problems. A multilevel octree grid algorithm (MLTG) of discrete ordinates method was developed to improve the modeling accuracy and simulation efficiency on 3-D Cartesian grids. The Balakovo-3 VVER-1000 neutron dosimetry benchmark is calculated to verify and validate this numerical technique. A simplified S₂ synthetic acceleration is used in the MLTG calculation method to improve the convergence of the source iterations. For the triangularly arranged fuel pins, we adopt a source projection algorithm to generate pin-by-pin source distributions of hexagonal assemblies. MLTG provides accurate geometric modeling and flexible fixed source description at a lower cost than traditional Cartesian grids. The total number of meshes is reduced to 1.9 million from the initial 9.5 million for the Balakovo-3 model. The numerical comparisons show that the MLTG results are in satisfactory agreement with the conventional S_N method and experimental data, within the root-mean-square errors of about 4% and 10%, respectively. Compared to uniform fine meshing, approximately 70% of the computational cost can be saved using the MLTG algorithm for the Balakovo-3 computational model.

• Proceedings of the Korean Institute of Building Construction Conference
• /
• 2018.11a
• /
• 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
• /
• v.56 no.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
• /
• 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.