A high-order deterministic dynamical low-rank method for proton transport in heterogeneous media

TL;DR

This paper introduces a dynamical low-rank approximation method for efficiently solving high-dimensional proton transport equations in heterogeneous media, significantly reducing computational costs while maintaining accuracy.

Contribution

It develops a novel low-rank model order reduction technique for proton transport, combining phase space discretizations and a mixture model, enabling high-resolution calculations at lower costs.

Findings

Reproduces full-rank results with lower rank and cost

Achieves high accuracy in heterogeneous media

Efficiently computes multiple beam sources with minimal additional cost

Abstract

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.