传热学大作业——二维物体热传导问题的数值解法
1.二维热传导问题的物理描述:
本次需要解决的问题是结合给定的边界条件,通过二维导热物体的数值解法,求解出某建筑物墙角稳态下的温度分布t以及单位长度壁面上的热流量φ。
关于边界条件和研究对象选取的物理描述:如图所示为本次作业需要求解的建筑物墙壁的截面。尺寸如图中所标注。
由于墙角的对称性,A-A,B-B截面都是绝热面,并且由于对称性,我们只需要研究墙角的1/4即可(图中阴影部分)。假设在垂直纸面方向上不存在热量的传递,我们只需要对墙角进行二维问题的研究即可。
1.3 关于导热量计算截面的物理描述:本次大作业需要解决对流边界条件和等温边界条件下两类边界条件的问题。由于对称性,我们只需研究1/4墙角外表面和内表面的导热量再乘4,即是墙壁的总导热量。
2.二维热传导问题的数学描写:
本次实验的墙角满足二维,稳态无内热源的条件,因此:
壁面内满足导热微分方程:
?2
在绝热面处,满足边界条件:
-λ(
在对流边界处满足边界条件:
-
3.二维热传导问题离散方程的建立:
本次作业中墙角的温度场是一个稳态的连续的场。本次作业中将1/4墙角的温度场离散化,划分成若干小的网格,每个网格的节点看成以它为中心的一个小区域的代表。通过这些节点,采用“热平衡法”,建立起相应的离散方程,通过高斯-赛德尔迭代法,得到最终收敛的温度场,从而完成对墙角温度场的数值解。
对1/4墙角的网格划分如下:
选取步长Δx=
x,y坐标轴的方向如图所示,x,y轴的单位长度为步长Δx, 取左下角点为(1,1)点,其他点的标号为其在x,y
建立离散方程,要对导热物体中的节点根据其边界条件进行分类(特殊节点用阴影标出):首先以对流边界条件下的墙角为例
1.外壁面上,平直边界节点:
建立离散方程:
λΔ
以(i,j)为中心节点,进一步整理得:
t
2.外部角点:
建立离散方程:
ho·Δx
以(i,j)为中心节点,进一步整理得:
t
3.绝热+对流边界角点:
建立离散方程:
ho·
以(i,j)为中心节点,进一步整理得:
t
4.内部角点:
建立离散方程:
h
以(i,j)为中心节点,进一步整理得:
t
5.绝热平直边界节点:
建立离散方程:
λ
以(i,j)为中心节点,进一步整理得:
t
6.对于普通内部节点:
建立离散方程:
λ
以(i,j)为中心节点,进一步整理得:
t
等温边界条件下:等温边界下内部节点和绝热边界下的节点离散方程与上述5,6式形式相同,在等温壁面处,节点方程只需写成ti,j
4.方程的求解:
由上图可知,本题中有16*12=192个节点,相应地,就会有192个待求解的离散方程。在如此高阶次的方程组下,根据目前的计算机发展水平,采用克莱姆法则求解是不现实的,因此,采用方便计算机求解的高斯—赛德尔迭代法进行迭代求解。
根据数学上的“主对角线占优”原则,在我们采用热平衡法导出差分方程时,如果每一个方程都选用导出该方程的中心节点的温度作为迭代变量,那么迭代一定收敛。
在计算过程中往往需要进行足够多的次数,迭代才能收敛。判断收敛的方法是在相邻两次迭代值之差(或相对偏差)的绝对值足够小时,称已达到迭代收敛,迭代计算终止。
本次计算中采用绝对残差判据:
max?
下面是本次作业所采用的程序框图:
输入被模拟对象的物性参数和几何尺寸,确定步长根据设定的步长确定变量个数和分布,设定初场矩阵
输入被模拟对象的物性参数和几何尺寸,确定步长
根据设定的步长确定变量个数和分布,设定初场矩阵
开始
列
列写边界角点和特殊点的离散方程(同时迭代)
列
列写一般内部角点的离散方程(同时迭代)
残差是否满足要求?
残差是否满足要求?
否
,
是
作图,并
作图,并根据结果计算热流量
4.方程的求解(续):
对于对流边界条件下,内壁面和外壁面的热流量可以根据对流换热公式:
?
在等温边界条件下,由于环境未知,无法直接在等温壁面上进行计算。但由于稳态导热,可以借助等温壁面附近的截面(下图中标红)进行计算。公式为:
?
5.计算程序源代码:
(请见附件)
等温边界条件:
运算结果:
计算获得的各网格节点温度(分歧点用红色标出):
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0
30.0