对流扩散方程的非振荡守恒特征差分法
Non Oscillation Conservative Characteristic Difference Method for Solving Convection Diffusion Equations
摘要: 当求解对流占优扩散问题时,采用传统的二次拉格朗日插值的特征差分法,会出现数值振荡,且不满足质量守恒。结合算子分裂,本质非震荡和MMOCAA质量校正,提出了求解对流扩散方程的非震荡的守恒特征差分法。首先采用局部一维(LOD)法把一个二维偏微分方程分裂成x方向和y方向的两个一维的偏微分方程组;其次在每个方向上利用二阶本质非振荡和MMOCAA格式进行数值计算。数值实验验证格式满足非震荡和质量守恒,能够有效地解决大型对流占优扩散问题。
Abstract: The characteristic difference method based on quadratic Lagrange interpolation is used to solve the convection dominated diffusion problem, but there will be a larger numerical oscillation and not conservative. Combining the operator splitting, non-oscillatory and mass correction methods, the non-oscillation conservative characteristic difference methods are proposed to solve the convection diffusion equations. Firstly, the partial differential equations in two dimensions are splitting into two one-dimensional partial differential equations along the x-direction and the y-direction, respectively. Secondly, the second-order essentially non-oscillation and MMOCAA schemes are presented to compute the equations. By the numerical results, it shows that the scheme not only meets non-oscillatory and mass conservation, but also effectively solves the convection-dominant diffusion problems.
文章引用:王钱钱, 李琳, 赵玉庆, 周忠国. 对流扩散方程的非振荡守恒特征差分法[J]. 应用数学进展, 2018, 7(11): 1446-1457. https://doi.org/10.12677/AAM.2018.711169

参考文献

[1] 童小红, 王兴, 苏李君, 等. 对流扩散方程的特征线EFG算法[J]. 工程数学学报, 2018, 35(2): 179-192.
[2] Douglas, J. and Thomas, F. (1982) Numerical Methods for Convention Dominated Diffusion Problems Based on Combining the Method of Characteristics with Finite Element or Finite Difference Procedures. SIAM Journal on Mathematic Analysis, 10, 871-885.
[3] Gabriel, R., Burman, E. and Karakatsani, F. (2017) Blending Low-Order Stabilised Finite Element Methods: A Positivity Preserving Local Projection Method for the Convection-Diffusion Equation. Computer Methods in Applied Mechanics and Engineering, 139, 23-26.
[4] 梁栋. 对流扩散方程的一类迎风格式[J]. 计算数学, 1991(2): 133-141.
[5] 袁益让, 梁栋, 芮洪兴, 海水入侵及防治工程的数值模拟[J]. 计算物理, 2001, 18(6): 556-562.
[6] Zhou, Z. and Liang, D. (2017) The Mass-Preserving and Modified-Upwind Splitting DDM Scheme for Time Dependent Convection Diffusion Equations. Journal of Computational and Applied Mathematics, 317, 247-273. [Google Scholar] [CrossRef
[7] 陆金普, 刘晓遇, 杜正平. 对流占优扩散问题的一种特征差分方法[J]. 清华大学学报, 2002, 42(8): 111-113.
[8] 由同顺. 对流扩散方程的本质非振荡特征差分方法[J]. 应用数学, 2000, 21(4): 89-94.
[9] 黄素珍, 张鲁明. 对流扩散方程的一种高精度特征差分格式[J]. 南京师大学报, 2005, 28(2): 38-41.
[10] Douglas, J., Huang, C. and Pereira, F. (1999) The Modified Method of Characteristics with Adjust Advection. Numerische Mathematik, 83, 353-369. [Google Scholar] [CrossRef
[11] Rui, H. and Tabata, M. (2010) A Mass-Conservative Characteristic Finite Element Scheme for Convection-Diffusion Problem. Journal of Scientific Computing, 43, 416-432. [Google Scholar] [CrossRef
[12] Fu, K. and Liang, D. (2016) The Conservative Characteristic FD Methods for Atmospheric Aerosol Transport Problems. Journal of Computational Physics, 305, 494-520. [Google Scholar] [CrossRef
[13] Zhou, Z. and Liang, D. (2016) The Mass-Preserving S-DDM Scheme for Two-Dimensional Parabolic Equations. Communications in Computational Physics, 19, 411-441. [Google Scholar] [CrossRef
[14] Zhou, Z., Liang, D. and Wong, Y. (2018) The New Mass-Conserving S-DDM Scheme for Two-Dimensional Parabolic Equations with Variable Coefficients. Applied Mathematics and Computation, 338, 882-902. [Google Scholar] [CrossRef
[15] 柳忠伟. 求解高维对流扩散方程的高效本质非振荡特征差分方法[D]: [学士学位论文]. 泰安: 山东农业大学, 2017.
[16] 李清扬, 王能超, 易大义. 数值分析[M]. 第5版. 北京: 清华大学出版社, 2008.