• 论文 •

任意四边形网格上扩散问题的一个稳定九点格式

1. 中国工程物理研究院研究生院, 北京 100088
• 收稿日期:2018-10-29 出版日期:2019-03-15 发布日期:2019-03-15
• 基金资助:

国家自然科学基金（11871009，11771052）.

Hong Qi, Su Shuai. A STABLE NINE-POINT SCHEME FOR DIFFUSION PROBLEMS ON ARBITRARY QUADRILATERAL MESHES[J]. Journal of Numerical Methods and Computer Applications, 2019, 40(1): 51-67.

A STABLE NINE-POINT SCHEME FOR DIFFUSION PROBLEMS ON ARBITRARY QUADRILATERAL MESHES

Hong Qi, Su Shuai

1. Graduate School of China Academy of Engineering Physics, Beijing 100088, China
• Received:2018-10-29 Online:2019-03-15 Published:2019-03-15

We propose and analyze a stable nine-point scheme for diffusion problems by applying a modified Q1 finite volume element method (mQ1-FVEM) to treat the vertex unknowns. Based on the theoretical results of the mQ1-FVEM and a simple discrete functional technique, the stability result and error estimate of the resulting nine-point scheme both in H1 norm are obtained under a standard and weak geometric assumption. Our proposed scheme does not possess the so-called numerical heat-barrier issue suffered by many existing cell-centered or hybrid schemes. But, we need to solve two linear algebraic systems. Numerical experiments are also presented to show the efficiency of the proposed scheme and validate the theoretical analysis.

MR(2010)主题分类:

()
 [1] Morel J, Roberts R, Shashkov M. A local support-operators diffusion discretization scheme for quadrilateral r-z meshes[J]. J. Comput. Phys., 1998, 144(1):17-51.[2] Shashkov M, Steinberg S, Support-operator finite-difference algorithms for general elliptic problems[J]. J. Comput. Phys., 1995, 118(1):131-151.[3] Aavatsmark I. An introduction to multipoint flux approximations for quadrilateral grids[J]. Comput. Geosci., 2002, 6(3):405-432.[4] Aavatsmark I, Barkve T, Boe O, Mannseth T. Discretization on unstructured grids. Part I:Derivation of the methods for inhomogeneous, anisotropic media[J] SIAM J. Sci. Comput., 1998, 19(5):1700-1716.[5] Domelevo K, Omnes P. A finite volume method for the Laplace equation on almost arbitrary two-dimensional grids[J]. ESAIM:M2AN, 2005, 39(6):1203-1249.[6] Hermeline F. A finite volume method for the approximation of diffusion operators on distorted meshes[J]. J. Comput. Phys., 2000, 160(2):481-499.[7] Li D, Shui H, Tang M. On the finite difference scheme of two-dimensional parabolic equation in a non-rectangular mesh[J]. J. Numer. Methods Comput. Appl., 1980, 4(1):217-224.[8] Lipnikov K, Shashkov M, Svyatskiy D. Monotone finite volume schemes for diffusion equations on unstructured triangular and shape-regular polygonal meshes[J]. J. Comput. Phys., 2007, 227(1):492-512.[9] Coudière Y, Vila J, Villedieu P. Convergence rate of a finite volume scheme for a two dimensional convection-diffusion problem[J]. Math. Model. Numer. Anal., 1999, 33(3):493-516.[10] Gao Z, Wu J. A linearity-preserving cell-centered scheme for the heterogeneous and anisotropic diffusion equations on general meshes[J] Int. J. Numer. Meth. Fluids, 2011, 67(12):2157-2183.[11] Sheng Z, Yuan G. A nine point scheme for the approximation of diffusion operators on distored quadrilateral meshes[J]. SIAM J. Sci. Comput., 2008, 30(3):1341-1361.[12] Hong Q, Wu J. Coercivity results of a modified Q1-finite volume element scheme for anisotropic diffusion problems[J]. Adv. Comput. Math., 2018, 44(3):897-922.[13] Pei W. The construction of simulation algorithms for laser fusion[J]. Commun. Comput. Phys., 2007, 2(2):255-270.[14] Ding L, Wu J, Yang Z, Dai Z. Simulation of Z-pinch implosion using MARED code[J]. High Power Laser Particle Beans, 2008, 20(2):212-218.[15] Wu J, Dai Z, Gao Z, Yuan G. Linearity preserving nine-point schemes for diffusion equation on distorted quadrilateral meshes[J]. J. Comput. Phys., 2010, 229(9):3382-3401.[16] Zel'dovich Y, Raizer Y. Benchmark on discretization schemes for anisotropic diffusion problems on general grids[M]. New York:Academic Press, 1967.
 [1] 李辰, 郭启龙, 孙东, 刘朋欣. 一种求解双曲守恒律方程的中心型WENO格式[J]. 数值计算与计算机应用, 2020, 41(3): 246-258. [2] 杨洋, 曾维新, 高振, 王保山. 求解双曲守恒律的紧-WENO杂交格式及RBF-FD间断检测方法[J]. 数值计算与计算机应用, 2020, 41(3): 232-245. [3] 刘新龙, 杨晓忠. 时间分数阶四阶扩散方程的显-隐和隐-显差分格式[J]. 数值计算与计算机应用, 2020, 41(3): 216-231. [4] 徐捷, 高振, 曾维新, 王保山. 高阶驻点上精度保持的六阶WENO有限差分格式[J]. 数值计算与计算机应用, 2020, 41(1): 68-82. [5] 赵心仪, 董明哲. 一类非线性抛物型方程的紧差分格式[J]. 数值计算与计算机应用, 2019, 40(3): 188-206. [6] 尹旭, 卢朓, 姜海燕. 数值求解含时Wigner方程的一种高阶算法[J]. 数值计算与计算机应用, 2019, 40(1): 21-33. [7] 澈力木格, 何斯日古楞, 李宏. 大气污染模型的POD基降维有限差分算法[J]. 数值计算与计算机应用, 2018, 39(3): 172-182. [8] 郭峰. 广义Zakharov-Kuznetsov方程的若干能量守恒算法[J]. 数值计算与计算机应用, 2018, 39(3): 183-194. [9] 吴立飞, 杨晓忠, 张帆. 非线性Leland方程的一种并行本性差分方法[J]. 数值计算与计算机应用, 2014, 35(1): 69-80. [10] 钟巍, 田宙. 五阶精度WENO差分型格子波尔兹曼算法求解单守恒方程[J]. 数值计算与计算机应用, 2012, 33(4): 242-250. [11] 张琎, 王震. 双曲方程基于分离算子的交替方向差分算法[J]. 数值计算与计算机应用, 2011, 32(3): 159-164. [12] 张在斌, 孙志忠. 一类非线性延迟抛物偏微分方程的Crank-Nicolson型差分格式[J]. 数值计算与计算机应用, 2010, 31(2): 131-140. [13] 洪振英,袁光伟,傅学东. 动态中子输运方程迭代初值的选取方法[J]. 数值计算与计算机应用, 2008, 29(4): 302-312.