传热学数值计算大作业(共7页).docx

《传热学数值计算大作业(共7页).docx》由会员分享，可在线阅读，更多相关《传热学数值计算大作业(共7页).docx（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上传热学数值计算大作业一 选题 传热学第四版P179页 例题 4-3二 相关数据及计算方法1厚2=0.06m的无限大平板受对称冷却，故按一半厚度作为模型进行计算2. =0.03m，初始温度t0=100,流体温度t=0；=40W/(m.K),h=1000W/(m2.K),Bi=h*x/=0.25；3.设定Fo=0.25和Fo=1两种情况通过C语言编程（源程序文件见附件）进行数值分析计算；当Fo=0.25时，Fo1/(2*(1+Bi)，Fo0.5,理论上温度分布出现振荡，与实际情况不符。三 网格划分 将无限大平面的一半划分为6个控制体，共7个节点。x=0.03/N=0.03

2、/6=0.005，即空间步长为0.005m四 节点离散方程绝热边界节点即i=1时，tij+1=2Foti+1j+(1-2Fo)tij内部节点即0i0.5；对于对流边界节点，Fo1/(2*(1+Bi)。3 备注 在Fo=0.25时，为了反映较长时间后温度的分布，取T=600,并选取了其中部分时刻的温度输出进行画图。图像显示，随着时间的增长，各点温度趋向一致。 而当Fo=1时由于结果会出现振荡，只取T=6观察即可。附录1 C语言源程序当步长为0.005，Fo=1的程序#include#includevoid main()int N=6;/*网格划分后计算区域共N个控制体*/int T=6;/*共计

3、算0T时刻的温度变化*/float t77,c77;/*将节点和时间对应构成二维数组*/ /*参数设定*/int i,j,tf,r,h,min;float Bi,x,Fo;Fo=1;x=0.03/N;tf=0;/*边界对流换热温度为0摄氏度*/r=40;h=1000;Bi=h*x/r; /*0时刻温度均为100摄氏度,即对各点设定初值*/for(i=0;i0时刻各节点的温度进行计算*/for(j=0;jT;j+) for(i=0;i0.00001) for(i=0;i=N;i+) if(i=0) tij+1=2*Fo*ti+1j+(1-2*Fo)*tij;/*绝热边界离散方程*/else if

4、 (i!=0)&(i!=N) tij+1=Fo*(ti+1j+ti-1j)+(1-2*Fo)*tij;/*中间所有节点的离散方程*/else if (i=N) tij+1=tij*(1-2*Fo*Bi-2*Fo)+2*Fo*ti-1j+2*Fo*Bi*tf;/*对流换热边界节点的离散方程*/ min=0; for(i=0;i=0;j-) for(i=0;i=N;i+) printf(%7.2f ,tij); n=n+1; if(n=(N+1) printf(n); n=0; getchar();当步长为0.005，Fo=0.25的程序#include#includevoid main()int

5、 N=6;/*网格划分后计算区域共N个控制体*/int T=600;/*共计算0T时刻的温度变化*/float t7601,c7601;/*将节点和时间对应构成二维数组*/ /*参数设定*/int i,j,tf,r,h,min;float Bi,x,Fo;Fo=0.25;x=0.03/N;tf=0;/*边界对流换热温度为0摄氏度*/r=40;h=1000;Bi=h*x/r; /*0时刻温度均为100摄氏度,即对各点设定初值*/for(i=0;i0时刻各节点的温度进行计算*/for(j=0;jT;j+) for(i=0;i0.00001) for(i=0;i=N;i+) if(i=0) tij+

6、1=2*Fo*ti+1j+(1-2*Fo)*tij;/*绝热边界离散方程*/else if (i!=0)&(i!=N) tij+1=Fo*(ti+1j+ti-1j)+(1-2*Fo)*tij;/*中间所有节点的离散方程*/else if (i=N) tij+1=tij*(1-2*Fo*Bi-2*Fo)+2*Fo*ti-1j+2*Fo*Bi*tf;/*对流换热边界节点的离散方程*/ min=0; for(i=0;i=0;j=(j-5) for(i=0;i=31;j=(j-30) for(i=0;i=N;i+) printf(%7.2f ,tij); n=n+1; if(n=(N+1) printf(n); n=0; getchar();专心-专注-专业