• Journal of computational fluids engineering
• /
• v.4 no.3
• /
• pp.63-70
• /
• 1999

A Direct simulation Monte-Carlo (DSMC) code is developed, which employs the Monte-Carlo statistical sampling technique to investigate hypersonic rarefied gas flows accompanying chemical reactions. The DSMC method is a numerical simulation technique for analyzing the Boltzmann equation by modeling a real gas flow using a representative set of molecules. Due to the limitations in computational requirements. the present method is applied to a flow around a simple two-dimensional object in exit velocity of 7.6 km/sec at an altitude of 90 km. For the calculation of chemical reactions an air model with five species (O₂, N₂, O, N, NO) and 19 chemical reactions is employed. The simulated result showed various rarefaction effects in the hypersonic flow with chemical reactions.

• Journal of the Korean Society of Radiology
• /
• v.11 no.4
• /
• pp.189-195
• /
• 2017

Radiation protection equipment has widely used to protect human body from radiations, for example X-ray and gamma ray. The material of the radiation protection equipment is mainly lead (Pb) which has brought out lead poisoning and pollution when the equipment is fallen into disuse. This problem makes research and development find new Pb-free materials for use of radiation protection. Manufacturing and evaluation processes for developing those material were carried out repletely until obtaining the performance of protection rate. In this study, combination possibility of shielding material was studied using Geant4 monte carlo simulation. X-ray tube under the same condition in the real measurement of the protection rate was simulated, and X-ray tube spectrum was obtained. The X-ray tube spectrum was applied to study on the protection rate and lead equivalent. The porosity effect was simulated, and was one of key factors to determine protection rate or lead equivalent in radiation protection sheet of Pb-free.

• 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

• 대한방사선방어학회:학술대회논문집
• /
• 2011.11a
• /
• pp.188-189
• /
• 2011

• Progress in Medical Physics
• /
• v.22 no.2
• /
• pp.79-84
• /
• 2011

In this study, Geant4 based Monte Carlo simulations were carried out for medical linear accelerator. Modified Medical Linac2 toolkit was used for calculation. The energy spectrum, most probable energy and the photon mean energy compared with the published results using the EGS4 code. The results well agreed with published results. The calculated results of photon fluence, energy fluence and mean energy according to the radius from the centre of the beam were analyzed. Monte Carlo simulation using Medical Linac2 code is considered to be useful for analysis of medical linear accelerator. Because the calculated results varies depending on Physics List model for same head structure. It it important to choose the right model for research purpose. Monte Carlo simulation using GEANT4 Medical Linac2 is a valuable for any novice to adopt this code to the study related to 6 MV photon fluence from medical linear accelerator.

• Progress in Medical Physics
• /
• v.24 no.4
• /
• pp.259-264
• /
• 2013

The DICOM converter program of the Geant4 Monte Carlo simulation code for the application of radiotherapy was developed. We analysis the header part of the DICOM file and find various parameters, such as matrix size, pixel size, stored data bits, high bit, and padding values. Especially we evaluate every pixel value of the DICOM files. To conform the exact convert of the pixel values, we developed the verify program. As a result, the DICOM formats generated from difference CT vendors can be converted and verified for Genat4 calculations.

• Progress in Medical Physics
• /
• v.23 no.2
• /
• pp.81-90
• /
• 2012

The effect of setup uncertainties on CTV dose and the correlation between setup uncertainties and setup margin were evaluated by Monte Carlo based numerical simulation. Patient specific information of IMRT treatment plan for rectal cancer designed on the VARIAN Eclipse planning system was utilized for the Monte Carlo simulation program including the planned dose distribution and tumor volume information of a rectal cancer patient. The simulation program was developed for the purpose of the study on Linux environment using open source packages, GNU C++ and ROOT data analysis framework. All misalignments of patient setup were assumed to follow the central limit theorem. Thus systematic and random errors were generated according to the gaussian statistics with a given standard deviation as simulation input parameter. After the setup error simulations, the change of dose in CTV volume was analyzed with the simulation result. In order to verify the conventional margin recipe, the correlation between setup error and setup margin was compared with the margin formula developed on three dimensional conformal radiation therapy. The simulation was performed total 2,000 times for each simulation input of systematic and random errors independently. The size of standard deviation for generating patient setup errors was changed from 1 mm to 10 mm with 1 mm step. In case for the systematic error the minimum dose on CTV $D_{min}^{stat{\cdot}}$ was decreased from 100.4 to 72.50% and the mean dose $\bar{D}_{syst{\cdot}}$ was decreased from 100.45% to 97.88%. However the standard deviation of dose distribution in CTV volume was increased from 0.02% to 3.33%. The effect of random error gave the same result of a reduction of mean and minimum dose to CTV volume. It was found that the minimum dose on CTV volume $D_{min}^{rand{\cdot}}$ was reduced from 100.45% to 94.80% and the mean dose to CTV $\bar{D}_{rand{\cdot}}$ was decreased from 100.46% to 97.87%. Like systematic error, the standard deviation of CTV dose ${\Delta}D_{rand}$ was increased from 0.01% to 0.63%. After calculating a size of margin for each systematic and random error the "population ratio" was introduced and applied to verify margin recipe. It was found that the conventional margin formula satisfy margin object on IMRT treatment for rectal cancer. It is considered that the developed Monte-carlo based simulation program might be useful to study for patient setup error and dose coverage in CTV volume due to variations of margin size and setup error.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.3
• /
• pp.203-211
• /
• 2023

A Markov chain Monte Carlo (MCMC) simulation is proposed for probabilistic full waveform inversion (FWI) in a layered half-space. Dynamic responses on the half-space surface are estimated using the thin-layer method when a harmonic vertical force is applied. Subsequently, a posterior probability distribution function and the corresponding objective function are formulated to minimize the difference between estimations and observed data as well as that of model parameters from prior information. Based on the gradient of the objective function, a proposal distribution and an acceptance probability for MCMC samples are proposed. The proposed MCMC simulation is applied to several layered half-space examples. It is demonstrated that the proposed MCMC simulation for probabilistic FWI can estimate probabilistic material properties such as the shear-wave velocities of a layered half-space.

• Journal of radiological science and technology
• /
• v.43 no.6
• /
• pp.443-451
• /
• 2020

This study is to evaluate absorbed dose from right lung for brachytherapy and to estimate the effects of tissue heterogeneities on dose distribution for Iridium-192 source using Monte Carlo simulation. The study employed Geant4 code as Monte Carlo simulation to calculate the dosimetry parameters. The dose distribution of Iridium-192 source in solid water equivalent phantom including aluminium plate or steel plate inserted was calculated and compared with the measured dose by the ion chamber at various distances. And the simulation was used to evaluate the dose of gamma radiation absorbed in the lung organ and other organs around it. The dose distribution embedded in right lung was calculated due to the presence of heart, thymus, spine, stomach as well as left lung. The geometry of the human body was made up of adult male MIRD type of the computational human phantom. The dosimetric characteristics obtained for aluminium plate inserted were in good agreement with experimental results within 4%. The simulation results of steel plate inserted agreed well with a maximum difference 2.75%. Target organ considered to receive a dose of 100%, the surrounding organs were left the left lung of 3.93%, heart of 10.04%, thymus of 11.19%, spine of 12.64% and stomach of 0.95%. When the statistical error is performed for the computational human phantom, the statistical error of value is under 1%.

• 대한방사선방어학회:학술대회논문집
• /
• 2009.11a
• /
• pp.184-185
• /
• 2009