A Positivity-Preserving Finite Volume Scheme with Least Square Interpolation for 3D Anisotropic Diffusion Equation

被引：0
作者：
Hui Xie
Xuejun Xu
Chuanlei Zhai
Heng Yong
机构：
[1] Institute of Applied Physics and Computational Mathematics,LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences
[2] Chinese Academy of Sciences,School of Mathematical Sciences
[3] Tongji University,undefined
来源：
Journal of Scientific Computing | 2021年 / 89卷
关键词：
Anisotropic diffusion; Positivity; Vertex unknowns; Least square;
D O I：
暂无
中图分类号：
学科分类号：
摘要：
In this paper, a new nonlinear finite volume scheme preserving positivity for 3D anisotropic diffusion equation is proposed. A distinct feature of our new scheme is that the auxiliary vertex unknowns appearing in the nonlinear two-point flux approximation are interpolated by the least-square method combined with the so-called correction function. The correction function is local, and only solved when the coefficient discontinuity occurs. Compared with other existing 3D interpolation formulas, our interpolation formula is designed for the general polyhedron mesh and coefficient discontinuity with lower and adaptive computational complexity. The scheme is proven to be positivity-preserving. Numerical examples are presented to demonstrate the second-order accuracy and positivity-preserving property for various anisotropic diffusion problems on the distorted meshes.
引用
收藏
相关论文
共 138 条
  • [1] Aavatsmark I(2006)Numerical convergence of the MPFA O-method and U-method for general quadrilateral grids Int. J. Numer. Meth. Fluids 51 939-961
  • [2] Eigestad GT(2009)A nine-point finite volume scheme for the simulation of diffusion in heterogeneous media C. R. Acad. Sci. Paris Ser. I 673-676
  • [3] Agelas L(2012)Monotonicity recovering and accuracy preserving optimization methods for postprocessing finite element solutions J. Comput. Phys. 231 3126-3142
  • [4] Eymard R(2017)Variational inequality approach to enforcing the non-negative constraint for advection-diffusion equations Comput. Methods Appl. Mech. Eng. 320 287-334
  • [5] Herbin R(2013)Monotone corrections for generic cell-centered finite volume approximations of anisotropic diffusion equations Numer. Math. 125 387-417
  • [6] Burdakov O(2003)Difference scheme by integral interpolation method for three dimensional diffusion equations Chin. J. Comput. Phys. 20 205-209
  • [7] Kapyrin I(1973)Maximum principle and uniform convergence for the finite element method Comput. Methods Appl. Mech. Eng. 2 17-31
  • [8] Vassilevski Y(2009)A monotone nonlinear finite volume method for diffusion equations on conformal polyhedral meshes Russ. J. Numer. Anal. Math. Modell. 24 207-227
  • [9] Chang J(2004)Failure of the discrete maximum pricinple for an elliptic finite element problems Math. Comput. 74 1-23
  • [10] Nakshatrala KB(2014)Finite volume schemes for diffusion equations: introduction to and review of modern methods Math. Models Methods Appl. Sci. (M3AS) 24 1575-1619
12345678910