Total: 1
The linear Boltzmann model for proton beams is a six-dimensional partial differential equation (PDE). We propose a deterministic solver for the linear Boltzmann model based on scattering decomposition and depth-splitting methods. The main idea is to first divide the protons into primary protons and scattering protons, whose equations are derived using the source iteration method. We then treat depth as the time variable in classical time-evolutionary problems and apply the depth-splitting method. In the depth-splitting method, the full operator is decomposed into three parts, with each subsystem being easily parallelizable, which is crucial for efficient simulations. The resulting discretization exhibits second-order convergence in both the depth and energy variables. The dose distributions obtained from our solver are compared with those from Monte Carlo simulations for various materials and heterogeneous cases.