传热学数值计算大作业

传热学数值计算大作业一 选题 《传热学》第四版 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 时，Fo<1/(2*(1+Bi)),理论上出现正确的计算结果；当 Fo=1 时，Fo>1/(2*(1+Bi))，Fo>0.5,理论上温度分布出现振荡，与实际情况不符。三 网格划分 将无限大平面的一半划分为 6 个控制体，共 7 个节点。△x=0.03/N=0.03/6=0.005，即空间步长为 0.005m四 节点离散方程绝热边界节点即 i=1 时，tij+1=2Fo△ti+1j+(1-2Fo△)tij内部节点即 00.5；对于对流边界节点，Fo<1/(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;/*共计算 0~T 时刻的温度变化*/float t[7][7],c[7][7];/*将节点和时间对应构成二维数组*/ /*参数设定*/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;i<=N;i++)t[i][0]=100;/*对 T>0 时刻各节点的温度进行计算*/for(j=0;j