matlab 空气粘度函数,方腔内自然对流MATLAB程序数值传热学.pdf

自然对流传热

问题描述：

一个二维矩形腔体，高H 1 ，宽L 1 ，物理模型如图1。

绝热壁面

热 冷

壁 壁

面 面

T=100 T=0

绝热壁面

图1 物理模型

上下壁面为绝热壁面，左壁面为热壁面T 100 ，右壁面为冷壁面T 0 。腔体内为初始

温度T 0 的空气，计算由此引起的自然对流传热。

0

数学模型：

控制方程：

使用Boussinesq 近似来考虑温差所引起的自然对流，在N-S 方程中添加一个体积力项。

u 1 2

uu p  u g T T (1)

 0 

t 

其中， 是初始温度， 是运动粘度， 是流体的体胀系数。

T  

0

能量守恒方程如下：

T 2

u T  T (2)

t

其中， 是热扩散系数。



边界条件：

由无滑移边界条件得，四个壁面上的速度均为零，即：un us uw ue 0 。

T

在热壁面上T 100 ，在冷壁面上T 0 ，在上下绝热壁面上处 0 。

y

数值处理：

区域离散化如图2 所示。

对于动量守恒方程，在不考虑压力的情况下先计算出一个临时速度

 n

u u n n 2 n

u u  u g T T (3)

 0 

dt

用SOR (超松弛迭代)法求压力场