第四章导热问题的数值方法.docVIP

5 热传导问题的数值方法 5.1一维稳态导热 一维稳态导热在直角坐标系下的控制方程可表示为： (5-1) 式中k为导热系数，T是温度，s是单位容积的热产生率。 首先选定控制和网格，如图5.1所示，并对方程（5-1）在所选定的控制体进行积分，即得： (5-2) 图5.1 控制和网格 然后进行离散化。如果用分线段性分布来计算方程（5-2）中的微商，那么最终的方程为: (5-3) 假设源项s在任一控制中之值可以表示为温度的线性函数，即，则导出的离散化方程为： (5-4) 式中 (5-5) 式(-4)就是一维稳态导热方程的离散形式，系数分别代表了节点P与E间及W与P间导热阻力的数，它们的大小反映了节点W和E处的温度对P点的影响程度。式中的是控制容积中的e和w界面上的当量导系数。计算时，物理参数值存储在节点的位置上。为了确定，还规定由节点上的物理量来计算相应界面上的量的方法。常用