計算傳熱學(xué)課程設(shè)計報告--正六面體的非穩(wěn)態(tài)導(dǎo)熱模擬

上傳人：奔*** IP屬地：河北更新時間：2024-03-01 格式：docx 頁數(shù)：24 大?。?75.69KB 人氣指數(shù)：12 舉報 版權(quán)申訴

已閱讀1頁，還剩23頁未讀，繼續(xù)免費閱讀

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、　　計算傳熱學(xué)課程設(shè)計報告　　論文題目：正六面體的非穩(wěn)態(tài)導(dǎo)熱模擬　　2014年 7 月 5 日　　目錄　　一、設(shè)計題目3<b

2、>　　二、控制方程4　　三、數(shù)學(xué)模型5　　四、數(shù)值方法及程序流程：5　　4.1隱式格式下離散過程5　　4.2顯式格式下離散過程6　　4.3程序流程：8&

3、lt;/p>　　五、集總參數(shù)法的應(yīng)用8　　六、導(dǎo)熱方程的無量綱化10　　七、結(jié)果的分析與處理11　　7.1 顯式格式處理結(jié)果11　　7.2 隱式格式處理結(jié)果15　　參考文獻(xiàn)

4、18　　附錄A：組員分工19　　附錄B：計算環(huán)境及源程序19　　正六面體的非穩(wěn)態(tài)導(dǎo)熱模擬　　摘要 　　本報告對正六面體的的非穩(wěn)態(tài)導(dǎo)熱問題進(jìn)行了模擬，運(yùn)用了數(shù)值模

5、擬的方法在三位直角坐標(biāo)系中對模型的導(dǎo)熱方程進(jìn)行了離散, 根據(jù)離散方程進(jìn)行了編程，最終模擬出了溫度隨時間變化的關(guān)系，得到了相關(guān)的圖像。　　關(guān)鍵詞：三維；數(shù)值模擬；非穩(wěn)態(tài)導(dǎo)熱 　　Regular hexahedron unsteady heatconduction simulation　　Abstr

6、act　　This report for regular hexahedron of unsteady heat conduction problem is simulated, using the method of numerical simulation in the three rectangular coordinate system to the discrete

7、 model of heat conduction equation based on the discrete equations programming, finally simulate the relationship between temperature change over time, obtained the relevant images.　　Keywords：Three d

8、imension；Numerical simulation；Unsteady heat conduction　　設(shè)計題目　　正六面體的非穩(wěn)態(tài)導(dǎo)熱模擬。如圖所示為一個正六面體，該正六面體由不銹鋼材料制成，三個邊的邊長分別為a=20mm，b=10mm，c=30mm，初始時刻溫度為0℃，突然將該六面體置于溫度為40℃的熱空

9、氣中，外表面的對流換熱系數(shù)為20W/（m²·K），不銹鋼的物性參數(shù)為16W/（m·K），ρ=7200kg/m³，c=460J/（kg·K）。　　請完成以下題目：　　建立該問題的數(shù)學(xué)方程，并給出定解條件；　　采用集總參數(shù)法對該問

10、題進(jìn)行分析，獲得溫度隨時間變化的函數(shù)關(guān)系，并計算該物體的時間常數(shù)和達(dá)到穩(wěn)態(tài)過程所需要的時間；　　選擇合適的參數(shù)將該問題無量綱化，得到無量綱導(dǎo)熱方程及邊界條件；　　分別用隱式和顯式格式對該問題進(jìn)行數(shù)學(xué)模擬，得到三維溫度隨時間的變化關(guān)系；　　將數(shù)值結(jié)果與集總參數(shù)法結(jié)果進(jìn)行對比，看看二者差距如何？

11、　　二、控制方程　　由導(dǎo)熱微分方程　　可得控制方程為：　?。?）物理特性：直角坐標(biāo)系，三維，非穩(wěn)態(tài)，無內(nèi)熱源，導(dǎo)熱　?。?）材料特性：常

12、物性　　（3）邊界條件：（六個面兩兩對稱）　　=0， 　　=0, 　　=0,

13、;　　三、數(shù)學(xué)模型　　數(shù)值方法及程序流程：　　4.1隱式格式下離散過程　　用內(nèi)節(jié)點法將所求的區(qū)域在均分網(wǎng)格的情況下離散化　　非穩(wěn)態(tài)項積分：　?。?擴(kuò)散項積

14、分：　　整理可得　　=+　　+ 　　對隱式格式進(jìn)行整理得：　　其中：</

15、b>　　4.2顯式格式下離散過程　　=+　　+ 　　(1)非穩(wěn)態(tài)項積分：　　(2)擴(kuò)散項積分：</p&g

16、t;　　整理得:　　=+　　+ 　　對顯式格式進(jìn)行整理得：　　其中：<

17、;p>　　4.3程序流程：　　五、集總參數(shù)法的應(yīng)用　　（1）根據(jù)能量平衡，在時間內(nèi)物體內(nèi)能的增量等于物體與環(huán)境的對流換熱量，即：　　式中，為在時間段內(nèi)物體的內(nèi)能的增量，　　式中，V為物體體積，和c分別是物體的密度和比熱，而T是物

18、體在任意時刻時的溫度。　　物體在時間內(nèi)的對流換熱量為，　　其中，，，，A=　　整理得：　　引入過余溫度：<

19、b>　　初始條件為：　　故，溫度隨時間變化的關(guān)系式：　?。?）時間常數(shù)：，帶入數(shù)據(jù)，　　取穩(wěn)態(tài)溫度T=39，則1666s　　六、導(dǎo)熱方程的無量綱化　?。?）控制方程：</b&

20、gt;　?。?）引入下述無量綱量　　特征長度：　　無量綱時間t：　　無量綱溫度：　　，， &

21、lt;/b>　　則有　　T=，　?。?）將上述無量綱數(shù)帶入控制方程中，　　經(jīng)整理化簡得：　　（4）無量綱化的邊

22、界條件：　　由無量綱溫度：，故，　　，，　　，，　　，，　　七、結(jié)果的分析與處理<p&

23、gt;　　7.1 顯式格式處理結(jié)果　　取中間某一截面繪圖　　4s時溫度分布　　8s時溫度分布　　12s時溫度分布<

24、;/p>　　16s時溫度分布　　20s時溫度分布　　24s時溫度分布　　256秒時的溫度分布　　7.2 隱式格式處理結(jié)果</

25、p>　　選取i=5截面進(jìn)行繪圖　　4s時的溫度場分布　　8s時的溫度場分布　　12s時的溫度場分布　　16s時的溫度場分布

26、20s時的溫度場分布　　24s時的溫度場分布　　148s時的溫度場分布　　由上面不同時刻的截圖可以得出如下結(jié)論：　?。?）在該模型的導(dǎo)熱過程中熱量由四周向中心傳遞　?。?）沿同一坐標(biāo)方向上，同一截面上不同溫度的分布區(qū)域呈橢圓形&l

27、t;/p>　?。?）隨著時間的增加，截面上低溫區(qū)域的面積越來越小　?。?）中心溫度越來越高，最終各不同溫度區(qū)域的溫差會越來越小　?。?）整個區(qū)域的溫度也會越來越接近環(huán)境溫度。　　(6)顯示格式和隱式格式下的溫度變化情況大致相同。　　(7)隨著時

28、間的推移，溫度升高的速度趨于平緩。　　(8)整個矩形截面的對稱軸也是溫度的對稱軸。　　將數(shù)值結(jié)果與集總參數(shù)法的結(jié)果進(jìn)行對比，集總參數(shù)法的時間大于數(shù)值方法的時間，造成這種結(jié)果的原因是，由于集總參數(shù)法應(yīng)用的前提是〈0.1，而題目中取，，所以誤差較大。　　參考文獻(xiàn)

29、;　　[1] 黃善波劉中良編計算傳熱學(xué)基礎(chǔ).青島：中國石油大學(xué)（華東）出版社，2009. 　　[2]楊世銘陶文銓編著傳熱學(xué)（第四版）高等教育出版社，2006 　　[3]陶文銓編著.數(shù)值傳熱學(xué)（第二版）.西安：西安交通大學(xué)出版社，2

30、001. 　　[4] 袁恩熙主編工程流體力學(xué)（第一版）.石油工業(yè)出版社，1986　　[5] 魏東平朱連章于廣斌編 C程序設(shè)計語言.北京：電子工業(yè)出版社,2009　　附錄A：組員分工　　組長：楊偉茂編寫程序

31、方程離散　　組員：蔡佳鑫設(shè)計報告方程離散　　組員：李培強(qiáng) 制作ppt 建立模型　　附錄B：計算環(huán)境及源程序　　計算環(huán)境：Microsoft Visual C++ 6.0　　隱式格式主

32、程序：　　#include<stdio.h>　　#include<math.h>　　# define NX 15　　# define NY 15　　# define NZ 15</

33、p>　　double T[NX][NY][NZ],b[NX][NY][NZ],T0[NX][NY][NZ],TT[NX][NY][NZ];　　double aP,aE,aW,aS,aN,aF,aB;　　double dx,dy,dz,dt,XL,YL,ZL;　　int N,M,W,p,

34、c,lamd,Tw,Tf,h;　　double Tao;　　void main()　　{　　FILE *fp;　　int i,j,k;

35、　　XL=0.02;//三邊距離　　YL=0.03;　　ZL=0.01;　　N=12;//節(jié)點數(shù)　　M=12;&

36、lt;p>　　W=12;　　dt=0.5;//時間步長　　Tw=273;//物體開始溫度　　p=7200;//物體密度　　c=460;　　Tf=31

37、3;//環(huán)境溫度　　h=20;//對流傳熱系數(shù)　　lamd=16;//導(dǎo)熱系數(shù)　　Tao=0;//時間　　for (i=1;i<N-1;i++)　　for (j=1;j<M-1;j++)<

38、p>　　for(k=1;k<W-1;k++)　　{　　T0[i][j][k]=Tw;　　T[i][j][k]=Tw;　　}

39、　　//步長 　　dx=1.0*XL/(N-1);　　dy=1.0*YL/(M-1);　　dz=1.0*ZL/(W-1);　　//賦初始值　　for (j=0;j<M

40、;j++)　　for(k=0;k<W;k++)　　{ 　　T0[0][j][k]=Tf;　　T[0][j][k]=Tf;　　T0[N-1][j][k]=Tf;&

41、lt;p>　　T[N-1][j][k]=Tf;　　}　　for (i=0;i<N;i++)　　for(k=0;k<W;k++)　　{

42、　　T0[i][0][k]=Tf;　　T[i][0][k]=Tf;　　T0[i][M-1][k]=Tf;　　T[i][M-1][k]=Tf;　　}　　for(i=0;i<N;i++)&l

43、t;/p>　　for (j=0;j<M;j++)　　{ 　　T0[i][j][0]=Tf;　　T[i][j][0]=Tf;　　T0[i][j][W-1]=Tf;<p&g

44、t;　　T[i][j][W-1]=Tf;　　}　　fp=fopen("result.txt","w");　　do　　{ </b&

45、gt;　　aE=aW=lamd*dz*dy/dx;　　aN=aS=lamd*dx*dz/dy;　　aF=aB=lamd*dx*dy/dz;　　aP=p*c*dx*dy*dz/dt+2*aE+2*aN+2*aF;

46、　//結(jié)果　　loop:　　for(i=1;i<N-1;i++)　　for(j=1;j<M-1;j++)　　for(k=1;k<W-1;k++)<

47、b>　　{ 　　b[i][j][k]=p*c*dx*dy*dz*T0[i][j][k]/dt;　　T[i][j][k]=aE*T[i+1][j][k]+aW*T[i-1][j][k]+aN*T[i][j+1][k]+aS*T[i][j-1][k]+aF*T[i][j][k+1]+aB*T[i][j][k-1];</p&

48、gt;　　T[i][j][k]=T[i][j][k]+b[i][j][k];　　T[i][j][k]=T[i][j][k]/aP;　　}　　for(i=0;i<N;i++)　　for(j=0;j<W;

49、j++)　　for(k=0;k<M;k++)　　T0[i][j][k]=T[i][j][k];　　Tao=Tao+dt;　　if(fmod(Tao,4)!=0) goto loop;　　//輸出</

50、b>　　if(fmod(Tao,4)==0)　　{　　fprintf(fp,"%lf時刻的溫度場為\n",Tao);　　for(i=0;i<N;i++)　　fo

51、r(j=0;j<W;j++)　　for(k=0;k<M;k++)　　fprintf(fp,"T[%2d][%2d][%2d]=%lf %lf %lf %lf\n",i,j,k,i*dx,j*dy,k*dz,T[i][j][k]);　　for(i=5;i<6;i++)</

52、p>　　for(j=0;j<W;j++)　　for(k=0;k<M;k++)　　fprintf(fp,"%lf %lf %lf\n",j*dy,k*dz,T[i][j][k]); 　　}

53、;　　}while(Tao<30);　　}　　顯式格式主程序：　　#include<stdio.h>　　#include<math.h><

54、/p>　　# define NX 15　　# define NY 15　　# define NZ 15　　double T[NX][NY][NZ],b[NX][NY][NZ],T0[NX][NY][NZ],TT[NX][NY][NZ];　　do

55、uble aP,aE,aW,aS,aN,aF,aB;　　double dx,dy,dz,dt,XL,YL,ZL;　　int N,M,W,p,c,lamd,Tw,Tf,h;　　double Tao;　　void main()<

56、b>　　{　　FILE *fp;　　int i,j,k;　　XL=0.02;//三邊距離　　YL=0.03;

57、　ZL=0.01;　　N=12;//節(jié)點數(shù)　　M=12;　　W=12;　　dt=0.5;//時間步長　　Tw=273;//物體開始溫度

58、　　p=7200;//物體密度　　c=460;　　Tf=313;//環(huán)境溫度　　h=20;//對流傳熱系數(shù)　　lamd=16;//導(dǎo)熱系數(shù)　　Tao=

59、0;//時間　　for (i=1;i<N-1;i++)　　for (j=1;j<M-1;j++)　　for(k=1;k<W-1;k++)　　{　　T0[i][j][k]=Tw

60、;　　T[i][j][k]=Tw;　　}　　//步長 　　dx=1.0*XL/(N-1);　　dy=1.0*YL/(M-1);&

61、lt;p>　　dz=1.0*ZL/(W-1);　　//賦初始值　　for (j=0;j<M;j++)　　for(k=0;k<W;k++)　　{ <

62、p>　　T0[0][j][k]=Tf;　　T[0][j][k]=Tf;　　T0[N-1][j][k]=Tf;　　T[N-1][j][k]=Tf;　　}　　for (i=0;i<N;

63、i++)　　for(k=0;k<W;k++)　　{ 　　T0[i][0][k]=Tf;　　T[i][0][k]=Tf;　　T0[i][M-1][k]=Tf;&l

64、t;p>　　T[i][M-1][k]=Tf;　　}　　for(i=0;i<N;i++)　　for (j=0;j<M;j++)　　{

65、　T0[i][j][0]=Tf;　　T[i][j][0]=Tf;　　T0[i][j][W-1]=Tf;　　T[i][j][W-1]=Tf;　　}　　fp=fopen("result.

66、txt","w");　　do　　{ 　　aE=aW=lamd*dz*dy/dx;　　aN=aS=lamd*dx*dz/dy;<p&g

67、t;　　aF=aB=lamd*dx*dy/dz;　　aP=p*c*dx*dy*dz/dt+2*aE+2*aN+2*aF;　　//結(jié)果　　loop:　　for(i=1;i<N-1;i++)&

68、lt;/p>　　for(j=1;j<M-1;j++)　　for(k=1;k<W-1;k++)　　{ 　　b[i][j][k]=p*c*dx*dy*dz*T0[i][j][k]/dt;　　T[

69、i][j][k]=aE*T[i+1][j][k]+aW*T[i-1][j][k]+aN*T[i][j+1][k]+aS*T[i][j-1][k]+aF*T[i][j][k+1]+aB*T[i][j][k-1];　　T[i][j][k]=T[i][j][k]+b[i][j][k];　　T[i][j][k]=T[i][j][k]/aP;

70、　　}　　for(i=0;i<N;i++)　　for(j=0;j<W;j++)　　for(k=0;k<M;k++)　　T0[i][j][k]=T[i][j][k];<

71、p>　　Tao=Tao+dt;　　if(fmod(Tao,4)!=0) goto loop;　　//輸出　　if(fmod(Tao,4)==0)　　{<p&g

72、t;　　fprintf(fp,"%lf時刻的溫度場為\n",Tao);　　for(i=0;i<N;i++)　　for(j=0;j<W;j++)　　for(k=0;k<M;k++)　　fprintf(fp,"T[%2d][%2

73、d][%2d]=%lf %lf %lf %lf\n",i,j,k,i*dx,j*dy,k*dz,T[i][j][k]);　　for(i=5;i<6;i++)　　for(j=0;j<W;j++)　　for(k=0;k<M;k++)　　fprint

眾賞文庫> 全部分類> 畢業(yè)設(shè)計

溫馨提示

1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 眾賞文庫僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

計算傳熱學(xué)課程設(shè)計報告--正六面體的非穩(wěn)態(tài)導(dǎo)熱模擬

文檔簡介

溫馨提示

最新文檔

評論

計算傳熱學(xué)課程設(shè)計報告--正六面體的非穩(wěn)態(tài)導(dǎo)熱模擬

文檔簡介

溫馨提示

最新文檔

評論

免費下載