基于二阶格式的有限体积保正格式

计算物理 ›› 2020, Vol. 37 ›› Issue (4): 379-392.

• • 下一篇

基于二阶格式的有限体积保正格式

赵菲^1,2, 盛志强², 袁光伟²

1. 北方工业大学理学院, 北京 100144;
2. 北京应用物理与计算数学研究所计算物理实验室, 北京 100088

收稿日期:2019-05-16 修回日期:2019-08-06 出版日期:2020-07-25 发布日期:2020-07-25
通讯作者: 袁光伟,E-mail:yuan_guangwei@iapcm.ac.cn
作者简介:赵菲(1991-),博士生,研究方向为偏微分方程数值解,E-mail:zhaofei2020@ncut.edu.cn
基金资助:
国家自然科学基金（11971069，11571048，11871112），NSAF（U1630249）和科学挑战专题（TZ2016002）资助项目

A Positivity-preserving Finite Volume Scheme Based on Second-order Scheme

ZHAO Fei^1,2, SHENG Zhiqiang², YUAN Guangwei²

1. College of Science, North China Univevsity of Technology, Beijing 100144, China;
2. Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China

Received:2019-05-16 Revised:2019-08-06 Online:2020-07-25 Published:2020-07-25

摘要/Abstract

摘要： 从二阶线性格式出发，通过对法向通量进行重构，得到非线性两点通量，获得四面体网格上的单元中心型有限体积保正格式.该格式适用于求解间断和各向异性扩散系数问题.无需假设辅助未知量非负，避免了辅助未知量计算出负时"遇负置零"的人为处理方式；并且证明该格式在每个非线性Picard迭代步具有强保正性，即当源项和边界条件非负时，线性化格式的非平凡解是严格大于零的.数值算例验证该格式具有二阶收敛性且是保正的.

关键词: 扩散方程, 四面体网格, 有限体积格式, 强保正性

Abstract: Based on a second-order accurate linear scheme, a normal flux is reconstructed to obtain a nonlinear finite volume scheme with a two-point flux discrete stencil on tetrahedral meshes. It is suitable for discontinuous and anisotropic diffusion coefficient problems, and can be generalized to general polyhedral meshes. It is unnecessary to assume that auxiliary unknowns are non-negative, and avoids artificial processing of "setting negative to be zero" in calculating auxiliary unknowns. Moreover, it is proved that the linearized scheme at each nonlinear iteration step satisfies strong positivity-preserving, i.e., as the source term and boundary condition are non-negative, non-zero solution of the scheme is strictly greater than zero. Numerical tests verify that the scheme has second-order accuracy and is strong positivity-preserving.

Key words: diffusion equations, tetrahedral meshes, finite volume scheme, strong positivity-preserving

中图分类号:

O241.82

赵菲, 盛志强, 袁光伟. 基于二阶格式的有限体积保正格式[J]. 计算物理, 2020, 37(4): 379-392.

ZHAO Fei, SHENG Zhiqiang, YUAN Guangwei. A Positivity-preserving Finite Volume Scheme Based on Second-order Scheme[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2020, 37(4): 379-392.

参考文献

[1] YU C X, FAN Z F, LIU J, et al. Modeling of shell-mixing into central hotspot in inertial confinement fusion implosion[J]. Chinese J Comput Phys, 2017, 34(4):379-386.
[2] SHEN Y S, LI K B, SHI X M, et al. Study on disposing high-level transuranic waste aste in a fusion fission reactor[J]. Chinese J Comput Phys, 2017, 34(2):142-148.
[3] YUAN G W, HANG X D, SHENG Z Q, et al. Progress in numerical methods for radiation diffusion equations[J]. Chinese J Comput Phys, 2009, 26(4):475-500.
[4] LI D Y, SHUI H S, TANG M J. On the finite difference scheme of two-dimensional parabolic equation in a non-rectangular mesh[J]. J Numer Methods Comput Appl, 1980, 1:217-224.
[5] SHASHKOV M, STEINBERG S. Support-operator finite-difference algorithms for general elliptic problems[J]. J Comput Phys, 1995, 118:131-151.
[6] AAVATSMARK I. An introduction to multipoint flux approximations for quadrilateral grids[J]. Comput Geosci, 2002, 6:405-432.
[7] CAO F J, YAO Y Z. Conservative positivity-preserving algorithm for Kershaw scheme of anisotropic diffusion problem[J]. Chinese J Comput Phys, 2017, 34(3):283-293.
[8] YAO Y Z, YUAN G W. Enforcing positivity with conservation for nine-point scheme of nonlinear diffusion equations[J]. Comput Methods Appl Mech Eng, 2012, 223-224(1):161-172.
[9] NORDBOTTEN J, AAVATSMARK I, EIGESTAD G. Monotonicity of control volume methods[J]. Numer Math, 2007, 106(2):255-288.
[10] LE POTIER C. Finite volume monotone scheme for highly anisotropic diffusion operators on unstructured triangular meshes[J]. C R Math Acad Sci Paris, 2005, 341(12):787-792.
[11] LIPNIKOV K, SHASHKOV M, SVYATSKIY D, et al. Monotone finite volume schemes for diffusion equations on unstructured triangular and shape-regular polygonal meshes[J]. J Comput Phys, 2007, 227(1):492-512.
[12] LIPNIKOV K, SVYATSKIY D, VASSILEVSKI Y. Interpolation-free monotone finite volume method for diffusion equations on polygonal meshes[J]. J Comput Phys, 2009, 228(3):703-716.
[13] LIPNIKOV K, SVYATSKIY D, VASSILEVSKI Y. A monotone finite volume method for advection-diffusion equations on unstructured polygonal meshes[J]. J Comput Phys, 2010, 229(11):4017-4032.
[14] YUAN G W, SHENG Z Q. Monotone finite volume schemes for diffusion equations on polygonal meshes[J]. J Comput Phys, 2008, 227(12):6288-6312.
[15] SHENG Z Q, YUE J Y, YUAN G W. Monotone finite volume schemes of nonequilibrium radiation diffusion equations on distorted meshes[J]. SIAM J Sci Comput, 2009, 31(4):2915-2934.
[16] SHENG Z Q, YUAN G W. An improved monotone finite volume scheme for diffusion equation on polygonal meshes[J]. J Comput Phys, 2012, 231(9):3739-3754.
[17] SHENG Z Q, YUAN G W. A new nonlinear finite volume scheme preserving positivity for diffusion equations[J]. J Comput Phys, 2016, 315:182-193.
[18] DANILOV A, VASSILEVSKI Y. A monotone nonlinear finite volume method for diffusion equations on conformal polyhedral meshes[J]. Russ J Numer Anal Math Modelling, 2009, 24(3):207-227.
[19] NIKITIN K, VASSILEVSKI Y. A monotone nonlinear finite volume method for advection-diffusion equations on unstructured polyhedral meshes in 3D[J]. Russ J Numer Anal Math Modelling, 2010, 25(4):335-358.
[20] HERMELINE F. Approximation of 2-D and 3-D diffusion operators with variable full tensor coefficients on arbitrary meshes[J]. Comput Methods Appl Mech Eng, 2007, 196(21):2497-2526.
[21] HERMELINE F. A finite volume method for approximating 3D diffusion operators on general meshes[J]. J Comput Phys, 2009, 228(16):5763-5786.
[22] LAI X, SHENG Z Q, YUAN G W. A finite volume scheme for three-dimensional diffusion equations[J]. Commun Comput Phys, 2015, 18(3):650-672.
[23] LAI X, SHENG Z Q, YUAN G W. Monotone finite volume scheme for three dimensional diffusion equation on tetrahedral meshes[J]. Commun Comput Phys, 2017, 21(1):162-181.
[24] WANG S, HANG X D, YUAN G W. A positivity-preserving finite volume scheme for diffusion equations on polyhedral meshes[J]. Math Numer Sin, 2015, 37(3):247-263.
[25] WANG S, HANG X D, YUAN G W. A pyramid scheme for three-dimensional diffusion equations on polyhedral meshes[J]. J Comput Phys, 2017, 350:590-606.
[26] WANG S, HANG X D, YUAN G W. A positivity-preserving pyramid scheme for anisotropic diffusion problems on general hexahedral meshes with nonplanar cell faces[J]. J Comput Phys, 2018, 371:152-167.
[27] LAN B, SHENG Z Q, YUAN G W. A new finite volume scheme preserving positivity for radionuclide transport calculations in radioactive waste repository[J]. Int J Heat Mass Tran, 2018, 121:736-746.
[28] GAO Y N, YUAN G W, WANG S, et al. A finite volume element scheme with a monotonicity correction for anisotropic diffusion problems on general quadrilateral meshes[J]. J Comput Phys, 2020, 407:109-143.
[29] 袁光伟, 等. 扩散方程计算方法[M]. 北京:科学出版社, 2015.

基于二阶格式的有限体积保正格式

A Positivity-preserving Finite Volume Scheme Based on Second-order Scheme

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	冯舒婷, 戴厚平, 宋通政. 二维空间分数阶反应扩散方程组的格子Boltzmann方法[J]. 计算物理, 2022, 39(6): 666-676.
[2]	胡晓燕, 范征锋. 惯性约束聚变内爆中基于多块结构网格的高效辐射扩散并行算法[J]. 计算物理, 2022, 39(3): 277-285.
[3]	张淑娴, 杭旭登. 二维扩散方程的局部二阶线性保极值节点计算方法及应用[J]. 计算物理, 2022, 39(1): 17-32.
[4]	魏雪丹, 戴厚平, 李梦军, 郑洲顺. 一维空间Riesz分数阶对流扩散方程的格子Boltzmann方法[J]. 计算物理, 2021, 38(6): 683-692.
[5]	徐金景, 袁光伟. 多流管网格上九点格式的节点插值方法[J]. 计算物理, 2021, 38(2): 153-164.
[6]	赵丰祥, 潘亮, 王双虎. 基于非结构四边形网格的WENO格式[J]. 计算物理, 2018, 35(5): 525-534.
[7]	朱晓钢, 聂玉峰, 王俊刚, 袁占斌. 分数阶对流扩散方程的特征有限元方法[J]. 计算物理, 2017, 34(4): 417-424.
[8]	郭少冬, 章明宇, 周海兵, 熊俊, 张树道. 三维非匹配网格上的扩散方程数值求解[J]. 计算物理, 2017, 34(1): 19-28.
[9]	吕桂霞, 孙顺凯. 扩散方程的有限方向差分方法[J]. 计算物理, 2015, 32(6): 649-661.
[10]	舒适, 岳孝强, 周志阳, 徐小文. SAMR网格上扩散方程有限体格式的逼近性与两层网格算法[J]. 计算物理, 2014, 31(4): 390-402.
[11]	邓志红, 孙玉良, 李富, Rizwan-uddin. 改进的节块展开方法求解圆柱几何对流扩散方程[J]. 计算物理, 2013, 30(4): 475-482.
[12]	曾凡海, 李常品. 时间分数阶亚扩散方程的高阶差分方法[J]. 计算物理, 2013, 30(4): 491-500.
[13]	杭旭登, 李敬宏, 袁光伟. 多群辐射扩散方程组的分裂迭代算法收敛分析[J]. 计算物理, 2013, 30(1): 111-119.
[14]	张荣培, 蔚喜军, 崔霞, 冯涛. 间断有限元方法求解-维非平衡辐射扩散方程[J]. 计算物理, 2012, 29(5): 641-646.
[15]	张荣培, 蔚喜军, 崔霞, 冯涛. 隐-显积分因子间断Galerkin方法求解二维辐射扩散方程[J]. 计算物理, 2012, 29(5): 647-653.