冶金傳輸原理B-2熱傳導(dǎo)-5-2012秋

穩(wěn)態(tài)導(dǎo)熱問(wèn)題的數(shù)值模擬例：一圓形金屬棒，長(zhǎng)L=0.5m，橫截面積為A=0.01m2，其導(dǎo)熱系數(shù)為常數(shù)1000W/m.℃，無(wú)內(nèi)熱源，金屬棒兩端溫度已給定，分別為100℃、500℃，且不隨時(shí)間變化，金屬棒徑向的溫度變化忽略不計(jì)。求該金屬棒內(nèi)的溫度分布。解：解析解數(shù)值解

穩(wěn)態(tài)導(dǎo)熱問(wèn)題的數(shù)值模擬解：例：一圓形金屬棒，長(zhǎng)L=0.5m，橫截面積為A=0.01m2，其導(dǎo)熱系數(shù)為常數(shù)1000W/m.℃，無(wú)內(nèi)熱源，金屬棒兩端溫度已給定，分別為100℃、500℃，且不隨時(shí)間變化，金屬棒徑向的溫度變化忽略不計(jì)。求該金屬棒內(nèi)的溫度分布。

穩(wěn)態(tài)導(dǎo)熱問(wèn)題的數(shù)值模擬解：

穩(wěn)態(tài)導(dǎo)熱問(wèn)題的數(shù)值模擬解：

穩(wěn)態(tài)導(dǎo)熱問(wèn)題的數(shù)值模擬解：課堂作業(yè):15分鐘按此種網(wǎng)格劃分方法對(duì)該問(wèn)題進(jìn)行數(shù)值求解

穩(wěn)態(tài)導(dǎo)熱問(wèn)題的數(shù)值模擬穩(wěn)態(tài)導(dǎo)熱的有限差分方法生成網(wǎng)格：將求解區(qū)域劃分為有限個(gè)網(wǎng)格單元建立差分方程差分方程的建立－熱平衡法差分方程的建立－邊界節(jié)點(diǎn)例：一矩形板，長(zhǎng)400mm，寬200mm，其導(dǎo)熱系數(shù)為常數(shù)，無(wú)內(nèi)熱源，假定該板各邊界上的溫度已給定，且不隨時(shí)間變化，在厚度方向的溫度變化可忽略不計(jì)。求該矩形板內(nèi)的溫度分布。解：第一類(lèi)邊界條件下的二維穩(wěn)態(tài)導(dǎo)熱問(wèn)題例：一矩形板，長(zhǎng)400mm，寬200mm，其導(dǎo)熱系數(shù)為常數(shù)，無(wú)內(nèi)熱源，假定該板各邊界上的溫度已給定，且不隨時(shí)間變化，在厚度方向的溫度變化可忽略不計(jì)。求該矩形板內(nèi)的溫度分布。解：第一類(lèi)邊界條件下的二維穩(wěn)態(tài)導(dǎo)熱問(wèn)題差分方程組的求解雅可比迭代：特點(diǎn)：計(jì)算第k+1次的值時(shí)，全部使用第k次迭代的值，收斂速度比較慢解

差分方程的計(jì)算求解差分方程組的求解雅可比迭代：解高斯－賽德?tīng)柕禾攸c(diǎn)：總是使用節(jié)點(diǎn)溫度的最新值

差分方程的計(jì)算求解顯式差分方程的穩(wěn)定性問(wèn)題有一無(wú)限大平板，初始溫度為0℃，開(kāi)始時(shí)平板兩側(cè)溫度突然升高到100℃，以后保持不變。假定平板分成10份，選取F=1，內(nèi)部節(jié)點(diǎn)方程及空間步長(zhǎng)、時(shí)間步長(zhǎng)如下：不穩(wěn)態(tài)導(dǎo)熱問(wèn)題的數(shù)值模擬顯式差分方程的穩(wěn)定性問(wèn)題有一無(wú)限大平板，初始溫度為0℃，開(kāi)始時(shí)平板兩側(cè)溫度突然升高到100℃，以后保持不變。假定平板分成10份，選取F=1，內(nèi)部節(jié)點(diǎn)方程及空間步長(zhǎng)、時(shí)間步長(zhǎng)如下：節(jié)點(diǎn)溫度，℃01234501000000012100100000024100010000036100200-1001000048100-100400-2001000用顯式差分方程求解不穩(wěn)態(tài)導(dǎo)熱時(shí)，解出現(xiàn)的波動(dòng)現(xiàn)象稱(chēng)為解的不穩(wěn)定性問(wèn)題不穩(wěn)態(tài)導(dǎo)熱問(wèn)題的數(shù)值模擬顯式差分方程的穩(wěn)定性問(wèn)題顯式差分方程的穩(wěn)定性判據(jù)－－F值的大小應(yīng)滿(mǎn)足使差分方程中