High Order Positivity-Preserving Nodal Discontinuous Galerkin Methods for Anisotropic Diffusion Problems and Their Applications


Xinyuan Liu, Xiamen University


In this talk, we consider second and third order accurate positivity-preserving (PP) nodal discontinuous Galerkin (DG) methods for one and two dimensional anisotropic diffusion problems. The key idea is to first represent the cell average of its numerical approximation as a weighted summation of Gaussian quadrature point values used in the updating of nodal DG methods, and then transform these Gaussian quadrature point values to some other special chosen point values. We prove that by taking parameters in the definition of nodal DG methods appropriately, together with a suitable time stability condition, the cell averages can be kept positive. A polynomial scaling limiter is then applied to obtain positive numerical approximations on the whole cell without sacrificing accuracy. Stability analysis without the PP limiter is also rigorously established. We then develop a provable third order bound-preserving (BP) nodal DG method for multi-component compressible miscible displacements in porous media by applying the third order PP nodal DG method. Numerical experiments are given to demonstrate desired accuracy, PP (or BP) and good performances of our proposed approach. This is a joint work with Prof. Tao Xiong and Prof. Yang Yang.