亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        兩邊空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的一種加權(quán)顯式有限差分方法

        2016-05-06 01:59:24馬亮亮劉冬兵

        馬亮亮, 劉冬兵

        (攀枝花學(xué)院 數(shù)學(xué)與計(jì)算機(jī)學(xué)院, 四川 攀枝花 617000)

        ?

        兩邊空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的一種加權(quán)顯式有限差分方法

        馬亮亮,劉冬兵

        (攀枝花學(xué)院 數(shù)學(xué)與計(jì)算機(jī)學(xué)院, 四川 攀枝花 617000)

        摘要:考慮兩邊空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的初邊值問題,基于Grünwald公式和移位Grünwald-Letnikov公式,提出一種加權(quán)顯式有限差分解法.利用傅里葉變換和特征值法,得到差分格式的穩(wěn)定性.然后使用最大模估計(jì)法證明在相同的條件下,所提出的差分格式是收斂的.最后通過數(shù)值例子說明所提出的差分格式是可靠和有效的,并對(duì)方程的數(shù)值解與精確解進(jìn)行比較,驗(yàn)證了文中的理論結(jié)果.

        關(guān)鍵詞:分?jǐn)?shù)階對(duì)流-擴(kuò)散方程; 空間分?jǐn)?shù)階導(dǎo)數(shù); 加權(quán)差分格式; 收斂性; 穩(wěn)定性; 有限差分法

        1預(yù)備知識(shí)

        分?jǐn)?shù)階微分方程中含有非整數(shù)階的導(dǎo)數(shù),能夠有效地描述各種各樣物質(zhì)的記憶和遺傳性質(zhì),在工程、物理、金融、水文、地理、材料等領(lǐng)域發(fā)揮了愈來愈重要的作用[1-14].遺憾的是,大多數(shù)分?jǐn)?shù)階微分方程的解析解都含有復(fù)雜的級(jí)數(shù)或者特殊函數(shù),不利于進(jìn)行近似計(jì)算,于是對(duì)分?jǐn)?shù)階微分方程進(jìn)行數(shù)值求解變得尤為重要.目前,許多學(xué)者致力于其中一類分?jǐn)?shù)階微分方程的-(時(shí)間、空間、空間-時(shí)間)分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的研究[15-23].

        本文主要研究下列兩邊空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的數(shù)值解法

        (1)

        (2)

        其中,0≤n-1<α

        兩邊空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程(1)是一種反常擴(kuò)散方程,在地下水溶質(zhì)運(yùn)移方面,該方程可用來描述含水層中溶質(zhì)運(yùn)移過程中的反常擴(kuò)散現(xiàn)象[24].關(guān)于這類問題的數(shù)值解法,前人已做了一些研究.蘇麗娟等[24]給出了雙邊空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的一種隱式有限差分解法.M. Meerschaert等[25-26]分別對(duì)單邊和雙邊對(duì)流-擴(kuò)散方程利用改進(jìn)型Grünwald-Letnikov差分方法進(jìn)行了求解,但結(jié)果表明用Grünwald得到的顯式差分格式是不穩(wěn)定的,因此數(shù)值解不會(huì)收斂于方程的精確解.夏源等[27]對(duì)單邊的時(shí)間和空間分?jǐn)?shù)階對(duì)流-彌散方程給出了類似的算法.F. Liu等[28]用此類算法對(duì)空間分?jǐn)?shù)階???普朗克方程進(jìn)行了離散計(jì)算.馬維元等[29]對(duì)空間-時(shí)間分?jǐn)?shù)階擴(kuò)散方程的初邊值問題提出了一種加權(quán)平均差分格式.張紅玉等[30]給出了一類空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的加權(quán)平均有限差分解法.本文在標(biāo)準(zhǔn)Grünwald公式的基礎(chǔ)上,利用移位Grünwald-Letnikov公式,構(gòu)造出一種新的加權(quán)顯式有限差分格式,并對(duì)其穩(wěn)定性和收斂性進(jìn)行理論分析,最后通過數(shù)值例子驗(yàn)證了該差分格式的有效性和可靠性.

        2加權(quán)顯式有限差分法

        (5)

        (6)

        (7)

        (9)

        (10)

        (9)~(10)式在時(shí)間和空間上都只有一階精度,為了提高空間上的精度,對(duì)Grünwald公式和移位Grünwald-Letnikov公式進(jìn)行加權(quán)處理得

        于是得到下面的加權(quán)顯式有限差分格式

        (11)

        其中,i=1,2,…,M-1,n=0,1,…,N-1;ε是權(quán)參數(shù),且0<ε<1.

        (12)

        (13)

        3穩(wěn)定性和收斂性分析

        定義 1對(duì)于任意的u∈R,(iu)α=sign(u)|u|αe(iπα)/2.

        借助于定義1和2,對(duì)(14)和(15)式做傅里葉變換得

        因此對(duì)任意的x∈R,|φ(x,h)|≤C3h2.

        為半徑的圓盤上,因此要使顯式有限差分格(11)式穩(wěn)定,需|λ-Ai,i|

        現(xiàn)要使顯式有限差分格(11)式穩(wěn)定,只需證明Ai,i-ri≥-1即可.而要使

        由泰勒展開式定理和定理1得

        令R=O(τ+h2)(1,1,…,1)T,結(jié)合(2)、(3)和(13)式,(16)式可寫成矩陣形式

        所以

        進(jìn)而

        4數(shù)值實(shí)例

        考慮如下兩邊空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程

        表1給出了當(dāng)t=0.5時(shí)不同時(shí)間和不同空間步長(zhǎng)下數(shù)值解與精確解之間的最大模誤差及其誤差比率,其中誤差比率是由前一行的最大模誤差與該行的最大模誤差相比而得到.

        表 1 t=0.5時(shí)的最大模誤差和收斂階

        表 2 t=0.5時(shí)方程數(shù)值解與精確解的比較

        5結(jié)語(yǔ)

        本文考慮了兩邊空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的數(shù)值逼近問題,首先采用Grünwald公式和移位Grünwald-Letnikov公式近似空間分?jǐn)?shù)階導(dǎo)數(shù),得到一種加權(quán)顯式有限差分解法;然后利用傅里葉變換和特征值法給出差分格式的穩(wěn)定性分析,借助于最大模估計(jì)方法驗(yàn)證差分格式的收斂性,最后通過數(shù)值實(shí)例驗(yàn)證了文中方法的精確性、可靠性和實(shí)用性.

        致謝攀枝花學(xué)院校級(jí)培育項(xiàng)目(2012PY08)、攀枝花學(xué)院校級(jí)科研項(xiàng)目(2013YB05和2012YB21)、攀枝花學(xué)院院級(jí)科研創(chuàng)新項(xiàng)目(Y2013-04)和攀枝花學(xué)院教育教學(xué)研究與改革青年項(xiàng)目(JJ1376)對(duì)本文給予了資助,謹(jǐn)致謝意.

        參考文獻(xiàn)

        [1] PODLUBNY I. Fractional Differential Equations[M]. New York, London:Academic Press,1999.

        [2] HILFER R. Applications of Fractional Calculus in Physics[M]. Singapore:Word Scientific,2000.

        [3] MORTON K W, MARYERS D F. Numerical Solution of Partial Differential Equation[M]. Cambridge:Cambridge University Press,2005.

        [4] LUCHKO Y. Initial-boundary-value problems for the generalized multi-term time-fractional diffusion equation[J]. J Math Anal Appl,2011,374(2):538-548.

        [5] 馬亮亮,劉冬兵. 二維變系數(shù)空間分?jǐn)?shù)階電報(bào)方程數(shù)值解[J]. 遼寧工程技術(shù)大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,33(3):429-432.

        [6] 馬亮亮,劉冬兵. 一類維空間Riesz分?jǐn)?shù)階擴(kuò)散方程的解析解[J]. 合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,37(4):506-509.

        [7] 馬亮亮. 時(shí)間分?jǐn)?shù)階擴(kuò)散方程的數(shù)值解法[J]. 數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2013,43(10):248-253.

        [8] 馬亮亮. 一種Caputo分?jǐn)?shù)階反應(yīng)-擴(kuò)散方程初邊值問題的隱式差分格式[J]. 貴州師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,31(1):58-61.

        [9] 馬亮亮,田富鵬. 空間分?jǐn)?shù)階Edwards-Wilkinson方程的顯式差分近似[J]. 沈陽(yáng)大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,25(3):250-252.

        [10] 馬亮亮,劉冬兵. 一類反常次擴(kuò)散方程N(yùn)eumann問題的有限差分格式收斂性分析[J]. 五邑大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,28(1):1-4.

        [11] 馬亮亮,劉冬兵. 變系數(shù)分?jǐn)?shù)階反應(yīng)-擴(kuò)散方程的數(shù)值解法[J]. 沈陽(yáng)大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,26(1):76-80.

        [12] 馬亮亮. 時(shí)間分?jǐn)?shù)階擴(kuò)散方程的隱式差分近似[J]. 貴州師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,32(2):79-82.

        [13] 劉冬兵,馬亮亮. 變時(shí)間分?jǐn)?shù)階反應(yīng)擴(kuò)散方程的數(shù)值分析[J]. 江南大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,13(1):109-112.

        [14] 馬亮亮,劉冬兵. 高維非齊次時(shí)間分?jǐn)?shù)階電報(bào)方程的基本解[J]. 四川師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,38(1):77-83.

        [15] LIU F, ZHUANG P, ANH V, et al. Stability and convergence of the difference methods for the space-time fractional advection-diffusion equation[J]. Appl Math Comput,2007,191(1):2-20.

        [16] JIANG H, LIU F, TURNER I, et al. Analytical solutions for the multi-term time-space Caputo-Riesz fractional advection-diffusion equations on a finite domain[J]. J Math Anal Appl,2012,389(2):1117-1127.

        [17] CARELLA A R, DORAO C A. Least-squares spectral method for the solution of a fractional advection-dispersion equation[J]. J Comput Phys,2013,232(1):33-45.

        [18] 馬亮亮,劉冬兵. 高維分?jǐn)?shù)階cable方程隱式差分逼近[J]. 遼寧工程技術(shù)大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,33(4):544-547.

        [19] 馬亮亮,田富鵬. 變系數(shù)空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的隱式差分近似[J]. 中北大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,35(1):11-14.

        [20] 馬亮亮. 一種時(shí)間分?jǐn)?shù)階對(duì)流擴(kuò)散方程的隱式差分近似[J]. 西北民族大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,34(1):7-12.

        [21] 馬亮亮. 變系數(shù)階空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的有限差分解法[J]. 沈陽(yáng)大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,25(4):341-344.

        [22] 馬亮亮. 變時(shí)間分?jǐn)?shù)階非定常對(duì)流擴(kuò)散方程的數(shù)值分析[J]. 遼東學(xué)院學(xué)報(bào)(自然科學(xué)版),2013,20(3):220-223.

        [23] 馬亮亮,劉冬兵. 一類變時(shí)間分?jǐn)?shù)階含源項(xiàng)非定常奇異攝動(dòng)對(duì)流擴(kuò)散方程的數(shù)值分析[J]. 沈陽(yáng)大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,25(5):424-427.

        [24] 蘇麗娟,王文洽. 雙邊空間分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的一種有限差分解法[J]. 山東大學(xué)學(xué)報(bào)(理學(xué)版),2009,44(10):26-29.

        [25] MEERSCHAERT M, TADJERAN C. Finite difference approximations for fractional advection-dispersion flow equations[J]. J Comput Appl Math,2004,172(1):65-77.

        [26] MEERSCHAET M, TADJERAN C. Finite difference approximations for two-sided space fractional partial differential equations[J]. Appl Num Math,2006,56(1):80-90.

        [27] 夏源,吳吉春. 分?jǐn)?shù)階對(duì)流-彌散方程的數(shù)值求解[J]. 南京大學(xué)學(xué)報(bào)(自然科學(xué)版),2007,43(4):441-446.

        [28] LIU F, ANH V, TURNER I. Numerical solution of the space fractional Fokker-Plank equation[J]. J Comput Appl Math,2004,166(1):209-219.

        [29] 馬維元,劉華. 兩邊空間-時(shí)間分?jǐn)?shù)階擴(kuò)散方程的加權(quán)有限差分格式[J]. 華東師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2012(3):41-48.

        [30] 張紅玉,崔明榮. 兩類分?jǐn)?shù)階對(duì)流-擴(kuò)散方程的有限差分方法[J]. 山東大學(xué)學(xué)報(bào)(理學(xué)版),2012,47(6):40-48.

        [31] 汪向艷,朱琳,芮洪興. 空間分?jǐn)?shù)階擴(kuò)散方程的一種加權(quán)顯式差分方法[J]. 寧夏大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,35(1):1-5.

        [32] 徐仲,張凱院,陸全,等. 矩陣簡(jiǎn)明教程[M]. 北京:科學(xué)出版社,2001.

        2010 MSC:35R11; 39A99

        (編輯鄭月蓉)

        A Weighted Explicit Finite Difference Method for Two-Sided Space Fractional Advection Diffusion Equation

        MA Liangliang,LIU Dongbing

        (CollegeofMathematicsandComputer,PanzhihuaUniversity,Panzhihua617000,Sichuan)

        Abstract:Based on Grünwald formula and the shifted Grünwald-Letnikov formula, a weighted explicit finite difference method is proposed to solve initial boundary value problems of two-sided space fractional advection diffusion equation. Their stability is analyzed by means of Fourier transform and eigenvalue analysis. Using the technique of maximum norm analysis, it is proved that the scheme convergent under the same condition. Illustrative example is included to demonstrate the validity and applicability of the scheme, and a comparison between the exact analytical and the numerical prediction is made to demonstrate the theoretical results.

        Key words:fractional advection diffusion equation; space fractional derivative; weighted difference scheme; convergence; stability; finite difference method

        doi:10.3969/j.issn.1001-8395.2016.01.013

        中圖分類號(hào):O241.82

        文獻(xiàn)標(biāo)志碼:A

        文章編號(hào):1001-8395(2016)01-0076-07

        作者簡(jiǎn)介:馬亮亮(1986—),男,講師,主要從事模型優(yōu)化和微分方程的研究,E-mail:mllpzh@126.com

        基金項(xiàng)目:國(guó)家自然科學(xué)基金(10671132和60673192)、四川省科技廳項(xiàng)目(2013JY0125)和攀枝花市市級(jí)應(yīng)用技術(shù)研究與開發(fā)資金項(xiàng)目

        收稿日期:2014-05-14

        (2014CY-G-22)

        久久狼精品一区二区三区| www国产无套内射com| 中文字幕在线播放| 欧美性videos高清精品| 浪荡少妇一区二区三区| 成人无码激情视频在线观看| 99久久免费精品色老| 日韩三级一区二区三区四区| 久久综合噜噜激激的五月天| 欧美性受xxxx白人性爽| 国产亚洲av片在线观看18女人| 亚洲九九夜夜| 亚洲国产精品色一区二区| 一区二区三区日韩蜜桃| 亚洲av综合一区二区在线观看| 国产午夜精品久久久久免费视| 亚洲国产AV无码男人的天堂| 国产美女在线精品亚洲二区| 日本在线中文字幕一区| 亚洲精品午夜久久久九九| 天堂国产一区二区三区| 老司机亚洲精品影院| 国产成人av 综合 亚洲| 欧美精品中文| 日韩精品欧美激情国产一区| 亚洲精品女同在线观看| 亚洲性无码av中文字幕 | 18禁裸男晨勃露j毛免费观看 | 精品国产福利一区二区在线| 91美女片黄在线观看| 成年人视频在线播放麻豆| 免费的小黄片在线观看视频| 亚洲av无码久久精品狠狠爱浪潮| 亚洲男女免费视频| 手机免费在线观看日韩av| 亚洲国产精品久久艾草| 国产精品久久久久乳精品爆| 有码精品一二区在线| 女同久久精品国产99国产精| 337p日本欧洲亚洲大胆色噜噜 | 国产精品久久久久久影视|