计算传热学程序设计.docxVIP

中国石油大学 （华东 ） 储运与建筑工程学院热能与动力工程系 《计算传热学程序设计》设计报告 学生姓名 : 学 号： 专业班级： 指导教师 2012 年 7 月 7 日 1、设计题目 有一房屋的砖墙厚 8= m , /= W/（mrC ）, p =x 10 j/（ m3K），室内温度Tfi保持20C不变，表面传热系数hi=6W/（m2 C ）。开始时墙的温度处于稳定状态， 内墙表面温度Twi为15C寒潮入 侵后，室外温度Tf2下降为-10C,外墙的表面传热系数为 35W/（m2 C）。试分析寒潮入侵后多少 时间内墙壁面方可感受到外界气温的变化。 已知参数 壁厚，墙壁导热系数，密度与比热容的乘积，室内和寒潮入侵后室外空气温度，室内空气 和外墙的表面传热系数，开始时稳定状态下的内墙表面温度。 求解 寒潮入侵多少时间后内墙壁面可感受到外界气温的变化 2物理与数学模型 物理模型 该墙面为常物性，可以假设：（1）其为无限大平面，（2）只有在厚度方向传热，没有纵 向传热，则该问题转化为一维常物性无限大平面非稳态导热问题。 数学模型 以墙外表面为坐标原点，沿厚度方向为坐标正方向，建立坐标系。基于上述模型，取其在 匚（x 匚（ (1a) 初始条件: (1b)1 It=c _ *1 * hi (I fi * 1 i*iX $ (1b) 在两侧相应的边界条件是第三类边界条件，分别由傅立叶定律可描述如下: 左边界:TX右边界：x 0 h2(Tx 0 Tf2)(1c) 左边界: T X 右边界： x 0 h2(Tx 0 Tf2) (1c) x hi(Tx Tfi) (1d) 3数值处理与程序设计 数值处理 采用外点法用均匀网格对求解区域进行离散化，得到的网格系统如图 2所示。一共使用了 0?N-1共N个节点。 节点间距Sx为: (2a) (2b) (2c) 此例中墙壁导热系数为常值，无源项。则可采用有限体积法对控制方程离散化，得到离散 方程为： apTp aETE awTw b 式中： 0 ap aE aw ap a a a0 C x aE , aW , ap x x b aPTp0 (2d) 其中的上标“ 0表示此为上一时刻的值，分别为节点所在控制容积左右边界上的导热系数， 由于墙壁导热系数不变，故都等于 入△ T为时间步长。由元体能量平衡法可以得知左右边界节 点的离散方程分别为： 左边界节点: _ - — 1 - (3) 右边界节点： (扌十 hi - 盒)Tn-i ^jfTN-2 + hiTfi 十扌丁』】 (4) 离散方程的详细推导过程见附录。 程序设计 由物理模型可以知道本问题为一维导热问题，一维导热问题的离散方程在取遍所有节点后 形成的是三对角的代数方程组，采用追赶法进行求解。 程序构成和方法：程序由主程序和一个子程序构成。主程序进行变量定义和各已知参数的输 入，以及左右边界节点和内部节点控制方程的输入； 子程序tdma实现追赶法用来计算每个节点 新的温度。Thomas算法求解过程分为两步：消元和回代。消元是从系数矩阵的第二行起，逐一 将每一行的非零元素消去一个，使原来的三元方程化为二元方程。消元进行到最后一行时，二 元方程就化为一元方程，直接得到最后一个未知数的值。然后逐一往前回代，由各二元方程求 出其它未知解。 程序特点：该程序有很强的适应性，一维常物性非稳态平壁导热问题都可以使用此程序， 只要适当更改边值条件即可。还可以进行修改解决非常物性问题。 程序中对输出节点，最大输出量都进行了控制， 对计算结果的分析有很大帮助。 而且Thoms 算法的优点需要内存小，工作量小，程序设计简单。 程序流程图：首先对变量赋值，然后由初始条件建立初始温度场，接着从左边界，内部节 点，到右边界进行迭代，直到满足精度要求为止，最后输出结果，程序结束。程序流程如下图 3。 4、模型与程序验证 模型 本题简化为厚度为2 =的一维非稳态模型如图 4所示，初始温度为15 C,在其中间建立 坐标系，左两边为对流换热，且换热系数相同都为 h=25 W/(m2 C )，且流体温度Tf=-1OC对于x 0，列出其导热微分方程式及定解条件： T 亍 ——a— (0 x , 0) (5) x (6)Xx (6) X x (6)Xx (6) X x 图3程序流程图 T(x,O) T°(0 T(x,) h T( , ) T T(x, (8) 引入过余温度: T(x, ) T (9) 0 0 5C To=1 右侧: 侧:h=25 侧: W/(m2 :C ) 0 图4 一维导热简化模型 直接根据公式得到解析解如下: Cnexp( 0. n 1 2 nFO)COS( n ) (10) 式中，Fox 式中，Fo —，系数Cn应该使上述无穷级数在 0是满足初始条件，由傅里叶 级