一维非稳态导热模型及数值解算.docx

《一维非稳态导热模型及数值解算.docx》由会员分享，可在线阅读，更多相关《一维非稳态导热模型及数值解算.docx（4页珍藏版）》请在第一文库网上搜索。

1、一维非稳态导热模型及数值解算1 .问题的提出将一块面积较大的高温金属薄板，放入恒温的自然流动的空气室内，高温块体会从内向外传递热量，逐渐冷却，那么金属薄板的表面温度随时间到底是怎么变化的呢？2 .问题的简化D由于金属薄板的长（/）和宽（相）远远大于厚度（”），故将金属薄板可以简化为无限大薄壁导热问题，即热量的传递沿厚度方向，为单项导热；2）室内空间较大，可将室内气体视为恒温条件；3）平板的厚度为d,密度为P,初始温度为To,导热系数为,比热为c,空气对流换热系数为h,室内温度为Ta；4）平板内部导热满足傅里叶定律，表面换热满足牛顿冷却公式。3 .数理方程描述由于平板换热关于中心线是对称的，仅对

2、平板一半区域进行计算即可，坐标X的原点选在平板中心线上，因而一般区域的非稳态导热的数学描述为：（1）数理方程（2）边界条件4 .数值离散建模1）求解域的离散空间步长为Ax,时间步长为,Ax、A大小的选择需要保证节点温度方程求解的稳定性，如图1。2）节点温度差分方程的建立运用热平衡法可以建立非稳态导热物体内部节点和边界节点温度差分方程。表示空间节点i在kA时刻（简称k时刻）的节点温度。a、内部节点温度差分方程：Ax=AAxpc1=Fo(1+1)+(1-2Fo)b、外边界节点温度差分方程:AJf一诸.,ZJ1IAAAXtO*1-tOa-77加侬T)=A万PCF-t1=2Fo(+/4)+(1-2BF

3、o-2Fo)C、中心线边界节点温度差分方程：5 .程序编写pub1icc1assstaticinti,k=0;staticint=20;staticdoub1e6=5;空气白燃对流换热系数staticdoub1eA=I09;黄铜的导热系数staticdoub1eAx=0.01;staticdoub1ec=3.9*100;铜的比热staticdoub1ep=8900;铜的密度staticdoub1eE；staticdoub1eac;staticdoub1eFoa*;staticdoub1eBih*;staticdoub1eTo=I00;staticdoub1e71=20;staticdoub1e=newdoub1e+1;pub1icstaticvoidmain(Stringargs)for(i=0;i=A/;i+)TiTo;System.out.printf(,%7.1fjT);System.out.printin();do/?+;for(i=1;i/V;i+)T=Fo*(-1+Ti+1)+(1-2*Fo)*;i=0;Ti=2*Fo*(+1-i)+Ti;i=N;T=2*Fo*(-i-1+*)+(1-2*B*Fo-2*Fo)*;if(e%7200=0)for(i=0;i=0.1);System,out.print1n(/?);6.结果120图.黄铜金属薄板表面温度随时间的变化曲线