传热学 matlab,数值传热-传热学作业-matlab.pdf

5. 厚度为250mm 无限大平壁，其导热系数λ=43+0.08t w/(m ·k)，平壁 侧温度为250℃，

另 侧温度为46 ℃，试用数值方法确定平壁内的温度分布，并确定通过该平壁的热流密度。

解：

(1) 建立模型

本题属于非常物性，无热源的 维稳态导热过程

将平壁沿厚度方向(x 方向)划分为N 个均匀相等的间距。

节点布置如图所示

1 2 3 4 N N+1

本题给出N=10。

(2 ) 通过热流法建立离散方程

a) 内节点离散方程

L P R

对节点 P(i)所代表的微元体，在x 方向上与节点 P 相邻的节点分别为 L(i-1)和

R(i+1)。由于节点之间的间距很小，可以认为相邻节点间的温度分布是线性的。

节点 P 所代表的网格单元与它周围各网格单元之间的导热量可根据傅里叶定

律直接写为:

t t t t

l i1 i 1 l i1 i 1

LP LP RP RP

x x

其中

(t t )

l 43 0.08 i i1

RP

2

(t t )

l 43 0.08 i1 i

LP

2

对节点P(i)所代表的微元写热平衡式，即可得节点P(i)温度的离散方程

0

LP RP

t t t t

l i1 i 1 l i1 i 1 0

LP RP

x x

l t l t

LP i1 LP i 1

t

i

l l

LP RP

b) 边界节点离散方程

由于本题的壁面温度属于第 类边界条件，因此

t 250 t 46

1 n1

c) 热流密度的计算

n

l (t t )

l RP i 1 i

q t i1

d d

(3 ) 编写流程图

计算机程序中使用的变量标识符： 计算机程序中输入数据：

i 节点的坐标变量 n 沿x 方向的网格划分数

t(i) 节点i 的温度 e 控制迭代过程终止的误差

tt 前 次算出的节点温度