方腔内自然对流MATLAB程序数值传热学
自然对流传热
问题描述:
一个二维矩形腔体,高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 (超松弛迭代)法求压力场
p 1 c p p 1 p