附录一一维不稳定导热的数值解法

1、附录一 一维不稳定导热的数值解法一维不稳定导热通常用有限差分法来求解。这种方法的基础是以不连续过程来代替空间里和时间上连续变化的过程。如将一维不稳定导热的微分方程式：用有限差量代替无限小量，则写为：这样解的精确性差，但简单得多，如果x与选择愈小，计算结果愈精确。采用有限差分法来计算平壁内温度分布的原理及方法叙述如下。如图1所示，一个厚度为x的平壁，将此平壁分为若干薄层，厚度各为x，在某一个时间内具有如图所示的温度分布。图 1 有限差分法图解 各平壁内部温度计算首先写出abcd区域的热平衡，这个区域的平均温度以t1表示，在平面ab上的温度降度为，在平面cd上的温度降度为.通过ab面的热流密度通过

2、cd面的热流密度，这个差额便是热流通过时留下来的部分，它用来加热物体，使abcd区域的平均温度由t1升至t1，其热平衡关系为:因为导温系数,即整理上式，则令傅立叶准数，为了保证差分方程正常收敛，必须遵守稳定性条件，即:取，则这就是说，只要知道前一个时间间隔相邻两层温度，就可用算术平均值的方法求出下一个时间间隔后该层温度。即该层温度为前一时间间隔相邻两层温度和的一半。若用n表示层数，k表示时间，写成通式为必须强调的是：因为取，当选定了x值后，就不能再任意选择，而是与被选定的x有如下关系： 各平壁之间的温度计算当平壁由两种或两种以上材料组成时，由于所取的时间间隔必须相同，因此当第一种材料xA确定后

3、，第二种材料xB不能任意选取，二者之间关系为故对于两种或两种以上材料组成的平壁，由于两种材料的性质不同，影响了交界面处的温度分布。如图1所示，由两种材料A和B组成的平壁，在同一时间间隔，ab表面的温度为tA，cd表面温度为tB,交界面的温度为t，经过一个时间间隔后交界面处的温度t可按如下方法求出： ， (a)其中： (b) （c）将(b)和(c)代入（a），并经整理得或式中：RA和RB分别为A、B两种材料的热阻。 外壁温度计算当已知气体温度和气体与壁面的热交换条件时，即第三类边界条件，固体壁面温度可按下述方法计算。由固体壁面向气体所传递的热流密度为：假定在最外一层平壁1x或(n-1) x中热流

4、不变，即从1x或（n-1）x处至表面的热流密度为经整理得式中例窑墙由耐火粘土砖、轻质砖和建筑砖砌成，厚度分别为0.23、0.115和0.24m，整个窑墙起始温度为20，外界空气温度为20且始终不变，窑内壁温度在8h内上升至820以后保温一个间隔时间，试计算此时窑墙的温度分布及蓄热量，窑墙材料物性参数见表：窑墙材料物性参数（W/（m）Cp（J/（kg）(kg/m3)耐火粘土砖=0.698+0.6410-3tCp=838+0.285t2000轻质砖=0.262+0.23210-3tCp=838+0.285t800建筑砖=0.465+0.5110-3tCp=838+0.285t1750解：假定内壁温

5、度从20上升至820以后保温一个间隔时间后墙内各材料的平均温度分别为240.9、41.95和21.2。则将耐火砖粘土砖分为5层.取n2=3取n3=6墙外表面的对流换热系数：工程计算中，常用下列简化公式计算平壁对无限大空间的自由对流换热系数式中K系数，垂直平壁K=2.56;热表面朝上K=3.26;热表面朝下K=1.63设外表面温度为20.01，则墙外表面的辐射换热系数：墙外表面的综合换热系数进行各节点温度计算:内壁温度经8h从20升至820，因此内壁升温速率内壁每时间间隔（0.625h）内温度升高1000.625=62.5。计算结果见表1。计算举例：（1）1 0x1时（2）3 1x1时4 1x1时（3）7 5x1（界面）时（4）8h的温度分布根据7.5h和8.125h的温度，采用内插法求解。（5）8h后保温1温度分布则以8h温度分布为基准再按有限差分法计算验算初始所假定的平均温度是否恰当：第一种材料耐火砖平均温度：相对误差第二种材料轻质砖平均温度：相对误差第三种材料建筑砖平均温度：相对误差相对误差E1