• 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.

• Nuclear Engineering and Technology
• /
• v.48 no.3
• /
• pp.751-757
• /
• 2016

The present work discusses two different models of boiling water reactor (BWR) bundle to compare the neutronic characteristics of uranium dioxide ($UO_2$) and uranium zirconium hydride ($UZrH_{1.6}$) fuel. Each bundle consists of four assemblies. The BWR assembly fueled with $UO_2$ contains $8{\times}8$ fuel rods while that fueled with $UZrH_{1.6}$ contains $9{\times}9$ fuel rods. The Monte Carlo N-Particle Transport code, based on the Mont Carlo method, is used to design three dimensional models for BWR fuel bundles at typical operating temperatures and pressure conditions. These models are used to determine the multiplication factor, pin-by-pin power distribution, axial power distribution, thermal neutron flux distribution, and axial thermal neutron flux. The moderator and coolant (water) are permitted to boil within the BWR core forming steam bubbles, so it is important to calculate the reactivity effect of voiding at different values. It is found that the hydride fuel bundle design can be simplified by eliminating water rods and replacing the control blade with control rods. $UZrH_{1.6}$ fuel improves the performance of the BWR in different ways such as increasing the energy extracted per fuel assembly, reducing the uranium ore, and reducing the plutonium accumulated in the BWR through burnup.

• Nuclear Engineering and Technology
• /
• v.50 no.7
• /
• pp.1051-1059
• /
• 2018

Fast periodic pulsed reactor is a type of reactor in which the fission bursts are formed entirely with external reactivity modulation with a specified time periodicity. This type of reactors could generate much larger intensity of neutron beams for experimental use, compared with the steady state reactors. In the design of fast periodic pulsed reactors, the time dependent simulation of the power pulse is majorly based on a point kinetic model, which is known to have limitations. A more accurate calculation method is desired for the design analyses of fast periodic pulsed reactors. Monte Carlo computer code MCNP6 is used for this task due to its three dimensional transport capability with a continuous energy library. Some new routines were added to simulate the rotation of the movable reflector parts in the time dependent calculation. Fast periodic pulsed reactor IBR-2M was utilized to validate the new routines. This reactor is periodically in prompt supercritical state, which lasts for ${\sim}400{\mu}s$, during the equilibrium state. This generates long neutron fission chains, which requires tremendously large amount of computation time during Monte Carlo simulations. Russian Roulette was applied for these very long neutron chains in MCNP6 calculation, combined with other approaches to improve the efficiency of the simulations. In the power pulse of the IBR-2M at equilibrium state, there is some discrepancy between the experimental measurements and the calculated results using the point kinetics model. MCNP6 results matches better the experimental measurements, which shows the merit of using MCNP6 calculation relative to the point kinetics model.

• Nuclear Engineering and Technology
• /
• v.56 no.2
• /
• pp.673-684
• /
• 2024

STREAM, developed by the Computational Reactor Physics and Experiment laboratory (CORE) of the Ulsan National Institute of Science and Technology (UNIST), is a deterministic neutron- and photon-transport code primarily designed for light water reactor (LWR) analysis. Initially, the photon module in STREAM did not account for fluorescence and bremsstrahlung photons. This article presents recent developments regarding the integration of atomic relaxation and bremsstrahlung models into the existing photon module, thus allowing for the transport of secondary photons. The photon flux and photon heating computed with the newly incorporated models is compared to results obtained with the Monte Carlo code MCS. The incorporation of secondary photons has substantially improved the accuracy of photon flux calculations, particularly in scenarios involving strong gamma emitters. However, it is essential to note that despite the consideration of secondary photon sources, there is no noticeable improvement in the photon heating for LWR problems when compared to the photon heating obtained with the previous version of STREAM.

• Journal of Radiation Protection and Research
• /
• v.38 no.1
• /
• pp.22-28
• /
• 2013

The skyshine effect is an essential and important phenomenon in the shielding design of the high energy accelerator. In this study, a new estimation method of neutron skyshine was proposed and was verified by comparison with existing methods. The effective dose of secondary neutrons and photons at the locations that was far away from high-energy electron accelerator was calculated using FLUKA and PHITS Monte Carlo code. The transport paths of secondary radiations to reach a long distance were classified as skyshine, direct, groundshine and multiple-shine. The contribution of each classified component to the total effective dose was evaluated. The neutrons produced from the thick copper target irradiated by 10 GeV electron beam was applied as a source term of this transport. In order to evaluate a groundshine effect, the composition of soil on the PAL-XFEL site was considered. At a relatively short distance less than 50 m from the accelerator tunnel, the direct and groundshine components mostly contributed to the total effective dose. The skyshine component was important at a long distance. The evaluated dose of neutron skyshine agreed better with the results using Rindi's formula, which was based on the experimental results at high energy electron accelerator. That also agreed with the estimated dose using the simple evaluation code, SHINE3, within about 20%. The total effective dose, including all components, was 10 times larger than the estimated doses using other methods for this comparison. The influence of multiple-shine path in this evaluation of the estimation method was investigated to be bigger than one of pure skyshine path.

• Nuclear Engineering and Technology
• /
• v.48 no.3
• /
• pp.765-772
• /
• 2016

In analyzing residual radiation, researchers generally use a two-step Monte Carlo (MC) simulation. The first step (MC1) simulates neutron transport, and the second step (MC2) transports the decay photons emitted from the activated materials. In this process, the stochastic uncertainty estimated by the MC2 appears only as a final result, but it is underestimated because the stochastic error generated in MC1 cannot be directly included in MC2. Hence, estimating the true stochastic uncertainty requires quantifying the propagation degree of the stochastic error in MC1. The brute force technique is a straightforward method to estimate the true uncertainty. However, it is a costly method to obtain reliable results. Another method, called the adjoint-based method, can reduce the computational time needed to evaluate the true uncertainty; however, there are limitations. To address those limitations, we propose a new strategy to estimate uncertainty propagation without any additional calculations in two-step MC simulations. To verify the proposed method, we applied it to activation benchmark problems and compared the results with those of previous methods. The results show that the proposed method increases the applicability and user-friendliness preserving accuracy in quantifying uncertainty propagation. We expect that the proposed strategy will contribute to efficient and accurate two-step MC calculations.

• Radiation Oncology Journal
• /
• v.19 no.1
• /
• pp.66-73
• /
• 2001

Purpose : For the research of Boron Neutron Capture Therapy (BNCT), fast neutrons generated from the MC-50 cyclotron with maximum energy of 34.4 MeV in Korea Cancer Center Hospital were moderated by 70 cm paraffin and then the dose characteristics were investigated. Using these results, we hope to establish the protocol about dose measurement of epi-thermal neutron, to make a basis of dose characteristic of epi-thermal neutron emitted from nuclear reactor, and to find feasibility about accelerator-based BNCT. Method and Materials : For measuring the absorbed dose and dose distribution of fast neutron beams, we used Unidos 10005 (PTW, Germany) electrometer and IC-17 (Far West, USA), IC-18, ElC-1 ion chambers manufactured by A-150 plastic and used IC-l7M ion chamber manufactured by magnesium for gamma dose. There chambers were flushed with tissue equivalent gas and argon gas and then the flow rate was S co per minute. Using Monte Carlo N-Particle (MCNP) code, transport program in mixed field with neutron, photon, electron, two dimensional dose and energy fluence distribution was calculated and there results were compared with measured results. Results : The absorbed dose of fast neutron beams was $6.47\times10^{-3}$ cGy per 1 MU at the 4 cm depth of the water phantom, which is assumed to be effective depth for BNCT. The magnitude of gamma contamination intermingled with fast neutron beams was $65.2{\pm}0.9\%$ at the same depth. In the dose distribution according to the depth of water, the neutron dose decreased linearly and the gamma dose decreased exponentially as the depth was deepened. The factor expressed energy level, $D_{20}/D_{10}$, of the total dose was 0.718. Conclusion : Through the direct measurement using the two ion chambers, which is made different wall materials, and computer calculation of isodose distribution using MCNP simulation method, we have found the dose characteristics of low fluence fast neutron beams. If the power supply and the target material, which generate high voltage and current, will be developed and gamma contamination was reduced by lead or bismuth, we think, it may be possible to accelerator-based BNCT.

• Bulletin of the Korean Space Science Society
• /
• 2011.04a
• /
• pp.28.1-28.1
• /
• 2011

Radiation-induced displacement damage which has caused the increase of the dark current in the focal plane adopted in the Ozone Monitoring Instrument (OMI) was studied in regards of the primary protons and the secondaries generated by the protons in the orbit. By using the Monte Carlo N-Particle Transport Code System (MCNPX) version 2.4.0 along with the Stopping and Range of Ions in Matter version 2010 (SRIM2010), effects of the primary protons as well as secondary particles including neutron, electron, and photon were investigated. After their doses and fluxes that reached onto the charge-coupled device (CCD) were examined, displacement damage induced by major sources was presented.

• Journal of Environmental Science International
• /
• v.13 no.2
• /
• pp.175-180
• /
• 2004

A prompt gamma-ray neutron activation (PGNA) system was simulated by the Monte Carlo N-Particle transport code (MCNP-4A) to estimate the level at which the scattered photon fluence rate, the absolute efficiency of the HPGe-detector, the volume of the concrete sample and the $^{35}$ /Cl(n, ${\gamma}$) reaction rate in this sample contribute to the count rate in the NaCl concentration measurement. The n- ${\gamma}$ fluence rates at the ST-2 beam tube exit of the HANARO reactor were used as input data, and the GAMMA-X type HPGe detector was modeled to tally 1.1649 MeV ${\gamma}$ -rays emitted from the $^{35}$ Cl(n, ${\gamma}$) reaction in the concrete sample. For three cylindrical concrete samples of 13.8, 46.8 and 157.1 ㎤ volumes, respectively, the relations between the NaCl weight fractions of 0.1, 1, 2 and 5 % in each of the concrete samples and the 1.1 649 MeV pulses created in the HPGe detector model were studied. As a result, it was found that the count rate at the same NaCl concentration nearly depends on the volume of the samples in a simulated condition of the same NaCl concentration samples, and that the linearities of the NaCl concentration calibration curves were reasonable in the narrow range of the NaCl weight fraction.

• Nuclear Engineering and Technology
• /
• v.46 no.3
• /
• pp.343-352
• /
• 2014

Uncertainty and sensitivity analyses with respect to nuclear data are performed with depletion calculations for BWR and PWR fuel assemblies specified in the framework of the UAM-LWR Benchmark Phase II. For this, the GRS sampling based tool XSUSA is employed together with the TRITON depletion sequences from the SCALE 6.1 code system. Uncertainties for multiplication factors and nuclide inventories are determined, as well as the main contributors to these result uncertainties by calculating importance indicators. The corresponding neutron transport calculations are performed with the deterministic discrete-ordinates code NEWT. In addition, the Monte Carlo code KENO in multi-group mode is used to demonstrate a method with which the number of neutron histories per calculation run can be substantially reduced as compared to that in a calculation for the nominal case without uncertainties, while uncertainties and sensitivities are obtained with almost the same accuracy.