• Nuclear Engineering and Technology
• /
• 제48권3호
• /
• pp.624-634
• /
• 2016

Nuclear criticality safety analyses (NCSAs) considering burnup credit were performed for the GBC-32 cask. The used nuclear fuel assemblies (UNFAs) discharged from Hanbit Nuclear Power Plant Unit 3 Cycle 6 were loaded into the cask. Their axial burnup distributions and average discharge burnups were evaluated using the DeCART and Multi-purpose Analyzer for Static and Transient Effects of Reactors (MASTER) codes, and NCSAs were performed using SCALE 6.1/STandardized Analysis of Reactivity for Burnup Credit using SCALE (STARBUCS) and Monte Carlo N-Particle transport code, version 6 (MCNP 6). The axial burnup distributions were determined for 20 UNFAs with various initial enrichments and burnups, which were applied to the criticality analysis for the cask system. The UNFAs for 20- and 30-year cooling times were assumed to be stored in the cask. The criticality analyses indicated that $k_{eff}$ values for UNFAs with nonuniform axial burnup distributions were larger than those with a uniform distribution, that is, the end effects were positive but much smaller than those with the reference distribution. The axial burnup distributions for 20 UNFAs had shapes that were more symmetrical with a less steep gradient in the upper region than the reference ones of the United States Department of Energy. These differences in the axial burnup distributions resulted in a significant reduction in end effects compared with the reference.

• Journal of Radiation Protection and Research
• /
• 제11권1호
• /
• pp.29-36
• /
• 1986

본 연구에서는 NaI검출기에서 감마스펙트럼 시뮤레이션 코드인 GESS를 개발하였다. 감마선에 의한 모든 상호작용이 시뮤레이션 과정에서 고려되었으며, 생성된 하전입자의 스펙트럼은 CSDA 모델에 기반을 두어 계산하였다. 매질내에서 입자 수송에 대한 해석수단으로는 몬테칼로 방법을 적용하였다. 코드의 검증을 위하여 1.33MeV의 입사 감마선에 대한 스펙트럼이 본 연구에서 개발된 코드에 의해 계산되었으며, 계산된 스펙트럼은 대체적으로 실험에서 얻은 스펙트럼과 거의 동일한 분포를 나타내고 있다.

• Nuclear Engineering and Technology
• /
• 제55권5호
• /
• pp.1885-1891
• /
• 2023

The present work aims to fabricate Na_1+xZr₂Si_xP_3-xO₁₂ compound at various calcination temperatures based on the zircon mineral. The fabricated compound was calcinated at 250, 500, and 1000℃. The effect of calcination temperature on the structure, crystal phase, and radiation shielding properties was studied for the fabricated compound. The X-ray diffraction diffractometer demonstrates that, the monoclinic crystal phase appeared at a calcination temperature of 250℃ and 500℃ is totally transformed to a high-symmetry hexagonal crystal phase under a calcination temperature of 1000℃. The radiation shielding capacity was also qualified for the fabricated compounds using the Monte Carlo N-Particle transport code in the g-photons energy interval between 15keV and 122keV. The impacts of calcination temperature on the g-ray shielding behavior were clarified in the present study, where the linear attenuation coefficient was enhanced by 218% at energy of 122keV, when the calcination temperature increased from 250 to 1000℃, respectively.

• 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
• /
• 제43권6호
• /
• pp.567-572
• /
• 2011

A sensitivity analysis of the methods used to evaluate the transport properties of a CdZnTe detector was performed using two different radiations (${\alpha}$ particle and gamma-ray) emitted from an $^{241}Am$ source. The mobility-lifetime products of the electron-hole pair in a planar CZT detector ($5{\times}5{\times}2\;mm^3$) were determined by fitting the peak position as a function of biased voltage data to the Hecht equation. To verify the accuracy of these products derived from ${\alpha}$ particles and low-energy gamma-rays, an energy spectrum considering the transport property of the CZT detector was simulated through a combination of the deposited energy and the charge collection efficiency at a specific position. It was found that the shaping time of the amplifier module significantly affects the determination of the (${\mu}{\tau}$) products; the ${\alpha}$ particle method was stabilized with an increase in the shaping time and was less sensitive to this change compared to when the gamma-ray method was used. In the case of the simulated energy spectrum with transport properties evaluated by the ${\alpha}$ particle method, the peak position and tail were slightly different from the measured result, whereas the energy spectrum derived from the low-energy gamma-ray was in good agreement with the experimental results. From these results, it was confirmed that low-energy gamma-rays are more useful when seeking to obtain the transport properties of carriers than ${\alpha}$ particles because the methods that use gamma-rays are less influenced by the surface condition of the CZT detector. Furthermore, the analysis system employed in this study, which was configured by a combination of Monte Carlo simulation and the Hecht model, is expected to be highly applicable to the study of the characteristics of CZT detectors.

• Journal of Radiation Protection and Research
• /
• 제45권2호
• /
• pp.69-75
• /
• 2020

Background: Dose conversion coefficients (DCCs) have been commonly used to estimate radiation-dose absorption by human organs based on physical measurements of fluence or kerma. The International Commission on Radiological Protection (ICRP) has reported a library of DCCs, but few studies have been conducted on their applicability to non-Caucasian populations. In the present study, we collected a total of 8 Korean pediatric and adult voxel phantoms to calculate the organ DCCs for idealized external photon-irradiation geometries. Materials and Methods: We adopted one pediatric female phantom (ETRI Child), two adult female phantoms (KORWOMAN and HDRK Female), and five adult male phantoms (KORMAN, ETRI Man, KTMAN1, KTMAN2, and HDRK Man). A general-purpose Monte Carlo radiation transport code, MCNPX2.7 (Monte Carlo N-Particle Transport extended version 2.7), was employed to calculate the DCCs for 13 major radiosensitive organs in six irradiation geometries (anteroposterior, posteroanterior, right lateral, left lateral, rotational, and isotropic) and 33 photon energy bins (0.01-20 MeV). Results and Discussion: The DCCs for major radiosensitive organs (e.g., lungs and colon) in anteroposterior geometry agreed reasonably well across the 8 Korean phantoms, whereas those for deep-seated organs (e.g., gonads) varied significantly. The DCCs of the child phantom were greater than those of the adult phantoms. A comparison with the ICRP Publication 116 data showed reasonable agreements with the Korean phantom-based data. The variations in organ DCCs were well explained using the distribution of organ depths from the phantom surface. Conclusion: A library of dose conversion coefficients for major radiosensitive organs in a series of pediatric and adult Korean voxel phantoms was established and compared with the reference data from the ICRP. This comparison showed that our Korean phantom-based data agrees reasonably with the ICRP reference data.

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

• 지구물리와물리탐사
• /
• 제15권4호
• /
• pp.219-226
• /
• 2012

이 연구에서는 국내에서 연구가 미약했었던 중성자검층 수치모델링을 이용하여 다양한 시추공 환경에서의 검출기 반응을 분석하였다. 이를 위해 중성자검층 환경을 MCNP 알고리듬으로 구현하여 시뮬레이션을 수행하였다. MCNP 알고리듬은 방사선 수송 시뮬레이션이 및 3차원 기하구조 표현이 가능하여 다양한 분야에서 전세계적으로 많이 이용되고 있다. 먼저 시뮬레이션 결과를 검증하기 위해, 기존 연구의 검출기반응 결과 그래프를 이용하여 비교 분석하였다. 중성자 검층 반응 분석이 가능한 중성자 검층기의 일반적인 특징에 기초하여 수학적으로 중성자검층기 모형을 구성하여 반응을 계산하였다. 먼저, 석회암, 사암, 돌로마이트 등과 같은 매질에서 공극률을 다양하게 변화시켜 가며 수치 계산함으로써, 이 연구에서 고려하고 있는 중성자검층기의 교정곡선(calibration chart)을 도출하였다. 이에 기초하여, 실제 중성자검층 시 고려해야 할 공내수 유무에 의한 반응 변화, 염수가 중성자검층에 미치는 영향 등을 분석함으로써 시추공 환경 변화에 따라 보다 정확하게 공극률을 해석할 수 있는 기틀을 마련하고자 한다.

• Journal of Radiation Protection and Research
• /
• 제41권4호
• /
• pp.344-349
• /
• 2016

Background: Carbon ion therapy has achieved satisfactory results. However, patients have a risk to get a secondary cancer. In order to estimate the risk, it is essential to understand particle transportation and nuclear reactions in the patient's body. The particle transport Monte Carlo simulation code is a useful tool to understand them. Since the code validation for heavy ion incident reactions is not enough, the experimental data of the elementary reaction processes are needed. Materials and Methods: We measured neutron production double-differential cross-sections (DDXs) on a carbon bombarded with 430 MeV/nucleon carbon beam at PH2 beam line of HIMAC facility in NIRS. Neutrons produced in the target were measured with NE213 liquid organic scintillators located at six angles of 15, 30, 45, 60, 75, and $90^{\circ}$. Results and Discussion: Neutron production double-differential cross-sections for carbon bombarded with 430 MeV/nucleon carbon ions were measured by the time-of-flight method with NE213 liquid organic scintillators at six angles of 15, 30, 45, 60, 75, and $90^{\circ}$. The cross sections were obtained from 1 MeV to several hundred MeV. The experimental data were compared with calculated results obtained by Monte Carlo simulation codes PHITS, Geant4, and FLUKA. Conclusion: PHITS was able to reproduce neutron production for elementary processes of carbon-carbon reaction precisely the best of three codes.

• 대한전자공학회논문지
• /
• 제27권4호
• /
• pp.534-540
• /
• 1990

In this paper the characteristics of submicron gate GaAs MESFET's have been studied using a particle model which takes into account the hot-electron transport phenomena, i.e., the velocity overshoot. \ulcornervalley(<000> direction), L valley (<111>direction), X valley (<100>direction) as the GaAs conduction energy band and optical phonon, acoustic phonon, equivalent intervalley, nonequivalent intervalley scattering as the scattering models, have been considered in this simulation. And the GaAs material and the device simulation have been done by determination of the free flight time, scattering mechanism and scattering angle according to Monte-Carlo algorithm which makes use of a particle model. As a result of the particle simulation, firstly the electron distribution, the potential energy distribution and the situation of electron displacement in 0.6 \ulcorner gate length device have been obtained. Secondly, the cutoff frequency, obtained by this method, is k47GHz which is in good agreement with the calculated result of theory. And the current-voltage characteristics curve which takes account of the buffer layer effect has been obtained. Lastly it has been verified that parasitic current at the buffer layer can be analyzed using channel depth modulation.