An Approximation of the M₂ Closure: Application to Radiotherapy Dose Simulation

Particle transport in radiation therapy can be modelled by a kinetic equation which must be solved numerically. Unfortunately, the numerical solution of such equations is generally too expensive for applications in medical centers. Moment methods provide a hierarchy of models used to reduce the numerical cost of these simulations while preserving basic properties of the solutions. Moment models require a closure because they have more unknowns than equations. The entropy-based closure is based on the physical description of the particle interactions and provides desirable properties. However, computing this closure is expensive. We propose an approximation of the closure for the first two models in the hierarchy, the M₁ and M₂ models valid in one, two or three dimensions of space. Compared to other approximate closures, our method works in multiple dimensions. We obtain the approximation by a careful study of the domain of realizability and by invariance properties of the entropy minimizer. The M₂ model is shown to provide significantly better accuracy than the M₁ model for the numerical simulation of a dose computation in radiotherapy. We propose a numerical solver using those approximated closures. Numerical experiments in dose computation test cases show that the new method is more efficient compared to numerical solution of the minimum entropy problem using standard software tools.