基于多項(xiàng)式插值逼近的分?jǐn)?shù)階偏微分方程高精度差分方法
本文關(guān)鍵詞:基于多項(xiàng)式插值逼近的分?jǐn)?shù)階偏微分方程高精度差分方法 出處:《東南大學(xué)》2016年博士論文 論文類型:學(xué)位論文
更多相關(guān)文章: 分?jǐn)?shù)階波方程 多項(xiàng)時(shí)間分?jǐn)?shù)階導(dǎo)數(shù) Bloch-Torrey方程 非線性時(shí)間分?jǐn)?shù)階反應(yīng)-擴(kuò)散方程 Caputo導(dǎo)數(shù) Riesz導(dǎo)數(shù) 超收斂點(diǎn) 穩(wěn)定性 收斂性
【摘要】:近幾十年來,由于分?jǐn)?shù)階導(dǎo)數(shù)具有非局部性質(zhì),比整數(shù)階導(dǎo)數(shù)更適合描述具有記憶和遺傳性質(zhì)的材料和過程.因此,分?jǐn)?shù)階微分方程更能準(zhǔn)確地刻畫許多自然界的現(xiàn)象,得到了越來越多的學(xué)者的關(guān)注.關(guān)于分?jǐn)?shù)階偏微分方程的數(shù)值解法無論對(duì)工程技術(shù)領(lǐng)域還是對(duì)數(shù)學(xué)本身都具有重要的價(jià)值.本文主要是在超收斂點(diǎn)處對(duì)時(shí)間分?jǐn)?shù)階波方程、多項(xiàng)時(shí)間分?jǐn)?shù)階波方程、時(shí)空分?jǐn)?shù)階Bloch-Torrey方程以及非線性時(shí)間分?jǐn)?shù)階四階反應(yīng)-擴(kuò)散方程等初邊值問題構(gòu)造數(shù)值解法,并給出相應(yīng)的理論分析.本文首先研究的是時(shí)間分?jǐn)?shù)階波方程的數(shù)值解法.對(duì)一維和二維時(shí)間分?jǐn)?shù)階波方程利用降階法得到等價(jià)的方程組,然后利用L2-1σ公式(Alikhanov,J.Comput.Phys.280(2015),424-438)對(duì)等價(jià)方程組建立時(shí)間方向二階精度,空間方向分別為二階和四階精度的有限差分格式.利用離散能量法,嚴(yán)格證明了格式在H1范數(shù)下的無條件穩(wěn)定性和收斂性.同時(shí)還給出了三維時(shí)間分?jǐn)?shù)階波方程的差分格式.數(shù)值算例驗(yàn)證了格式的計(jì)算精度和有效性.其次,對(duì)多項(xiàng)時(shí)間分?jǐn)?shù)階波方程建立時(shí)間二階精度的有限差分格式.利用降階法得到等價(jià)方程組,再對(duì)方程組中的多項(xiàng)時(shí)間分?jǐn)?shù)階導(dǎo)數(shù)在其超收斂點(diǎn)處離散,從而對(duì)多項(xiàng)時(shí)間分?jǐn)?shù)階波方程分別建立時(shí)間和空間方向都為二階精度的差分格式和時(shí)間二階、空間四階精度的差分格式.我們證明了兩個(gè)格式是唯一可解的,且在最大模下是無條件穩(wěn)定的和收斂的,收斂階分別為O(τ2 + h2)和O(τ2 + h4).數(shù)值實(shí)驗(yàn)表明格式的有效性,驗(yàn)證了差分格式的理論分析精度.隨后,討論了一維和二維時(shí)空分?jǐn)?shù)階Bloch-Torrey方程的差分方法.利用L2-1σ公式來離散時(shí)間分?jǐn)?shù)階Caputo導(dǎo)數(shù),分別應(yīng)用分?jǐn)?shù)階二階中心差分格式(C.Celik,M.Duman,J.Comput.Phys.231(2012),1743-1750.)和四階緊算子(X.Zhao,Z.Z.Sun,Z.P.Zhao,SIAM J.Sci.Comput.36(2014),A2865-A2886.)對(duì)空間分?jǐn)?shù)階Riesz導(dǎo)數(shù)進(jìn)行逼近,從而對(duì)一維和二維Bloch-Torrey方程構(gòu)造有限差分格式.同時(shí)我們給出了分?jǐn)?shù)階二階中心差分算子的權(quán)系數(shù)和的下界的一個(gè)估計(jì)式.利用離散能量法以及權(quán)系數(shù)的下界估計(jì)式,我們對(duì)格式的穩(wěn)定性和收斂性給出了嚴(yán)格的理論證明.對(duì)二維問題,我們還給出了兩個(gè)ADI格式來求解方程.數(shù)值算例驗(yàn)證了差分格式的有效性.最后一部分考慮了非線性時(shí)間分?jǐn)?shù)階四階反應(yīng)-擴(kuò)散方程的數(shù)值逼近.首先利用降階法,得到一個(gè)等價(jià)方程組,運(yùn)用L2-1σ公式對(duì)時(shí)間Caputo分?jǐn)?shù)階導(dǎo)數(shù)進(jìn)行離散,對(duì)空間整數(shù)階導(dǎo)數(shù)采用二階格式離散,進(jìn)而構(gòu)造一個(gè)三層線性化的有限差分格式.利用離散能量法,我們給出了格式在L2模下的無條件穩(wěn)定性和收斂性的嚴(yán)格的理論證明,收斂階為O(τ2 +h12+h22).對(duì)于差分格式的收斂性證明是理論分析的一個(gè)難點(diǎn),我們主要應(yīng)用二維網(wǎng)格函數(shù)空間的一個(gè)嵌入定理給出了格式的收斂性分析.數(shù)值實(shí)驗(yàn)驗(yàn)證了格式的理論分析的精度.
[Abstract]:In this paper , the finite difference scheme of time fractional order wave equation is established by means of discrete energy method . The numerical experiments show the validity of the format and verify the theoretical analysis accuracy of the difference scheme . Then , we discuss the difference method of the two - dimensional space - time fractional order ' s ( C . Celik , M . Duman , J . Comput . Phys . 231 ( 2012 ) , 1743 - 1750 . ) and the fourth - order compact operator ( X . Zhao , Z . Z . Sun , Z . P . Zhao , SIAM J . Sci . Comput . 36 ( 2014 ) , A2865 - A2886 . ) . In this paper , we give a rigorous theoretical proof for the stability and convergence of a two - dimensional finite difference scheme by using the discrete energy method and the lower bound estimation formula of the weight coefficients . The convergence of the difference scheme is a difficult problem for the theoretical analysis . We mainly apply one embedding theorem of the two - dimensional grid function space to give the convergence analysis of the format . The numerical experiments verify the accuracy of the theoretical analysis of the format .
【學(xué)位授予單位】:東南大學(xué)
【學(xué)位級(jí)別】:博士
【學(xué)位授予年份】:2016
【分類號(hào)】:O241.82
【相似文獻(xiàn)】
相關(guān)期刊論文 前10條
1 趙羽,蔡平,周敏東;分?jǐn)?shù)階Fourier變換的數(shù)值計(jì)算[J];哈爾濱工程大學(xué)學(xué)報(bào);2002年06期
2 王德金;鄭永愛;;分?jǐn)?shù)階混沌系統(tǒng)的延遲同步[J];動(dòng)力學(xué)與控制學(xué)報(bào);2010年04期
3 楊晨航,劉發(fā)旺;分?jǐn)?shù)階Relaxation-Oscillation方程的一種分?jǐn)?shù)階預(yù)估-校正方法[J];廈門大學(xué)學(xué)報(bào)(自然科學(xué)版);2005年06期
4 王發(fā)強(qiáng);劉崇新;;分?jǐn)?shù)階臨界混沌系統(tǒng)及電路實(shí)驗(yàn)的研究[J];物理學(xué)報(bào);2006年08期
5 夏源;吳吉春;;分?jǐn)?shù)階對(duì)流——彌散方程的數(shù)值求解[J];南京大學(xué)學(xué)報(bào)(自然科學(xué)版);2007年04期
6 張隆閣;;一類參數(shù)不確定混沌系統(tǒng)的分?jǐn)?shù)階自適應(yīng)同步[J];中國(guó)科技信息;2009年15期
7 陳世平;劉發(fā)旺;;一維分?jǐn)?shù)階滲透方程的數(shù)值模擬[J];高等學(xué)校計(jì)算數(shù)學(xué)學(xué)報(bào);2010年04期
8 辛寶貴;陳通;劉艷芹;;一類分?jǐn)?shù)階混沌金融系統(tǒng)的復(fù)雜性演化研究[J];物理學(xué)報(bào);2011年04期
9 黃睿暉;;分?jǐn)?shù)階微方程的迭代方法研究[J];長(zhǎng)春理工大學(xué)學(xué)報(bào);2011年06期
10 蔣曉蕓,徐明瑜;分形介質(zhì)分?jǐn)?shù)階反常守恒擴(kuò)散模型及其解析解[J];山東大學(xué)學(xué)報(bào)(理學(xué)版);2003年05期
相關(guān)會(huì)議論文 前10條
1 李西成;;經(jīng)皮吸收的分?jǐn)?shù)階藥物動(dòng)力學(xué)模型[A];中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)'2009論文摘要集[C];2009年
2 謝勇;;分?jǐn)?shù)階模型神經(jīng)元的動(dòng)力學(xué)行為及其同步[A];第四屆全國(guó)動(dòng)力學(xué)與控制青年學(xué)者研討會(huì)論文摘要集[C];2010年
3 張碩;于永光;王亞;;帶有時(shí)滯和隨機(jī)擾動(dòng)的不確定分?jǐn)?shù)階混沌系統(tǒng)準(zhǔn)同步[A];中國(guó)力學(xué)大會(huì)——2013論文摘要集[C];2013年
4 李常品;;分?jǐn)?shù)階動(dòng)力學(xué)的若干關(guān)鍵問題及研究進(jìn)展[A];中國(guó)力學(xué)大會(huì)——2013論文摘要集[C];2013年
5 李常品;;分?jǐn)?shù)階動(dòng)力學(xué)簡(jiǎn)介[A];第三屆海峽兩岸動(dòng)力學(xué)、振動(dòng)與控制學(xué)術(shù)會(huì)議論文摘要集[C];2013年
6 蔣曉蕓;徐明瑜;;時(shí)間依靠分?jǐn)?shù)階Schr銉dinger方程中的可動(dòng)邊界問題[A];中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)'2009論文摘要集[C];2009年
7 王花;;分?jǐn)?shù)階混沌系統(tǒng)的同步在圖像加密中的應(yīng)用[A];第二屆全國(guó)隨機(jī)動(dòng)力學(xué)學(xué)術(shù)會(huì)議摘要集與會(huì)議議程[C];2013年
8 王在華;;分?jǐn)?shù)階動(dòng)力系統(tǒng)的若干問題[A];第三屆全國(guó)動(dòng)力學(xué)與控制青年學(xué)者研討會(huì)論文摘要集[C];2009年
9 張碩;于永光;王莎;;帶有時(shí)滯和隨機(jī)擾動(dòng)的分?jǐn)?shù)階混沌系統(tǒng)同步[A];第十四屆全國(guó)非線性振動(dòng)暨第十一屆全國(guó)非線性動(dòng)力學(xué)和運(yùn)動(dòng)穩(wěn)定性學(xué)術(shù)會(huì)議摘要集與會(huì)議議程[C];2013年
10 李西成;;一個(gè)具有糊狀區(qū)的分?jǐn)?shù)階可動(dòng)邊界問題的相似解研究[A];中國(guó)力學(xué)大會(huì)——2013論文摘要集[C];2013年
相關(guān)博士學(xué)位論文 前10條
1 陳善鎮(zhèn);兩類空間分?jǐn)?shù)階偏微分方程模型有限差分逼近的若干研究[D];山東大學(xué);2015年
2 任永強(qiáng);油藏與二氧化碳埋存問題的數(shù)值模擬與不確定性量化分析以及分?jǐn)?shù)階微分方程的數(shù)值方法[D];山東大學(xué);2015年
3 蔣敏;分?jǐn)?shù)階微分方程理論分析與應(yīng)用問題的研究[D];電子科技大學(xué);2015年
4 卜紅霞;基于分?jǐn)?shù)階傅里葉域稀疏表征的CS-SAR成像理論與算法研究[D];北京理工大學(xué);2015年
5 楊變霞;分?jǐn)?shù)階Laplace算子的譜理論及其在微分方程中的應(yīng)用[D];蘭州大學(xué);2015年
6 邵晶;幾類微分系統(tǒng)的定性理論及其應(yīng)用[D];曲阜師范大學(xué);2015年
7 方益;分?jǐn)?shù)階Yamabe問題的一些緊性結(jié)果[D];中國(guó)科學(xué)技術(shù)大學(xué);2015年
8 王國(guó)濤;幾類分?jǐn)?shù)階非線性微分方程解的存在理論及應(yīng)用[D];西安電子科技大學(xué);2014年
9 陳明華;分?jǐn)?shù)階微分方程的高階算法及理論分析[D];蘭州大學(xué);2015年
10 尹學(xué)輝;基于分?jǐn)?shù)階PDE的圖像結(jié)構(gòu)保持型去噪算法研究[D];重慶大學(xué);2015年
相關(guān)碩士學(xué)位論文 前10條
1 黃志穎;非線性時(shí)間分?jǐn)?shù)階微分方程的數(shù)值解法[D];華南理工大學(xué);2015年
2 趙九龍;基于分?jǐn)?shù)階微積分的三維圖像去噪增強(qiáng)算法研究[D];寧夏大學(xué);2015年
3 楚彩虹;單載波分?jǐn)?shù)階傅里葉域均衡系統(tǒng)及關(guān)鍵技術(shù)研究[D];鄭州大學(xué);2015年
4 全曉靜;非線性分?jǐn)?shù)階積分方程的Adomian解法[D];寧夏大學(xué);2015年
5 黃潔;非線性分?jǐn)?shù)階Volterra積分微分方程的小波數(shù)值解法[D];寧夏大學(xué);2015年
6 莊嶠;復(fù)合介質(zhì)中時(shí)間分?jǐn)?shù)階熱傳導(dǎo)正逆問題及其應(yīng)用研究[D];山東大學(xué);2015年
7 高素娟;分?jǐn)?shù)階延遲偏微分方程的緊致有限差分方法[D];山東大學(xué);2015年
8 趙珊珊;時(shí)—空分?jǐn)?shù)階擴(kuò)散方程的快速算法以及MT-TSCR-FDE的快速數(shù)值解法[D];山東大學(xué);2015年
9 王珍;分?jǐn)?shù)階奇異邊值問題的研究[D];山東師范大學(xué);2015年
10 馮靜;一類分?jǐn)?shù)階奇異脈沖邊值問題正解的存在性研究[D];山東師范大學(xué);2015年
,本文編號(hào):1364287
本文鏈接:http://www.sikaile.net/shoufeilunwen/jckxbs/1364287.html