1、第四章第四章 導(dǎo)熱問題數(shù)值解法基礎(chǔ)導(dǎo)熱問題數(shù)值解法基礎(chǔ)第一節(jié) 建立離散方程的方法第二節(jié) 穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算第三節(jié) 非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算第四章第四章 導(dǎo)熱問題數(shù)值解法基礎(chǔ)導(dǎo)熱問題數(shù)值解法基礎(chǔ) 數(shù)值解法為一種近似解法，但其近似精度可人為控制，故又是一種十分精確的解法。數(shù)值解法的主要思路和步驟：數(shù)值解法的主要思路和步驟：1.1.將物體從空間上或時(shí)間上劃分成若干個(gè)網(wǎng)格單元（即離散化），以單元上節(jié)點(diǎn)的值替代單元值；2.2.用各種方法建立起各節(jié)點(diǎn)與相鄰節(jié)點(diǎn)的關(guān)系方程式，構(gòu)成節(jié)點(diǎn)方程組；3.3.解節(jié)點(diǎn)方程組，得各節(jié)點(diǎn)數(shù)值（如：ti,j，傳熱中多為溫度值，此時(shí)即認(rèn)為得溫度分布）。第一節(jié)第一節(jié) 建立離散

2、方程的方法建立離散方程的方法 如圖，i、j分別為沿x、y方向節(jié)點(diǎn)的序號(hào)，x、y為步長(zhǎng)，網(wǎng)絡(luò)線與邊界的交點(diǎn)為邊界節(jié)點(diǎn)。 圖中每一個(gè)節(jié)點(diǎn)都代表以它為中心的一個(gè)小區(qū)域（的溫度值等）。 節(jié)點(diǎn)越多，結(jié)果越精確。 同樣可把時(shí)間也分割成許多間隔，同樣越小，結(jié)果越精確。一、區(qū)域和時(shí)間的離散化 空間上，把物體劃分為有限數(shù)目的網(wǎng)格單元，將原來在空間上連續(xù)的物理量，轉(zhuǎn)變?yōu)橛邢迋€(gè)離散的網(wǎng)格單元節(jié)點(diǎn)（又稱結(jié)點(diǎn)）。ix（ i +1）x（ j +1）y jy（ j -1）y（ i -1）x第一節(jié)第一節(jié) 建立離散方程的方法建立離散方程的方法二、泰勒級(jí)數(shù)展開法（有限差分法）二、泰勒級(jí)數(shù)展開法（有限差分法） 對(duì)于連續(xù)函數(shù)t，相鄰

3、兩節(jié)點(diǎn)間關(guān)系可用泰勒級(jí)數(shù)描述。 1! 3,!2, 1333222xjixtxjixtjixtjijixtt 2! 3,!2, 1333222xjixtxjixtjixtjijixtt xxttjixtjiji0, 1, xxttjixtjiji0, 1,1.1.一階導(dǎo)級(jí)的一階導(dǎo)級(jí)的向前差分向前差分表達(dá)式：表達(dá)式：舍去式x2后各項(xiàng)，則有：2.2.一階導(dǎo)級(jí)的一階導(dǎo)級(jí)的向后差分向后差分表達(dá)式：表達(dá)式：舍去式x2后各項(xiàng)，則有：第一節(jié)第一節(jié) 建立離散方程的方法建立離散方程的方法二、泰勒級(jí)數(shù)展開法（有限差分法）二、泰勒級(jí)數(shù)展開法（有限差分法） 對(duì)于連續(xù)函數(shù)t，相鄰兩節(jié)點(diǎn)間關(guān)系可用泰勒級(jí)數(shù)描述。 1! 3,

4、!2, 1333222xjixtxjixtjixtjijixtt 2! 3,!2, 1333222xjixtxjixtjixtjijixtt 22,0, 1, 1xxttjixtjiji 3022,2, 1, 122xxtttjixtjijiji 4022,2,1,1,22yytttjiytjijiji3.3.一階導(dǎo)級(jí)的一階導(dǎo)級(jí)的中心差分中心差分表達(dá)式：表達(dá)式：-式且忽略后項(xiàng)，則有：4.4.二階導(dǎo)級(jí)的二階導(dǎo)級(jí)的中心差分中心差分表達(dá)式：表達(dá)式：+式且忽略后項(xiàng)，則有：第一節(jié)第一節(jié) 建立離散方程的方法建立離散方程的方法三、熱平衡法三、熱平衡法 基于能量守恒定律，對(duì)微元體內(nèi)，認(rèn)為相鄰節(jié)點(diǎn)間溫度分布呈線

5、性，如圖，節(jié)點(diǎn)P與周圍節(jié)點(diǎn)的導(dǎo)熱量由傅里葉定律有： 1, 1yxttLPjiji1, 1yxttRPjiji1,1,xyttTPjiji1,1,xyttBPjiji0221,1, 1, 1jijijiyxjijijixytttttt若無熱源、常物性、穩(wěn)態(tài)：LP+RP+TP+BP=0則有：第二節(jié)第二節(jié) 穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算一、內(nèi)節(jié)點(diǎn)離散（節(jié)點(diǎn)）方程的建立一、內(nèi)節(jié)點(diǎn)離散（節(jié)點(diǎn)）方程的建立 以常物性、無熱源的二維穩(wěn)態(tài)導(dǎo)熱為例，一般為使計(jì)算簡(jiǎn)便，常使x=y，于是由式、的有限差分法和熱平衡法均可得到各內(nèi)節(jié)點(diǎn)間的關(guān)系為： ti+1,j+ti-1,j+ti,j+1+ti,j-1-4t

6、i,j=0 或?qū)懽鳎?ti,j=（ti+1,j+ti-1,j+ti,j+1+ti,j-1）/4 值得注意的是：上節(jié)點(diǎn)方程中五個(gè)節(jié)點(diǎn)均應(yīng)為內(nèi)節(jié)點(diǎn)，否則，上方程就可能不再適用。 對(duì)于內(nèi)節(jié)點(diǎn)方程，當(dāng)導(dǎo)熱過程為常物性、無熱源或熱源均勻分布的穩(wěn)態(tài)導(dǎo)熱時(shí)可用有限差分法得到，否則，采用有限差分法會(huì)產(chǎn)生麻煩，而熱平衡法則可克服這些麻煩，這是熱平衡法的優(yōu)點(diǎn)。但此法需對(duì)每個(gè)節(jié)點(diǎn)均作熱平衡分析。 一般邊界節(jié)點(diǎn)與內(nèi)節(jié)點(diǎn)間的聯(lián)系方程多用熱平衡法建立一般邊界節(jié)點(diǎn)與內(nèi)節(jié)點(diǎn)間的聯(lián)系方程多用熱平衡法建立。第二節(jié)第二節(jié) 穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算二、邊界節(jié)點(diǎn)離散方程的建立二、邊界節(jié)點(diǎn)離散方程的建立 第一類邊界條

7、件直接給出了邊界溫度值，可直接以數(shù)值形式參加到邊界節(jié)點(diǎn)與相鄰內(nèi)節(jié)點(diǎn)的離散方程中，第二、三類邊界條件則應(yīng)視具體情況，建立具體的離散方程。 下面以對(duì)流邊界條件下的拐角節(jié)點(diǎn)為例，介紹建立方法。第二節(jié)第二節(jié) 穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算二、邊界節(jié)點(diǎn)離散方程的建立二、邊界節(jié)點(diǎn)離散方程的建立yx12, 1yxttLPjiji12,1,xyttBPjiji12,xjifEPtth12,yjifFPtth0222,1, 1fxhjixhjijitttt令：x=y 且有：LP+RP+TP+BP=0代入整理得：第二節(jié)第二節(jié) 穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算三、節(jié)點(diǎn)離散方程組的求解三、節(jié)

8、點(diǎn)離散方程組的求解 對(duì)n個(gè)未知溫度的節(jié)點(diǎn)，用前述方法即可得到n個(gè)方程，將節(jié)點(diǎn)溫度ti,j按順序號(hào)1、2、n編號(hào)，且整理后總可得到如下方程組： t1=a11t1+a12t2+a1ntn+c1 t2=a21t1+a22t2+a2ntn+c2 tn=an1t1+an2t2+anntn+cn 采用計(jì)算機(jī)求解上述方程組時(shí)，常用迭代法和高斯消元法，這里僅介紹經(jīng)常使用的迭代法。第二節(jié)第二節(jié) 穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算三、節(jié)點(diǎn)離散方程組的求解三、節(jié)點(diǎn)離散方程組的求解1.1.簡(jiǎn)單迭代法簡(jiǎn)單迭代法A.先任設(shè)一組節(jié)點(diǎn)溫度的初始值：t10 t20 tn0；B.將t10 t20 tn0代入方程組右側(cè)，

9、得一組新值：t11t21tn1；C.將各項(xiàng)分別與各項(xiàng)相比較，看其中Max|ti1-ti0|是否小于規(guī)定的精度值，若小于時(shí)則退出循環(huán)，打印結(jié)果，否則繼續(xù)迭代；D.將t11 t21 tn1代入，得t12 t22 tn2，比較Max|ti2-ti1|是否小于規(guī)定的精度值，若小于時(shí)則退出循環(huán)，打印結(jié)果，否則繼續(xù)迭代；E.將t1k t2k tnk代入方程組，得t1k+1 t2k+1 tnk+1，直至Max|tik+1-tik|，或Max|（tik+1-tik）/tik|（此時(shí)為相對(duì)誤差），退出循環(huán)，打印結(jié)果，t1k+1 t2k+1 tnk+1這組值即為方程組的解。第二節(jié)第二節(jié) 穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算穩(wěn)態(tài)

10、導(dǎo)熱問題的數(shù)值計(jì)算三、節(jié)點(diǎn)離散方程組的求解三、節(jié)點(diǎn)離散方程組的求解2.2.高斯高斯- -賽德爾迭代法賽德爾迭代法 在簡(jiǎn)單迭代法基礎(chǔ)上作改進(jìn)，每次迭代時(shí)，總是使用節(jié)點(diǎn)溫度的最新值，例如第k次迭代后的值代入t1的表達(dá)式后得t1第k+1次迭代值t1k+1 ，則t2k+1 t2k+1 tnk+1應(yīng)為： t2k+1=a21t1k+1+a22t2k+a2ntnk+c2 t3k+1=a31t1k+1+a32t2k+1+a33t3k+a3ntnk+1+c3 tnk+1=an1t1k+1+an2t2k+1+ann-1tn-1k+1+anntnk+cn 這種迭代方法能使結(jié)果很快收斂，且可減少計(jì)算時(shí)間和計(jì)算機(jī)存貯量

11、，很快達(dá)到精度要求。第三節(jié)第三節(jié) 非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算一、顯式差分格式一、顯式差分格式 對(duì)一維、常物性、無內(nèi)熱源的非穩(wěn)態(tài)導(dǎo)熱沿x方向分割成n段，時(shí)間從 =0開始，接間隔分離，對(duì)其導(dǎo)熱微分方程式 中的右側(cè)用中心差分，左側(cè)用向前差分，則有：22xttax21112xtttttkikikikikiakixakikixakitttt2221111kikikikitFottFot211112122021xaxa即：整理得：或?qū)懽鳎涸谶x擇及x時(shí)，應(yīng)遵循：或： Fo1/2同理得：二維非穩(wěn)態(tài)導(dǎo)熱數(shù)值解的穩(wěn)定性條件為：Fo1/4第三節(jié)第三節(jié) 非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)

12、值計(jì)算二、隱式差分格式二、隱式差分格式 若將溫度對(duì)時(shí)間的導(dǎo)數(shù)采用向后差分，則有：21112xtttttkikikikikia21111112xtttttkikikikikiakikikixakixatttt111112221kikikikitttFotFo1111121整理后得：或： 上式與顯式不同的是：不能據(jù)k時(shí)刻的溫度分布直接得到（k+1）時(shí)刻的溫度分布，僅得到此時(shí)刻各節(jié)點(diǎn)間的方程式，要得到溫度分布，還需解節(jié)點(diǎn)方程，故此種格式是無條件穩(wěn)定的，即 、x可據(jù)情況任意選取。令k=k+1，上式可等價(jià)地寫成：第三節(jié)第三節(jié) 非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算三、邊界節(jié)點(diǎn)方程式的建立三、邊

13、界節(jié)點(diǎn)方程式的建立 對(duì)二、三類邊界，常用熱平衡的方法建立節(jié)點(diǎn)方程，據(jù)采用時(shí)間節(jié)點(diǎn)的前后不同，依然有顯式與隱式之分。xx/2=（k+1）=kt1k+1t1kt2k+1t2kAxttkk1221Atthkkf1對(duì)流/1112kkxttAcu或?qū)懽鳎簁kxckkfxhkktttttt111211122顯式隱式1.1.顯式差分格式顯式差分格式如圖，對(duì)于邊界節(jié)點(diǎn)1、t1k：A.內(nèi)節(jié)點(diǎn)2導(dǎo)入1的熱量：B.表面對(duì)流換熱量：C.內(nèi)能增量u：據(jù)熱平衡A+B=C并整理得：第三節(jié)第三節(jié) 非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算三、邊界節(jié)點(diǎn)方程式的建立三、邊界節(jié)點(diǎn)方程式的建立1.1.顯式差分格式顯式差分格式

14、t1k+1=2Fo（t2k+Bitfk）+（1-2BiFo-2Fo）t1k 顯然要使上式獲得穩(wěn)定解的條件仍然是t1k項(xiàng)前系數(shù)： 1-2BiFo-2Fo0 即： Fo1/（2Bi+2）討論：討論：當(dāng)問題為第三類邊界時(shí)，內(nèi)外節(jié)點(diǎn)均需滿足上述穩(wěn)定性條件，因Bi0，只要滿足上式條件，則必有：Fo1/2；對(duì)于絕熱邊界，因h=0，Bi=0，仍可用Fo1/2條件；對(duì)于二維非穩(wěn)態(tài)導(dǎo)熱，其內(nèi)外節(jié)點(diǎn)方程的穩(wěn)定性條件為： 421BiFo第三節(jié)第三節(jié) 非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算三、邊界節(jié)點(diǎn)方程式的建立三、邊界節(jié)點(diǎn)方程式的建立2.2.隱式差分格式隱式差分格式 對(duì)前圖，若以點(diǎn)t1k+1進(jìn)行熱平衡分析

15、，經(jīng)整理后可得： （1+2BiFo+2Fo）t1k+1=2Fo（t2k+1+Bitfk+1）+t1k四、顯式與隱式差分格式的區(qū)別：四、顯式與隱式差分格式的區(qū)別：1.顯式差分格式有穩(wěn)定性條件，即： x、 的選取要受Fo1/（2Bi+2）（一維）或Fo1/（2Bi+4）（二維）條件的制約，而隱式格式則無此約束，x、可任意選?。?.顯式只要知道了k 時(shí)刻的溫度分布，則（k+1）時(shí)刻的溫度分布可直接簡(jiǎn)單地解得，而隱式則需解出（k+1） 時(shí)刻的節(jié)點(diǎn)方程組，才能得此時(shí)的溫度分布。 第三節(jié)第三節(jié) 非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算五、節(jié)點(diǎn)離散方程組的求解五、節(jié)點(diǎn)離散方程組的求解1.1.顯式差分格式的數(shù)值解法步驟顯式差分格式的數(shù)值解法步驟在滿足數(shù)值解穩(wěn)定性條件的基礎(chǔ)上，選取適當(dāng)?shù)?、x，并對(duì)各節(jié)點(diǎn)進(jìn)行編號(hào)；用顯式格式寫出各內(nèi)、外節(jié)點(diǎn)方程式；據(jù)初始條件逐個(gè)計(jì)算各節(jié)點(diǎn)在 時(shí)刻的節(jié)點(diǎn)溫度，再據(jù) 時(shí)刻計(jì)算2 （n-1） 時(shí)刻計(jì)算n ，即得溫度分布。第三節(jié)第三節(jié) 非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算非穩(wěn)態(tài)導(dǎo)熱問題的數(shù)值計(jì)算五、節(jié)點(diǎn)離散方程組的求解五、節(jié)點(diǎn)離散方程組的求解2.2.隱式差分格式的數(shù)值解法步驟隱式差分格式的數(shù)值解法步驟選取適當(dāng)