Operations Research and Fuzziolgy 运筹与模糊学, 2013, 3, 7-14 http://dx.doi.org/10.12677/orf.2013.32002 Published Online May 2013 (http://www.hanspub.org/journal/orf.html) An Implicit Finite Difference Scheme for Space-Time Fractional Partial Differential Equation* Yang Zhang1, Ruiyi Wang2 1School of Mathematical Science, Nankai University, Tianjin 2School of Science, Hebei University of Technology, Tianjin Email: zhangyang@nankai.edu.cn Received: Apr. 6th, 2013; revised: Apr. 10th, 2013; a cc e p te d : M ay 1 0th, 2013 Copyright © 2013 Yang Zhang, Ruiyi Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract: Fractional order differential equations are generalizations of classical differential equations. They are widely used in the fields of diffusive transport, finance, nonlinear dynamics, signal processing and others. In this paper, an implicit finite difference method for a class of initial-boundary value space-time fractional two-sided space partial differential equatio ns with variable coefficients on a finite domain is established. The stability and convergence order are analyzed for the resulted implicit scheme. With mathematical induction skills, the scheme is proved to be unconditionally stable and convergent. Keywords: Fractional Derivative; Implicit Methods; Stability; Convergence; Error Estimate 一类分数阶导数微分方程的隐式差分解法* 张 阳1,王瑞怡 2 1南开大学数学科学学院,天津 2河北工业大学理学院,天津 Email: zhangyang@nankai.edu.cn 收稿日期:2013 年4月6日;修回日期:2013 年4月10 日;录用日期:2013 年5月10 日 摘 要:分数阶导数微分方程作为通常微分方程的推广,被广泛地应用于工程,物理,信息处理,金 融,水文等领域。本文给出了数值求解一类时间空间分数阶导数的双边空间微分方程的一种隐式差分 格式,并对其稳定性和收敛性进行了理论分析,证明了格式的无条件稳定性并给出了收敛阶估计。 关键词:分数阶导数;隐格式;稳定性;收敛性;误差估计 1. 引言 分数阶导数方程是指方程中含有非整数阶的导数,在反常扩散,粘弹性本构建模,信号处理,控制,流体 力学,图像处理,软物质研究、金融、水文等领域中利用它们来建立数学模型,可以精确地描述这些现象。近 年来,越来越多的学者关注分数阶导数。分数阶导数微分方程,是传统的整数阶微分方程的推广[1-8]。分数阶导 数微积分的研究与兴起,伴随着数字计算机计算技术的提高而迅速发展。双边空间分数阶对流–扩散方程是一 种反常扩散,在地下水溶质运移方面,它可描述含水层中溶质运移过程中的反常扩散现象。关于这类问题的数 值解法 Meerschaert 等人分别对单边对流–扩散方程和双边扩散方程,进行差分求解,如文献[9]针对变系数反应 扩散方程 *资助信息:该工作由中国国家自然科学基金资助(11171168,11071132)。 Copyright © 2013 Hanspub 7 张阳,王瑞怡 一类分数阶导数微分方程的隐式差分解法 ,,, , crtcrt crt vrdrf rt tr r 建立显式和隐式两种差分解法,文献[10,11]建立了双边空间分数阶导数微分方程 ,,, ,, uxtuxtuxt cxt csxt txx 数值解法并进行了稳定性分析,文献[12]针对二维分数阶扩散方程 ,,,, ,, ,, uxytuxyt uxyt dxyexyqxyt txy ,, 建立了交替方向隐式差分解法;文献[13]对一维时间分数阶对流–弥散方程 2 2 ,,cxt cxt cxt vD x tx , 给出了差分算法;文献[14,15]用另一种方法对这个问题进行离散计算。文献[16]中针对二维分数阶对流–弥散方程 22 22 ,,,, ,,,,,,c xytc xytc xytc xytc xyt vD xy tx y 建立了数值解法。文献[17]则对时间空间分数阶导数微分方程 ,,,, ,, uxtuxtuxt uxt b xtcxtcs xt txxx , 的显式差分解法进行了研究。目前,关于带时间分数阶导数的双边空间分数阶导数微分方程的隐式解法尚未见 到,本文针对文献[17]中的问题(1 情形)建立了相应的隐式数值解法并进行理论分析,得到了很好的理论结果。 采用隐式差分解法可以有效地增加稳定性能,本文得到的计算格式是无条件稳定的,并且给出了收敛阶的估计, 这些结论拓广了该领域的研究成果。 本文结构如下:第二节针对研究的对象进行了描述,并给出了两个常见的引理;第三节建立了所描述方程 的隐式差分解法并进行了方法的稳定性分析;第四节给出了方法的收敛性分析和收敛阶的估计;第五节总结了 结论与展望。 2. 问题描述 本文针对带分数阶时间空间导数的双边空间微分方程建立隐式差分解法并进行理论分析,研究如下方程: uxtuxtuxt uxt bxtc xtcsxt x tx x (1) 这里,0。其中LxR tT 01 ,12 , ,bxt, ,cxt , ,cxt 均为非负有界函数,并且假定 初始条件: ,以及边界条件 ,0t f xuxL x R ,,Rt 0ux Ltux 。 方程(1)右端的左边(+)以及右边(–)的分数阶导数为 Riemann-Liouville 分数阶导数。按如下方式定义[10]: 1 d1d d dd nx n L fx f n xx xn (2) 和 1 d1 dd dd nnR n x fx f n xx x n (3) Copyright © 2013 Hanspub 8 张阳,王瑞怡 一类分数阶导数微分方程的隐式差分解法 这里 是正整数,并且满足 n1nn 。 方程(1)左端的 ,uxt t 由Caputo 分数阶导数给出定义[18]: 0 ,, 1d 1 t uxt ux t t (4) 这里是 Gamma 函数。 根据 Grünwald 估计,为了得到稳定的数值方法[9],定义左右两边分数阶导数分别按如下形式近似(Shifted Grünwald): 0 d1 lim 1 d M k Mk fx g fx kh xh (5) 和 0 d() 1 lim 1 d M k Mk fx g fx kh xh (6) 这里 M , M 均为正整数, x L hM ,Rx hM , 0 11 1,1, 1,2, ! k k k gg k k 。 引理 1. 上式定义的 k g 满足如下条件: 1) ,2) 01g10g ,3) 11 10 ! k k k gk k , 2 2 。 对于引理 1,结论(1)、(2)是显然的,下面只说明结论(3)。 因为1 ,所以 0 ,10 ,0k , 。从而有 2k 2 22 11 121 111 !! 12 1 10 ! kkk k k kk gkk k k 0 引理成立。对于 ,根据二项式定理,有 0 1, m m zz m 1z 1z。令 ,有 0 0 j j g 。由引 理1知 10, 1 , 1,2,, i j jj g gi m (7) 3. 隐格式的建立及稳定性分析 先建立方程(1)在前文所述初边值条件下的隐格式并通过圆盘定理证明该差分格式是无条件稳定的,并将采 用另一种方法来证明格式的稳定性及收敛性,并给出收敛阶的估计。 引入步长 h,时间步长 ,则有 , 0; ,0, 0,1,,, 0,1,, jn x LjhhtnjmnK 。其中 ,0 TRL h Km 。因此 i Lx R 。记 为逼近 n i u in uxt 的近似解。类似地,记 nn iinii bbxtc cxt n n 及 n ii s sxt。对于 , a a uxt x , , a a ux x t,用 shifted Gr ünwald形式的逼近,有如下表达式[9]: 1 0 ,11, i nki n k uxt g ux khth xh (8) Copyright © 2013 Hanspub 9 张阳,王瑞怡 一类分数阶导数微分方程的隐式差分解法 1 0 ,11, mi nki n k uxt g ux khth xh (9) 1111 O jnjnj n uxtuxtux th xh (10) 对于时间分数阶导数有如下估计[18]: 1 11 0 1 1 01 1111 0 ,, 1d 1 ,, 1d 1 ,, 1 2 n t in i n nj ij ij j jn ninjinj j uxt ux t t uxt uxt t uxtuxt jj (11) 故可对方程(1)建立如下的隐格式: 1 111 0 11 11 1111 1,1 ,1 00 1 2 1 nj nj nii j nn imi nnnn ii ikiikki kk uu jj uu bgcugcu hh 11 nn iki s (12) 易知,该方法是相容的,其舍入误差为 h 。若记 11 1 22 2 nn n ii i iii cc b rh hh 2 11 i 则有 11 1 11111 11 1 000 1 nimi nj njnnnnn iiiiii kiikkiik jkk uujjrurugugu s (13) 记 11 1 j cj j ,则当 时有 0n 11 021 1120 11 1101 11 33 1 iii iiiiiiii imi kiikkiikii kk 1 1 g gurggu rggu guguu s (14) 当 时有 0n 11 021 11201 11 1 111 11 1 33 1 1 22 nn iii iiiiiiii imi n nn nnj kiikkiikij inii kk j ggurggu rggu 1 01 n g ugu uducu s 1 (15) 其中 11 1 21 j dj jj , 。方程(14),(15)又可改写成 12j n 1 10 1 1110 12 0 0 nnn n nn DU US DUdUdUdUcUSn Uf (16) 其中 ,, T 01 ,,, nnn n m Uuuu T 11 0, ,,,0 nnn m Sss T 11 0,, ,,0 m ffxfx 。 Copyright © 2013 Hanspub 10 张阳,王瑞怡 一类分数阶导数微分方程的隐式差分解法 而 具有以下形式: 1, ,1;1, ,1; ij Dim jm 11 20 02 1 1 1 1 1 1 1 iii iii iji i iji ji i rg gji rg gji Dggji g ji g ji 当 当 当 当 当 其中 ,,; 00 1D00 j D1,2,,jm1 mm D ,0 mj D ,0,1, ,1jm 。 定理 1. 逼近方程(1)的隐格式(16)在01 ,12 时是无条件稳定的。 证明. 根据Gerschgorin定理,以及引理1,知当 10g 。时,矩阵的特征值位于 个圆 的并集,圆心 为 0, 1 k gk Dm ,ii D 1 11 iiiiiii i Dr gr 而半径 11 11 0,0, 10, 1 mimi iikikikiiiiii kkikkkk Dggrrggr i 利用(7),矩阵 的特征值满足 Diii ii DD i ,因 此 112 2 ii r i , 的谱半径不小于 1,因而 的谱半径不超过1,从而稳定性得证。 D 1 D 下面用另一种方法给出隐格式关于最大模的稳定性分析。 引理 2. 上面所定义的 11 1 j cj j 满足以下条件 1) , 0, 0,1, j cj 2) , 1,0,1, jj cc j 3) ,即 。 1 11 1 22 1 k jj k j cc c 1 11 1 22 1 k jk j dc 对于(1),(3)是显然的,(2)利用单调性也可以得到。令 11 1 10211120 11 11 11 33 1 nn n iiii iiiiiii imi nn kiikkiik kk Luggurgg urggu gu gu 1 1 n i 0 1 210 1 1 0 22 0 i nn innj ijini j un Lu uducun 当 当 1 (14),(15)可以改写为 11 12 ,0,1,2,, nnn iii LuLusnK K 。 设 是由(14),(15)计算的近似值,则误差(扰动) 满足 121 12 j iimj u 121 12 jj j i ii uimj K u 1 12 0121 121 nn ii LLnKi m 记 ,则可导出下面的定理 T 12 1 kkk k m E 定理 2. 误差满足 k E1001 ,1 n EEn K ,其中 11 max kk i im E 。 证明. 用数学归纳法。 1) 当 时,0n1 11 max l im 1 i ,利用引理1,我们得到 Copyright © 2013 Hanspub 11 张阳,王瑞怡 一类分数阶导数微分方程的隐式差分解法 11 111 11 021112011 33 11 1111 021 1120111 33 10 12 1 1 lml lllllllllllllk lklk lk kk lml llll lllllllkllkkllk kk ll ggrggrgggg ggrgg rgggg LL 1 1 1 00 lE 2) 假设 0012 k EEk n成立,则 11 11 max nn il im E 1n ,利用引理 1,有 11 11111 021112011 33 1 1111 021 1120111 3 1 1 lml nnnnn lllllllllllllklklk lk kk l nnnn lllllllllll kllkkllk k ggrgg rgggg ggrgg rgggg 1 1 1 n n 1 3 1 11 0 12 1 1 100 1 22 22 ml k n nnn nj ll ljlnl j n jn j LL dc dcE E 定理得证。 由上述定理知格式(16)无条件稳定。 4. 隐格式的收敛性分析与收敛阶估计 下面进行收敛性分析。令 是方程(1)在相应初边值条件下于点 ,1,2,,1;1,2, ik uxtim kK , ik x t的 精确解。定义 ,1,2,,1;1,2, kk iiki euxtuimkK 并且T 12 1 ,,, kkk k m Yee e ,则有如下表达式 11 12 0 10 kk ii LeLe R e k 其中。记 为局部截断误差,则存在常数,使得 1, 2,,1;0,1, 2,imkK 1k i R0C 1, 1,2,,1;0,1,2,1 k i RC himkK 定理 3. 误差 满足 k Y 11kk YCch0,1,2, ,1k K, 。这里 11 max kk e i im Y ,C为正常数。 证明. 用数学归纳法。 1˚ 当 时,令 0k1 11 max l im e 1 i e,利用引理1,有 11 111 11 021 112011 33 11 1111 021 1120111 33 10 11 1 1 lml llllllllllll lklklklk kk lml lllllllllll kllkkllk kk ll eggergg erggegege ggerggerggegege Le Le 1 1 1 11 ll RRC h 注意到 ,因此 1 01c 11 0 YCc h。 2˚ 假设 1 0,1,2, n YCc hn k成立,下面估计 11 11 max kk il im Ye 1k e 利用引理1,有 Copyright © 2013 Hanspub 12 张阳,王瑞怡 一类分数阶导数微分方程的隐式差分解法 11 11 111 021 1120111 33 1 11111 021 1120111 3 1 1 lml kk kkk llll llll llllslslsls ss l kkkkk llllllll lllsllsslls s gerggerggege ge ggergge rggegege 1k l eg 1 3 11 12 1 10 1 1 1 1 22 22 ml s kkk lll k kkj ljlkl j k kkj j j LeLe R edeceC h YdYCh 利用引理2,,因此 11 1 kk cc 1,1,2, jk YCc hj k,由上面表达式可得 11k ed Cch 1 1 1 1 22 k ijkkk j cCch 故 11kk YCc h定理得证。 根据定理 3,并注意到 11 11 1 1 lim lim k ckk lim 1 11 11 kk k kkk k ,因此存在常数 C ,使得 11kk YCc h ,因 为kT ,进而 k YCT h,从而有下述定理。 定理 方程(1)在相应初边值条件下在点4. 假设 ,1;1,2,m kK 是 ,1,2, ik uxt i , ik x t处的精确解。 是 k i u隐格式(16)所得到的数值解,则存在正常数 C ,使得 ,, 1, k ik i uxtuCh i 2,, 1; 1,2,m kK 注:由定理4知该数值方法的误差估计达到最优阶。 5. 结 分数阶导数微分方程的时间导数推广到分数阶,建立了一种隐式差分求解格式并进行了相 应的 论与展望 本文将双边空间 理论分析,得到了格式的无条件稳定性及丰满的一阶时间和空间的收敛阶逼近。所得结论亦可推广到方程 右端项包含更为复杂的情形,如 xt uxt uxt cxtcxtdxy uxtsxt x txx 所采用的方法也适用于高维空间情形,至于高维问题如何进行算子分裂尚有待探讨。最后,作者对审稿人对稿 参考文献 (References) urant at the end of the random walk: Recent developments in the description of anomalous transport by frac- ics, 2004, A37: 161-208. ce of the difference methods for the space-time fractional advec- uxt u b xt 件提出的修改意见表示感谢。 [1] R. Metzler, J. Klafter. The resta tional dynamics. Journal of Phys [2] Z. Deng, V. P. Singh and L. Bengtsson. Numerical solution of fractional advection-dispersion equation. Journal of Hydraulic Engineering, 2004, 130(5): 422-431. [3] R. Gorenflo, F. Mainardi, E. Scalas and M. Raberto. Fractional calculus and continuous-time finance III: The diffusion limit. In: Kolhmann, S. Tang, Eds. Mathematical Finance. Basel: Birkhauser Verlag, 2001: 171-180. [4] 苏丽娟, 王文洽. 双边分数阶对流–扩散方程的一种有限差分解法[J]. 山东大学学报(理学版), 2009, 10: 29-32. [5] F. Liu, P. Zhuang, V. Anh, I. Turner and K. Burrage. Stability and convergen Copyright © 2013 Hanspub 13 张阳,王瑞怡 一类分数阶导数微分方程的隐式差分解法 Copyright © 2013 Hanspub 14 on and Hong Kong: World Scientific Publication Com- C. Tadjeran. Finite difference approximations for fractional advection-dispersion flow equations. Journal of Computational 6, 211: 249-261. tion of the space fractional Fokker-Planck equation. Journal of Computational & Applied Mathe- 李宁平, 陈璐. 分数阶导数双边空间微分方程的显式差分解法[J]. 运筹与模糊学, 2012, 2, 1-7. onal diffusion equations. Physical tion-diffusion equation. Applied Mathematics and Computation, 2007, 191(1): 12-21. [6] 周激流, 蒲亦非, 廖科. 分数阶微积分原理及其在现代信号分析与处理中的应用[M]. 北京: 科学出版社, 2010. [7] R. Hilfer. Application of fractional calculus in physics. Singapore, New Jersey, Lond pany, 2000. [8] A. A. Kibas, H. M. Srivastava and J. J. Trujillo. Theory and application of fractional differential equation s. Amsterdam: Elservier, 2006. [9] M. M. Meerschaert, & Applied Mathematics, 2004, 172(1): 65-77. [10] M. M. Meerschaert, C. Tadjeran. Finite difference approximations for two-sided space-fractional partial differential equations. Applied Nu- merical Mathematics, 2006, 56(1): 80-90. [11] V. K. Tuan, R. Gorenflo. Extrapolation to the limit for numerical fractional differentiation. Zeitschrift für Angewandte Mathematik und Me- chanik, 1995, 75: 646-648. [12] M. M. Meerschaert, H. P. Scheffler and C. Tadjeran. Finite difference method for two dimensional fractional dispersion equations. Journal of Computational Physics, 200 [13] 夏源, 吴吉春. 分数阶对流——弥散方程的数值求解[J]. 南京大学学报(自然科学版), 2007, 43(4): 44-446. [14] F. Liu, V. Ahn and I. Turner. Numerical solu matics, 2004, 166(1): 209-219. [15] Y. Zhang. A finite difference method for fractional partial differential equation. Applied Mathematics and Computation, 2009, 215(2): 524- 529. [16] 周璐莹, 吴吉春, 夏源. 二维分数阶对流–弥散方程的数值解[J]. 高校地质学报, 2009, 15(4): 569-575. [17] 张阳, [18] M. M. Meerschaert, D. A. Benson, H. P. Scheffler and B. Baeumer. Stochastic solution of space-time fracti Review, 2002, E65: 1103-1106. |