A high-order deterministic dynamical low-rank method for proton transport in heterogeneous media
By: Pia Stammer , Niklas Wahl , Jonas Kusch and more
Potential Business Impact:
Speeds up cancer treatment planning with faster math.
Dose calculations in proton therapy require the fast and accurate solution of a high-dimensional transport equation for a large number of (pencil) beams with different energies and directions. Deterministically solving this transport problem at a sufficient resolution can however be prohibitively expensive, especially due to highly forward peaked scattering of the protons. We propose using a model order reduction approach, the dynamical low-rank approximation (DLRA), which evolves the solution on the manifold of low-rank matrices in (pseudo-)time. For this, we compare a collided-uncollided split of the linear Boltzmann equation and its Fokker-Planck approximation. We treat the uncollided part using a ray-tracer and combine high-order phase space discretizations and a mixture model for materials with DLRA for the collided equation. Our method reproduces the results of a full-rank reference code at significantly lower rank, and thus computational cost and memory, and further makes computations feasible at much higher resolutions. At higher resolutions, we also achieve good accuracy with respect to TOPAS MC in homogeneous as well as heterogeneous materials. Finally, we demonstrate that several beam sources with different angles can be computed with little cost increase compared to individual beams.
Similar Papers
Discontinuous Galerkin discretization of conservative dynamical low-rank approximation schemes for the Vlasov-Poisson equation
Numerical Analysis
Simulates how charged particles move and interact.
Interpolatory Dynamical Low-Rank Approximation: Theoretical Foundations and Algorithms
Numerical Analysis
Makes computer models run much faster.
A Unified Framework from Boltzmann Transport to Proton Treatment Planning
Probability
Improves cancer treatment planning with math.