• Progress in Medical Physics
• /
• v.16 no.3
• /
• pp.138-147
• /
• 2005

The Monte Carlo method cannot have been used for routine treatment planning because of heavy time consumption for the acceptable accuracy. Since calculation time is proportional to particle histories, we can save time by decreasing the number of histories. However, a small number of histories can cause serious uncertainties. In this study, we proposed Monte Carlo dose computation time and uncertainty reduction method using specially designed filters and adaptive denoising process. Proposed algorithm was applied to 6 MV photon and 21 MeV electron dose calculations in homogeneous and heterogeneous phantoms. Filtering time was negligible comparing to Monte Carlo simulation time. The accuracy was improved dramatically in all situations and the simulation of 1 $\%$ to 10$\%$ number of histories of benchmark in photon and electron dose calculation showed the most beneficial result. The empirical reduction of necessary histories was about a factor of ten to fifty from the result.

• Journal of Radiation Protection and Research
• /
• v.39 no.1
• /
• pp.30-37
• /
• 2014

Recently High-Definition Reference Korean-Man (HDRK-Man) and High-Definition Reference Korean-Woman (HDRK-Woman) were constructed in Korea. The HDRK phantoms were designed to represent respectively reference Korean male and female to calculate effective doses for Korean by performing Monte Carlo dose calculation. However, the Monte Carlo dose calculation requires detailed knowledge on computational human phantoms and Monte Carlo simulation technique which regular researchers in radiation protection dosimetry and practicing health physicists do not have. Recently the UFPE (Federal University of Pernambuco) research group has developed, and opened to public, an online Monte Carlo dose calculation system called CALDOSE_X(www.caldose.org). By using the CALDOSE_X, one can easily perform Monte Carlo dose calculations. However, the CALDOSE_X used caucasian phantoms to calculate organ doses or effective doses which are limited for Korean. The present study developed an online reference Korean dose calculation system which can be used to calculate effective doses for Korean.

• Progress in Medical Physics
• /
• v.14 no.4
• /
• pp.240-248
• /
• 2003

The Monte Carlo calculation is the most accurate means of predicting radiation dose, but its accuracy is accompanied by an increase in the amount of time required to produce a statistically meaningful dose distribution. In this study, the effects on calculation time by introducing variance reduction techniques and increasing computing power, respectively, in the Monte Carlo dose calculation for a 6 MV photon beam from the Varian 600 C/D were estimated when maintaining accuracy of the Monte Carlo calculation results. The EGSnrc­based BEAMnrc code was used to simulate the beam and the EGSnrc­based DOSXYZnrc code to calculate dose distributions. Variance reduction techniques in the codes were used to describe reduced­physics, and a computer cluster consisting of ten PCs was built to execute parallel computing. As a result, time was more reduced by the use of variance reduction techniques than that by the increase of computing power. Because the use of the Monte Carlo dose calculation in clinical practice is yet limited by reducing the computational time only through improvements in computing power, introduction of reduced­physics into the Monte Carlo calculation is inevitable at this point. Therefore, a more active investigation of existing or new reduced­physics approaches is required.

• Communications for Statistical Applications and Methods
• /
• v.3 no.2
• /
• pp.291-298
• /
• 1996

통계학분야 그리고 또 다른 많은 분야에서 수치적 계산을 다루는 문제가 자주 발생한다. 적당한 컴퓨터 컴퓨터 시간안에 상당한 정도의 정확성을 줄 수 있고 또한 보다 광범위하게 사용 가능한 유용한 알고리즘의 필요성을 느낀다. 이러한 문제에 가능한 하나의 몬테칼로 알고리즘인 주표본 알고리즘을 소개하였다. 그리고 특히 본 눈문에서는 축차확률비검정의 오차확률을 계산하는 곳에 주표본 알고리즘을 적용하고 결과를 비교분석하였다.

• Progress in Medical Physics
• /
• v.10 no.1
• /
• pp.55-62
• /
• 1999

In an effort to study the characteristics of x-rays utilized in radiation therapy, we calculated the energy distribution and the mean energy of x-rays generated from a tungsten target bombarded by 6, 10, and 15 MeV electron beams, using a Monte Carlo technique. The average photon energies calculated as a function of the beam radius lied in 1.4 ∼ 1.6, 2.1 ∼ 2.5 and 2.8 ∼ 3.3 MeV ranges for 4, 10, and 15 MV electron beams, respectively, which turned out to have no strong dependence on the radius. Using the energy distributions of 6,10, and 15 MV x-rays obtained for the target distance of 100 cm, percentage depth doses were determined using Monte Carlo calculations. For the case 10 MV, a comparison was made between our calculation and measurement performed by others. The calculated percentage depth dose appeared somewhat smaller than the measured one except in the surface region. We conclude that this is due to the fact that the beam hardening effect resulting from the flattening filter was not properly allowed for in our Monte Carlo calculations.

• Journal of radiological science and technology
• /
• v.31 no.1
• /
• pp.89-95
• /
• 2008

The quality correction factor for used beam and qualities is strongly required for clinical dosimetry by TRS-398 protocol of IAEA. In this study the quality correction factors for a commercial plane-parallel ionization chamber in high energy electron beams were calculated by Monte Carlo code(DOSRZnrc/EGSnrc). In comparison of quality correction factor, the difference between this study and TRS-398 were within 1% in 5-20 MeV. In case of 4MeV the difference was 1.9%. As an independent method of determination of quality correction factor this study can be applied to evaluate values in the protocol or calculate the factor for a new chamber.

• Progress in Medical Physics
• /
• v.21 no.3
• /
• pp.253-260
• /
• 2010

Most brachytherapy treatment planning systems employ a dosimetry formalism based on the AAPM TG-43 report which does not appropriately consider tissue heterogeneity. In this study we aimed to set up a simple Monte Carlo-based intracavitary high-dose-rate brachytherapy (IC-HDRB) plan verification platform, focusing particularly on the robustness of the direct Monte Carlo dose calculation using material and density information derived from CT images. CT images of slab phantoms and a uterine cervical cancer patient were used for brachytherapy plans based on the Plato (Nucletron, Netherlands) brachytherapy planning system. Monte Carlo simulations were implemented using the parameters from the Plato system and compared with the EBT film dosimetry and conventional dose computations. EGSnrc based DOSXYZnrc code was used for Monte Carlo simulations. Each $^{192}Ir$ source of the afterloader was approximately modeled as a parallel-piped shape inside the converted CT data set whose voxel size was $2{\times}2{\times}2\;mm^3$. Bracytherapy dose calculations based on the TG-43 showed good agreement with the Monte Carlo results in a homogeneous media whose density was close to water, but there were significant errors in high-density materials. For a patient case, A and B point dose differences were less than 3%, while the mean dose discrepancy was as much as 5%. Conventional dose computation methods might underdose the targets by not accounting for the effects of high-density materials. The proposed platform was shown to be feasible and to have good dose calculation accuracy. One should be careful when confirming the plan using a conventional brachytherapy dose computation method, and moreover, an independent dose verification system as developed in this study might be helpful.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.10a
• /
• pp.115-120
• /
• 1997

영광3호기 방사선관리구역에 대한 중성자선량률을 정확히 평가하기 위하여 MCNP4A 전산코드를 이용, 방사선관리구역에서의 중성자 스펙트럼 계산을 수행하였다. 영광3호기에 대한 보다 정확하고 정밀한 3차원 몬테칼로 모델을 구축하기 위하여 핵연료집합체 구성요소 및 원자로심을 둘러싸고 있는 baffle, barrel,압력용기 등을 정확하게 묘사하였으며, 특히 방사선관리구역 주위의 구조물에 대해서도 3자원 MCNP 모델을 구축함으로써 원자로심부터 방사선관리구역까지 완전한 몬테칼로 모사(full-scope Monte Carlo simulation)를 이용한 계산을 수행하였다. 계산결과는 에너지 구간에 따른 중성자속 스펙트럼으로 나타내었으며 이 결과를 바탕으로 중성자속에 대한 선량률 환산인자를 고려하여 중성자선량률을 계산할 수 있다.

• Progress in Medical Physics
• /
• v.14 no.4
• /
• pp.234-239
• /
• 2003

We examined the variation of percent depth dose (PDD) curves for 10 MV X-rays in the presence of magnetic fields. The EGS4 Monte Carlo code was applied and modified to take account of the effect of electron deflection under magnetic field was used. We defined and tested DI (dose improvement) and DR (dose reduction) to describe variation of PDD curves under various magnetic fields. For a magnetic field of 3 T applied at the depth region of 5-10 cm and field size of 10${\times}$10 $\textrm{cm}^2$, the DI is 1.56 (56% improvement) and DR is 0.68 (32% reduction). We explained the results from the Lorentz law and the concept of electron equilibrium. We suggested that the dose optimization in radiotherapy can be achieved from using the characteristics of dose distributions under magnetic fields.

• Proceedings of the Korean Statistical Society Conference
• /
• 2000.11a
• /
• pp.135-140
• /
• 2000

베이지안 신경망 모형(Bayesian Neural Networks Models)에서 주어진 입력값(input)은 블랙 박스(Black-Box)와 같은 신경망 구조의 각 층(layer)을 거쳐서 출력값(output)으로 계산된다. 새로운 입력 데이터에 대한 예측값은 사후분포(posterior distribution)의 기대값(mean)에 의해 계산된다. 주어진 사전분포(prior distribution)와 학습데이터에 의한 가능도함수(likelihood functions)를 통해 계산되어진 사후분포는 매우 복잡한 구조를 갖게 됨으로서 기대값의 적분계산에 대한 어려움이 발생한다. 이때 확률적 추정에 의한 근사 방법인 몬테칼로 적분을 이용한다. 이러한 방법으로서 Hybrid Monte Carlo 알고리즘은 우수한 결과를 제공하여준다(Neal 1996). 본 논문에서는 Hybrid Monte Carlo 알고리즘과 기존에 많이 사용되고 있는 Gibbs sampling, Metropolis algorithm, 그리고 Slice Sampling등의 몬테칼로 방법들을 비교한다.