传热传质上机实习题(参考资料c语言)

《传热传质上机实习题(参考资料c语言)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、练习题三：一维非稳态导热的数值计算非稳态导热问题由于有时间变量，其数值计算出现了一些新的特点。在非稳态导热微分方程中，与时间因素相关的非稳态项是温度对时间的一阶导数，这给差分离散带来了新的特点。由于这个特点，可以采用不同的方法构造差分方程，从而得到几种不同的差分格式，即所谓的显式、隐式和半显式。我们仍从一个具体问题出发来研究非稳态导热问题的数值计算。3.1问题一块无限大平板（如图3所示），其一半厚度为L=0.1m，初始温度T0=1000，突然将其插入温度T=20的流体介质中。平板的导热系数=34.89W/（m），密度=7800kg/m3,比热c=712J/（kg），平板与介质的对流换热系数为

2、h=233W/（m2），求平板内各点的温度分布。程序#include#include#defineN10#defineK11main(){inti,j,l;floatcha;floata,x,y,Fo,Bi;floatt[N][K],b[N][K];/*打印出题目*/printf("ttt一维非稳态导热问题tt");printf("tttttt----");printf("题目：练习题三");y=1;/*y代表Δτ*/x=0.05/(N-1);a=34.89/(7800*712);Fo=(a*y)/(x*x);Bi=23

3、3*x/34.89;printf("显示格式条件:");printf("1、Fo=%3.1f<0.5t",Fo);printf("t2、1-2Fo*Bi-2Fo=%4.2f>0",1-2*Fo*Bi-2*Fo);/*时刻为零时，赋予初场温度*/for(i=0;i0.001){for(i=

4、0;i

5、i][j]);cha=cha/N;}}/*输出温度分布，其中l控制输出值的排列;这个结果是横轴为x，纵轴为τ的直角坐标下从左上角开始依次的*/printf("经数值离散计算的温度分布为：");l=0;for(j=K-1;j>=0;j--)for(i=0;i999.99)printf("%6.1f",t[i][j]);elseprintf("%6.2f",t[i][j]);l=l+1;if(l==N){printf("");l=0;}}getchar();/*为了是生成的exe文件结果算的后不会立即退出，方便观看*/}