Abstract
Purpose:The Boltzmann‐Fokker‐Planck equation (BFPE) accurately models the migration of photons/charged particles in tissues. While the Monte Carlo (MC) method is popular for solving BFPE in a statistical manner, we aim to develop a deterministic BFPE solver based on various state‐of‐art numerical acceleration techniques for rapid and accurate dose calculation.Methods:Our BFPE solver is based on the structured grid that is maximally parallelizable, with the discretization in energy, angle and space, and its cross section coefficients are derived or directly imported from the Geant4 database. The physical processes that are taken into account are Compton scattering, photoelectric effect, pair production for photons, and elastic scattering, ionization and bremsstrahlung for charged particles.While the spatial discretization is based on the diamond scheme, the angular discretization synergizes finite element method (FEM) and spherical harmonics (SH). Thus, SH is used to globally expand the scattering kernel and FFM is used to locally discretize the angular sphere. As a Result, this hybrid method (FEM‐SH) is both accurate in dealing with forward‐peaking scattering via FEM, and efficient for multi‐energy‐group computation via SH. In addition, FEM‐SH enables the analytical integration in energy variable of delta scattering kernel for elastic scattering with reduced truncation error from the numerical integration based on the classic SH‐based multi‐energy‐group method.Results:The accuracy of the proposed BFPE solver was benchmarked against Geant4 for photon dose calculation. In particular, FEM‐SH had improved accuracy compared to FEM, while both were within 2% of the results obtained with Geant4.Conclusion:A deterministic solver of the Boltzmann‐Fokker‐Planck equation is developed for dose calculation, and benchmarked against Geant4.Xiang Hong and Hao Gao were partially supported by the NSFC (#11405105), the 973 Program (#2015CB856000) and the Shanghai Pujiang Talent Program (#14PJ1404500).
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have