Discontinuous Galerkin FEMs for Radiation Transport Problems
The linear Boltzmann transport equation (LBTE), a high dimensional partial integrodifferential equation, is used to model radiation transport. Radiation transport is an area of physics that is concerned with the propagation and distribution of radiative particles, such as photons and electrons withi...
| Main Author: | |
|---|---|
| Format: | Thesis (University of Nottingham only) |
| Language: | English |
| Published: |
2023
|
| Subjects: | |
| Online Access: | https://eprints.nottingham.ac.uk/76685/ |
| _version_ | 1848800927152078848 |
|---|---|
| author | Widdowson, Richard |
| author_facet | Widdowson, Richard |
| author_sort | Widdowson, Richard |
| building | Nottingham Research Data Repository |
| collection | Online Access |
| description | The linear Boltzmann transport equation (LBTE), a high dimensional partial integrodifferential equation, is used to model radiation transport. Radiation transport is an area of physics that is concerned with the propagation and distribution of radiative particles, such as photons and electrons within a material medium.
In this thesis, we present a high-order discontinuous Galerkin finite element method (DGFEM) discretisations of the steady state linear Boltzmann transport equation in the spatial, angular, and energy domains. Comparisons between the tensor discontinuous Galerkin finite element method and discrete ordinates method show that the former is higher order than the latter. A method of block diagonalising the resulting matrix into a sequence of transport equations coupled by the right hand side while retaining high order convergence, is demonstrated for both the angular and energy domains. This new method offers the arbitrary order convergence rates of the discontinuous Galerkin finite element method, but it can be implemented in an almost identical form to standard multigroup discrete ordinates methods.
The assembly of the matrix for the resulting transport equations for a variety of different type of elements is discussed. The generation of meshes formed of general polytopes is discussed, and a comparison between the time to solve the transport equation on meshes formed of the different element types follows.
An efficient implementation of the discontinuous Galerkin finite element method for transport equations is then presented. This algorithm exploits the fixed wind direction of the transport problems resulting from the discontinuous Galerkin finite element method of the LBTE, to solve the transport problem while never forming the matrix. We then compare this algorithm to a direct matrix solver for both convex and non-convex polytopes.
An a posteriori error bound for the discontinuous Galerkin finite element method of the LBTE and the transport problem is then derived. We use this error bound to develop an adaptive framework for the LBTE. Three different adaptive algorithms for the LBTE are then presented and compared. The h-refinement algorithm, which marks the element with the error of each tensor element it is part of, shows a clear advantage over the other methods. |
| first_indexed | 2025-11-14T20:59:20Z |
| format | Thesis (University of Nottingham only) |
| id | nottingham-76685 |
| institution | University of Nottingham Malaysia Campus |
| institution_category | Local University |
| language | English |
| last_indexed | 2025-11-14T20:59:20Z |
| publishDate | 2023 |
| recordtype | eprints |
| repository_type | Digital Repository |
| spelling | nottingham-766852023-12-13T16:12:30Z https://eprints.nottingham.ac.uk/76685/ Discontinuous Galerkin FEMs for Radiation Transport Problems Widdowson, Richard The linear Boltzmann transport equation (LBTE), a high dimensional partial integrodifferential equation, is used to model radiation transport. Radiation transport is an area of physics that is concerned with the propagation and distribution of radiative particles, such as photons and electrons within a material medium. In this thesis, we present a high-order discontinuous Galerkin finite element method (DGFEM) discretisations of the steady state linear Boltzmann transport equation in the spatial, angular, and energy domains. Comparisons between the tensor discontinuous Galerkin finite element method and discrete ordinates method show that the former is higher order than the latter. A method of block diagonalising the resulting matrix into a sequence of transport equations coupled by the right hand side while retaining high order convergence, is demonstrated for both the angular and energy domains. This new method offers the arbitrary order convergence rates of the discontinuous Galerkin finite element method, but it can be implemented in an almost identical form to standard multigroup discrete ordinates methods. The assembly of the matrix for the resulting transport equations for a variety of different type of elements is discussed. The generation of meshes formed of general polytopes is discussed, and a comparison between the time to solve the transport equation on meshes formed of the different element types follows. An efficient implementation of the discontinuous Galerkin finite element method for transport equations is then presented. This algorithm exploits the fixed wind direction of the transport problems resulting from the discontinuous Galerkin finite element method of the LBTE, to solve the transport problem while never forming the matrix. We then compare this algorithm to a direct matrix solver for both convex and non-convex polytopes. An a posteriori error bound for the discontinuous Galerkin finite element method of the LBTE and the transport problem is then derived. We use this error bound to develop an adaptive framework for the LBTE. Three different adaptive algorithms for the LBTE are then presented and compared. The h-refinement algorithm, which marks the element with the error of each tensor element it is part of, shows a clear advantage over the other methods. 2023-12-12 Thesis (University of Nottingham only) NonPeerReviewed application/pdf en cc_by https://eprints.nottingham.ac.uk/76685/1/RSJWiddowson_14342838_thesis%20%2824%29.pdf Widdowson, Richard (2023) Discontinuous Galerkin FEMs for Radiation Transport Problems. PhD thesis, University of Nottingham. linear Boltzmann transport equations Galerkin finite element methods radiation transport |
| spellingShingle | linear Boltzmann transport equations Galerkin finite element methods radiation transport Widdowson, Richard Discontinuous Galerkin FEMs for Radiation Transport Problems |
| title | Discontinuous Galerkin FEMs for Radiation Transport Problems |
| title_full | Discontinuous Galerkin FEMs for Radiation Transport Problems |
| title_fullStr | Discontinuous Galerkin FEMs for Radiation Transport Problems |
| title_full_unstemmed | Discontinuous Galerkin FEMs for Radiation Transport Problems |
| title_short | Discontinuous Galerkin FEMs for Radiation Transport Problems |
| title_sort | discontinuous galerkin fems for radiation transport problems |
| topic | linear Boltzmann transport equations Galerkin finite element methods radiation transport |
| url | https://eprints.nottingham.ac.uk/76685/ |