• Nuclear Engineering and Technology
• /
• v.54 no.5
• /
• pp.1769-1780
• /
• 2022

A new method of dose calculation algorithm, called GPU-accelerated Monte Carlo and collapsed cone Convolution (GMCC) was developed to improve the calculation speed of BNCT treatment planning system. The GPU-accelerated Monte Carlo routine in GMCC is used to simulate the neutron transport over whole energy range and the Collapsed Cone Convolution method is to calculate the gamma dose. Other dose components due to alpha particles and protons, are calculated using the calculated neutron flux and reaction data. The mathematical principle and the algorithm architecture are introduced. The accuracy and performance of the GMCC were verified by comparing with the FLUKA results. A water phantom and a head CT voxel model were simulated. The neutron flux and the absorbed dose obtained by the GMCC were consistent well with the FLUKA results. In the case of head CT voxel model, the mean absolute percentage error for the neutron flux and the absorbed dose were 3.98% and 3.91%, respectively. The calculation speed of the absorbed dose by the GMCC was 56 times faster than the FLUKA code. It was verified that the GMCC could be a good candidate tool instead of the Monte Carlo method in the BNCT dose calculations.

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2002.09a
• /
• pp.116-118
• /
• 2002

It has been noted that Monte Carlo simulations are the most accurate method to calculate dose distributions in any material and geometry. Monte Carlo transport algorithms determine the absorbed dose by following the path of representative particles as they travel through the medium. Accurate Monte Carlo dose calculations rely on detailed modeling of the radiation source. We modeled the effects of beam modifiers such as collimators, blocks, wedges, etc. of our accelerator, Varian Clinac 600C/D to ensure accurate representation of the radiation source using the EGSnrc based BEAM code. These were used in the EGSnrc based DOSXYZ code for the simulation of particles transport through a voxel based Cartesian coordinate system. Because Monte Carlo methods use particle-by-particle methods to simulate a radiation transport, more particle histories yield the better representation of the actual dose. But the prohibitively long time required to get high resolution and accuracy calculations has prevented the use of Monte Carlo methods in the actual clinical spots. Our ultimate aim is to develop a Monte Carlo dose calculation system designed specifically for radiation therapy planning, which is distinguished from current dose calculation methods. The purpose of this study in the present phase was to get dose calculation results corresponding to measurements within practical time limit. We used parallel processing and some variance reduction techniques, therefore reduced the computational time, preserving a good agreement between calculations of depth dose distributions and measurements within 5% deviations.

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

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2002.09a
• /
• pp.106-108
• /
• 2002

Recent advances in radiation transport algorithms, computer hardware performance, and parallel computing make the clinical use of Monte Carlo based dose calculations possible. Monte Carlo treatment planning requires accurate beam information as input to generate accurate dose distributions. The procedures to obtain this accurate beam information are called "commissioning", which includes accelerator head modeling. In this study, we would like to investigate how much accurately Monte Carlo based dose calculations can predict the measured beam data in various conditions. The Siemens 6MV photon beam and the BEAM Monte Carlo code were used. The comparisons including the percentage depth doses and off-axis profiles of open fields and wedges, output factors will be presented.

• Proceedings of the Korean Society of Medical Physics Conference
• /
• 2002.09a
• /
• pp.291-293
• /
• 2002

Accurate dose calculation in radiation treatment planning is most important for successful treatment. Since human body is composed of various materials and not an ideal shape, it is not easy to calculate the accurate effective dose in the patients. Many methods have been proposed to solve the inhomogeneity and surface contour problems. Monte Carlo simulations are regarded as the most accurate method, but it is not appropriate for routine planning because it takes so much time. Pencil beam kernel based convolution/superposition methods were also proposed to correct those effects. Nowadays, many commercial treatment planning systems, including Pinnacle and Helax-TMS, have adopted this algorithm as a dose calculation engine. The purpose of this study is to verify the accuracy of the dose calculated from pencil beam kernel based treatment planning system Helax-TMS comparing to Monte Carlo simulations and measurements especially in inhomogeneous region. Home-made inhomogeneous phantom, Helax-TMS ver. 6.0 and Monte Carlo code BEAMnrc and DOSXYZnrc were used in this study. Dose calculation results from TPS and Monte Carlo simulation were verified by measurements. In homogeneous media, the accuracy was acceptable but in inhomogeneous media, the errors were more significant.

• Journal of Radiation Protection and Research
• /
• v.26 no.2
• /
• pp.93-99
• /
• 2001

Voxel head phantom for overcoming the limitation of mathematical phantom in depleting anatomical details was constructed and example dose calculation for BNCT was performed. The repeated structure algorithm of the general purpose Monte Carlo code, MCNP4B was applied for yokel Monte Carlo calculation. Simple binary yokel phantom and combinatorial geometry phantom composed of two materials were constructed for validating the voxel Monte Carlo calculation system. The tomographic images of VHP man provided by NLM(National Library of Medicine) were segmented and indexed to construct yokel head phantom. Comparison of doses for broad parallel gamma and neutron beams in AP and PA directions showed decrease of brain dose due to the attenuation of neutron in eye balls in case of yokel head phantom. The spherical tumor volume with diameter, 5cm was defined in the center of brain for BNCT dose calculation in which accurate 3 dimensional dose calculation is essential. As a result of BNCT dose calculation for downward neutron beam of 10keV and 40keV, the tumor dose is about doubled when boron concentration ratio between the tumor to the normal tissue is $30{\mu}g/g$ to $3{\mu}g/g$. This study established the voxel Monte Carlo calculation system and suggested the feasibility of precise dose calculation in therapeutic radiology.

• 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.11 no.2
• /
• pp.91-99
• /
• 2000

This is a preliminary study for developing the method of the dose reconstruction in the patients, irradiated by mega-voltage photon beams from the linear accelerator, using the transit dose distributions. In this study we present the method of three-dimensional dose reconstruction and evaluate the method by computer simulation. To acquire the dose distributions in the patients (or phantoms) we first calculate the differences between the doses at the arbitrary points in the patients and the doses at the corresponding points where the transit doses are measured. Then, we can get the dose in the patients from the measured transit dose and the calculated value of the difference. The dose differences are calculated by applying the inverse square law and using the linear attenuation coefficient. The scatter to primary dose ratios, which are calculated by the Monte Carlo program using the CT data of the patient (or phantoms), are also used in the calculations. For the evaluation of this method we used various kinds of homogeneous and inhomogeneous phantoms and calculated the transit dose distributions with the Monte Carlo program. From the distributions we reconstructed the dose distributions in the phantom. We used mono-energy Photon beam of 1.5MeV and Monte Carlo program EGS4. The comparison between the dose distributions reconstructed using the method and the distributions calculated by the Monte Carlo program was done. They agreed within errors of -4%∼+2%. This method can be used to predict the dose distributions in the patient

• 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.29 no.2
• /
• pp.105-115
• /
• 2004

The Monte Carlo simulation requires very much time to obtain a result of acceptable accuracy. Therefore we should know the optimum number of history not to sacrifice time as well as the accuracy. In this study, we have investigated the effects of statistical uncertainties of the photon dose calculation. BEAMnrc and DOSXYZnrc systems were used for the Monte Carlo dose calculation and the case of mediastinum was simulated. The several dose calculation result from various number of histories had been obtained and analyzed using the criteria of isodose curve comparison, dose volume histogram comparison(DVH) and root mean-square differences(RMSD). Statistical uncertainties were observed most evidently in isodose curve comparison and RMSD while DVHs were less sensitive. The acceptable uncertainties $(\bar{{\Delta}D})$ of the Monte Carlo photon dose calculation for the radiation therapy were estimated within total 9% error or 1% error for over than $D_{max}/2$ voxels or voxels at maximum dose.