有限體積法(FVM)是一種重要的求解偏微分方程的數(shù)值離散方法.由于其能夠保持物理量局部積分的守恒性,該方法在流體力學(xué)、地質(zhì)學(xué)等領(lǐng)域被廣泛的應(yīng)用.本文主要研究一維三次元以及二維雙三次元有限體積法的收斂性規(guī)律,從理論(一維)和數(shù)值實(shí)驗(yàn)(一維和二維)上,給出了當(dāng)對(duì)偶剖分滿足對(duì)稱性變化時(shí),相應(yīng)格式按L2模的收斂階及超收斂的結(jié)果.首先,考慮兩點(diǎn)邊值問(wèn)題:其中 I =[a,b],p ∈ C1(I),pmax ≥ p(x)≥pmin0,f ∈ L2(I).對(duì)區(qū)間I =[a,b]作剖分Tb,節(jié)點(diǎn)取為a = x0x1……xn,=b,[xi-1,xi]為任一單元(i= 1,2,…,n).試探函數(shù)空間Uh為相應(yīng)于Th的Lagrange型三次元有限元空間.相應(yīng)的計(jì)算節(jié)點(diǎn)(即Lagrange插值節(jié)點(diǎn))取為每個(gè)單元的三等分點(diǎn)(含端點(diǎn)),記作 xi-k/3(k= 1,2).相應(yīng)于Th的對(duì)偶剖分Th*,其節(jié)點(diǎn)取為a =x0xa*x1/2*x1-α*x1……xixi-α*xi+1/2*xi-1……xn=b.其中xi+α*=α(xi+1-xi)+xi,xi+1-α*=(1-α)(xi+1-xi)+ xi,xi+1/2*=1/2(xi+1xi),i=1,2,…,n.且0α1/3為用于調(diào)節(jié)對(duì)偶剖分位置的參數(shù).檢驗(yàn)函數(shù)空間Vh為相應(yīng)于Th*的分片常數(shù)函數(shù)空間.與試探函數(shù)空間Uh的計(jì)算節(jié)點(diǎn)xj,xj-1/3,xj-2/3,(j=1,…,n)相應(yīng)的對(duì)偶單元分別為「xj-α*,xj+α*」,「xj-1/2*,xj-α*」,「xj-1+α*,xj-1/2*」.一維三次元有限體積法格式為:求uh∈Uh,使得ah(uh,vh)=(f,vh),(?)vh∈ Vh......(1.2)相應(yīng)的單元形式為其中Φj(x),Φj-1/3(x),Φj-2/3(x)為與xj,xj-1/3,xj-2/3相應(yīng)對(duì)偶單元上的特征函數(shù),{Φj(x),Φj-1/3(x),Φj-2/3(x):1 ≤j≤n為檢驗(yàn)函數(shù)空間Vh的一組基底.定義П*為Uh到Vh的投影算子(詳細(xì)定義見(jiàn)2.3節(jié)),則(1.2)式等價(jià)于αh(uh,Пhuh)=(f,]Пh*,uh),(?)uh∈Uh,進(jìn)而我們得到了格式的穩(wěn)定性及H1誤差估計(jì).定理1(穩(wěn)定性)當(dāng) 0α1、3且h充分小時(shí),離散雙線性形式ah(uh,Пh*uh)是正定的,即存在與子空間Uh無(wú)關(guān)的常數(shù)γ0,使得ah(uh,Пh*uh)≥γ‖uh‖12(?)uh ∈ Uh.定理2(H1誤差估計(jì))設(shè)問(wèn)題(1.1)的解u∈H4(I),uh是三次元有限體積格式的解,則0α1/3且h充分小時(shí)有如下誤差估計(jì):‖u-uh‖≤Ch3|u|4.我們擴(kuò)展了正交性條件[13]在一維的限制,并證明了所有具有對(duì)稱對(duì)偶結(jié)構(gòu)的一維三次有限體積格式均有最佳的按L2模收斂階.定理3(L2誤差估計(jì))設(shè)問(wèn)題(1.1)的解u∈HE1(I)∩H5(I),uh是三次元有限體積格式的解,則當(dāng)0α1/3且h充分小時(shí)有如下誤差估計(jì):‖u-uh‖0≤Ch4 |u|5.我們通過(guò)數(shù)值實(shí)驗(yàn),驗(yàn)證了上述收斂性結(jié)果.同時(shí)獲得了如下超收斂的結(jié)論.結(jié)論1:在均勻的網(wǎng)格剖分Th上,對(duì)偶剖分隨α對(duì)稱性變化時(shí),單元端點(diǎn)及中點(diǎn)處(不含邊界單元)數(shù)值解導(dǎo)數(shù)的平均值與精確解導(dǎo)數(shù)有如下超收斂結(jié)果.結(jié)論2:在均勻的網(wǎng)格剖分Th上,整體超逼近性質(zhì)僅在對(duì)偶剖分節(jié)點(diǎn)取為應(yīng)力佳點(diǎn)時(shí)成立,即‖uI-uh‖=O(h4).其中uI為u的分片三次Lagrange插值.再考慮二維Poisson方程:其中 Ω={(x,y)|a≤x≤b,c≤y≤d},αΩ為Ω的邊界.Th為Ω上的均勻矩形網(wǎng)格剖分,試探函數(shù)空間Uh為相應(yīng)于Th的Lagrange型雙三次元有限元空間.對(duì)于每個(gè)矩形單元,將每條邊如一維情形按α進(jìn)行剖分,連結(jié)對(duì)邊上的對(duì)偶剖分節(jié)點(diǎn),我們得到了每個(gè)矩形單元內(nèi)的對(duì)偶剖分.即與Th相應(yīng)的對(duì)偶剖分Th*,檢驗(yàn)函數(shù)空間Vh為相應(yīng)于Th*的分片常數(shù)函數(shù)空間.二維雙三次元有限體積法格式為:求uh∈Uh,使得αh(uh,vh)=(f,vh),(?)vh∈Vh.相應(yīng)的單元形式為其中N為所有計(jì)算節(jié)點(diǎn),ΦP為相應(yīng)于P的對(duì)偶單元KP*上的特征函數(shù).我們通過(guò)數(shù)值實(shí)驗(yàn),獲得了如下結(jié)論.結(jié)論1:對(duì)偶剖分隨α對(duì)稱性變化時(shí),H1與L2誤差估計(jì)均能達(dá)到最佳收斂階,即‖u-uh‖1= O(h3),‖u-uh‖n=O(h4).結(jié)論2:在均勻的網(wǎng)格剖分Th上,對(duì)偶剖分隨α對(duì)稱性變化時(shí),在對(duì)稱點(diǎn)處(各單元的頂點(diǎn)、邊中點(diǎn)以及中心點(diǎn))的數(shù)值解梯度的平均值與精確解梯度有如下超收斂結(jié)果.結(jié)論3:在均勻的網(wǎng)格剖分Th上,整體超逼近性質(zhì)僅在對(duì)偶剖分節(jié)點(diǎn)取為應(yīng)力佳點(diǎn)時(shí)成立,即‖uI-uh‖=O(h4).其中u1為u的雙三次Lagrange插值.
【學(xué)位單位】：吉林大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位年份】：2018
【中圖分類】：O241.82
【文章目錄】：
摘要
Abstract
第1章 緒論
第2章 一維Lagrange型三次元
    2.1 試探函數(shù)空間和檢驗(yàn)函數(shù)空間
    2.2 三次元有限體積法的一般格式
1誤差估計(jì)'>    2.3 穩(wěn)定性和H¹誤差估計(jì)
    2.4 正交性條件
2誤差估計(jì)'>    2.5 L²誤差估計(jì)
    2.6 數(shù)值實(shí)驗(yàn)
第3章 二維Lagrange型雙三次元
    3.1 試探函數(shù)空間和檢驗(yàn)函數(shù)空間
    3.2 雙三次元有限體積法的一般格式
    3.3 數(shù)值實(shí)驗(yàn)
第4章 結(jié)論和展望
參考文獻(xiàn)
致謝

【參考文獻(xiàn)】

相關(guān)期刊論文 前5條

1 于長(zhǎng)華;王曉玲;李永海;;解兩點(diǎn)邊值問(wèn)題的一類修改的三次有限體積元法[J];計(jì)算數(shù)學(xué);2010年04期

2 于長(zhǎng)華;李永海;;解兩點(diǎn)邊值問(wèn)題的基于應(yīng)力佳點(diǎn)的二次有限體積元法[J];吉林大學(xué)學(xué)報(bào)(理學(xué)版);2009年04期

3 倪平,吳微;廣義Galerkin方法的超收斂估計(jì)[J];高等學(xué)校計(jì)算數(shù)學(xué)學(xué)報(bào);1986年02期

4 吳微,李榮華;解一維二階橢圓和拋物型微分方程的廣義差分法[J];數(shù)學(xué)年刊A輯(中文版);1984年03期

5 向新民;解兩點(diǎn)邊值問(wèn)題的廣義差分法——Lagrange二次元[J];黑龍江大學(xué)自然科學(xué)學(xué)報(bào);1982年02期

相關(guān)博士學(xué)位論文 前1條

1 陳艷利;二階橢圓型方程有限體積法的若干研究[D];吉林大學(xué);2014年

本文編號(hào)：2862033