计算传热学程序设计

《计算传热学程序设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、中国石油大学(华东)储运与建筑工程学院热能与动力工程系《计算传热学程序设计》设计报告学生姓名:学号：专业班级：指导教师2012年7月7日211、设计题目有一房屋的砖墙厚δ=0.3m，λ=0.85W/(m·℃)，ρc=1.05×106J/(m3·K)，室内温度Tf1保持20℃不变，表面传热系数h1=6W/(m2·℃)。开始时墙的温度处于稳定状态，内墙表面温度Tw1为15℃寒潮入侵后，室外温度Tf2下降为-10℃，外墙的表面传热系数为35W/(m2·℃)。试分析寒潮入侵后多少时间内墙壁面方可感受到外界气温的变化。室外寒

2、流入侵室内0x图1墙壁简化图1.1已知参数壁厚，墙壁导热系数，密度与比热容的乘积，室内和寒潮入侵后室外空气温度，室内空气和外墙的表面传热系数，开始时稳定状态下的内墙表面温度。1.2求解寒潮入侵多少时间后内墙壁面可感受到外界气温的变化？2物理与数学模型2.1物理模型该墙面为常物性，可以假设：（1）其为无限大平面，（2）只有在厚度方向传热，没有纵向传热，则该问题转化为一维常物性无限大平面非稳态导热问题。2.2数学模型以墙外表面为坐标原点，沿厚度方向为坐标正方向，建立坐标系。基于上述模型，取其在x方向上的微元作为研究对象

3、，则该问题的数学模型可描述如下：（1a）初始条件：T

4、τ=0=Tw1-h1(Tf1-Tw1)(δ-x)/λ（1b）在两侧相应的边界条件是第三类边界条件，分别由傅立叶定律可描述如下：左边界：21（1c）右边界：（1d）3数值处理与程序设计3.1数值处理采用外点法用均匀网格对求解区域进行离散化，得到的网格系统如图2所示。一共使用了0~N-1共N个节点。节点间距δx为:δx=δN-1图2墙壁内的网格划分此例中墙壁导热系数为常值，无源项。则可采用有限体积法对控制方程离散化，得到离散方程为：（2a）式中：（2b），，（2c）

5、（2d）其中的上标“0”表示此为上一时刻的值,分别为节点所在控制容积左右边界上的导热系数，由于墙壁导热系数不变，故都等于λ，△τ为时间步长。由元体能量平衡法可以得知左右边界节点的离散方程分别为：左边界节点：21ap02+h2+λδxT0=λδxT1+h2Tf2+ap02T00（3）右边界节点：ap02+h1+λδxTN-1=λδxTN-2+h1Tf1+ap02TN-10（4）离散方程的详细推导过程见附录。3.2程序设计由物理模型可以知道本问题为一维导热问题，一维导热问题的离散方程在取遍所有节点后形成的是三对角的代数

6、方程组，采用追赶法进行求解。程序构成和方法:程序由主程序和一个子程序构成。主程序进行变量定义和各已知参数的输入，以及左右边界节点和内部节点控制方程的输入；子程序tdma实现追赶法用来计算每个节点新的温度。Thomas算法求解过程分为两步:消元和回代。消元是从系数矩阵的第二行起，逐一将每一行的非零元素消去一个，使原来的三元方程化为二元方程。消元进行到最后一行时，二元方程就化为一元方程，直接得到最后一个未知数的值。然后逐一往前回代，由各二元方程求出其它未知解。程序特点：该程序有很强的适应性，一维常物性非稳态平壁导热问题

7、都可以使用此程序，只要适当更改边值条件即可。还可以进行修改解决非常物性问题。程序中对输出节点，最大输出量都进行了控制，对计算结果的分析有很大帮助。而且Thoms算法的优点需要内存小，工作量小，程序设计简单。程序流程图：首先对变量赋值，然后由初始条件建立初始温度场，接着从左边界，内部节点，到右边界进行迭代，直到满足精度要求为止，最后输出结果，程序结束。程序流程如下图3。4、模型与程序验证4.1模型本题简化为厚度为2=0.3m的一维非稳态模型如图4所示，初始温度为15℃，在其中间建立坐标系，左两边为对流换热，且换热系数

8、相同都为h=25W/(m2·℃)，且流体温度Tf=-10℃对于x0，列出其导热微分方程式及定解条件：（5）21开始进行参数赋值用迭代法求解温度场计数并判断是否完成温度场计算输入至屏幕输入至文件程序结束NY图3程序流程图（6）（7）（8）引入过余温度：（9）21图4一维导热简化模型0左侧：hTf右侧：h=25W/(m2·℃)Tf=-10℃δx-δT0=15℃直接根据公式得到解析解如下：（10）式中，，系数应该使上述无穷级数在是满足初始条件，由傅里叶级数理论可得：（11）是超越方程的根，称为特征根。…（12）其中。4.

9、2程序验证（1）由模型可以得到相关信息然后进行编程，同等时间下计算出中心处温度的解析解和数值解进行比较，数据记录在表1。然后计算出相对误差，作图5，观察数值解与分析解的比较曲线。由图表中可以发现，平壁中心不同时刻温度值的分析解和数值解相差不是很大，二者吻合的比较好，可以说明所编制的数值解法的程序是正确的。相对误差先增大后减小，增大的原因是此时温度接近零度，相