Abstract
Purpose: To apply deterministic techniques to provide a numerical solution for the 3‐D integro‐differential form of the Boltzmann transport equation using the discrete ordinates (SN) method for a radiotherapy unit. Method and Materials: Comparisons were made between deterministic and Monte Carlo simulations considering angular quadrature, spatial discretization (also spatial differencing), and cross section library selection. In order to assess these issues, we simulated the 60Co beam from a standard Eldorado radiotherapy unit. In addition, we simulated the gantry, air gaps, and a water phantom. A 3‐D spatial distribution was generated to yield a 24 ‘z‐level’ model of 150,000 fine mesh cells. To accurately compare our deterministic results with Monte Carlo, we considered an equivalent MCNP5 model with ENDF/B‐VI cross‐section data libraries and volumetric flux (F4) tallies along the central axis of the model. Results: Overall, deterministic SN results yielded good agreement with Monte Carlo results within the Monte Carlo stochastic error for the reference case. At the same time, the parallel SN model provided a full 3D‐scalar flux distribution over the entire geometry within an acceptable running time. To accurately represent the radiation transport, a quadrature level of at least S24 was required with an added ordinate splitting refinement of 6 directions of interest. These results where obtained when the Directional Theta‐Weighted (DTW) differencing scheme was locked for the regions with low density material (air regions), and the adaptive differencing scheme with upgrades up to the Exponential‐Directional Weighted (EDW) differencing scheme was prescribed for higher density regions. Conclusions: The discrete ordinates (SN) method provides a viable alternative to the Monte Carlo method for solving the Boltzmann transport equation. With the proper discretization, an adaptive differencing scheme, and suitable quadratures with ordinate splitting, the SN method is a useful tool in determining the resulting radiation transport throughout a radiotherapy unit.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.