Pure Mathematics
Vol.07 No.01(2017), Article ID:19461,9 pages
10.12677/PM.2017.71001

A Compact Finite Difference Scheme for Sobolev Equations

Xin Jing*, Luming Zhang

School of Mathematics, Nanjing University of Aeronautics and Astronautics, Nanjing Jiangsu

Received: Dec. 17th, 2016; accepted: Jan. 2nd, 2017; published: Jan. 5th, 2017

Copyright © 2017 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

ABSTRACT

A compact finite difference scheme is presented for Sobolev equations. It is proved by the discrete energy method that the compact scheme is unconditionally stable and convergent in norm, and the order of convergence is. The numerical experiment results show that the theory is accurate.

Keywords:Sobolev Equations, Compact Finite Difference Scheme, Convergence

Sobolev方程的一个紧致差分格式

经鑫*,张鲁明

南京航空航天大学理学院,江苏 南京

收稿日期:2016年12月17日;录用日期:2017年1月2日;发布日期:2017年1月5日

摘 要

本文对Sobolev方程提出了一个紧致差分格式。并用能量方法证明了该差分格式是以无穷模范数无条件收敛和稳定的,收敛阶为。数值实验结果验证了理论分析的正确性。

关键词 :Sobolev方程,紧致差分格式,收敛性

1. 引言

考虑常系数Sobolev方程的定解问题:

(1)

(2)

(3)

其中为未知函数,为终结时间,为正常数,分别为为源函数、边值与初值,且都是已知函数,

Sobolev方程是一个重要的数学物理方程,被用来描述流体穿过岩石或土壤及不同介质中的流体运动,土壤中湿气的迁移等 [1] [2] [3] 。由于实际问题中Sobolev方程的源函数、初边值或计算域会很复杂,故有效的求解方法是求其数值解。

关于二维常系数Sobolev方程定解问题(1)~(3)的数值求解,有很多这方面的研究。 [4] - [10] 提出了H1-Galerkin混合有限元方法、最小二乘Galerkin有限元方法、时间间断的Galerkin有限元方法、最低阶的半离散和全离散混合有限元方法以及有限体积元方法等来求解。对于求解Sobolev方程定解问题最有效方法之一的有限差分方法,罗振东等 [11] 提出基于POD (奇异值分解和特征投影分解)降阶外推差分算法。先建立全二阶Crank-Nicolson隐格式,再基于POD方法对C-N格式的降阶外推得到新的差分格式。此格式有效的减少了计算过程中的误差积累,提高了计算效率。但这些处理方法的时间及空间精度均不超过2阶。显然,对问题(1)~(3)研究高精度的数值方法有着重要的实际意义,但就我们的阅读所及,还未见这方面的研究。本文对问题(1)~(3)提出了一个紧致差分格式,证明了该差分格式的解以无穷模范数无条件收敛,且收敛阶是。数值实验结果表明理论结果是正确的。在数值计算时,系数矩阵为三对角块矩阵,我们采用块三对角追赶法求解以提高计算效率。

2. 差分格式的建立

2.1. 一些记号

首先对区域上进行网格剖分。引入

符号。其中分别表示空间步长与时间步长,取整数。设指标集

上的网格函数,记,类似地定义网格函数

引入如下记号:

上的网格函数,定义范数:

类似的可以定义:

我们还用到:

2.2. 差分格式

在点处考虑微分方程(1)

(4)

由Taylor展开式

(5)

(6)

(7)

(8)

(9)

(10)

(11)

将(5)~(9)式代入(4)式得到:

(12)

对(12)式左乘算子AB,再利用(10)和(11),因算子A、B可交换,故有:

(13)

即:

(14)

其中

由(2)~(3)初边值条件,我们可以得到:

(15)

(16)

略去小量项,将替换,得到定解问题(1)~(3)的紧致差分格式:

(17)

(18)

(19)

3. 数值解的先验估计

引理1:设上的网格函数,且,则有:

证:

类似的方向也有相同的结论。

引理2 [12] :设上的网格函数,则

引理3 [12] :对于时间序列,任意的。有:

引理4 [12] :假设是一个非负序列,,并且非负序列满足

那么满足

定理1:设为差分方程

的解,其中,那么

(20)

证:在差分方程两边分别与作内积,得:

(21)

其中计算分别如下,

由引理1,有:

(22)

由引理1,有:

(23)

(24)

(25)

将(22)~(25)代入(21),然后不等式两边同乘以,得:

(26)

。则(26)式有:

递推之即得:

由引理3,取,有:

(27)

又因为。由引理1

故:

由上式及(27)式,有

移项,不等式两边再同乘以,有:

又由引理4,有:

4. 差分格式的收敛性与稳定性

定理2:设为定解问题(1)~(3)的解,为差分格式(17)~(19)的解,记

那么

证:式(14)~(16)减去(17)~(19)得到误差方程

(28)

(29)

(30)

由定理1的(20)式和引理2,则存在常数,有:

又由故:

类似地,可以证明数值解的稳定性。

定理3:在定理2的条件下差分格式(17)~(19)的数值解依无穷模范数无条件稳定。

5. 数值算例

其中,定解问题的精确解

表1给出了取不同步长时数值解的误差

其中。为了验证此格式收敛阶为,利用公式

Table 1. Errors of numerical solution and convergence rate of example at different step sizes

表1. 取不同步长时数值解的误差及收敛比率

Figure 1. The surface of exact solution and numerical solution when t = 1 (h = 1/20, τ = 1/400)

图1. t = 1时的精确解曲面(h = 1/20, τ = 1/400)与数值解曲面

Figure 2. The surface of absolute error when t = 1 (h = 1/20, τ = 1/400)

图2. t = 1时的误差曲面(h = 1/20, τ = 1/400)与数值解曲面

计算出相邻两次取不同空间步长时的比率。从表1中可以看出比率接近4,从而验证了此格式的收敛阶为图1为算例时的精确解曲面与数值解曲面。可以看出此格式的数值解与精确解是吻合的。图2为算例时的误差曲面。进一步说明了此格式求解的准确性。

文章引用

经 鑫,张鲁明. Sobolev方程的一个紧致差分格式
A Compact Finite Difference Scheme for Sobolev Equations[J]. 理论数学, 2017, 07(01): 1-9. http://dx.doi.org/10.12677/PM.2017.71001

参考文献 (References)

  1. 1. Barenblett, G.I., Zheltov, I.P. and Kochian, I.N. (1990) Basic Concepts in the Theory of Homogeneous Liquids in Fissured Rocks. Journal of Applied Mathematics and Mechanics, 24, 1286-1303. https://doi.org/10.1016/0021-8928(60)90107-6

  2. 2. 施德明. 非线性湿气迁移方程的边值问题[J]. 应用数学学报, 1990, 13(1): 33-40.

  3. 3. Ting, T.W. (1974) A Cooling Process According to Two-Temperature Theory of Heat Conduction. Journal of Mathematical Analysis and Applications, 45, 23-31. https://doi.org/10.1016/0022-247X(74)90116-4

  4. 4. 郭会, 芮洪兴. Sobolev方程的最小二乘Galerkin有限元方法[J]. 应用数学学报, 2006, 29(4): 610-618.

  5. 5. 郭玲, 陈焕贞. Sobolev方程的H1-Galerkin混合有限元方法[J]. 系统科学与数学, 2006, 26(3): 302-314.

  6. 6. 何斯日古楞, 李宏. Sobolev方程的时间间断Galerkin有限元方法[J]. 高校应用数学学报, 2011, 26(4): 467-473.

  7. 7. 李宏, 罗振东, 安静, 孙萍. Sobolev方程的全离散有限体积元格式及数值模拟[J]. 计算数学, 2012, 34(2): 164- 172.

  8. 8. 李宏, 周文文, 方志朝. Sobolev方程的CN全离散化有限元格式[J]. 计算数学, 2013, 35(1): 41-47.

  9. 9. 史艳华, 石东洋. Sobolev方程新混合元方法的高精度分析[J]. 系统科学与数学, 2014, 34(4): 452-463.

  10. 10. Zhao, Z.H., Li, H. and Luo, Z.D. (2016) A New Space-Time Continuous Galerkin Method with Mesh Modification for Sobolev Equations. Journal of Mathematical Analysis and Applications, 400, 86-105. https://doi.org/10.1016/j.jmaa.2016.03.035

  11. 11. 罗振东, 张博. Sobolev方程基于POD的降阶外推差分算法[J]. 应用数学和力学, 2016, 37(1): 108-116.

  12. 12. Liao, H.L. and Sun, Z.Z. (2010) Maximum Norm Error Bounds of ADI and Compact ADI Methods for Solving Parabolic Equations. Numerical Methods Partial Differential Equations, 26, 37-60. https://doi.org/10.1002/num.20414

期刊菜单