传热学数值计算大作业 一 选题 《传热学》第四版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 #include void 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...