一维非稳态导热模型及数值解算1.问题的提出将一块面积较大的高温金属薄板,放入恒温的自然流动的空气室内,高温块体会从内向外传递热量,逐渐冷却,那么金属薄板的表面温度随时间到底是怎么变化的呢?2.问题的简化1)由于金属薄板的长(l)和宽(m)远远大于厚度(d),故将金属薄板可以简化为无限大薄壁导热问题,即热量的传递沿厚度方向,为单项导热;2)室内空间较大,可将室内气体视为恒温条件;3)平板的厚度为d,密度为ρ,初始温度为To,导热系数为λ,比热为c,空气对流换热系数为h,室内温度为Ta;4)平板内部导热满足傅里叶定律,表面换热满足牛顿冷却公式。3.数理方程描述由于平板换热关于中心线是对称的,仅对平板一半区域进行计算即可,坐标x的原点选在平板中心线上,因而一般区域的非稳态导热的数学描述为:(1)数理方程(2)边界条件4.数值离散建模1)求解域的离散空间步长为Dx,时间步长为Dt,Dx、Dt大小的选择需要保证节点温度方程求解的稳定性,如图1。2)节点温度差分方程的建立运用热平衡法可以建立非稳态导热物体内部节点和边界节点温度差分方程。表示空间节点i在kDt时刻(简称k时刻)的节点温度。a、内部节点温度差分方程:b、外边界节点温度差分方程:c、中心线边界节点温度差分方程:5.程序编写publicclass{staticinti,k=0;staticintN=20;staticdoubleh=5;//空气自燃对流换热系数staticdoubleλ=109;//黄铜的导热系数staticdoubleΔx=0.01;staticdoublec=3.9*100;//铜的比热staticdoubleρ=8900;//铜的密度staticdoubleτ=1;staticdoublea=λ/ρ/c;staticdoubleFo=a*τ/Δx/Δx;staticdoubleBi=h*Δx/λ;staticdoubleTo=100;staticdoubleT1=20;staticdouble[]T=newdouble[N+1];publicstaticvoidmain(String[]args){for(i=0;i<=N;i++){T[i]=To;System.out.printf("%7.1f",T[i]);}System.out.println();do{k++;for(i=1;i=0.1);System.out.println(k);}}6.结果图.黄铜金属薄板表面温度随时间的变化曲线