传热学大作业.docxVIP

传热学大作业报告院 系：机电工程学院班 级：机械1112班学 号：0806111734学生姓名：高 辉 洲任课老师： 郑 煜 一厚度为0.1m的无限大平壁，两侧均为对流换热边界条件，初始时两侧流体温度与壁内温度一致，tf1=tf2=t0=5℃,已知两侧对流换热系数分别为h1=11w/（m2?K），h2=23w/(m2?K),壁的导热系数?=0.43w/(m?K)，导温系数?=0.3437*10-6m2/s.如果一侧的环境温度tf1,突然升高为50℃，，并维持不变，计算在其他参数不变的情况下，平壁内温度的分布及两侧对流密度随时间的变化规律（用图形表示）。要求：将全部计算内容（包括网格的划分、节点离散方程组的建立、计算框图、程序及计算结果）用A4纸打印。解：假设：本题是平壁一维非稳态导热过程，常物性，无内热源，初始条件为定温，两侧边界条件均为第三类边界条件。计算：（1）数学描述：导热微分方程和定解条件构成问题的完整描述：（2）时间-空间区域的离散化：将平壁10等分，共11个节点，步长?x=0.01m，?t=60s(经过验算满足稳定性条件??如下图所示：一维非稳态导热问题时间-空间区域离散化（3）微分方程的离散，即节点离散方程组：1)内节点（n=2~10)2)边界节点（n=1和n=11)由上最终可得到如下离散方程组：（4）计算框图:（5）程序及计算结果：#includestdio.h#includemath.hvoid main(){int N=10;int T=1000;float t[11][1002],q1[1000],q2[1000];inti,j;floatx,r=0.43,tf1,tf2,h1=11.0,h2=23.0; x=0.01;tf1=50.0;tf2=5.0;for(i=0;i=10;i++)t[i][0]=5.0;for(j=0;jT;j++){ for(i=0;i=10;i++){ if(i==0)t[i][j+1]=0.4821*t[i][j]+0.4124*t[i+1][j]+0.1055*tf1;if(i0i10) t[i][j+1]=0.2062*(t[i+1][j]+t[i-1][j])+0.5876*t[i][j];if(i==10)t[i][j+1]=0.3670*t[i][j]+0.4124*t[i-1][j]+0.2206*tf2;}q1[j+1]=h1*(50-t[0][j+1]);q2[j+1]=h2*(t[10][j+1]-5);}printf(\n当步长为0.01，时间步长为60s，间隔20个时间步长时计算的温度和热流密度分布为：\n);for(j=0;jT;j++) {if(j%20==1){for(i=0;i=10;i++)printf(%7.2f,t[i][j]);printf(\n);printf(%7.2f,%7.2f,q1[j],q2[j]);printf(\n);} }（6）运行结果：9.75 5.00 5.00 5.00 5.00 5.00 5.00 5.00 5.00 5.00 5.00442.78 0.0023.10 16.89 12.21 9.01 7.04 5.94 5.39 5.15 5.05 5.01 5.01295.87 0.1427.15 21.67 17.05 13.36 10.54 8.51 7.12 6.22 5.68 5.37 5.21251.31 4.7929.61 24.64 20.26 16.53 13.47 11.03 9.17 7.80 6.82 6.16 5.71224.25 16.2531.35 26.76 22.62 18.97 15.84 13.22 11.10 9.41 8.10 7.09 6.32205.11 30.2532.68 28.39 24.45 20.91 17.78 15.08 12.79 10.88 9.30 8.00 6.91190.50 43.9733.74 29.69 25.93 22.49 19.39 16.64 14.24 12.16 10.37 8.82 7.45178.88 56.3034.60 30.75 27.14 23.79 20.72 17.95 15.47 13.26 11.29 9.52 7.91169.42 66.9735.31 31.63 28.14 24.87 21.83 19.04 16.50 1