• Title/Summary/Keyword: Monte Carlo 계산

Search Result 515, Processing Time 0.022 seconds

Measurement and Monte Carlo Simulation of 6 MV X-rays for Small Radiation Fields (선형가속기의 6 MV X-선에 대한 소형 조사면 측정과 몬테 카를로 시뮬레이션)

  • Jeong Dong Hyeok;Lee Jeong Ok;Kang Jeong Ku;Kim Soo Kon;Kim Seung Kon;Moon Sun Rock
    • Radiation Oncology Journal
    • /
    • v.16 no.2
    • /
    • pp.195-202
    • /
    • 1998
  • Purpose : In order to obtain basic data for treatment plan in radiosurgery, we measured small fields of 6 MV X-rays and compared the measured data with our Monte Carlo simulations for the small fields. Materials and Methods : The small fields of 1.0, 2.0 and 3.0 cm in diameter were used in this study. Percentage depth dose (PDD) and beam Profiles of those fields were measured and calculated. A small semiconductor detector, water phantoms, and a remote control system were used for the measurement Monte Carlo simulations were Performed using the EGS4 code with the input data prepared for the energy distribution of 6 MV X-rays, beam divergence, circular fields and the geometry of the water phantoms. Results : In the case of PDD values, the calculated values were lower than the measured values for all fields and depths, with the differences being 0.3 to 5.7% at the depths of 20 to 20.0 cm and 0.0 to 8.9% at the surface regions. As a result of the analysis of beam profiles for all field sizes at a depth of loom in water phantom, the measured 90% dose widths were in good agreement with the calculated values, however, the calculated Penumbra radii were 0.1 cm shorter than measured values. Conclusion : The measured PDDs and beam profiles agreement with the Monte Carlo calculations approximately. However, it is different when it comes to calculations in the area of phantom surface and penumbra because the Monte Carlo calculations were performed under the simplified geometries. Therefore, we have to study how to include the actual geometries and more precise data for the field area in Monte Carlo calculations. The Monte Carlo calculations will be used as a useful tool for the very complicated conditions in measurement and verification.

  • PDF

Analysis of Electromagnetic Wave Scattering From a Perfectly Conducting One Dimensional Fractal Surface Using the Monte-Carlo Moment Method (몬테칼로 모멘트 방법을 이용한 1차원 프랙탈 완전도체 표면에서의 전자파 산란 해석)

  • 최동묵;김채영
    • Journal of the Institute of Electronics Engineers of Korea TC
    • /
    • v.39 no.12
    • /
    • pp.566-574
    • /
    • 2002
  • In this paper, the scattered field from a perfectly conducting fractal surface by the Monte-Carlo moment method was computed. An one-dimensional fractal surface was generated by using the fractional Brownian motion model. Back scattering coefficients are calculated with different values of the spectral parameter(S$\_$0/), and fractal dimension(D) which determine characteristics of the fractal surface. The number of surface realization for the computed field, the point number, and the width of surface realization are set to be 80, 2048, and 64L, respectively. In order to verify the computed results these results are compared with those of small perturbation methods, which show good agreement between them.

Evaluation of Factors Used in AAPM TG-43 Formalism Using Segmented Sources Integration Method and Monte Carlo Simulation: Implementation of microSelectron HDR Ir-192 Source (미소선원 적분법과 몬테칼로 방법을 이용한 AAPM TG-43 선량계산 인자 평가: microSelectron HDR Ir-192 선원에 대한 적용)

  • Ahn, Woo-Sang;Jang, Won-Woo;Park, Sung-Ho;Jung, Sang-Hoon;Cho, Woon-Kap;Kim, Young-Seok;Ahn, Seung-Do
    • Progress in Medical Physics
    • /
    • v.22 no.4
    • /
    • pp.190-197
    • /
    • 2011
  • Currently, the dose distribution calculation used by commercial treatment planning systems (TPSs) for high-dose rate (HDR) brachytherapy is derived from point and line source approximation method recommended by AAPM Task Group 43 (TG-43). However, the study of Monte Carlo (MC) simulation is required in order to assess the accuracy of dose calculation around three-dimensional Ir-192 source. In this study, geometry factor was calculated using segmented sources integration method by dividing microSelectron HDR Ir-192 source into smaller parts. The Monte Carlo code (MCNPX 2.5.0) was used to calculate the dose rate $\dot{D}(r,\theta)$ at a point ($r,\theta$) away from a HDR Ir-192 source in spherical water phantom with 30 cm diameter. Finally, anisotropy function and radial dose function were calculated from obtained results. The obtained geometry factor was compared with that calculated from line source approximation. Similarly, obtained anisotropy function and radial dose function were compared with those derived from MCPT results by Williamson. The geometry factor calculated from segmented sources integration method and line source approximation was within 0.2% for $r{\geq}0.5$ cm and 1.33% for r=0.1 cm, respectively. The relative-root mean square error (R-RMSE) of anisotropy function obtained by this study and Williamson was 2.33% for r=0.25 cm and within 1% for r>0.5 cm, respectively. The R-RMSE of radial dose function was 0.46% at radial distance from 0.1 to 14.0 cm. The geometry factor acquired from segmented sources integration method and line source approximation was in good agreement for $r{\geq}0.1$ cm. However, application of segmented sources integration method seems to be valid, since this method using three-dimensional Ir-192 source provides more realistic geometry factor. The anisotropy function and radial dose function estimated from MCNPX in this study and MCPT by Williamson are in good agreement within uncertainty of Monte Carlo codes except at radial distance of r=0.25 cm. It is expected that Monte Carlo code used in this study could be applied to other sources utilized for brachytherapy.

Monte-Carlo Based Interference Analysis of 900MHz Passive RFID systems (몬테카를로(Monte-Carlo) 방법을 적용한 수동형 900MHz 대역의 RFID 간섭 분석)

  • Kim Young-Hwan;Eo Pill-Seon;Yang Hoon-Gee;Park Seung-Keun;Kang Bongsoon;Kim Young-Soo;Yook Jong-Gwan
    • Journal of the Institute of Electronics Engineers of Korea TC
    • /
    • v.43 no.1 s.343
    • /
    • pp.9-18
    • /
    • 2006
  • This paper presents a Monte-Carlo based method to obtain a probability of interference in the 900MHz passive RFID systems. We show an efficient algorithm to calculate not only in-band interference due to unwanted emission of interfering sources for a given emission mask, but out-of-band interference, which depends on the blocking performances of a victim receiver filter. We next apply the proposed method to two types of RFID systems, one is to use simple FHSS and the other to adopt a hybrid use of FHSS and LBT, which senses the channel before transmission. Simulation is first performed with a SEAMCAT. We next make simulation with a MATLAB software which implements SEAMCAT algorithm and show both results from two softwares are similar.

Comparative Study on the Applicability of Point Estimate Methods in Combination with Numerical Analysis for the Probabilistic Reliability Assessment of Underground Structures (수치해석과 연계한 지하구조물의 확률론적 신뢰성 평가를 위한 점추정법의 적용성에 관한 비교 연구)

  • Park, Do-Hyun;Kim, Hyung-Mok;Ryu, Dong-Woo;Choi, Byung-Hee;Han, Kong-Chang
    • Tunnel and Underground Space
    • /
    • v.22 no.2
    • /
    • pp.86-92
    • /
    • 2012
  • Point estimate method has a less accuracy than Monte Carlo simulation that is usually considered as an exact probabilistic method, but this method still remains popular in probability-based reliability assessment in geotechnical and rock engineering, because it significantly reduce the number of sampling points and produces the statistical moments of a performance function in a reasonable accuracy. In the present study, we investigated the accuracy and applicability of point estimate methods proposed by Rosenblueth and Zhou & Nowak by comparing the results of these two methods with those of Monte Carlo simulations. The comparison was carried out for the problem of a lined circular tunnel in an elastic medium where an closed-form analytical solution is given. The comparison results showed that despite the non-linearity of the analytical solution, the statistical moments calculated by the point estimate methods and the Monte Carlo simulations agreed well with an average error of roughly 1-2%. This average error demonstrates the applicability of the two point estimate methods for the probabilistic reliability assessment of underground structures in combination with numerical analysis.

Dose Verification Study of Brachytherapy Plans Using Monte Carlo Methods and CT Images (CT 영상 및 몬테칼로 계산에 기반한 근접 방사선치료계획의 선량분포 평가 방법 연구)

  • Cheong, Kwang-Ho;Lee, Me-Yeon;Kang, Sei-Kwon;Bae, Hoon-Sik;Park, So-Ah;Kim, Kyoung-Joo;Hwang, Tae-Jin;Oh, Do-Hoon
    • 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.

Benchmark Results of a Monte Carlo Treatment Planning system (몬데카를로 기반 치료계획시스템의 성능평가)

  • Cho, Byung-Chul
    • Progress in Medical Physics
    • /
    • v.13 no.3
    • /
    • pp.149-155
    • /
    • 2002
  • Recent advances in radiation transport algorithms, computer hardware performance, and parallel computing make the clinical use of Monte Carlo based dose calculations possible. To compare the speed and accuracies of dose calculations between different developed codes, a benchmark tests were proposed at the XIIth ICCR (International Conference on the use of Computers in Radiation Therapy, Heidelberg, Germany 2000). A Monte Carlo treatment planning comprised of 28 various Intel Pentium CPUs was implemented for routine clinical use. The purpose of this study was to evaluate the performance of our system using the above benchmark tests. The benchmark procedures are comprised of three parts. a) speed of photon beams dose calculation inside a given phantom of 30.5 cm$\times$39.5 cm $\times$ 30 cm deep and filled with 5 ㎣ voxels within 2% statistical uncertainty. b) speed of electron beams dose calculation inside the same phantom as that of the photon beams. c) accuracy of photon and electron beam calculation inside heterogeneous slab phantom compared with the reference results of EGS4/PRESTA calculation. As results of the speed benchmark tests, it took 5.5 minutes to achieve less than 2% statistical uncertainty for 18 MV photon beams. Though the net calculation for electron beams was an order of faster than the photon beam, the overall calculation time was similar to that of photon beam case due to the overhead time to maintain parallel processing. Since our Monte Carlo code is EGSnrc, which is an improved version of EGS4, the accuracy tests of our system showed, as expected, very good agreement with the reference data. In conclusion, our Monte Carlo treatment planning system shows clinically meaningful results. Though other more efficient codes are developed such like MCDOSE and VMC++, BEAMnrc based on EGSnrc code system may be used for routine clinical Monte Carlo treatment planning in conjunction with clustering technique.

  • PDF

Reliability Assessment for Corroded Pipelines by Separable Monte Carlo Method (Separable Monte Carlo 방법을 적용한 부식배관 신뢰도평가)

  • Lee, Jin-Han;Jo, Young-Do;Kim, Lae Hyun
    • Journal of the Korean Institute of Gas
    • /
    • v.19 no.5
    • /
    • pp.81-86
    • /
    • 2015
  • A deterministic stress-based methodology has traditionally been applied in pipeline design. Meanwhile, reliability based design and assessment (RBDA) methodology has been extensively applied in offshore or nuclear structures. Lately, the release of ISO standard on reliability based limit state methods for pipelines ISO16708 indicates that the RBDA methodology is one of the newest directions of natural gas pipeline design method. This paper presents a case study of the RBDA procedure for predicting the time-dependent failure probability of pipelines with corrosion defects, where separable Monte Carlo (SMC) method is applied in the reliability estimation for corroded pipeline instead of traditional, crude Monte Carlo(CMC) Method. The result shows the SMC method take advantage of improving accuracy in reliability calculation.

Assessment of Dose Distribution using the MIRD Phantom at Uterine Cervix and Surrounding Organs in High Doserate Brachytheraphy (자궁주위 방사선 근접치료시 MIRD 팬텀을 이용한 주변장기의 피폭환경평가)

  • Lee, Yun-Jong;Nho, Young-Chang;Lee, Jai-Ki
    • Korean Journal of Environmental Biology
    • /
    • v.24 no.4
    • /
    • pp.387-391
    • /
    • 2006
  • Computational and experimental dosimetry of Henschke applicator with respect to high dose rate brachytherapy using the MIRD phantom and a remote control afterloader were performed. A comparison of computational dosimetry was made between the simulated Monte Carlo dosimetry and GAMMADOT brachytherapy Planning system's dosimetry. Dose measurements was performed using ion chamber in a water phantom. Dose rates are calculated using Monte Carlo code MCNP4B and the GAMMADOT. Thecomputational models include the detailed geometry of Ir-192 source, tandem tube, and shielded ovoids for accurate estimation. And transit dose delivered during source extension to and retraction from a given dwell position was estimated by Monte Carlo simulations. Point doses at ICRU bladder/rectal pointswhich have been recommened by ICRU 38 was assessed. Calculated and measured dose distribution data agreed within 4% each other. The shielding effect of ovoids leads to 19% and 20% dose reduction at bladder surface and rectal points.

Calculation of the Detection Range for a Given Cumulative Probability in Airborne Surveillance Radars (탐색 레이다에서 누적확률에 기인한 탐지거리 계산에 관한 연구)

  • Kim, Eun Hee;Roh, Ji-Eun
    • The Journal of Korean Institute of Electromagnetic Engineering and Science
    • /
    • v.29 no.1
    • /
    • pp.24-27
    • /
    • 2018
  • The performance measure of airborne radars is the range at which the cumulative probability of detection has some specified value, because the per-scan detection probability is an oscillatory function of the target range in airborne radars operating with the dynamic clutter environment. As a result, no one range, at which the per-scan detection probability has a given value, can give a meaningful description of the range performance. In this paper, we provide the equation to calculate the cumulative detection probability and show that the result of Monte Carlo simulation is same as the calculated value in a simple scenario. This verified Monte Carlo model will be used to evaluate the performance of airborne radars in various operating scenarios, at which the numerical calculation is difficult.