参考2.多维导问题的数值解原理.docVIP

第二章 多维导热问题 2.1 二维非稳态导热全隐格式的通用离散方程 三种二维坐标系中的网格系统见下图2-1。采用控制容积积分法导出的离散方程以二维直角坐标系下的为例，根据二维非稳态导热方程： （2.1） 取全隐格式，假设节点之间温度线性分布，界面上热流密度均匀分布。非稳态项积分： 扩散项积分： 源项积分： 上述结果整理成： （2.2） 其中各系数为： （2.2） ， ，， （2.3a） （2.3b） （2.3c） 仍然需要记住，式（2.3a）表示的是各节点之间的热导（热阻的倒数），分子上的、代表的是各控制容积面上的面积；在二维问题中，的乘积是控制容积的体积。代表的是控制容积的热惯性。由此可见，利用上述系数计算式的物理含义，很容易写出三维导热问题的离散化方程及它的系数。 对于圆柱轴对称坐标和极坐标，同样可以利用系数的物理含义写出各系数计算式，离散方程与式（2.2）相同。不过要注意，在圆柱轴对称坐标中，选用一个弧度角的范围，极坐标取垂直于纸面一个单位长度（1m）。这样三种坐标系下的离散方程的系