@Article{NMTMA-13-220, author = {Su , Shuai and Wu , Jiming}, title = {A Vertex-Centered and Positivity-Preserving Finite Volume Scheme for Two-Dimensional Three-Temperature Radiation Diffusion Equations on General Polygonal Meshes}, journal = {Numerical Mathematics: Theory, Methods and Applications}, year = {2019}, volume = {13}, number = {1}, pages = {220--252}, abstract = {

Two-dimensional three-temperature (2-D 3-T) radiation diffusion equations are widely used to approximately describe the evolution of radiation energy within a multimaterial system and explain the exchange of energy among electrons, ions and photons. In this paper, we suggest a new positivity-preserving finite volume scheme for 2-D 3-T radiation diffusion equations  on general polygonal meshes. The vertex unknowns are treated as primary ones  for which the finite volume equations are constructed. The edge-midpoint and cell-centered unknowns are used as auxiliary ones and interpolated by the primary unknowns, which makes the final scheme a pure vertex-centered one. By comparison, most existing positivity-preserving finite volume schemes are cell-centered and based on the convex decomposition of the co-normal. Here, the co-normal decomposition is  not convex in general, leading to  a fixed stencil of  the flux  approximation  and avoiding a certain search algorithm on complex grids. Moreover, the new scheme effectively alleviates the numerical heat-barrier issue suffered by most existing cell-centered or hybrid schemes in solving strongly nonlinear radiation diffusion equations. Numerical experiments demonstrate the second-order accuracy and the positivity of the solution on various distorted grids. For the problem without analytic solution, the contours of the numerical solutions obtained by our scheme on distorted meshes accord with those on smooth quadrilateral meshes.

}, issn = {2079-7338}, doi = {https://doi.org/10.4208/nmtma.OA-2018-0121}, url = {http://global-sci.org/intro/article_detail/nmtma/13438.html} }