• Nuclear Engineering and Technology
• /
• v.19 no.2
• /
• pp.85-98
• /
• 1987

To analyze the aerosol dynamics in severe accidents of LMFBR, a new computer code entitled MCAD (Multicomponent Aerosol Dynamics) has been developed. The code can treat two component aerosol system using relative collision probability of each particles as sequences of accident scenarios. Coagulation and removal mechanisms incorporating Brownian diffusion and gravitational sedimentation are included in this model. In order to see the effect of particle geometry, the code makes use of the concept of density correction factor and shape factors. The code is verified using the experimental result of NSPP-300 series and compared to other code. At present, it fits the result of experiment well and agrees to the existing code. The input variables included are very uncertain. Hence, it requires uncertainty and sensitivity analysis as a supplement to code development. In this analysis, 14 variables are selected to analyze. The input variables are compounded by experimental design method and Latin hypercube sampling. Their results are applied to Response surface method to see the degree of regression. The stepwise regression method gives an insight to which variables are significant as time elapse and their reasonable ranges. Using Monte Carlo Method to the regression model of LHS, the confidence level of the results of MCAD and their variables is improved.

• Journal of the Korean Society of Radiology
• /
• v.15 no.6
• /
• pp.849-853
• /
• 2021

PET images used in medical diagnoses are created using positron emitting radionuclides. The radiation used for imaging is generated at 0.511 MeV by p-annihilation. The CSDA range is the distance the particle radiation flew physically, and it is different from the range shown in PET images. This study proposes a novel method that uses radiological criteria to measure this range. The experiment was conducted by applying the MCNP6 simulation to positron emitting nuclides ¹⁸F, ¹¹C, ¹³N, and ¹⁵O. Radiological criteria were based on the location of the p-annihilation event, which is also the image signal. Results showed the radiological range of positrons to be 2.3, 3.9, 5.0, and 7.9 mm for ¹⁸F, ¹¹C, ¹³N, and ¹⁵O, respectively. The higher the positron energy, the larger its difference from the CSDA range. Positron emitting nuclide is being developed and studied as a nuclide for dosimetry or radiotherapy. Further research needs to be conducted into various positron ranges.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.217-217
• /
• 2015

The objective of this study was mainly to evaluate the water resources potential of Lake Tana Basin (LTB) by using Soil and Water Assessment Tool (SWAT). From SWAT simulation of LTB, about 5236 km2 area of LTB is gauged watershed and the remaining 9878 km2 area is ungauged watershed. For calibration of model parameters, four gauged stations were considered namely: Gilgel Abay, Gummera, Rib, and Megech. The SWAT-CUP built-in techniques, particle swarm optimization (PSO) and generalized likelihood uncertainty estimation (GLUE) method was used for calibration of model parameters and PSO method were selected for the study based on its performance results in four gauging stations. However the level of sensitivity of flow parameters differ from catchment to catchment, the curve number (CN2) has been found the most sensitive parameters in all gauged catchments. To facilitate the transfer of data from gauged catchments to ungauged catchments, clustering of hydrologic response units (HRUs) were done based on physical similarity measured between gauged and ungauged catchment attributes. From SWAT land use/ soil use/slope reclassification of LTB, a total of 142 HRUs were identified and these HRUs are clustered in to 39 similar hydrologic groups. In order to transfer the optimized model parameters from gauged to ungauged catchments based on these clustered hydrologic groups, this study evaluates three parameter transfer schemes: parameters transfer based on homogeneous regions (PT-I), parameter transfer based on global averaging (PT-II), and parameter transfer by considering Gilgel Abay catchment as a representative catchment (PT-III) since its model performance values are better than the other three gauged catchments. The performance of these parameter transfer approach was evaluated based on values of Nash-Sutcliffe efficiency (NSE) and coefficient of determination (R2). The computed NSE values was found to be 0.71, 0.58, and 0.31 for PT-I, PT-II and PT-III respectively and the computed R2 values was found to be 0.93, 0.82, and 0.95 for PT-I, PT-II, and PT-III respectively. Based on the performance evaluation criteria, PT-I were selected for modelling ungauged catchments by transferring optimized model parameters from gauged catchment. From the model result, yearly average stream flow for all homogeneous regions was found 29.54 m3/s, 112.92 m3/s, and 130.10 m3/s for time period (1989 - 2005) for region-I, region-II, and region-III respectively.

• Nuclear Engineering and Technology
• /
• v.54 no.1
• /
• pp.162-176
• /
• 2022

In Lead-based reactor (LBR) severe accident, the meltdown and migration inside the reactor core will lead to fuel fragment concentration, which may further cause re-criticality and even core disintegration. Accurately predicting the migration and solidification behavior of melt in LBR severe accidents is of prime importance for safety analysis of LBR. In this study, the Moving Particle Semi-implicit (MPS) method is validated and used to simulate the migration and solidification behavior. Two main surface tension models are validated and compared. Meanwhile, the MPS method is validated by the L-plate solidification test. Based on the improved MPS method, the migration and solidification behavior of melt in LBR severe accident was studied furthermore. In the Pb-Bi coolant, the melt flows upward due to density difference. The migration and solidification behavior are greatly affected by the surface tension and viscous resistance varying with enthalpy. The whole movement process can be divided into three stages depending on the change in velocity. The heat transfer of core melt is determined jointly by two heat transfer modes: flow heat transfer and solid conductivity. Generally, the research results indicate that the MPS method has unique advantage in studying the migration and solidification behavior in LBR severe accident.

• 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 Radiation Protection and Research
• /
• v.48 no.4
• /
• pp.197-203
• /
• 2023

Background: The thermoluminescent dosimeter (TLD) and Monte Carlo (MC) dosimetry are carried out to determine the occupational dose for personnel in the handling of ¹²⁵I seed sources. Materials and Methods: TLDs were placed in different layers of the Alderson-Rando phantom in the thyroid, lung and also eyes and skin surface. An ¹²⁵I seed source was prepared and its activity was measured using a dose calibrator and was placed at two distances of 20 and 50 cm from the Alderson-Rando phantom. In addition, the Monte Carlo N-Particle Extended (MCNPX 2.6.0) code and a computational phantom with a lattice-based geometry were used for organ dose calculations. Results and Discussion: The comparison of TLD and MC results in the thyroid and lung is consistent. Although the relative difference of MC dosimetry to TLD for the eyes was between 4% and 13% and for the skin between 19% and 23%, because of the existence of a higher uncertainty regarding TLD positioning in the eye and skin, these inaccuracies can also be acceptable. The isodose distribution was calculated in the cross-section of the head phantom when the ¹²⁵I seed was at two distances of 20 and 50 cm and it showed that the greatest dose reduction was observed for the eyes, skin, thyroid, and lungs, respectively. The results of MC dosimetry indicated that for near the head positions (distance of 20 cm) the absorbed dose rates for the eye lens, eye and skin were 78.1±2.3, 59.0±1.8, and 10.7±0.7 µGy/mCi/hr, respectively. Furthermore, we found that a 30 cm displacement for the ¹²⁵I seed reduced the eye and skin doses by at least 3- and 2-fold, respectively. Conclusion: Using a computational phantom to monitor the dose to the sensitive organs (eye and skin) for personnel involved in the handling of ¹²⁵I seed sources can be an accurate and inexpensive method.

• Journal of Radiation Protection and Research
• /
• v.41 no.1
• /
• pp.71-79
• /
• 2016

Background: It is necessary to analyze radioactivity of naturally occurring radioactive materials (NORM) in products to ensure radiological safety required by Natural Radiation Safety Management Act. The pretreatments for the existing analysis methods require high technology and time. Such destructive pretreatments including grinding and dissolution of samples make impossible to reuse products. We developed a rapid and simple procedure of radioactivity analysis for thorium series in the products containing NORM. Materials and Methods: The developed method requires non-destructive or minimized pretreatment. Radioactivity of the product without pretreatment is initially measured using gamma spectroscopy and then the measured radioactivity is adjusted by considering material composition, mass density, and geometrical shape of the product. The radioactivity adjustment can be made using scaling factors, which is derived by radiation transport Monte Carlo simulation. Necklace, bracelet, male health care product, and tile for health mat were selected as representative products for this study. The products are commonly used by the public and directly contacted with human body and thus resulting in high radiation exposure to the user. Results and Discussion: The scaling factors were derived using MCNPX code and the values ranged from 0.31 to 0.47. If radioactivity of the products is measured without pretreatment, the thorium series may be overestimated by up to 2.8 times. If scaling factors are applied, the difference in radioactivity estimates are reduced to 3-24%. Conclusion : The developed procedure in this study can be used for other products with various materials and shapes and thus ensuring radiological safety.

• Progress in Medical Physics
• /
• v.27 no.1
• /
• pp.31-36
• /
• 2016

The paper discusses radiation dose of dual energy CT on which copper modulation layer, is mounted in order to improve diagnostic performance of the dual energy CT. The radiation dose is estimated using MCNPX and its results are compared with that of the conventional dual energy CT system. CT X-ray spectra of 80 and 120 kVp, which are usually used for thorax, abdominal, head, and neck CT scans, were generated by the SPEC78 code and were used for the source specification 'SDEF' card for MCNPX dose modeling. The copper modulation layer was located 20 cm away from a source covering half of the X-ray window. The radiation dose was measured as changing its thickness from 0.5 to 2.0 mm at intervals of 0.5 mm. Since the MCNPX tally provides only normalized values to a single particle, the dose conversion coefficients of F6 tally for the modulation layer-based dual energy CBCT should be calculated for matching the modeling results into the actual dose. The dose conversion coefficient is $7.2*10^4cGy/output$ that is obtained from dose calibration curve between F6 tally and experimental results in which GAFCHORMIC EBT3 films were exposed by an already known source. Consequently, the dose of the modulation layer-based dual energy cone beam CT is 33~40% less than that of the single energy CT system. On the basis of the results, it is considered that scattered dose produced by the copper modulation layer is very small. It shows that the modulation layer-based dual energy CBCT system can effectively reduce radiation dose, which is the major disadvantage of established dual energy CT.

• Journal of Radiation Protection and Research
• /
• v.39 no.4
• /
• pp.199-205
• /
• 2014

In the field of neutron detections, $^3He$ gas, the so-called "the gold standard," is the most widely used material for neutron detections because of its high efficiency in neutron capturing. However, from variable causes since early 2009, $^3He$ is being depleted, which has maintained an upward pressure on its cost. For this reason, the demands for $^3He$ replacements are rising sharply. Research into neutron converting materials, which has not been used well due to a neutron detection efficiency lower than the efficiency of $^3He$, although it can be chosen for use in a neutron detector, has been highlighted again. $^{10}B$, which is one of the $^3He$ replacements, such as $BF_3$, $^6Li$, $^{10}B$, $Gd_2O_2S$, is being researched by various detector development groups owing to a number of advantages such as easy gamma-ray discrimination, non-toxicity, low cost, etc. One of the possible techniques for the detection is an indirect neutron detection method measuring secondary radiation generated by interactions between neutrons and $^{10}B$. Because of the mean free path of alpha particle from interactions that are very short in a solid material, the thickness of $^{10}B$ should be thin. Therefore, to increase the neutron detection efficiency, it is important to make a $^{10}B$ thin film. In this study, we fabricated a $^{10}B$ thin film that is about 60 um in thickness for neutron detection using well-known technology for the manufacturing of a thin electrode for use in lithium ion batteries. In addition, by performing simple physical tests on the conductivity, dispersion, adhesion, and flexibility, we confirmed that the physical characteristics of the fabricated $^{10}B$ thin film are good. Using the fabricated $^{10}B$ thin film, we made a proportional counter for neutron monitoring and measured the neutron pulse height spectrum at a neutron facility at KAERI. Furthermore, we calculated using the Monte Carlo simulation the change of neutron detection efficiency according to the number of thin film layers. In conclusion, we suggest a fabrication method of a $^{10}B$ thin film using the technology used in making a thin electrode of lithium ion batteries and made the $^{10}B$ thin film for neutron detection using suggested method.