Journal of Radiation Protection and Research

The Korean Association for Radiation Protection (대한방사선방어학회)
권호 표지
DOI QR코드

DOI QR Code

TET2MCNP: A Conversion Program to Implement Tetrahedral-mesh Models in MCNP

  • Han, Min Cheol (Department of Nuclear Engineering, Hanyang University) ;
  • Yeom, Yeon Soo (Department of Nuclear Engineering, Hanyang University) ;
  • Nguyen, Thang Tat (Department of Nuclear Engineering, Hanyang University) ;
  • Choi, Chansoo (Department of Nuclear Engineering, Hanyang University) ;
  • Lee, Hyun Su (Department of Nuclear Engineering, Hanyang University) ;
  • Kim, Chan Hyeong (Department of Nuclear Engineering, Hanyang University)
  • Received : 2016.08.22
  • Accepted : 2016.09.19
  • Published : 2016.12.31
Download PDF
⟨ Previous Next ⟩

Abstract

Background: Tetrahedral-mesh geometries can be used in the MCNP code, but the MCNP code accepts only the geometry in the Abaqus input file format; hence, the existing tetrahedral-mesh models first need to be converted to the Abacus input file format to be used in the MCNP code. In the present study, we developed a simple but useful computer program, TET2MCNP, for converting TetGen-generated tetrahedral-mesh models to the Abacus input file format. Materials and Methods: TET2MCNP is written in C++ and contains two components: one for converting a TetGen output file to the Abacus input file and the other for the reverse conversion process. The TET2MCP program also produces an MCNP input file. Further, the program provides some MCNP-specific functions: the maximum number of elements (i.e., tetrahedrons) per part can be limited, and the material density of each element can be transferred to the MCNP input file. Results and Discussion: To test the developed program, two tetrahedral-mesh models were generated using TetGen and converted to the Abaqus input file format using TET2MCNP. Subsequently, the converted files were used in the MCNP code to calculate the object- and organ-averaged absorbed dose in the sphere and phantom, respectively. The results show that the converted models provide, within statistical uncertainties, identical dose values to those obtained using the PHITS code, which uses the original tetrahedral-mesh models produced by the TetGen program. The results show that the developed program can successfully convert TetGen tetrahedral-mesh models to Abacus input files. Conclusion: In the present study, we have developed a computer program, TET2MCNP, which can be used to convert TetGen-generated tetrahedral-mesh models to the Abaqus input file format for use in the MCNP code. We believe this program will be used by many MCNP users for implementing complex tetrahedral-mesh models, including computational human phantoms, in the MCNP code.

Keywords

References

  1. Goorley JT, et al. Initial MCNP6 release overview-MCNP6 version1.0. Los Alamos National Laboratory. LA-UR-13-22934. 2013;1.
  2. Sato T, et al. Particle and heavy ion transport code system, PHITS, version 2.52. J Nucl Sci Technol. 2013;50:913-923. https://doi.org/10.1080/00223131.2013.814553
  3. Agostinelli S, et al. GEANT4-a simulation toolkit. Nucl. Instrum. Meth. A. 2003;506:250-303. https://doi.org/10.1016/S0168-9002(03)01368-8
  4. Cowan P, Shuttleworth E, Bird A, Cooper A. The launch of MCBEND 10. Radiat. Prot. Dosim. 2005;115:403-7. https://doi.org/10.1093/rpd/nci179
  5. Badal A, Kyprianou I, Banh DP, Badano A, Sempau J. penMesh-Monte Carlo radiation transport simulation in a triangle mesh geometry. IEEE Trans. Med. Imaging. 2009;28:1894-1901. https://doi.org/10.1109/TMI.2009.2021615
  6. Yeom YS, Jeong JH, Han MC, Kim CH. Tetrahedral-mesh-based computational human phantom for fast Monte Carlo dose calculations. Phys. Med. Biol. 2014;59:3173-3185. https://doi.org/10.1088/0031-9155/59/12/3173
  7. Han MC, Yeom YS, Kim CH, Kim S, Sohn JW. New approach based on tetrahedral-mesh geometry for accurate 4D Monte Carlo patient-dose calculation. Phys Med Biol. 2015;60:1601-1612. https://doi.org/10.1088/0031-9155/60/4/1601
  8. Nguyen TT, et al. Incorporation of detailed eye model into polygon-mesh versions of ICRP-110 reference phantoms. Phys. Med. Biol. 2015;60:8695-8707. https://doi.org/10.1088/0031-9155/60/22/8695
  9. Dogdas B, Stout D, Chatziioannou AF, Leahy RM. Digimouse: a 3D whole body mouse atlas from CT and cryosection data. Phys. Med. Biol. 2007;52:577-587. https://doi.org/10.1088/0031-9155/52/3/003
  10. Poole CM, Cornelius I, Trapp JV, Langton CM. Fast tessellated solid navigation in GEANT4. Ieee Trans. Nucl. Sci. 2012;59:1695-1701. https://doi.org/10.1109/TNS.2012.2197415
  11. Fonseca GP, Landry G, White S, D'Amours M, Yoriyaz H, Beaulieu L, Reniers B, Verhaegen F. The use of tetrahedral mesh geometries in Monte Carlo simulation of applicator based brachytherapy dose distributions. Phys. Med. Biol. 2014;59:5921-5935. https://doi.org/10.1088/0031-9155/59/19/5921
  12. Fang QQ. Mesh-based Monte Carlo method using fast ray-tracing in Plucker coordinates. Biomed. Opt. Express. 2010;1:165-175. https://doi.org/10.1364/BOE.1.000165
  13. Si H. TetGen, a Delaunay-based quality tetrahedral mesh generator. Acm. Trans. Math. Software. 2015;41(2):1-36.
  14. Martz R. The MCNP6 book on unstructured mesh geometry: User's guide. Los Alamos National Laboratory, LA-UR-11-05 668-rev8. 2014;19-24.
  15. Kim CH, Jeong JH, Bolch WE, Cho KW, Hwang SB. A polygonsurface reference Korean male phantom (PSRK-Man) and its direct implementation in Geant4 Monte Carlo simulation. Phys. Med. Biol. 2011;56:3137-3161. https://doi.org/10.1088/0031-9155/56/10/016

Cited by

  1. Multi-threading performance of Geant4, MCNP6, and PHITS Monte Carlo codes for tetrahedral-mesh geometry vol.63, pp.9, 2018, https://doi.org/10.1088/1361-6560/aabd20
  2. CAD 모델 기반의 4D CT 데이터 제작 의용공학 융합 프로그램 개발 vol.8, pp.4, 2016, https://doi.org/10.15207/jkcs.2017.8.4.131
  3. Mesh-type reference Korean phantoms (MRKPs) for adult male and female for use in radiation protection dosimetry vol.64, pp.8, 2016, https://doi.org/10.1088/1361-6560/ab0b59
  4. Computation Speeds and Memory Requirements of Mesh-Type ICRP Reference Computational Phantoms in Geant4, MCNP6, and PHITS vol.116, pp.5, 2016, https://doi.org/10.1097/hp.0000000000000999
  5. POLY2TET: a computer program for conversion of computational human phantoms from polygonal mesh to tetrahedral mesh vol.40, pp.4, 2020, https://doi.org/10.1088/1361-6498/abb360
  6. Optimization of target, moderator, and collimator in the accelerator-based boron neutron capture therapy system: A Monte Carlo study vol.53, pp.6, 2016, https://doi.org/10.1016/j.net.2020.12.006