传热学实验—-墙角matlab导热问题

二维导热物体温度场的数值模拟姓名小明学号111111班级能动学院能动一、问题描述有一墙角模型，尺寸如图1所示，导热系数0.53W/(m·K)，墙角内外壁为第一类边界条件。求解该模型的温度分布及导热量。图1二、计算原理根据热平衡法列出节点方程，各方向导入单元体的热量之和为零。内节点和绝热边界点（图1点划线上的点）的方程形式不同。图2图2所示的内节点和绝热边界节点方程如下：内节点：q=00)()()()(1,,1,,1,1,,1,xyttxyttyxttyxttjijijijijijijijiWESN绝热边界点：0)(02)(2)(1,,1,1,,1,xyttyxttyxttjijijijijijiWESN三、计算过程用Matlab7.1语言编写计算程序，初取网格步长myx1.0运行结果：图1：各个点的温度数值图2：分层设色等温线分布图3：等温线分布（每两条线间隔为三度）四、小结本次数值模拟是运用matlab程序用于数值计算。小组成员共同讨论并复习了热传导问题的数学描述和热平衡法；从模拟过程中练习了不同节点迭代方程的建立；并简单学习了matlab语言的使用。这次大作业对于我们以后的学习和可能的研究来说是一个很好的锻炼机会。