有限體積法介紹_第1頁(yè)
有限體積法介紹_第2頁(yè)
有限體積法介紹_第3頁(yè)
有限體積法介紹_第4頁(yè)
有限體積法介紹_第5頁(yè)
已閱讀5頁(yè),還剩3頁(yè)未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、有限體積法1 有限體積法基本原理上一章講到的有限差分法將數(shù)值網(wǎng)格的節(jié)點(diǎn)上定義為計(jì)算節(jié)點(diǎn),并在網(wǎng)格節(jié)點(diǎn)上對(duì)微分形式的流體基本方程進(jìn)行離散,用網(wǎng)格節(jié)點(diǎn)上的物理量的代數(shù)方程作為原PDE的近似。在本章所要學(xué)習(xí)的有限體積法則采用了不同的離散形式。首先,有限體積法離散的是積分形式的流體力學(xué)基本方程:(1)計(jì)算域用數(shù)值網(wǎng)格劃分成若干小控制體。和有限差分法不同的是,有限體積法的網(wǎng)格定義了控制體的邊界,而不是計(jì)算節(jié)點(diǎn)。有限體積法的計(jì)算節(jié)點(diǎn)定義在小控制體內(nèi)部。一般有限體積法的計(jì)算節(jié)點(diǎn)有兩種定義方法,一種是將網(wǎng)格節(jié)點(diǎn)定義在控制體的中心,另一種方法中,相鄰兩個(gè)控制體的計(jì)算節(jié)點(diǎn)到公共邊界的距離相等。第一種方法的優(yōu)點(diǎn)在于

2、用計(jì)算節(jié)點(diǎn)的值作為控制體上物理量的平均值具有二階的精度;第二種方法的好處是在控制體邊界上的中心差分格式具有較高的精度。積分形式的守恒方程在小控制體和計(jì)算域上都是成立的。為了獲得每一個(gè)控制體上的代數(shù)方程,面積分和體積分需要用求面積公式來(lái)近似。2 面積分的近似采用結(jié)構(gòu)化網(wǎng)格,在二維情況下,每一個(gè)控制體有4個(gè)面,二維情況,每一個(gè)控制體有6個(gè)表面。計(jì)算節(jié)點(diǎn)用大寫(xiě)字母表示,控制體邊界和節(jié)點(diǎn)用小寫(xiě)字母表示。為了保證守恒性,控制體不能重疊,每一個(gè)面都是相鄰兩個(gè)控制體的唯一公共邊界??刂企w邊界上的積分等于控制體個(gè)表面的積分的和:(2)上式中,f可以表示或。PENWSNNEEWWSSNWNESWSEewnsne

3、nwsesw顯然,為了獲得邊界上的積分,必須知道f在邊界上的詳細(xì)分布情況,這是不可能實(shí)現(xiàn)的,由于只是計(jì)算節(jié)點(diǎn)上的函數(shù)值,因此必須采用近似的方法來(lái)計(jì)算積分。整個(gè)近似過(guò)程分成兩步第一步:用邊界上幾個(gè)點(diǎn)的近似積分公式第二步:邊界點(diǎn)上的函數(shù)值用計(jì)算節(jié)點(diǎn)函數(shù)值的插值函數(shù)近似面積分可采用以下不同精度的積分公式:二階精度積分:(3)上式中為邊界中點(diǎn)出的函數(shù)值。近似為方格中心點(diǎn)的值乘以方格的面積。三階精度積分:(4)四階精度積分:(5)應(yīng)該注意的是,采用不同精度的積分公式,在相應(yīng)的邊界點(diǎn)的插值時(shí)也應(yīng)采用相應(yīng)精度的插值函數(shù)。積分公式的精度越高,近似公式就越復(fù)雜。3 體積分的近似和面積分相似,體積分也有不同精度的

4、近似公式二階精度積分公式(6)采用雙二次樣條函數(shù)(7)可以得到四階精度的積分公式:(8)4 函數(shù)的插值在上節(jié)講到的積分的近似公式中用到了非計(jì)算節(jié)點(diǎn)上的函數(shù)值,被積函數(shù)f中包含了多個(gè)物理量及其偏微分,如對(duì)流項(xiàng),擴(kuò)散項(xiàng),在源項(xiàng)中也有類(lèi)似情況,這里假定流場(chǎng)和流體的物性參數(shù)是已知的,物理量及其偏導(dǎo)數(shù)在控制面上的值需要通過(guò)計(jì)算節(jié)點(diǎn)上物理量的插值得到。下面已e面為例進(jìn)行討論。4.1 迎風(fēng)插值(UDS)用上游計(jì)算節(jié)點(diǎn)的函數(shù)值近似相當(dāng)于對(duì)一階偏導(dǎo)數(shù)采用迎風(fēng)格式,因此用UDS來(lái)表示這種近似方法,在UDS中:(9)UDS是唯一無(wú)條件滿(mǎn)足有界性要求的近似格式,在數(shù)值過(guò)程中不會(huì)產(chǎn)生數(shù)值振蕩。UDS存在數(shù)值粘性。根據(jù)T

5、aylor公式,該格式具有一階精度,并具有數(shù)值粘性:(10)在多維問(wèn)題中,如果流動(dòng)方向和網(wǎng)格是斜交的,截?cái)嗾`差會(huì)在垂直于流動(dòng)方向以及流線(xiàn)方向產(chǎn)生擴(kuò)散,這是一種非常嚴(yán)重的誤差,函數(shù)的峰值或函數(shù)值的快速變化會(huì)被抹平,為了得到高精度結(jié)果需要采用非常精細(xì)的網(wǎng)格。4.2 線(xiàn)性插值(CDS)(11)(12)線(xiàn)性插值具有二階精度,線(xiàn)性插值相當(dāng)于FDM中的CDS格式,因此用CDS表示。CDS格式會(huì)產(chǎn)生數(shù)值振蕩。對(duì)于擴(kuò)散項(xiàng)(13)4.3 三階迎風(fēng)格式(QUICK)和UDS類(lèi)似,QUICK格式也和流動(dòng)方向有關(guān)(14)其中:;(15a);(15b)4.4 高階格式(4階精度CDS)采用三次曲線(xiàn)可擬合出四階精度的中心

6、插值公式,在均勻網(wǎng)格中,四階公式為:(16)(17)5 邊界的處理對(duì)于對(duì)流項(xiàng),在入口處一般給出了流量或函數(shù)值,在邊界和對(duì)稱(chēng)面上流量為零,在出口處假設(shè)和出口的法向坐標(biāo)無(wú)關(guān),因此可采用迎風(fēng)格式。對(duì)于擴(kuò)散項(xiàng)則可能需要采用偏心格式。6 有限體積法應(yīng)用舉例例:考慮一標(biāo)量在已知流場(chǎng)中的輸運(yùn)過(guò)程(如圖4.4所示),輸運(yùn)方程為:(18)邊界條件:;北部入口邊界;西部壁面邊界對(duì)稱(chēng)條件;南部邊界梯度為0;東部出口條件,流線(xiàn)方程f(y)壁面對(duì)稱(chēng)邊界入口,f0出口流線(xiàn),xy=c對(duì)流項(xiàng):(19)為質(zhì)量通量。(20)若采用UDS格式,代數(shù)方程組中各項(xiàng)系數(shù)為:;(21)若采用CDS格式,代數(shù)方程組中各項(xiàng)系數(shù)為;(22)根據(jù)

7、連續(xù)性方程:(23)相鄰CV之間的關(guān)系:;(24)其余相鄰CV有類(lèi)似關(guān)系擴(kuò)散項(xiàng)采用CDS格式(25)代數(shù)方程組中擴(kuò)散項(xiàng)系數(shù)為:;(26)對(duì)于任意控制體(27),l為任意指標(biāo)P,E,W,S,N。(28)邊界條件的處理:對(duì)于西部和北部邊界,由于給定了函數(shù)值,對(duì)流項(xiàng)可直接代入函數(shù)值而無(wú)需插值,擴(kuò)散項(xiàng)則采用一側(cè)差分(29)這里,W點(diǎn)和P的w邊中點(diǎn)重合。南邊和西邊的梯度為零,以南邊為例,由于梯度為零,代數(shù)方程變?yōu)椋?30)6 SIMPLE方法考慮定常不可壓流動(dòng)問(wèn)題,控制方程為:連續(xù)性方程:(31)動(dòng)量方程:(32)不可壓縮問(wèn)題求解的困難在于壓力場(chǎng)的求解。主要原因在于壓力p沒(méi)有獨(dú)立的方程組。先考慮一維問(wèn)題

8、:對(duì)于動(dòng)量方程:(33)若采用CDS格式簡(jiǎn)化后得:(34)根據(jù)連續(xù)性方程,則有,由于相鄰節(jié)點(diǎn)之間的壓力沒(méi)有聯(lián)系方程,容易造成壓力交錯(cuò)現(xiàn)象。為了解決這一問(wèn)題,可采用交錯(cuò)網(wǎng)格技術(shù),即速度場(chǎng)和壓力場(chǎng)采用不同的網(wǎng)格。以二維問(wèn)題為例,交錯(cuò)網(wǎng)格的布置如下圖所示:PEWNSueuwvnvs主控制體為壓力控制體(黑色實(shí)線(xiàn)網(wǎng)格),u的控制體(紅色虛線(xiàn)網(wǎng)格)的計(jì)算節(jié)點(diǎn)在主控制體的e邊,控制體的e,w邊界通過(guò)主控制體的計(jì)算節(jié)點(diǎn),v控制體(藍(lán)色雙點(diǎn)劃線(xiàn)網(wǎng)格)的計(jì)算節(jié)點(diǎn)在主控制體的n邊,該控制體的n,s面經(jīng)過(guò)主控制體的計(jì)算節(jié)點(diǎn)。在u的控制體中,采用有限體積法離散可得u的代數(shù)方程:(35)壓力場(chǎng)的求解采用壓力校正方法。即采用預(yù)估的壓力場(chǎng)求速度,再用連續(xù)性方程校正壓力場(chǎng)。當(dāng)連續(xù)性方程得到滿(mǎn)足時(shí),壓力場(chǎng)就是真實(shí)的壓力場(chǎng)。具體步驟如下1 預(yù)測(cè)壓力場(chǎng)2 將預(yù)測(cè)壓力場(chǎng)代入動(dòng)量方程,分別求解速度場(chǎng)(36)3 用連續(xù)性方程校正壓力設(shè)方程的精確解為u,v,p;(37)其

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論