SUET558: Monte Carlo Photon Transport Simulations On GPU with Quadric Geometry
Abstract
Purpose: Monte Carlo simulation on GPU has experienced rapid advancements over the past a few years and tremendous accelerations have been achieved. Yet existing packages were developed only in voxelized geometry. In some applications, e.g. radioactive seed modeling, simulations in more complicated geometry are needed. This abstract reports our initial efforts towards developing a quadric geometry module aiming at expanding the application scope of GPUbased MC simulations. Methods: We defined the simulation geometry consisting of a number of homogeneous bodies, each specified by its material composition and limiting surfaces characterized by quadric functions. A tree data structure was utilized to define geometric relationship between different bodies. We modified our GPUbased photon MC transport package to incorporate this geometry. Specifically, geometry parameters were loaded into GPU’s shared memory for fast access. Geometry functions were rewritten to enable the identification of the body that contains the current particle location via a fast searching algorithm based on the tree data structure. Results: We tested our package in an example problem of HDRbrachytherapy dose calculation for shielded cylinder. The dose under the quadric geometry and that under the voxelized geometry agreed in 94.2% of total voxels within 20% isodose line based on a statisticalmore »
 Authors:
 The University of Texas Southwestern Medical Ctr, Dallas, TX (United States)
 Publication Date:
 OSTI Identifier:
 22496273
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Medical Physics; Journal Volume: 42; Journal Issue: 6; Other Information: (c) 2015 American Association of Physicists in Medicine; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 61 RADIATION PROTECTION AND DOSIMETRY; 62 RADIOLOGY AND NUCLEAR MEDICINE; ALGORITHMS; BRACHYTHERAPY; CHARGEDPARTICLE TRANSPORT; COMPUTERIZED SIMULATION; GEOMETRY; MONTE CARLO METHOD; PHOTON TRANSPORT; RADIATION DOSES; SHIELDING
Citation Formats
Chi, Y, Tian, Z, Jiang, S, and Jia, X. SUET558: Monte Carlo Photon Transport Simulations On GPU with Quadric Geometry. United States: N. p., 2015.
Web. doi:10.1118/1.4924920.
Chi, Y, Tian, Z, Jiang, S, & Jia, X. SUET558: Monte Carlo Photon Transport Simulations On GPU with Quadric Geometry. United States. doi:10.1118/1.4924920.
Chi, Y, Tian, Z, Jiang, S, and Jia, X. Mon .
"SUET558: Monte Carlo Photon Transport Simulations On GPU with Quadric Geometry". United States.
doi:10.1118/1.4924920.
@article{osti_22496273,
title = {SUET558: Monte Carlo Photon Transport Simulations On GPU with Quadric Geometry},
author = {Chi, Y and Tian, Z and Jiang, S and Jia, X},
abstractNote = {Purpose: Monte Carlo simulation on GPU has experienced rapid advancements over the past a few years and tremendous accelerations have been achieved. Yet existing packages were developed only in voxelized geometry. In some applications, e.g. radioactive seed modeling, simulations in more complicated geometry are needed. This abstract reports our initial efforts towards developing a quadric geometry module aiming at expanding the application scope of GPUbased MC simulations. Methods: We defined the simulation geometry consisting of a number of homogeneous bodies, each specified by its material composition and limiting surfaces characterized by quadric functions. A tree data structure was utilized to define geometric relationship between different bodies. We modified our GPUbased photon MC transport package to incorporate this geometry. Specifically, geometry parameters were loaded into GPU’s shared memory for fast access. Geometry functions were rewritten to enable the identification of the body that contains the current particle location via a fast searching algorithm based on the tree data structure. Results: We tested our package in an example problem of HDRbrachytherapy dose calculation for shielded cylinder. The dose under the quadric geometry and that under the voxelized geometry agreed in 94.2% of total voxels within 20% isodose line based on a statistical ttest (95% confidence level), where the reference dose was defined to be the one at 0.5cm away from the cylinder surface. It took 243sec to transport 100million source photons under this quadric geometry on an NVidia Titan GPU card. Compared with simulation time of 99.6sec in the voxelized geometry, including quadric geometry reduced efficiency due to the complicated geometryrelated computations. Conclusion: Our GPUbased MC package has been extended to support photon transport simulation in quadric geometry. Satisfactory accuracy was observed with a reduced efficiency. Developments for charged particle transport in this geometry are currently in progress.},
doi = {10.1118/1.4924920},
journal = {Medical Physics},
number = 6,
volume = 42,
place = {United States},
year = {Mon Jun 15 00:00:00 EDT 2015},
month = {Mon Jun 15 00:00:00 EDT 2015}
}

Purpose: In proton dose calculation, clinically compatible speeds are now achieved with Monte Carlo codes (MC) that combine 1) adequate simplifications in the physics of transport and 2) the use of hardware architectures enabling massive parallel computing (like GPUs). However, the uncertainties related to the transport algorithms used in these codes must be kept minimal. Such algorithms can be checked with the socalled “Fano cavity test”. We implemented the test in two codes that run on specific hardware: gPMC on an nVidia GPU and MCsquare on an Intel Xeon Phi (60 cores). Methods: gPMC and MCsquare are designed for transportingmore »

SUET313: Dosimetric Deviation of Misaligned Beams for a 6 MV Photon Linear Accelerator Using Monte Carlo Simulations
Purpose: To quantify the dosimetric variations of misaligned beams for a linear accelerator by using Monte Carlo (MC) simulations. Method and Materials: Misaligned beams of a Varian 21EX Clinac were simulated to estimate the dosimetric effects. All the linac head components for a 6 MV photon beam were implemented in BEAMnrc/EGSnrc system. For incident electron beam parameters, 6 MeV with 0.1 cm fullwidthhalfmax Gaussian beam was used. A phase space file was obtained below the jaw per each misalignment condition of the incident electron beam: (1) The incident electron beams were tilted by 0.5, 1.0 and 1.5 degrees on themore » 
SUET58: A Novel Monte Carlo Photon Transport Simulation Scheme and Its Application in Cone Beam CT Projection Simulation
Purpose: Monte Carlo (MC) simulation is an important tool to solve radiotherapy and medical imaging problems. Low computational efficiency hinders its wide applications. Conventionally, MC is performed in a particleby particle fashion. The lack of control on particle trajectory is a main cause of low efficiency in some applications. Take cone beam CT (CBCT) projection simulation as an example, significant amount of computations were wasted on transporting photons that do not reach the detector. To solve this problem, we propose an innovative MC simulation scheme with a pathbypath sampling method. Methods: Consider a photon path starting at the xray source.more » 
SUET112: An OpenCLBased CrossPlatform Monte Carlo Dose Engine (oclMC) for Coupled PhotonElectron Transport
Purpose: Low computational efficiency of Monte Carlo (MC) dose calculation impedes its clinical applications. Although a number of MC dose packages have been developed over the past few years, enabling fast MC dose calculations, most of these packages were developed under NVidia’s CUDA environment. This limited their code portability to other platforms, hindering the introduction of GPUbased MC dose engines to clinical practice. To solve this problem, we developed a crossplatform fast MC dose engine named oclMC under OpenCL environment for external photon and electron radiotherapy. Methods: Coupled photonelectron simulation was implemented with standard analogue simulation scheme for photon transportmore » 
SUET29: A Web Application for GPUBased Monte Carlo IMRT/VMAT QA with Delivered Dose Verification
Purpose: To enable an existing web application for GPUbased Monte Carlo (MC) 3D dosimetry quality assurance (QA) to compute “delivered dose” from linac logfile data. Methods: We added significant features to an IMRT/VMAT QA web application which is based on existing technologies (HTML5, Python, and Django). This tool interfaces with python, ccode libraries, and command linebased GPU applications to perform a MCbased IMRT/VMAT QA. The web app automates many complicated aspects of interfacing clinical DICOM and logfile data with cuttingedge GPU software to run a MC dose calculation. The resultant web app is powerful, easy to use, and is ablemore »