危國華

(福建廣播電視大學(xué)三明分校，福建 三明 365000)

0 引言

近年來，分?jǐn)?shù)階微積分以及分?jǐn)?shù)階微分方程在模擬自然界的各種現(xiàn)象中得到了廣泛應(yīng)用[1-4]，隨后分?jǐn)?shù)階微分方程的數(shù)值求解也成為研究者們的一個(gè)研究熱點(diǎn)[5-11]。針對(duì)空間分?jǐn)?shù)階微分方程的數(shù)值方法，已有文獻(xiàn)采用有限差分方法[12]、有限元方法[13]和譜方法[14]進(jìn)行處理。在離散微分方程之前，有限差分方法需要對(duì)區(qū)域利用規(guī)則均勻點(diǎn)進(jìn)行劃分，有限元方法和譜方法需要對(duì)區(qū)域進(jìn)行劃分并生成網(wǎng)格，這些方法使得在處理二維情形時(shí)僅限于規(guī)則區(qū)域上的問題。鑒于無網(wǎng)格方法不需要生成網(wǎng)格，可以方便處理不規(guī)則區(qū)域，本文首次嘗試將其應(yīng)用到空間分?jǐn)?shù)階微分方程中，采用多項(xiàng)式基點(diǎn)插值配置法處理帶有雙側(cè)分?jǐn)?shù)階導(dǎo)數(shù)的空間分?jǐn)?shù)階微分方程，在數(shù)值例子中分別采用等間距節(jié)點(diǎn)以及不規(guī)則散點(diǎn)離散空間變量，均得到了較好的結(jié)果。

本文討論如下帶有雙側(cè)分?jǐn)?shù)階導(dǎo)數(shù)的空間分?jǐn)?shù)階微分方程：

(1)

本文假設(shè)空間分?jǐn)?shù)階微分方程的解滿足u(x,·)∈C2(0,T)。

1 離散格式

1.1 多項(xiàng)式點(diǎn)插值

首先，在問題域[a,b]內(nèi)生成場節(jié)點(diǎn)，節(jié)點(diǎn)可以隨機(jī)生成，也可以人工加入節(jié)點(diǎn)。設(shè)有內(nèi)部節(jié)點(diǎn)x1,x2,…,xd和邊界節(jié)點(diǎn)x0,xd+1。

對(duì)每一個(gè)場節(jié)點(diǎn)，都做一個(gè)包含該節(jié)點(diǎn)的支持域，在實(shí)際計(jì)算中，一般取節(jié)點(diǎn)平均間距的2～3倍大小[15]。設(shè)第i個(gè)節(jié)點(diǎn)的支持域?yàn)棣竔，支持域Ωi內(nèi)包含ni個(gè)節(jié)點(diǎn)xi1,xi2,…,xini。為書寫方便，記Di={i1,i2,…,ini}，這也就意味著，當(dāng)l∈Di時(shí)，xl∈Ωi。

接下來在節(jié)點(diǎn)的支持域內(nèi)對(duì)函數(shù)u(x)及其左側(cè)分?jǐn)?shù)階導(dǎo)數(shù)利用基函數(shù)進(jìn)行插值逼近。為了描述方便，假設(shè)計(jì)算點(diǎn)x的支持域內(nèi)包含n個(gè)場節(jié)點(diǎn)x1,x2,…,xn，連續(xù)函數(shù)u(x)對(duì)應(yīng)于這些場節(jié)點(diǎn)上的函數(shù)值分別為u1,u2,…,un。因此，函數(shù)u(x)可以由這組場節(jié)點(diǎn)近似表示為：

(2)

其中pj(x)(j=1,2,…,m)為空間坐標(biāo)變量的單項(xiàng)式，被稱為多項(xiàng)式基，m是多項(xiàng)式基的個(gè)數(shù)，aj(j=1,2,…,m)為一組待定常數(shù)。通常使用的一維k次完備多項(xiàng)式基為

(3)

其中基函數(shù)的個(gè)數(shù)滿足m=k+1。

Us=Pm·a，

(4)

由于n=m，矩陣Pm是n×n維的方陣。由式(4)，有：

(5)

(6)

1.2 全離散格式

接下來離散時(shí)間變量，記tn=nτ,n=0,1,…,N，其中τ=T/N為時(shí)間步長。全離散格式為：

(7)

2 數(shù)值結(jié)果

為了避免在計(jì)算中產(chǎn)生Runge現(xiàn)象，在每個(gè)場節(jié)點(diǎn)xi處，取n=6，也就是選擇6個(gè)臨近節(jié)點(diǎn)作為其支持域{xl,l=i1,i2,…,i6}，于是k=5，基函數(shù)為5次完備多項(xiàng)式基。

先將問題域[0，1]用規(guī)則節(jié)點(diǎn)xi=ih(i=0,1,2,…,50)劃分，采用全離散格式式(7)進(jìn)行求解(計(jì)算中取m=5)，從圖1中可以看出，精確解和數(shù)值解吻合得很好。表1中列出時(shí)間步長變化時(shí)，數(shù)值解與精確解的誤差以及誤差階。

接下來將問題域[0，1]用不規(guī)則節(jié)點(diǎn)劃分，采用全離散格式式(7)進(jìn)行求解。圖2給出此時(shí)精確解與數(shù)值解，表2列出針對(duì)不規(guī)則點(diǎn)劃分空間變量時(shí)間步長變化時(shí)，數(shù)值解與精確解的誤差以及誤差階?？梢钥闯?，本文所提出的數(shù)值方法仍然適用于不規(guī)則點(diǎn)劃分。

表1 采用規(guī)則點(diǎn)劃分求解具有單側(cè)導(dǎo)數(shù)的方程的數(shù)值解與精確解的誤差以及誤差階

τε2order2ε∞order∞1/102.214 9e-002—3.708 8e-002—1/201.131 1e-0020.969 51.894 0e-0020.969 51/405.715 4e-0030.984 99.569 8e-0030.984 91/802.872 7e-0030.992 54.809 9e-0030.992 51/1601.440 1e-0030.996 22.411 3e-0030.996 21/3207.210 1e-0040.998 11.207 2e-0030.998 1

取α= 1.8，分別將問題域[0，1]用規(guī)則節(jié)點(diǎn)xi=ih(i=0,1,2,…,50)和不規(guī)則節(jié)點(diǎn)劃分，采用全離散格式式(7)進(jìn)行求解(計(jì)算中取m=5)，從圖3和圖4中可以看出，無論是規(guī)則點(diǎn)，還是不等距節(jié)點(diǎn)，精確解和數(shù)值解吻合得很好。表3和表4給出規(guī)則點(diǎn)和不規(guī)則點(diǎn)劃分空間變量，時(shí)間步長變化時(shí)的誤差以及誤差階。

表2 采用不規(guī)則點(diǎn)劃分求解具有單側(cè)導(dǎo)數(shù)方程的數(shù)值解與精確解的誤差以及誤差階

τε2order2ε∞order∞1/102.184 1e-002—3.702 5e-002—1/201.115 4e-0020.969 51.890 8e-0020.969 51/405.635 9e-0030.984 99.553 5e-0030.984 91/802.832 7e-0030.992 54.801 8e-0030.992 51/1601.420 1e-0030.996 22.407 1e-0030.996 21/3207.109 6e-0040.998 11.205 1e-0030.998 1

表3 采用規(guī)則點(diǎn)劃分求解具有雙側(cè)導(dǎo)數(shù)的方程的數(shù)值解與精確解的誤差以及誤差階

τε2order2ε∞order∞1/103.867 4e-003—2.838 5e-004—1/201.933 7e-0031.000 01.419 2e-0041.000 01/409.668 6e-0041.000 07.096 2e-0051.000 01/804.834 3e-0041.000 03.548 1e-0051.000 01/1602.417 2e-0041.000 01.774 1e-0051.000 01/3201.208 7e-0040.999 98.870 7e-0061.000 0

表4 采用規(guī)則點(diǎn)劃分求解具有雙側(cè)導(dǎo)數(shù)的方程的數(shù)值解與精確解的誤差以及誤差階

τε2order2ε∞orde∞1/104.883 4e-003—3.073 9e-004—1/202.441 7e-0031.000 01.537 0e-0041.000 01/401.220 8e-0031.000 17.684 3e-0051.000 11/806.102 8e-0041.000 23.841 6e-0051.000 21/1603.050 4e-0041.000 51.920 2e-0051.000 51/3201.524 2e-0041.000 99.594 8e-0061.000 9

3 結(jié)論

本文將基于多項(xiàng)式基點(diǎn)插值配置法處理帶有雙側(cè)分?jǐn)?shù)階導(dǎo)數(shù)的空間分?jǐn)?shù)階微分方程，在數(shù)值例子中分別采用等間距節(jié)點(diǎn)以及不規(guī)則散點(diǎn)離散空間變量，均得到了較好的逼近結(jié)果。但是，由于無網(wǎng)格方法理論上的匱乏，即使在討論傳統(tǒng)的整數(shù)階微分方程時(shí)，此方法在空間上的逼近階數(shù)也不能得到保證。由于傳統(tǒng)方法在不規(guī)則區(qū)域上的局限性，在后續(xù)工作中，筆者將繼續(xù)探討將其應(yīng)用到二維不規(guī)則區(qū)域上的空間分?jǐn)?shù)階微分方程。

[]

[1]HAVLIN S，BEN-AVRAHAM D.Diffusion in disordered media：advances in physics[J].2002，51:187-292.DOI:10.1016/0169-7439(91)80040-W.

[2]ATHAR M，F(xiàn)ETECAU C.Unsteady flow of a generalized Maxwell fluid with fractional derivative due to a constantly accelerating plate[J].Computers and Mathematics with Applications， 2009，57:596-603.DOI:10.1016/j.camwa.2008.09.052.

[3]YANG D，ZHU K Q.Start-up flow of a viscoelastic fluid in a pipe with the fractional Maxwell’s model[J].Computers & Mathematics with Applications，2010，60(8):2231-2238.DOI:10.1016/ j.camwa.2010.08.013.

[4]王建宏，殷姝.一類分?jǐn)?shù)階混沌系統(tǒng)的滑膜控制[J].機(jī)械制造與自動(dòng)化，2016，45(3):180-183.DOI:10.19344/j.cnki.issn1671-5276.2016.03.054.

[5]LIU F，ANH V，TURNER I,et al.Numerical simulation for solute transport in fractal porous media[J].Australian and New Zealand Industrial and Applied Mathematics Journal，2004，45:461- 473.DOI:10.21914/ANZIAMJ.V45I0.901.

[6]ROOP J P.Computational aspects of FEM approximation of fractional advection dispersion equation on bounded domains inR2[J].Journal of Computational and Applied Mathematics，2006，193(1):243-268.DOI:10.1016/j.cam.2005.06.005.

[7]LIN Y M，XU C J.Finite difference/spectral approximations for the time-fractional diffusion equation[J].Journal of Computational Physics，2007，225(2):1079-1095.DOI:10.1016/j.jcp.2007.02.001.

[8]LANGLANDS T，HENRY B，WEARNE S.Solutions of a fractional cable equation:finite case[R].Sydney:University of New South Wales，2005.

[9]LIU Q，GU Y T，ZHUANG P，et al.An implicit RBF meshless approach for time fractional diffusion equations[J].Computational Mechanics，2011，48:1-12.DOI:10.1007/s00466-011-0573-x.

[10]LI C，DENG W H，SHEN X Q.Exact solutions and their asymptotic behaviors for the averaged generalized fractional elastic models[J].Communications in Theoretical Physics，2014，62(4):443- 450.

[11]ZHUANG P，LIU F，TURNER I，et al.Galerkin finite element method and error analysis for the fractional cable equation[J].Numerical Algorithms，2016，72:447- 466.DOI:10.1007/s11075-015-0055-x.

[12]ZENG F H，LIU F W，LI C P，et al. A Crank-Nicolson ADI spectral method for a two-dimensional Riesz space fractional nonlinear reaction-diffusion equation[J].SIAM Journal on Numerical Analysis，2014，52(6):2599-2622.DOI:10.1137/130934192.

[13]ZHAO Y Z，BU W B，HUANG J F，et al.Finite element method for two-dimensional space-fractional advection-dispersion equations[J].Applied Mathematics and Computation，2015，257:553-565.DOI:10.1016/j.amc.2015.01.016.

[14]ZHANG Y X，DING H F.High-order algorithm for the two-dimension Riesz space-fractional diffusion equation[J].International Journal of Computer Mathematics，2017，94(10):2063-2073.DOI:10.1080/00207160.2016.1274746.

[15]LIU G R.Mesh free methods:moving beyond the finite element method[M].Boca Raton:CRC Press，2005.

[16]PODLUBNY I.Fractional differential equations[M].NewYork:Academic Press，1999.