数值传热学4.pdfVIP

数值传热学 第四章 扩散方程的数值解及其应用(2) 主讲 陶文铨 西安交通大学能源与动力工程学院 热流中心 CFD-NHT-EHT CENTER 2011年10月10 日, 西安 1/60 第4 章教学目录 4.1 一维导热问题 4.2 多维非稳态导热的全隐格式 4.3 源项与边界条件的处理 4.4 求解代数方程的TDMA及ADI 方法 4.5 管道内充分发展对流换热概说 4.6 圆管内充分发展的对流换热 4.7长方形截面通道内充分发展的对流换热 2/60 4.4 求解代数方程的TDMA及ADI 方法 4.4.1 求解一维导热问题代数方程的三对角阵算法 1.一维导热问题代数方程通用形式 2.Thomas算法 3.第一类边界条件的处理 4.4.2 求解多维非稳态导热全隐格式的ADI方法 1.求解方法概述 2. Peaceman-Rachford的ADI迭代 3/60 4.4 求解代数方程的TDMA及ADI 方法 4.4.1 求解一维导热问题代数方程的三对角阵算法 1. 一维导热问题代数方程通用形式 每行三个未知数 稳态及非稳态隐式 （f0)都要联立求解一 组代数方程： a T a T a T b P P E E W W 其系数矩阵是一三对角阵 (Tri-diagonal matrix )。 4/60 2. Thomas算法的一般形式 将上式改写为： AT BT C T D , i 1,2,M 1 (a) i i i i1 i i1 i 端点条件：i=1, C =0; i＝M1, B ＝0 i i (1) 消元过程－把每行的未知数由三个减少为二个。 设消元后方程形式为： T P T Q (b) i1 i1 i i1 系数除到等号后的项 5/60 所谓消元就是要找出系数P ,Q 与A ,B ,C ,D 间的关系。 i i i i i i 将(b)乘以Ci,并与(a)式相加： AT BT C T D (a) i i i i1 i i1 i C T C P T C Q (b) i i 1 i i1 i i i1 AT C P T BT D C Q i i