A deterministic solver for the linear Boltzmann model of a single mono-directional proton beam

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Zhang, Xiaojiang, Bai, Xuemin, Tang, Min
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866908292672913408
author Zhang, Xiaojiang
Bai, Xuemin
Tang, Min
author_facet Zhang, Xiaojiang
Bai, Xuemin
Tang, Min
contents 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.
format Preprint
id arxiv_https___arxiv_org_abs_2504_00340
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle A deterministic solver for the linear Boltzmann model of a single mono-directional proton beam
Zhang, Xiaojiang
Bai, Xuemin
Tang, Min
Numerical Analysis
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.
title A deterministic solver for the linear Boltzmann model of a single mono-directional proton beam
topic Numerical Analysis
url https://arxiv.org/abs/2504.00340