c语言求解热传导方程,二维稳态导热问题的数值解法.docx

核科学与技术学院

《传热学》

二维稳态导热问题的

数值解法作业

姓名：罗晓

学号： 2014151214

班级：任课教师：李磊，张智刚

哈尔滨工程大学

核科学与技术学院

2016 年 11月 28 日

1

问题重述：

第一题：

如图所示，一个无限长矩形柱体，其横截面的边长分别为

L1 和

L2，常物性。该问题可视为二维稳态导热问题， 边界条件如图中所示，

其中 L1

， 2

，

w1

℃， w2

℃， λ

。

=0.6m L =0.4m

T =60

T

=20

=200 W/(mK)

1)

编写程序求解二维导热方程。

2)

绘制 x=L1/2 和 y=L2/2 处的温度场，并与解析解进行比较。已

知矩形内的温度场的解析解为：

t x, y

t w1

tw2 sin x L1

sh y

L1 。

sh L2

L1

第二题

将第一题中 y

L2 处的边界条件变为 t tw 2 ，其他条件不变。

1) 编写程序求解二维导热方程并