工程流體力學(xué)的計算方法cfd基礎(chǔ)

1、第13章 工程流體力學(xué)的計算方法（CFD基礎(chǔ)）,§13.1代數(shù)方程的牛頓迭代法,牛頓迭代法用于求解超越方程,的根，在曲線,上取一點,求,顯然,是方程,的一個比,更精確的解，,重復(fù)以上計算可以得到任意精確的解。,例：水從池中經(jīng)管道流出，已知管長,沿程阻力損失系數(shù),局部阻力損失系數(shù),,水徑,設(shè)計流量,試求管徑d,解：列水面和管道出口截面的伯努利方程：,代入數(shù)據(jù)化簡得：,令,則上式化為：,選,作為初值,經(jīng)3次迭代后得,誤差小于,

2、因此取,退出,§13.2差分法,解析函數(shù),可以在點,領(lǐng)域展開成泰勒級數(shù),設(shè)有,三個差分節(jié)點，,其坐標(biāo)為,設(shè)函數(shù)在這三個節(jié)點的值為：,設(shè)節(jié)點間距為,則有泰勒展開式,①,②,退出,則有,③,一階導(dǎo)數(shù)向后差分式,④,一階導(dǎo)數(shù)向前差分式,可見,具有,的一階精度,①②上述兩式相減則有：,⑤,①②上述兩式相加則有：,⑥,對于形如,,的微分方程也可以求出y的泰勒展開式,,退出,兩式相減得：,可見：,具有三階精度。,退出,在平面勢流中，流函數(shù)

3、和速度勢函數(shù)均滿足拉普拉斯方程：,現(xiàn)將計算區(qū)域分成若干網(wǎng)格，每個,網(wǎng)格的邊長都是,，節(jié)點,簡記為,其二階導(dǎo)數(shù)可以用式⑥近似表示，則拉普拉斯的差分式為：,令,則：,⑦,對每個網(wǎng)格節(jié)點都建立形如上式的差分方程，就得到各節(jié)點的流函數(shù)的代數(shù)方程組，給出邊界條件，用迭代法可求出其數(shù)值解。,例：水在兩平板間流動，上板壁的滲透速度v0=1m/s下壁不可滲透，入口和出口速度均勻分布，分別為u1=3m/s和u2=1m/s 和，設(shè)板長h

4、=3m 寬h=1.5m,將長和寬分成3等份：,退出,①下壁面是一條流線，取其流函數(shù)為零，即,②左邊入口處：,因而,同理,③右邊出口處：,④在上壁面，,退出,各節(jié)點的代數(shù)方程，,由⑦式得：,點（2，2）,點（2，3）,點（3，2）,點（3，3）,退出,利用高斯法解此線性方程組得：,于是各節(jié)點的流函數(shù)的數(shù)值為：,4.50 3.50 2.50 1.50 3.00 2.33 1.67 1.001.50 1

5、.17 0.83 0.500.00 0.00 0.00 0.00,※ 流函數(shù)的物理意義：平面流動中，流過兩條流線間任一曲線的體積流量（單位厚度）等于兩條流函數(shù)之差。,退出,§13.3特征線法,特征線法用于求解一維非定?？蓧嚎s流動問題的數(shù)值解，水擊壓力波在管道內(nèi)的傳播，高速列車進入隧道時所產(chǎn)生的壓力波的傳播，都屬于這種流動。,圓截面中可壓縮粘性的非定常流動的運動微分方程：,①,管壁粘性摩擦應(yīng)力,可以用沿程阻

6、力損失系數(shù),表示：,②,②代入①得：,③,退出,連續(xù)性方程：,④,壓力波的傳播速度：,音波（微壓波）,水擊波：,D：管道直徑,E：流體體積彈性系數(shù),E固：管壁材料的彈性模量,δ：管壁厚度,ρ：流體密度,水擊波的傳播速度C=1200～1400m/s,退出,這樣連續(xù)性方程可改寫成:,③，⑤兩式相加，減得：,如果：,則：,如果：,則：,退出,表示一條x–t的關(guān)系曲線，記作C+，稱為特征曲線C+。,同樣,稱為特征線C－ 。,這樣⑥⑦式可以寫成：

7、,沿特征線C+：,⑧,沿特征線C－：,⑨,退出,下面用特征線法研究水擊波和隧道壓力波問題,1.水擊波,水擊波傳播速度C=1200～1400m/s，遠大于水速度u，因此，特征線方程可以寫成：,沿C+：,⑽,沿C－：,⑾,設(shè)管長為L，管壁為D，截面,管內(nèi)流速為u，閥門過流截面A隨時間而變,，過流速為V，V>u,退出,在初始時刻（t=0），管道定常出流，根據(jù)伯努利方程：,此時，管流速度為定常出流速度，記作u0,距管道入口x的截面上壓強水

8、頭,的伯努利方程：,即：壓強水頭h沿x的分布為：,退出,當(dāng)t>0時，管道入口的壓強水頭恒為h0即,在閥門處，V和h的關(guān)系：,此時管內(nèi)流速u：,取L為特征長度，h0為特征水頭，u0取為特征速度，L0/C特征時間,將式⑧⑨⑩⑾⑿寫成無量綱的形式：又,沿C+：,⒀,沿C－：,⒁,式中：沿,（邊界條件：）,退出,初始條件：,邊界條件：,⒂,（初始條件：）,目的是要計算出管道每一截面處的流動參數(shù)從時刻t=0到t=T任一時間里的變化值。,將管

9、長L分為N等分，,即無量綱的空間步長,時間步長與空間步長相等,退出,計算x–t圖上的網(wǎng)絡(luò)如上圖所示。,如果,時刻管軸上每一個節(jié)點,的流動參數(shù)u，h已經(jīng)算出，,利用特征線方程就可以計算下一時刻,各個節(jié)點,的流動參數(shù),將時刻,的節(jié)點（i+1，j）記為P，過P作特征線C+和C－必,通過節(jié)點和（i，j-1）（i-1，j）和（i，j-1）分別記為W和E,退出,根據(jù)特征線方程⒀⒁有：,或,這樣就可以求出P點的h和u，即：,其中：,退出,在x=0處，

10、只有一條特征線C－，但由于壓強水頭恒為1，因此：,在x=1處，即j=N處，只有一條特征線C+但u 可以由式⒂給出：,由以上兩式聯(lián)立可解出,和,逐個時刻進行計算，就可以得到各時刻管道上個節(jié)點的流動參數(shù)，u和h的值。,退出,2.隧道壓力波,波速C是小擾動波的傳播，,氣流速度u和C相比不一定,是微量，不能忽略，另外C是變化的，小擾動波的傳播可視為等熵,過程。,由,得到：,退出,于是特征線方程可以寫成：,沿C+：,⒃,沿C－：,⒄,式中：,特征

11、線不經(jīng)過節(jié)點，因為波動,方程的差分要滿足：,因空間步長△x要大于時間步長△t,所以w點落在節(jié)點（i,j-1）的右側(cè),E點落在節(jié)點（i,j+1）的左側(cè),退出,特征方程⒃⒄的差分形式為：,的位置未知，要反復(fù)迭代才能確定其位置，當(dāng),確定后，點的參數(shù)用節(jié)點,的參數(shù)內(nèi)插得到，符號表,示兩點間的平均值。,當(dāng)各節(jié)點的u和c求出后，由等熵關(guān)系求出節(jié)點的ρ ，R，T,其初始條件和邊界條件的確定比較復(fù)雜。,退出,§13.4有限元的插值函數(shù),一．線

12、性插值,如果已知曲線上的幾個點坐標(biāo),求曲線,方程，可以近似用折線表示：,如果將坐標(biāo)原點放在x1處，即x1=0且令x2–x1=L,則有：,令,退出,∴,顯然：,稱為樣條函數(shù)或插值函數(shù),已知n個點曲線可以表示為：,一.加權(quán)余量法，設(shè)有流動問題的方程和邊界條件是：,退出,一般化為：,設(shè)u近似為：,定義余量,為：,要使,最小，則存在加權(quán)函數(shù),與,在區(qū)域Ω內(nèi)E正,交，即內(nèi)積為零：,如果選,為加權(quán)函數(shù)，則有,即：,解此方程可求出系數(shù)u1,…,um這