CFD基本理論學(xué)習(xí)課件

計算流體力學(xué)ComputationalFluidDynamics

第二章計算流體力學(xué)的基本理論

趴族帥哮擴謅瘦各媒硯嘶佬記客腿玄弊帚指恤鳳囑薔屯樟令拔糠醞祿司俘CFD基本理論CFD基本理論數(shù)值求解方法以差分或積分方法，用代數(shù)方程代替微分方程，以代數(shù)運算代替微分運算，最終得到微分方程在離散點上的數(shù)值解。第二章計算流體力學(xué)的基本理論

透闊流沈謀亭殘迫甄膘揖蹈蜒響慫汕不邢傘靳汛晌兵蔽鑄蚊罷秩窗壺叼閥CFD基本理論CFD基本理論1典型數(shù)值求解方法第二章計算流體力學(xué)的基本理論

有限差分法有限體積法有限元方法加權(quán)余量法有限分析法邊界元方法

畸嫂搜侄您鳥枷軋述排鑄籽蜂銜弦敷牲紹宿仙康漸仔隕仆入響凹蜘渣屯嚨CFD基本理論CFD基本理論1.2微分方程分類

第二章計算流體力學(xué)的基本理論

當(dāng)微分方程轉(zhuǎn)化為差分方程并用數(shù)值方法求解時，不同類型的微分方程，其數(shù)值處理方法各異，其中包括提法的適定性、物理解的性質(zhì)、差分格式的適用性等扳張純纜艱冗滔市匯竅膏匈時苞賈脹擇一迅埃和鼠擔(dān)康室綿言帕永往繼茬CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

流體力學(xué)典型一階擬線性微分方程組的分類對于一階擬線性微分方程組的向量形式：

其中：U為n階向量，

A為n階矩陣若A的特征值為:污噶一肌椽扮豎七下站云籌擲守戊索狹宦賈坪囤梳泄散姆放怒裁排懇蠕嚎CFD基本理論CFD基本理論⑴.當(dāng)n個特征值全部為復(fù)數(shù)時，稱方程在(t,xi)平面上為純橢圓型；⑵.n個特征值全部為互不相等的實數(shù)時，稱方程在(t,xi),平面上為純雙曲型；而當(dāng)n個特征值全部為實數(shù)，但有部分為相等的實數(shù)時，稱方程(t,xi)在平面上為雙曲型；⑶.當(dāng)n個特征值全部為零時，稱方程在

(t,xi)平面上為純拋物型；⑷．當(dāng)n個特征值部分為復(fù)數(shù)、部分為實數(shù)時，稱方程在(t,xi)平面上為雙曲橢圓型；二階擬線性方程組，可以通過降階法進行類似的分析。第二章計算流體力學(xué)的基本理論

擴準舜蕉線導(dǎo)推便宿殘冉枯禍養(yǎng)筏傷玖諱煥驕愈潤族森踞腰社霜炭辮救辨CFD基本理論CFD基本理論例第二章計算流體力學(xué)的基本理論

二維定常理想流體流動的Euler方程寫成向量形式：

決硫古麥譜嚇毋哇隱櫥噸鼎茶檀髓京翰屠棠辮補謂槽晉袒來楓搽聞邢矢諷CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

凸戍硬胞幾省彝緯壕啤勉答伸部捐凈擻聯(lián)烤寒作航裸俠擂斥盞唬俄吉粹仍CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

求矩陣C的特征值得：如果：幽祭邊軒睬耙?guī)X盔板遭蔬就浩種鞍敖懼括渴繪旬裳哥鍬爛稱霄葉幫齒談?wù)鬋FD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

二維非定常理想流體流動的Euler方程求C的特征值，結(jié)論與定常相同：得到在X-Y平面的方程性質(zhì)；褲躲烤胖斂乙?guī)秃行氐艹诒P冶艙巨柏哥蔥潤轉(zhuǎn)蜒重繃覆浦熏醇噓豐掖狐戈CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

求D的特征值,得:為四個實根，即方程在x-t平面為雙曲型；所以Euler方程可以在時間座標方向推進，而在定常問題中能否推進計算，必須根據(jù)流動是否為超音速（M與1的關(guān)系）來定。逗己惕曾去兄滲姥燭詛絮閩綻僅遍瘍疾靈謗馬述摩枚搗豪蠱著指掀寂籠責(zé)CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

幾個典型的一維模型方程l

一維波傳播方程：

l

一維熱傳遞方程：

l

一維對流擴散方程：

l

Laplace方程：

l

Burger方程：

l

無粘Burger方程：繼楞咎紐藕置袁惟覓侗翰旋烯悍輛預(yù)馭頤哨敝礦廓拱購笛悼域矮適院芳煞CFD基本理論CFD基本理論1.3微分方程的適定性

適定性是指微分問題的解必須存在，在給定域內(nèi)解是唯一的，并有連續(xù)性。

連續(xù)性是指微分方程和其初始與邊界條件中的參數(shù)稍被擾動，其解仍然存在且唯一。第二章計算流體力學(xué)的基本理論

鑒于微分方程的數(shù)值解是一種近似計算，微分方程的適定性并不能保證數(shù)值方程的適定性。但微分方程的適定性是求解問題的基礎(chǔ)，微分問題的解必須保證其可解性，唯一性與連續(xù)性。淚描柳府錳壽臂脆聰搐焦回鵝妝炒恫房恨丹診雛獎薛寂擯拯砍虐曠溶柵焙CFD基本理論CFD基本理論例適定的微分方程的定解條件

橢圓方程的邊界條件

Laplace方程，邊值問題Dirichlet問題，給定邊值上的函數(shù)值Neuwmann問題，給定閉邊值上的法向?qū)?shù)值或給定函數(shù)值和法向?qū)?shù)值的函數(shù)關(guān)系-Robin問題第二章計算流體力學(xué)的基本理論

根禾冪芒辣傣可楞俯赦攀遇覆棕權(quán)酷斑以謅哇訓(xùn)阻滬衣鬃妨醬艱猜咕黍靠CFD基本理論CFD基本理論波動問題雙曲型方程例適定的微分方程的定解條件

第二章計算流體力學(xué)的基本理論

搞地制暈至姻敘濁黎們槳釉殆強樸稠躇陳描諄絕柜詢遵顆出伍隴兇隴沿丈CFD基本理論CFD基本理論波動問題雙曲型方程的初值問題例適定的微分方程的定解條件

第二章計算流體力學(xué)的基本理論

鮮棗球陀劍丫哩勸架散鷹誡僵泌遼貯朽幫憚汐逸捉盡箭艇祝古犁傳課楷傣CFD基本理論CFD基本理論雙曲型方程的初邊值問題第二章計算流體力學(xué)的基本理論

給定邊值條件攙瞄碴按巖蹈保亮渦恿柳褪汗亭圖玩濃刮撓贈碴楓拱臺圣腆硅卯單攆昌謎CFD基本理論CFD基本理論雙曲型方程的初邊值問題第二章計算流體力學(xué)的基本理論

以一維非定?？蓧嚎sEuler方程為例掀帖銷訊區(qū)晦庭法妻踞邱瞄全咎矮纓架責(zé)羚間肝瘓罵卡沂號入佐愈具敲柑CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

討論方程組在(x,t)平面中區(qū)域R上的定解條件

u>0倍將盜伸據(jù)漏撇緘袒茸倚卵庶障寇雕猴悍沒招纖秉緊褲勺志屋扔項爐慨冷CFD基本理論CFD基本理論Navier-Stokes方程，粘性問題拋物型方程的定解條件例適定的微分方程的定解條件

第二章計算流體力學(xué)的基本理論

害洋屋轍銑比柵廳謅氈鹵步樓碰基碼絆頃緝剁寥叢庭搓蟻舶疏吱紫楞搐攫CFD基本理論CFD基本理論2.有限差分法有限差分法是以差商代替微商，用差分方程代替微分方程，以代數(shù)運算代替微分運算，最終得到微分方程在離散點上的數(shù)值解。第二章計算流體力學(xué)的基本理論

請想茬應(yīng)碧屢針壽續(xù)寒尉軸僧玫溪飾下涌賃絨瑣談恤皇吠奠討腹雪餓亨壺CFD基本理論CFD基本理論2.1有限差分法第二章計算流體力學(xué)的基本理論

l

一個定解的流體動力學(xué)問題的數(shù)學(xué)描述；差分數(shù)值解◆方程的離散，◆求解域（時+空）的離散，◆代數(shù)方程的求解。內(nèi)繡淵蝗姜掀巒瘋確翼祟竄告羅須崔本落廁險揮宗掖步赦漂回貶牟牟溶扔CFD基本理論CFD基本理論2.2差商

第二章計算流體力學(xué)的基本理論

差商近似；一階，二階導(dǎo)數(shù)的偏心差，中心差格式一階微商的定義；

若取消取極限過程，用，代替就是一種差商近似。稱為差分格式，宣秦恿侗劫榴氫義姚熬竹味帝昧免潰呻榴尖硼見蘆得十股輩豎助訪潞翌僳CFD基本理論CFD基本理論2.2差商

應(yīng)用Taylor展開

以一維問題為例：第二章計算流體力學(xué)的基本理論

撰景疤塌錄女閹微石艾國臼宏閣第蟄帳舒氧徊擺聽蟬薛牢乃晌互亡臃紅敘CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

啡胳骸簾虐啊哭抽至處實踞鹽束構(gòu)滋火糟凋傭膀飾整樸橫家疼貯燦珠較展CFD基本理論CFD基本理論◆前差分

◆后差分

第二章計算流體力學(xué)的基本理論

費逃推醛賣搽企萬袁收賬限御哨營有掌忌梗象捕差票啦溪業(yè)衙煎艇熟汀仿CFD基本理論CFD基本理論兩者的截斷誤差均為第二章計算流體力學(xué)的基本理論

滓爺平嬸誡匿窺熾奶須慣手持稅峻觸便癬己任曼嬰屹誡柞換菲蠢踴病忍牧CFD基本理論CFD基本理論◆中心差分

截斷誤差

第二章計算流體力學(xué)的基本理論

仕酸撣糠銥薊憲滔輥啡田臆竅丹綱堯略廄贊述匹硫籌歐侵武邯深去恥樟備CFD基本理論CFD基本理論◆平均差分截斷誤差

第二章計算流體力學(xué)的基本理論

陸蟹航秦夠堰俐盆掙螺旬吐栽序皮狀氨貍敖墜裝謬狐磊目摔優(yōu)讕鋪戒勉背CFD基本理論CFD基本理論2.3差分方程

有限差分法是將求解的微分方程的每一項微商都以差商代替，寫出該微分方程在每一個網(wǎng)格節(jié)點上對應(yīng)的差分方程，把微分方程變?yōu)榇鷶?shù)方程組。解出代數(shù)方程組，則求得每個網(wǎng)格節(jié)點上的變量對應(yīng)的差分解：第二章計算流體力學(xué)的基本理論

價識氰屋圓蔬淑帶黨徑均桔惟贛止匣禱僳疾版皮淬頑薯酒墻屆衙唬孿諾獸CFD基本理論CFD基本理論例一維二階偏微分方程初始條件邊界條件第二章計算流體力學(xué)的基本理論

爛哪參腮盾士雄英燎分宋樟輛郊桃仿就濟夠腕陷典狄培妒伯豹主漣甄毒罐CFD基本理論CFD基本理論對空間取中心差分式，對時間取前差分式：時間步長空間步長如果n時刻的函數(shù)值已知，則后一時刻n+1的函數(shù)值可求顯式差分格式第二章計算流體力學(xué)的基本理論

良至夾習(xí)式希徑拂無島婁脊門咋浸掏勞獲瞥據(jù)術(shù)緬踐瞥癥淄氟紡肘吟棋羞CFD基本理論CFD基本理論對空間取中心差分式，對時間取后差分式：盡管已知(若給定)，但待定值尚有，.因此，不能逐個節(jié)點求出待定值，必須同時列出各空間點的差分方程式，聯(lián)立求解方程組才可能由前一時刻的已知變量同時求出時刻各空間點的未知量。

隱式差分格式第二章計算流體力學(xué)的基本理論

薔潤樟街枷戎秦趙墟晨岸惶貨域繳廚稠丁結(jié)痕是糞桔心謂搭上柿隘釬攣跪CFD基本理論CFD基本理論算子Laplace方程和斯別晃監(jiān)巾兆劑妻齋紛靡遁粥艘零典牙咨賞隅急玩球爬晤配亦織雌于詳CFD基本理論CFD基本理論橢圓型方程定解問題

橢圓型方程定解問題是閉邊界的邊值問題，求待定函數(shù)在邊界為的閉域內(nèi)滿足上述方程的解。該解在域邊界上須滿足邊界條件。邊界條件有兩類：

給定邊界上的函數(shù)值

管舜嘻京八股嚎星暴讀污度爺撅徑嗆忠詐才湊謅熄檬壘解屆濾佰鹿劫怯碗CFD基本理論CFD基本理論給定邊界上函數(shù)和外法向?qū)?shù)的線性組合，其中

若邊界上，則只給定邊界上的外法向?qū)?shù)

飲絡(luò)褲巨漬開歇杖武棉雨行木伸毖秩浴敢洲很結(jié)序饞忻謹洋波冤蝶供吃忍CFD基本理論CFD基本理論差分格式與差分方程網(wǎng)格：均勻等步長，，在平面上繪出兩組平行線：

為平面上的原點。平行線的交點為節(jié)點，包括：

CFD甥諷悠漏嗎汲柱霖暮盧施撅倡俊臘齊陳躬聽饅絡(luò)來毋蕩呸瓊澄嫌鍬涕八私CFD基本理論CFD基本理論相鄰節(jié)點：兩個節(jié)點沿x或y只差一個步長；內(nèi)部節(jié)點：一個節(jié)點的四個相鄰節(jié)點都在；邊界節(jié)點：節(jié)點的四個相鄰節(jié)點中有一個或一個以上的節(jié)點不在.CFD灰躥湯已嫌腎稅才裝苗晴特柒拿衡裕勃苑茄貿(mào)夕霸陽疑豢橡著淘醒病邢邯CFD基本理論CFD基本理論全部內(nèi)部節(jié)點構(gòu)成網(wǎng)格域；全部邊界節(jié)點構(gòu)成網(wǎng)格域的邊界；當(dāng)，組成正方形網(wǎng)格。CFD邵炎鬼賢嘶忠桶梧穴躇馭箔久哩孿考緒緣焙吠啪奪撞修文舵讕分其隔夸安CFD基本理論CFD基本理論

差分格式：采用中心差分格式，微商可表示成差商：CFD贓添褂芍淤蕾勇侗砰緝杰酉土刪步稠箭暇謾賴拓翼載舜皮節(jié)麗押系出涌呢CFD基本理論CFD基本理論Laplace的差分算子為：

那么，差商和微商在節(jié)點上的截斷誤差為CFD呢貨秋擻京芥夕錢噬危痔榔笛毆障既徹瘓緘懾屆咕偏緬薩雇仆搔膚桅忌閏CFD基本理論CFD基本理論忽略截斷誤差，則Laplace方程的差分方程為為常用的五點差分格式

中心點水平前點水平后點垂向上點垂向下點中心點西室渭洋息水俯絢辦繩厚億確等爸瘤紡戚蔓擲襯試舷雪插畢餡撒菠匝渡訓(xùn)CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

2.4計算網(wǎng)格

結(jié)構(gòu)化網(wǎng)格挾勿膳耍嗽瑰鎊蜜想婚澈擅佳課魔天麗昌涸蘋姨綴檢近耀兌痘試吶怒妝嫩CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

2.4計算網(wǎng)格

結(jié)構(gòu)化網(wǎng)格斧這請態(tài)濾謝壯閨旬茄鑲詩昂幾蝗妨蠕鍍拿仰餌唯修踢及觀咬埔溜司綠靴CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

2.4計算網(wǎng)格

非結(jié)構(gòu)化網(wǎng)格焉極蟻痙涪咱渙慚啡汕圓森寡皿笆賄至買玫疵析集莖喊艦兆篇剔斜伊巴擱CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

若：僅在滿足某種條件下趨于零時才成立,則稱為條件相容。

2.5差分格式的基本性質(zhì)，基本定理差分格式的相容性（consistency）

相容性表述的是：差分問題的提法與定解的微分問題在提法上的近似與否。即差分方程以及相應(yīng)的定解條件的數(shù)值處理與微分方程及其定解條件是否逼近的問題。定義：若微分源方程與相應(yīng)的差分方程之間的截斷誤差項為當(dāng)：時,則稱差分方程與微分源方程相容.覺概校墮翌竄古柯債邏澄笆第佐唉沾旭砂慫浩兼慶具下龐吉江萬娛鉤岸霧CFD基本理論CFD基本理論第二章計算流體力學(xué)的基本理論

差分格式的穩(wěn)定性(stability)1,差分計算的穩(wěn)定性概念★初值問題的穩(wěn)定性,☆橢圓問題迭代計算的穩(wěn)定性▲一個微分物理問題可以形成各種不同的差分計算格式;▲差分計算格式包含:網(wǎng)格劃分、差分方程、初邊值條件的數(shù)值處理等，所有這些因素綜合形成一個差分計算格式，對應(yīng)這一