【摘要】：本文首先介紹非匹配網(wǎng)格上橢圓問題的混合有限元方法,在此基礎(chǔ)上我們提出了非匹配網(wǎng)格上的混合有限體積元方法,并給出相應(yīng)的格式.橢圓方程的Neumann邊值問題可以表示為：其中u和p分別表示流和壓力.考慮到介質(zhì)在Ω內(nèi)不同區(qū)域的性質(zhì)不同,在求解上述橢圓問題時(shí),我們有必要對(duì)區(qū)域Ω進(jìn)行分塊處理.假設(shè)開集Q可分割為若干非重疊的子區(qū)域Ωi,i=1,2,…,n,即Q是閉集∪i=1nΩi(?)R2的內(nèi)部,其中Ωi為開集Ωi的閉包.Γi是子區(qū)域Ωi的內(nèi)部邊界,即Γi=aΩi\aΩ.定義區(qū)域Ωi與Ωj的公共邊界為Γij,即Γij=Γi∩Γj,其中記Γii=0進(jìn)而有Γ=∪i,jΓij.為了保證非匹配邊界上p和u的連續(xù)性,在交界面上需要滿足如下連續(xù)性條件利用該條件,非匹配網(wǎng)格上的橢圓問題可以表示為：定義函數(shù)空間則上述橢圓方程的變分形式可以表示為：其中L2(Ωi)或(L2(Ωi))2上內(nèi)積記為(·,·)i,L22(Γi)上內(nèi)積記為的內(nèi)積記為(·,·ij在Ω上定義有限元空間：則非匹配網(wǎng)格上的混合有限元方法可以表示為：求uh∈vh,ph∈Wh,λh∈Λh使得對(duì)于i=1,2,…,n有其中uh,i=uh|Ωi,ph,i=ph|Ωi,λh,i=ph|Γi,這種非匹配網(wǎng)格上的混合有限元方法解的存在性和唯一性以及收斂性都已經(jīng)得到了證明(文獻(xiàn)[30]).類似于混合有限元方法,利用邊界上的連續(xù)性條件(4),我們給出了非匹配網(wǎng)格上的有限體積元方法.在有限體積元方法中,試探函數(shù)空間選為式(13)-(15)定義的函數(shù)空間Uh×Wh×Γh,其中Uh選取最低階Raviart-Thomas元,Wh,Γh選取分片常數(shù)空間,即可以看出如果(j,l)∈Iik,那么對(duì)偶單元Tjl就與Tik有公共邊.相應(yīng)的檢驗(yàn)函數(shù)空間為Vh×Wh×∧h,其中定義從試探函數(shù)空間到檢驗(yàn)函數(shù)空間的投影算子γh:Uh→Vh,表示為其中|Γik|代表三角形邊Γik的長度.顯然投影算子γh建立了Uh和Vh空間的一一對(duì)應(yīng)關(guān)系.空間Vh=γhUh同構(gòu)與空間Uh,進(jìn)而算子γh：Uh→Vh是有界的,一對(duì)一的.因此非匹配有限體積元方法可以表示為：求uh∈Uh,ph∈Wh,經(jīng)過數(shù)值實(shí)驗(yàn),結(jié)果顯示,非匹配三角形網(wǎng)格上混合有限體積元和混合有限體積元方法的收斂速度一致.實(shí)驗(yàn)結(jié)果同時(shí)說明非匹配網(wǎng)格上的有限體積元方法的格式的正確性.
【學(xué)位授予單位】：吉林大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2015
【分類號(hào)】：O241.82

【參考文獻(xiàn)】

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

1 ;GENERALIZED DIFFERENCE METHODS ON ARBITRARYQUADRILATERAL NETWORKS[J];Journal of Computational Mathematics;1999年06期

2 田明忠,陳仲英;橢圓型方程的廣義差分法(二次元)[J];高等學(xué)校計(jì)算數(shù)學(xué)學(xué)報(bào);1991年02期

3 陳仲英;;熱傳導(dǎo)方程的一種廣義差分方法[J];中山大學(xué)學(xué)報(bào)(自然科學(xué)版);1990年01期

4 李潛;一類非線性雙曲型方程的廣義Galerkin方法[J];計(jì)算數(shù)學(xué);1986年02期

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

6 祝丕琦,李榮華;二階橢圓偏微分方程的廣義差分法(Ⅱ)——四邊形網(wǎng)情形[J];高等學(xué)校計(jì)算數(shù)學(xué)學(xué)報(bào);1982年04期

本文編號(hào)：2763720