An adaptive moving mesh discontinuous Galerkin method for the radiative transfer equation
For researchers solving the high-dimensional radiative transfer equation, this method offers a mesh adaptation strategy to handle discontinuities and sharp layers, but the results are incremental as they only show performance on selected examples without SOTA comparisons.
The paper develops an adaptive moving mesh discontinuous Galerkin method for the radiative transfer equation, combining dynamic mesh adaptation with discrete ordinate angular discretization. Numerical examples in 1D and 2D demonstrate improved accuracy and efficiency.
The radiative transfer equation models the interaction of radiation with scattering and absorbing media and has important applications in various fields in science and engineering. It is an integro-differential equation involving time, space and angular variables and contains an integral term in angular directions while being hyperbolic in space. The challenges for its numerical solution include the needs to handle with its high dimensionality, the presence of the integral term, and the development of discontinuities and sharp layers in its solution along spatial directions. Its numerical solution is studied in this paper using an adaptive moving mesh discontinuous Galerkin method for spatial discretization together with the discrete ordinate method for angular discretization. The former employs a dynamic mesh adaptation strategy based on moving mesh partial differential equations to improve computational accuracy and efficiency. Its mesh adaptation ability, accuracy, and efficiency are demonstrated in a selection of one- and two-dimensional numerical examples.