SU‐FF‐I‐50: Whole Body and Distal Organ‐Specific Dosimetry Using Parallel SN Methods

• Published in: Medical physics (Lancaster) Vol. 34; no. 6; p. 2349
• Main Authors: Al‐Basheer, A, Ghita, M, Sjoden, G, Bolch, W, Lee, C
• Format: Journal Article
• Language: English
• Published: American Association of Physicists in Medicine 01-06-2007
• Subjects: Anatomy; Boltzmann equations; Dosimetry; Medical physics; Monte Carlo methods
• ISSN: 0094-2405; 2473-4209
• DOI: 10.1118/1.2760427

Purpose: To create a new methodology for dose computations applicable to general medical physics applications, based upon a direct deterministic solution of the (3‐D) Boltzmann transport equation (BTE). Method and Materials: Using the multigroup discrete ordinates (SN) deterministic method in the PENTRAN‐MP (Parallel Environment Neutral‐particle TRANsport‐Medical Physics) code system, fluxes and corresponding doses were determined for a clinical test case using a UF Series B whole‐body voxelized pediatric patient phantom. With a 90 keV planar x‐ray source, we performed two sets of transport calculations based on the same phantom model, employing both the Monte Carlo (MC) and SN methods. SN calculations were performed using PENTRAN, and the MCNP5 code was used for MC calculations. Post processing in the PENTRAN‐MP code system includes seamless parallel data extraction using the PENDATA code, followed by application of the 3D‐DOSE code to compute dose in each phantom voxel, with dose‐volume histograms for critical organs of interest. Results: We demonstrate that the deterministic SN and MC solutions are, in the vicinity of the source region, completely in agreement; the largest statistical error of the MC simulation was <5%, typically averaging <2%. Moreover, the deterministic SN method also provided globally converged results in each voxel at sites distal to the source (in ∼3 hrs on 16 processors), where the MC results on the same computing platform had significant out of field error, and would require a significantly longer running time (estimated at >2000 hrs) to produce meaningful results with an acceptable stochastic error to enable comparison with the well‐converged SN results. Conclusions: A new methodology for 3D dose calculations has been developed based on a parallel discrete ordinates (SN) solution of the BTE. With the proper discretization applied, the SN method presents an accurate, fast solution yielding a complete 3‐D dose distribution without stochastic error.