High-order conservative positivity-preserving DG-interpolation for deforming meshes and application to moving mesh DG simulation of radiative transfer

Abstract

Solution interpolation between deforming meshes is an important component for several applications in scientific computing, including indirect arbitrary-Lagrangian-Eulerian and rezoning moving mesh methods in numerical solution of partial differential equations. In this paper, a high-order, conservative, and positivity-preserving interpolation scheme is developed based on the discontinuous Galerkin solution of a linear time-dependent equation on deforming meshes. The scheme works for bounded but otherwise arbitrary mesh deformation from the old mesh to the new one. The cost and positivity preservation (with a linear scaling limiter) of the DG-interpolation are investigated. Numerical examples are presented to demonstrate the properties of the interpolation scheme. The DG-interpolation is applied to the rezoning moving mesh DG solution of the radiative transfer equation, an integro-differential equation modeling the conservation of photons and involving time, space, and angular variables. Numerical results obtained for examples in one and two spatial dimensions with various settings show that the resulting rezoning moving mesh DG method maintains the same convergence order as the standard DG method, is more efficient than the method with a fixed uniform mesh, and is able to preserve the positivity of the radiative intensity.

arXiv comment28 pages

References in corpus (5)

• A geometric discretization and a simple implementation for variational mesh generation and adaptation
• On the mesh nonsingularity of the moving mesh PDE method
• A quasi-Lagrangian moving mesh discontinuous Galerkin method for hyperbolic conservation laws
• A comparative numerical study of meshing functionals for variational mesh adaptation
• An adaptive moving mesh discontinuous Galerkin method for the radiative transfer equation

Cited by in corpus (3)

• A high-order well-balanced positivity-preserving moving mesh DG method for the shallow water equations with non-flat bottom topography
• A well-balanced moving mesh discontinuous Galerkin method for the Ripa model on triangular meshes
• A Metric Tensor Approach to Data Assimilation with Adaptive Moving Meshes