Gespeichert in:
| Hauptverfasser: | , , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2025
|
| Schlagworte: | |
| Online-Zugang: | https://arxiv.org/abs/2509.07689 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866916942805204992 |
|---|---|
| author | Moujaes, Paul Kuzmin, Dmitri Bäumer, Christian |
| author_facet | Moujaes, Paul Kuzmin, Dmitri Bäumer, Christian |
| contents | We discretize the $M_1$ model of radiative transfer using continuous finite elements and propose a tailor-made monolithic convex limiting (MCL) procedure for enforcing physical realizability. The $M_1$ system of nonlinear balance laws for the zeroth and first moments of a probability distribution function is derived from the linear Boltzmann equation and equipped with an entropy-based closure for the second moment. To ensure hyperbolicity and physical admissibility, evolving moments must stay in an invariant domain representing a convex set of realizable states. We first construct a low-order method that is provably invariant domain preserving (IDP). Introducing intermediate states that represent spatially averaged exact solutions of homogeneous Riemann problems, we prove that these so-called bar states are realizable in any number of space dimensions. This key auxiliary result enables us to show the IDP property of a fully discrete scheme with a diagonally implicit treatment of reactive terms. To achieve high resolution, we add nonlinear correction terms that are constrained using a two-step MCL algorithm. In the first limiting step, local bounds are imposed on each conserved variable to avoid spurious oscillations and maintain positivity of the scalar-valued zeroth moment (particle density). The second limiting step constrains the magnitude of the vector-valued first moment to be realizable. The flux-corrected finite element scheme is provably IDP. Its ability to prevent nonphysical behavior while attaining high-order accuracy in smooth regions is verified in a series of numerical tests. The developed methodology provides a robust simulation tool for dose calculation in radiotherapy. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2509_07689 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Realizability-preserving monolithic convex limiting in continuous Galerkin discretizations of the M1 model of radiative transfer Moujaes, Paul Kuzmin, Dmitri Bäumer, Christian Numerical Analysis We discretize the $M_1$ model of radiative transfer using continuous finite elements and propose a tailor-made monolithic convex limiting (MCL) procedure for enforcing physical realizability. The $M_1$ system of nonlinear balance laws for the zeroth and first moments of a probability distribution function is derived from the linear Boltzmann equation and equipped with an entropy-based closure for the second moment. To ensure hyperbolicity and physical admissibility, evolving moments must stay in an invariant domain representing a convex set of realizable states. We first construct a low-order method that is provably invariant domain preserving (IDP). Introducing intermediate states that represent spatially averaged exact solutions of homogeneous Riemann problems, we prove that these so-called bar states are realizable in any number of space dimensions. This key auxiliary result enables us to show the IDP property of a fully discrete scheme with a diagonally implicit treatment of reactive terms. To achieve high resolution, we add nonlinear correction terms that are constrained using a two-step MCL algorithm. In the first limiting step, local bounds are imposed on each conserved variable to avoid spurious oscillations and maintain positivity of the scalar-valued zeroth moment (particle density). The second limiting step constrains the magnitude of the vector-valued first moment to be realizable. The flux-corrected finite element scheme is provably IDP. Its ability to prevent nonphysical behavior while attaining high-order accuracy in smooth regions is verified in a series of numerical tests. The developed methodology provides a robust simulation tool for dose calculation in radiotherapy. |
| title | Realizability-preserving monolithic convex limiting in continuous Galerkin discretizations of the M1 model of radiative transfer |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2509.07689 |