• Journal of the Korean Physical Society
• /
• v.73 no.11
• /
• pp.1764-1773
• /
• 2018

A-PET is a quad-head PET scanner developed for use in small-animal imaging. The dimensions of its volumetric field of view (FOV) are $46.1{\times}46.1{\times}46.1mm^3$ and the gap between the detector modules has been minimized in order to provide a highly sensitive system. However, such a small FOV together with the quad-head geometry causes image quality degradation. The main factor related to image degradation for the quad-head PET is the mispositioning of events caused by the penetration effect in the detector. In this paper, we propose a precise method for modelling the system at the high spatial resolution of the A-PET using a LOR (line of response) based ML-EM (maximum likelihood expectation maximization) that allows for penetration effects. The proposed system model provides the detection probability of every possible ray-path via crystal sampling methods. For the ray-path sampling, the sub-LORs are defined by connecting the sampling points of the crystal pair. We incorporate the detection probability of each sub-LOR into the model by calculating the penetration effect. For comparison, we used a standard LOR-based model and a Monte Carlo-based modeling approach, and evaluated the reconstructed images using both the National Electrical Manufacturers Association NU 4-2008 standards and the Geant4 Application for Tomographic Emission simulation toolkit (GATE). An average full width at half maximum (FWHM) at different locations of 1.77 mm and 1.79 mm are obtained using the proposed system model and standard LOR system model, which does not include penetration effects, respectively. The standard deviation of the uniform region in the NEMA image quality phantom is 2.14% for the proposed method and 14.3% for the LOR system model, indicating that the proposed model out-performs the standard LOR-based model.

• Journal of Radiation Protection and Research
• /
• v.46 no.2
• /
• pp.58-65
• /
• 2021

Background: The evaluation of skyshine distribution, release of airborne radioactive nuclides, and soil activation and groundwater migration were required for radiological assessment of the impact on the environment surrounding In-Flight (IF)-system facility of the RAON (Rare isotope Accelerator complex for ON-line experiment) accelerator complex. Materials and Methods: Monte Carlo simulation by MCNPX code was used for evaluation of skyshine and activation analysis for air and soil. The concentration model was applied in the estimation of the groundwater migration of radionuclides in soil. Results and Discussion: The skyshine dose rates at 1 km from the facility were evaluated as 1.62 × 10^-3 μSv·hr^-1. The annual releases of ³H and ¹⁴C were calculated as 9.62 × 10^-5 mg and 1.19 × 10^-1 mg, respectively. The concentrations of ³H and ²²Na in drinking water were estimated as 1.22 × 10^-1 Bq·cm^-3 and 8.25 × 10^-3 Bq·cm^-3, respectively. Conclusion: Radiological assessment of environmental impact on the IF-facility of RAON was performed through evaluation of skyshine dose distribution, evaluation of annual emission of long-lived radionuclides in the air and estimation of soil activation and groundwater migration of radionuclides. As a result, much lower exposure than the limit value for the public, 1 mSv·yr^-1, is expected during operation of the IF-facility.

• Journal of the Korean Society of Radiology
• /
• v.17 no.1
• /
• pp.141-148
• /
• 2023

The algebraic reconstruction technique is one of the reconstruction methods in CT and shows good image quality against noise-dominant conditions. The number of iteration is one of the key factors determining the execution time for the algebraic reconstruction technique. However, there are some rules for determining the number of iterations that result in more than a few hundred iterations. Thus, the rules are difficult to apply in practice. In this study, we proposed a method to determine the number of iterations for practical applications. The reconstructed image quality shows slow convergence as the number of iterations increases. Image quality 𝜖 < 0.001 was used to determine the optimal number of iteration. The Shepp-Logan head phantom was used to obtain noise-free projection and projections with noise for 360, 720, and 1440 views were obtained using Geant4 Monte Carlo simulation that has the same geometry dimension as a clinic CT system. Images reconstructed by around 10 iterations within the stop condition showed good quality. The method for determining the iteration number is an efficient way of replacing the best image-quality-based method, which brings over a few hundred iterations.

• Smart Structures and Systems
• /
• v.31 no.2
• /
• pp.101-111
• /
• 2023

Structural health monitoring (SHM) covers a range of damage detection strategies for buildings. In real-time, SHM provides a basis for rapid decision making to optimise the speed and economic efficiency of post-event response. Previous work introduced an SHM method based on identifying structural nonlinear hysteretic parameters and their evolution from structural force-deformation hysteresis loops in real-time. This research extends and generalises this method to investigate the impact of a wide range of flag-shaped or pinching shape nonlinear hysteretic response and its impact on the SHM accuracy. A particular focus is plastic stiffness (K_p), where accurate identification of this parameter enables accurate identification of net and total plastic deformation and plastic energy dissipated, all of which are directly related to damage and infrequently assessed in SHM. A sensitivity study using a realistic seismic case study with known ground truth values investigates the impact of hysteresis loop shape, as well as added noise, on SHM accuracy using a suite of 20 ground motions from the PEER database. Monte Carlo analysis over 22,000 simulations with different hysteresis loops and added noise resulted in absolute percentage identification error (median, (IQR)) in K_p of 1.88% (0.79, 4.94)%. Errors were larger where five events (Earthquakes #1, 6, 9, 14) have very large errors over 100% for resulted K_p as an almost entirely linear response yielded only negligible plastic response, increasing identification error. The sensitivity analysis shows accuracy is reduces to within 3% when plastic drift is induced. This method shows clear potential to provide accurate, real-time metrics of non-linear stiffness and deformation to assist rapid damage assessment and decision making, utilising algorithms significantly simpler than previous non-linear structural model-based parameter identification SHM methods.

• Journal of Radiation Protection and Research
• /
• v.48 no.2
• /
• pp.90-99
• /
• 2023

Background: Quantitative risk assessments should be accompanied by uncertainty analyses of the risk models employed in the calculations. In this study, we aim to develop a computational code named SUMRAY for use in cancer risk projections from radiation exposure taking into account uncertainties. We also aim to make SUMRAY publicly available as a resource for further improvement of risk projection. Materials and Methods: SUMRAY has two versions of code written in R and Python. The risk models used in SUMRAY for all-solid-cancer mortality and incidence were those published in the Life Span Study of a cohort of the atomic bomb survivors in Hiroshima and Nagasaki. The confidence intervals associated with the evaluated risks were derived by propagating the statistical uncertainties in the risk model parameter estimates by the Monte Carlo method. Results and Discussion: SUMRAY was used to calculate the lifetime or time-integrated attributable risks of cancer under an exposure scenario (baseline rates, dose[s], age[s] at exposure, age at the end of follow-up, sex) specified by the user. The results were compared with those calculated using another well-known web-based tool, Radiation Risk Assessment Tool (RadRAT; National Institutes of Health), and showed a reasonable agreement within the estimated confidential interval. Compared with RadRAT, SUMRAY can be used for a wide range of applications, as it allows the risk projection with arbitrarily specified risk models and/or population reference data. Conclusion: The reliabilities of SUMRAY with the present risk-model parameters and their variance-covariance matrices were verified by comparing them with those of the other codes. The SUMRAY code is distributed to the public as an open-source code under the Massachusetts Institute of Technology license.

• Progress in Medical Physics
• /
• v.33 no.4
• /
• pp.80-87
• /
• 2022

Purpose: Proton therapy has different relative biological effectiveness (RBE) compared with X-ray treatment, which is the standard in radiation therapy, and the fixed RBE value of 1.1 is widely used. However, RBE depends on a charged particle's linear energy transfer (LET); therefore, measuring LET is important. We have developed a LET measurement method using the inefficiency characteristic of an EBT3 film on a proton beam's Bragg peak (BP) region. Methods: A Gafchromic EBT3 film was used to measure the proton beam LET. It measured the dose at a 10-cm pristine BP proton beam in water to determine the quenching factor of the EBT3 film as a reference beam condition. Monte Carlo (MC) calculations of dose-averaged LET (LET_d) were used to determine the quenching factor and validation. The dose-averaged LETs at the 12-, 16-, and 20-cm pristine BP proton beam in water were calculated with the quenching factor. Results: Using the passive scattering proton beam nozzle of the National Cancer Center in Korea, the LET_d was measured for each beam range. The quenching factor was determined to be 26.15 with 0.3% uncertainty under the reference beam condition. The dose-averaged LETs were measured for each test beam condition. Conclusions: We developed a method for measuring the proton beam LET using an EBT3 film. This study showed that the magnitude of the quenching effect can be estimated using only one beam range, and the quenching factor determined under the reference condition can be applied to any therapeutic proton beam range.

• Journal of Radiation Protection and Research
• /
• v.48 no.4
• /
• pp.184-196
• /
• 2023

Background: Radiation protection is crucial in various fields due to the harmful effects of radiation. Shielding is used to reduce radiation exposure, but gamma radiation poses challenges due to its high energy and penetration capabilities. Materials and Methods: This work investigates the radiation shielding properties of polyvinylidene fluoride (PVDF) samples containing different weight fraction of tungsten carbide (WC), tungsten trioxide (WO₃), and tungsten disulfide (WS₂). Parameters such as the mass attenuation coefficient (MAC), half-value layer (HVL), mean free path (MFP), effective atomic number (Z_eff), and macroscopic effective removal cross-section for fast neutrons (Σ_R) were calculated using the Phy-X/PSD software. EpiXS simulations were conducted for MAC validation. Results and Discussion: Increasing the weight fraction of the additives resulted in higher MAC values, indicating improved radiation shielding. PVDF-xWC showed the highest percentage increase in MAC values. MFP results indicated that PVDF-0.20WC has the lowest values, suggesting superior shielding properties compared to PVDF-0.20WO₃ and PVDF-0.20WS₂. PVDF-0.20WC also exhibited the highest Z_eff values, while PVDF-0.20WS₂ showed a slightly higher increase in Z_eff at energies of 0.662 and 1.333 MeV. PVDF-0.20WC has demonstrated the highest Σ_R value, indicating effective shielding against fast neutrons, while PVDF-0.20WS₂ had the lowest Σ_R value. The Monte Carlo N-Particle Transport version 5 (MCNP5) simulations showed that PVDF-xWC attenuates gamma radiation more than pure PVDF, significantly decreasing the dose equivalent rate. Conclusion: Overall, this research provides insights into the radiation shielding properties of PVDF mixtures, with PVDF-xWC showing the most promising results.

• Journal of Radiation Protection and Research
• /
• v.48 no.3
• /
• pp.117-123
• /
• 2023

Background: We investigated the impact of 0.35 T magnetic field on dose calculation for non-small cell lung cancer (NSCLC) stereotactic ablative radiotherapy (SABR) in the ViewRay system (ViewRay Inc.), which features a simultaneous use of magnetic resonance imaging (MRI) to guide radiotherapy for an improved targeting of tumors. Materials and Methods: Here, we present a comprehensive analysis of the effects induced by the 0.35 T magnetic field on various characteristics of SABR plans including the plan qualities and dose calculation for the planning target volume, organs at risk, and outer/inner shells. Therefore, two SABR plans were set up, one with a 0.35 T magnetic field applied during radiotherapy and another in the absence of the field. The dosimetric parameters were calculated in both cases, and the plan quality indices were evaluated using a Monte Carlo algorithm based on a treatment planning system. Results and Discussion: Our findings showed no significant impact on dose calculation under the 0.35 T magnetic field for all analyzed parameters. Nonetheless, a significant enhancement in the dose was calculated on the skin surrounding the tumor when the 0.35 T magnetic field was applied during the radiotherapy. This was attributed to the electron return effect, which results from the deviation of the electrons ejected from tissues upon radiation due to Lorentz forces. These returned electrons re-enter the tissues, causing a local dose increase in the calculated dose. Conclusion: The present study highlights the impact of the 0.35 T magnetic field used for MRI in the ViewRay system for NSCLC SABR treatment, especially on the skin surrounding the tumors.

• Journal of Radiation Protection and Research
• /
• v.48 no.3
• /
• pp.131-143
• /
• 2023

Background: This study examined the detection limit of thyroid screening monitoring conducted at the time of the Fukushima Daiichi Nuclear Power Plant (FDNPP) accident in 2011 using a Monte Carlo simulation. Materials and Methods: We calculated the detection limit of a NaI(Tl) survey meter to measure ¹³¹I accumulation in the thyroid gland of children. Mathematical phantoms of 1- and 5-year-old children were developed in the simulation of the Particle and Heavy Ion Transport code System code. Contamination of the body surface with eight radionuclides found after the FDNPP accident was assumed to have been deposited on the neck and shoulder area. Results and Discussion: The detection limit was calculated as a function of ambient dose rate. In the case of 40 Bq/cm² contamination on the body surface of the neck, the present simulations showed that residual thyroid radioactivity corresponding to thyroid dose of 100 mSv can be detected within 21 days after intake at the ambient dose rate of 0.2 µSv/hr and within 11 days in the case of 2.0 µSv/hr. When a time constant of 10 seconds was used at the dose rate of 0.2 µSv/hr, the estimated survey meter output error was 5%. Evaluation of the effect of individual differences in the location of the thyroid gland confirmed that the measured value would decrease by approximately 6% for a height difference of ±1 cm and increase by approximately 65% for a depth of 1 cm. Conclusion: In the event of a nuclear disaster, simple measurements carried out using a NaI(Tl) scintillation survey meter remain effective for assessing ¹³¹I intake. However, it should be noted that the presence of short-half-life radioactive materials on the body surface affects the detection limit.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.16 no.12
• /
• pp.8128-8139
• /
• 2015

Recently, the necessity of developing the load and resistance factor design(LRFD) for soft ground improvement method has been raised, since the limit state design is requested as international technical standard for the foundation of structures. In this study, to develop LRFD codes for foundation structures in Korea, target reliability index and resistance factor for static bearing capacity of driven steel pipe piles were calibrated in the framework of reliability theory. The 16 data(in Gwangyang) and the 57 data(Korea Institute of Construction Technology, 2008) sets of static load test and soil property tests conducted in the whole domestic area were collected along with available subsurface investigation results. The resistance bias factors were evaluated for the tow static design methods by comparing the representative measured bearing capacities with the expected design values. Reliability analysis was performed by two types of advanced methods : the First Order Reliability Method (FORM), and the Monte Carlo Simulation (MCS) method using resistance bias factor statistics. As a result, when target reliability indices of the driven pipe pile were selected as 2.0, 2.33, 2.5, resistance factor of two design methods for SPT N at pile tip less than 50 were evaluated as 0.611~0.684, 0.537~0.821 respectively, and STP N at pile tip more than 50 were evaluated as 0.545~0.608, 0.643~0.749 respectively. The result from this research will be useful for developing various foundations and soil structures under LRFD.