若干四階非線(xiàn)性偏微分方程的數(shù)值解法
本文關(guān)鍵詞: 四階非線(xiàn)性?huà)佄锓匠?差分法 擬譜方法 Fourier譜方法 長(zhǎng)時(shí)間行為 出處:《吉林大學(xué)》2017年博士論文 論文類(lèi)型:學(xué)位論文
【摘要】:高階非線(xiàn)性?huà)佄锓匠套鳛閿?shù)學(xué)模型可描述物理、化學(xué)、地理、環(huán)境科學(xué)等很多領(lǐng)域出現(xiàn)的現(xiàn)象,應(yīng)用廣泛,是非線(xiàn)性學(xué)科的一個(gè)重要組成部分.本文研究了三類(lèi)具有應(yīng)用背景的四階非線(xiàn)性偏微分方程的數(shù)值解法.理論上主要研究數(shù)值方法的先驗(yàn)估計(jì)和收斂性,特別是數(shù)值解的H~2模先驗(yàn)估計(jì),這是四階方程收斂性分析的基礎(chǔ).此外,文中也給出了數(shù)值實(shí)驗(yàn)以驗(yàn)證理論上的結(jié)果.首先,本文研究了描述晶體表面增長(zhǎng)的方程的差分法.關(guān)于該方程數(shù)值解法的研究不多,筆者曾經(jīng)和其他作者合作研究過(guò)該問(wèn)題的有限元法,但由于三次有限元形成剛度矩陣的計(jì)算量過(guò)大,因而又進(jìn)一步來(lái)研究方程的差分法.本文針對(duì)一維模型和二維模型各給出了一種差分格式,分別證明了差分解的存在唯一性、先驗(yàn)估計(jì)和收斂性,并給出了數(shù)值算例.該方程的非線(xiàn)性項(xiàng)是有理式形式,直接處理比較復(fù)雜,但通過(guò)一次分部積分可以大幅降低非線(xiàn)性項(xiàng)的復(fù)雜度.在有限元法的研究中是用連續(xù)內(nèi)積進(jìn)行分析,從而可以直接用分部積分處理非線(xiàn)性項(xiàng).但在構(gòu)造差分格式時(shí),如果只是簡(jiǎn)單移植,那么非線(xiàn)性項(xiàng)在離散內(nèi)積下不滿(mǎn)足分部積分.在不影響計(jì)算復(fù)雜度的前提下,為了能采取與有限元法相似的處理方法,本文針對(duì)該方程的一維模型構(gòu)造了一種半隱半顯型線(xiàn)性化差分格式,使得非線(xiàn)性項(xiàng)在相應(yīng)的離散內(nèi)積下可以進(jìn)行一次分部積分,從而降低非線(xiàn)性項(xiàng)的處理難度.在理論上證明了該數(shù)值格式在時(shí)間方向上為一階收斂,在空間方向上為二階收斂,數(shù)值算例驗(yàn)證了這一結(jié)果.針對(duì)二維模型本文構(gòu)造了 Crank-Nicolson型的差分格式,并定義了與之匹配的離散內(nèi)積,使得非線(xiàn)性項(xiàng)也可以進(jìn)行類(lèi)似的處理.所構(gòu)造的Crank-Nicolson型的差分格式雖然比線(xiàn)性化差分格式的計(jì)算量大,但在時(shí)間上達(dá)到二階收斂,綜合效果較好.數(shù)值結(jié)果也與該理論結(jié)果相吻合.其次,本文研究了二維擴(kuò)展的Fisher-Kolmogorov(EFK)方程u_t + γ△~2u-△u + f(u)= 0的擬譜方法.關(guān)于該方程的數(shù)值解法有較多的研究成果,方法涉及各類(lèi)有限元法和差分法.本文先對(duì)空間進(jìn)行離散,構(gòu)造了恰當(dāng)?shù)那蠼饪臻g和離散內(nèi)積,使得求解空間SN以及離散內(nèi)積的定義不僅與方程的初邊值條件相匹配,同時(shí)也滿(mǎn)足SN空間中基函數(shù)的正交性.在此基礎(chǔ)上,構(gòu)造了半離散擬譜近似格式,并證明了半離散解的存在性、先驗(yàn)估計(jì)及收斂性.進(jìn)而又建立了Crank-Nicolson型的全離散擬譜格式,證明了全離散解的先驗(yàn)估計(jì),并給出了按L2范數(shù)的最優(yōu)階誤差估計(jì).針對(duì)EFK方程的非線(xiàn)性項(xiàng)f(u)= u~3-u,在H~2半模先驗(yàn)估計(jì)的分析中引入能量泛函E_h(t)=γ/2||△u_h||~2+1/2||%絬_h||~2+(H(u_h),1)_h,u_h∈sN,其中H(u_h)=1/4(1-u_h~2)~2,并且H'(u_h)=f(u_h).利用能量泛函隨時(shí)間遞減的性質(zhì)證得H~2半模的先驗(yàn)估計(jì).全離散解的先驗(yàn)估計(jì)證明中也構(gòu)造和借助了相應(yīng)的離散能量泛函.給出了 Fourier擬譜方法的收斂性結(jié)果,并用數(shù)值算例驗(yàn)證了相應(yīng)的理論結(jié)果.最后,本文研究了一類(lèi)對(duì)流Cahn-Hilliard方程u_t+γ△~2u-△f(u)%絞(u)=0的Fourier譜近似解的長(zhǎng)時(shí)間行為,其f(u)為關(guān)于u的四次多項(xiàng)式,g(u)為關(guān)于u的五次多項(xiàng)式.文中給出了譜近似解的H~2模長(zhǎng)時(shí)間先驗(yàn)估計(jì)的證明,并對(duì)有限時(shí)間段分析了數(shù)值解的收斂性.數(shù)值解的先驗(yàn)估計(jì)是理論分析的基礎(chǔ),但由于該方程非線(xiàn)性項(xiàng)的特殊性,得到其數(shù)值解的H0模先驗(yàn)估計(jì)是難點(diǎn).為此本文針對(duì)非線(xiàn)性項(xiàng)%絞(u)給出了一種反對(duì)稱(chēng)格式,并利用該格式構(gòu)造了半離散譜近似格式,進(jìn)而建立了全離散格式.由于全離散格式中g(shù)(u)利用了反對(duì)稱(chēng)格式,所以有效地排除了該非線(xiàn)性項(xiàng)的影響,得到了 H0模的長(zhǎng)時(shí)間先驗(yàn)估計(jì).數(shù)值解H~2半模的先驗(yàn)估計(jì)中通常會(huì)借助于能量泛函,但該方程和H~2半模對(duì)應(yīng)的能量泛函難以找到.由于有了H0模的先驗(yàn)估計(jì),所以本文在H~2半模長(zhǎng)時(shí)間先驗(yàn)估計(jì)中,對(duì)非線(xiàn)性項(xiàng)的影響,采用充分地利用四階導(dǎo)數(shù)項(xiàng)和H0模的估計(jì),并多次使用Nirenberge不等式的分析方法,得到了較好的理論結(jié)果.
[Abstract]:In this paper , we study the finite element method for solving the nonlinear term . In this paper , we study the finite element method for solving the nonlinear term by using the finite element method . In this paper , we have given a priori estimate of the H ~ 2 half - norm by means of the energy functional function .
【學(xué)位授予單位】:吉林大學(xué)
【學(xué)位級(jí)別】:博士
【學(xué)位授予年份】:2017
【分類(lèi)號(hào)】:O241.82
【相似文獻(xiàn)】
相關(guān)期刊論文 前10條
1 林志瑗;場(chǎng)波問(wèn)題的一些數(shù)值解法[J];電子器件;1980年03期
2 路新玲;張啟敏;;固定資產(chǎn)投資系統(tǒng)數(shù)值解的收斂性[J];寧夏大學(xué)學(xué)報(bào)(自然科學(xué)版);2007年04期
3 何劍;;隨機(jī)Kuramoto-Sivashinsky方程數(shù)值解的收斂性[J];太原師范學(xué)院學(xué)報(bào)(自然科學(xué)版);2009年04期
4 林鵬程;關(guān)于二維薄片振動(dòng)問(wèn)題的數(shù)值解[J];福州大學(xué)學(xué)報(bào);1961年00期
5 耿兆銓;非恒定水流數(shù)值解空間失穩(wěn)問(wèn)題分析[J];水利學(xué)報(bào);1985年07期
6 崔耀東 ,林兆榮;鈑料壓延數(shù)值解[J];南京航空航天大學(xué)學(xué)報(bào);1986年01期
7 楊年甫;張質(zhì)文;;光貼片實(shí)驗(yàn)全場(chǎng)應(yīng)力的數(shù)值解法[J];實(shí)驗(yàn)力學(xué);1991年01期
8 徐永權(quán);微分方程的數(shù)值解法[J];天津理工學(xué)院學(xué)報(bào);1996年02期
9 閔弘,章光,王水林,王元漢;一種提高數(shù)值解精度的方法探討[J];巖石力學(xué)與工程學(xué)報(bào);2002年01期
10 項(xiàng)林川;一個(gè)追線(xiàn)問(wèn)題的數(shù)值解法及推廣[J];大學(xué)物理;2004年08期
相關(guān)會(huì)議論文 前10條
1 顧本立;顧寧;;二維逆散射成象數(shù)值解[A];1987年全國(guó)微波會(huì)議論文集(上)[C];1987年
2 石文靜;孟慶春;;含穿透分層雙金屬梁數(shù)值解與片條合成能量解法的比較[A];中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)'2009論文摘要集[C];2009年
3 楊吉生;徐立軍;孫天宇;葉聲華;;非常規(guī)形體散射近場(chǎng)的數(shù)值解法[A];1995年全國(guó)微波會(huì)議論文集(上冊(cè))[C];1995年
4 吳強(qiáng);;GSM(1,1)模型的結(jié)構(gòu)及數(shù)值解[A];數(shù)學(xué)及其應(yīng)用文集——中南模糊數(shù)學(xué)和系統(tǒng)分會(huì)第三屆年會(huì)論文集(上卷)[C];1995年
5 張智勇;陳予恕;;VRKF方法在求解某類(lèi)非線(xiàn)性常微分方程數(shù)值解中的應(yīng)用[A];第九屆全國(guó)動(dòng)力學(xué)與控制學(xué)術(shù)會(huì)議會(huì)議手冊(cè)[C];2012年
6 譚建國(guó);王洪禮;;泊松白噪聲激勵(lì)的隨機(jī)動(dòng)力學(xué)系統(tǒng)的數(shù)值解法[A];第八屆全國(guó)動(dòng)力學(xué)與控制學(xué)術(shù)會(huì)議論文集[C];2008年
7 張?zhí)煸?莫俊永;;核三體Faddeev方程數(shù)值解[A];第八屆全國(guó)核物理會(huì)議文摘集(上冊(cè))[C];1991年
8 劉心鑄;吳宣志;董英君;;直流穩(wěn)定電場(chǎng)的數(shù)值解法及其應(yīng)用[A];中國(guó)地質(zhì)科學(xué)院礦床地質(zhì)研究所文集(2)[C];1981年
9 翟俊輝;陳忠安;;平面彈性接觸問(wèn)題的一種數(shù)值解法[A];第14屆全國(guó)結(jié)構(gòu)工程學(xué)術(shù)會(huì)議論文集(第一冊(cè))[C];2005年
10 盧文全;;關(guān)于光線(xiàn)方程數(shù)值解的精度問(wèn)題[A];面向21世紀(jì)的科技進(jìn)步與社會(huì)經(jīng)濟(jì)發(fā)展(上冊(cè))[C];1999年
相關(guān)博士學(xué)位論文 前9條
1 劉鳳楠;若干四階非線(xiàn)性偏微分方程的數(shù)值解法[D];吉林大學(xué);2017年
2 梁芬;一些反問(wèn)題的數(shù)值解法研究[D];汕頭大學(xué);2010年
3 曾光;橢圓型方程的數(shù)值解法與穩(wěn)定性分析[D];電子科技大學(xué);2011年
4 邵和助;Schr(?)dinger方程的數(shù)值解法[D];上海大學(xué);2010年
5 鄒戰(zhàn)勇;一類(lèi)離散HJB方程的數(shù)值解法[D];湖南大學(xué);2008年
6 趙曉朋;高階非線(xiàn)性?huà)佄锓匠探獾男再|(zhì)及其數(shù)值解法[D];吉林大學(xué);2013年
7 孫越;廣義橢球函數(shù)解析解與數(shù)值解的研究[D];北京郵電大學(xué);2015年
8 陳宇;若干反問(wèn)題模型的數(shù)值求解[D];上海交通大學(xué);2009年
9 白冬梅;非線(xiàn)性Schr(?)dinger方程及其相關(guān)耦合方程的數(shù)值解法[D];南京航空航天大學(xué);2010年
相關(guān)碩士學(xué)位論文 前10條
1 左一凡;無(wú)限時(shí)區(qū)正倒向隨機(jī)微分方程及永續(xù)債券定價(jià)的數(shù)值解法[D];復(fù)旦大學(xué);2014年
2 謝夢(mèng)玲;一類(lèi)非線(xiàn)性多延遲偏微分方程的數(shù)值解法[D];華中科技大學(xué);2014年
3 穆錦榮;改進(jìn)的數(shù)值方法及簡(jiǎn)支梁?jiǎn)栴}的數(shù)值解[D];青海民族大學(xué);2016年
4 吳曉婷;木材熱壓及干燥過(guò)程的兩種數(shù)值方法[D];東北林業(yè)大學(xué);2016年
5 李盤(pán)潤(rùn);隨機(jī)固定資產(chǎn)模型數(shù)值解[D];北方民族大學(xué);2016年
6 李運(yùn)利;基于Bernstein多項(xiàng)式和樣條函數(shù)的高階常微分方程數(shù)值解法研究[D];合肥工業(yè)大學(xué);2010年
7 白雪;關(guān)于Gilson-Pickering方程數(shù)值解及分支研究[D];哈爾濱工業(yè)大學(xué);2012年
8 王倩;年齡相關(guān)的單一種群動(dòng)態(tài)模型數(shù)值解法[D];新疆大學(xué);2005年
9 李希文;馬爾可夫調(diào)制隨機(jī)脈沖微分方程數(shù)值解的收斂性與穩(wěn)定性[D];哈爾濱工業(yè)大學(xué);2009年
10 周曼;一類(lèi)熱導(dǎo)率的反演數(shù)值解法[D];浙江大學(xué);2007年
,本文編號(hào):1464793
本文鏈接:http://sikaile.net/shoufeilunwen/jckxbs/1464793.html