Score: 0

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

Published: August 6, 2025 | arXiv ID: 2508.04484v1

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.

Country of Origin
🇳🇱 Netherlands

Page Count
26 pages

Category
Mathematics:
Numerical Analysis (Math)