三維穩(wěn)態(tài)導(dǎo)熱問(wèn)題數(shù)值求解 實(shí)驗(yàn)內(nèi)容的ppt講解

1、實(shí)驗(yàn)一、三維穩(wěn)態(tài)導(dǎo)熱問(wèn)題數(shù)值求解實(shí)驗(yàn)一、三維穩(wěn)態(tài)導(dǎo)熱問(wèn)題數(shù)值求解 一正方體金屬塊，其長(zhǎng)寬高均為0.1m,上邊界面溫度為200，其他5個(gè)面溫度均為50，利用C語(yǔ)言在101010三維網(wǎng)格上編寫(xiě)該三維穩(wěn)態(tài)導(dǎo)熱問(wèn)題計(jì)算程序，并求出圖2中所示中間面的溫度分布。 實(shí)驗(yàn)一、三維穩(wěn)態(tài)導(dǎo)熱問(wèn)題數(shù)值求解實(shí)驗(yàn)一、三維穩(wěn)態(tài)導(dǎo)熱問(wèn)題數(shù)值求解n一正方體金屬塊，其長(zhǎng)寬高均為0.1m,上邊界面溫度為200，其他5個(gè)面溫度均為50，利用C語(yǔ)言在101010三維網(wǎng)格上編寫(xiě)該三維穩(wěn)態(tài)導(dǎo)熱問(wèn)題計(jì)算程序，并求出圖中所示中間面的溫度分布。 其余邊界50上邊界200導(dǎo)熱系數(shù)為常數(shù)、穩(wěn)態(tài)、無(wú)內(nèi)熱源時(shí)的導(dǎo)導(dǎo)熱系數(shù)為常數(shù)、穩(wěn)態(tài)、無(wú)內(nèi)熱源時(shí)的導(dǎo)熱

2、微分方程式（控制方程）為熱微分方程式（控制方程）為 2222220tttxyz第一類邊界條件第一類邊界條件 ttb=200； tbob=50； tlb=50； trb=50； tfb=50； tbab=50；劃分為劃分為101010101010的三維網(wǎng)格后，的三維網(wǎng)格后，x=y=z 邊界點(diǎn)方程為邊界點(diǎn)方程為1,1, 1, 1, 1, 1,61kjikjikjikjikjikjikjittttttt內(nèi)節(jié)點(diǎn)離散方程內(nèi)節(jié)點(diǎn)離散方程 第一類邊界條件第一類邊界條件 ti,j,n=200； ti,j,0 =50；ti,0,k =50 ； ti,m,k=50； tl,j,k=50；t0,j,k=50； 實(shí)驗(yàn)

3、二、二維非穩(wěn)態(tài)導(dǎo)熱問(wèn)題數(shù)值求解實(shí)驗(yàn)二、二維非穩(wěn)態(tài)導(dǎo)熱問(wèn)題數(shù)值求解 n 已知二維平板長(zhǎng)寬分別為0.2m、0.1m，平板的熱擴(kuò)散系數(shù)a=0.510-4，平板初始溫度為50，AB邊界溫度為200，BC邊界為100，AD及CD邊均維持50恒溫，編制該二維非穩(wěn)態(tài)導(dǎo)熱問(wèn)題的計(jì)算程序，并求解時(shí)刻二維平板內(nèi)溫度分布。 導(dǎo)熱系數(shù)為常數(shù)、無(wú)內(nèi)熱源，導(dǎo)熱微分方程式導(dǎo)熱系數(shù)為常數(shù)、無(wú)內(nèi)熱源，導(dǎo)熱微分方程式（控制方程）（控制方程）222222()ttttaxyz式中，式中， ，稱為熱擴(kuò)散率。，稱為熱擴(kuò)散率。)/( ca對(duì)于二維問(wèn)題，導(dǎo)熱微分方程式為對(duì)于二維問(wèn)題，導(dǎo)熱微分方程式為 2222ytxtat第一類邊界條件第一類邊界條件 tAB=200； tBC=100； tCD=50； tDA=50；初始條件（時(shí)間條件）初始條件（時(shí)間條件） tC0500tt劃分為劃分為10105 5的二維網(wǎng)格后，的二維網(wǎng)格后，x=y邊界點(diǎn)方程為邊界點(diǎn)方程為內(nèi)節(jié)點(diǎn)離散方程（顯式差分格式）內(nèi)節(jié)點(diǎn)離散方程（顯式差分格式） 第一類邊界條件第一類邊界條件 tAB=ti,M=200； tBC=tL,j=100 ； tCD=ti,0=50； tDA=t0,j=50