2023年传热学大作业.pdf

传热学大作业一一二维物体热传导问题的数值解法1 .二维热传导问题的物理描述:本次需要解决的问题是结合给定的边界条件，通过二维导热物体的数值解法，求解出某建筑物墙角稳态下的温度分布t 以及单位长度壁面上的热流量中1.1 关于边界条件和研究对象选取的物理描述：如图所示为本次作业需规定解的建筑物墙壁的截面尺寸如图中所标注1.2 由于墙角的对称性,AA,B-B 截面都是绝热面，并且由于对称性，我们只需要研究墙角的1/4 即可（图中阴影部分）假设在垂直纸面方向上不存在热量的传递，我们只需要对墙角进行二维问题的研究即可1.3 关于导热量计算截面的物理描述:本次大作业需要解决对流边界条件和等温边界条件下两类边界条件的问题由于对称性，我们只需研究1/4墙角外表面和内表面的导热量再乘4,即是墙壁的总导热量2.二维热传导问题的数学描写：本次实验的墙角满足二维，稳态无内热源的条件，因此:壁面内满足导热微分方程:9dxz+w=dyz在绝热面处，满足边界条件:嗑=在对流边界处满足边界条件:一 人（氯”w-t f）3.二维热传导问题离散方程的建立:本次作业中墙角的温度场是一个稳态的连续的场本次作业中将1/4墙角的温度场离散化，划提成若干小的网格,每个网格的节点当作以它为中心的一个社区域的代表。

通过这些节点，采 用“热平衡法”，建立起相应的离散方程，通过高斯-赛德尔迭代法，得到最终收敛的温度场,从而完毕对墙角温度场的数值解对 1/4墙角的网格划分如下：(h i,t f 1/t w i)选取步长Ax=Ay=0.1m,为了方便研究，对导热物体的网格节点进行编码，编码规则如下：0 x,y坐标轴的方向如图所示，x,y 轴的单位长度为步长A x,取左下角点为（1,1）点，其他点的标号为其在x,y 轴上的坐标以此进行编码，进行离散方程的建立0建立离散方程，要对导热物体中的节点根据其边界条件进行分类（特殊节点用阴影标出）：一方面以对流边界条件下的墙角为例1 .外壁面上，平直边界节点:建立离散方程:(i,j-1)L,tfo a,j)入A y 7-+入 c A+入 c A+h o A x(t f 0 -&j)=0y A x 2 A y 2 A y v /0以(i,j)为中心节点，进一步整理得:_ ,(S ij-i+公+1)+入.o+i j +九tf02 2 +/i o%团建立离散方程:前 加 回 一3)+入 学 4+浮.=0,0 lJJ 2 Lx 2 A y以(i,j)为中心节点，进一步整理得：2.(t i+i j +G j-i)+h。

tfo“J X+ho-Lxo3.绝热+对流边界角点:建立离散方程：心 功 小，、工；）%叼+1 1%工 3G+IJ-以（i j）为中心节点，进一步整理得：_ /，g j+l +t i+l j）+八竽%1 x hA+ho,-24.内部角点:建立离散方程:同-处)+心.如 守+处.+冷.巧 卢以(i,j)为中心节点，进一步整理得:42，(G+i j +t/j-i)+t i j)+田 A x tf i3A+hi-A x5.绝热平直边界节点:建立离散方程:竺 知+1 1%竺 ttj_1-Gj1 Ay1 X jT+A A y ,=0以（i ,j）为中心节点，进一步整理得:6.对于普通内部节点:建立离散方程:A A x -乂+A A x -+A A y +=0 y A y%A x以（i,j）为中心节点，进一步整理得:幻_ +&-1J+L+lj)4A等温边界条件下：等温边界下内部节点和绝热边界下的节点离散方程与上述5,6式形式相同，在等温壁面处，节点方程只需写成%=tw即可4.方程的求解:由上图可知，本题中有16*1 2=1 92个节点，相应地,就会有1 9 2 个待求解的离散方程在如此高阶次的方程组下，根据目前的计算机发展水平，采用克莱姆法则求解是不现实的，因此，采用方便计算机求解的高斯一赛德尔迭代法进行迭代求解。

根据数学上的“主对角线占优”原则，在我们采用热平衡法导出差分方程时，假如每一个方程都选用导出该方程的中心节点的温度作为迭代变量，那么迭代一定收敛在计算过程中往往需要进行足够多的次数,迭代才干收敛判断收敛的方法是在相邻两次迭代值之差（或相对偏差）的绝对值足够小时，称已达成迭代收敛，迭代计算终止3本次计算中采用绝对残差判据：max|铲）-E0下面是本次作业所采用的程序框图:V输入被模拟对象的物性参数和V根据设定的步长拟定变量个数列写边界角点和特殊点的离散列写一般内部角点的离散方程a 0是V作图，并根据结果计算热流量4.方程的求解（续）：回对于对流边界条件下，内壁面和外壁面的热流量可以根据对流换热公式:巾=h Ay At+W h Ax At团在等温边界条件下，由于环境未知，无法直接在等温壁面上进行计算但由于稳态导热，可以借助等温壁面附近的截面（下图中标红）进行计算公式为:ZAt v AtA%,F/A Ay Ay 乙5.计算程序源代码:（请见附件）等温边界条件:运算结果：计算获得的各网格节点温度（分歧点用红色标出）：热流量:30.030.030.030.030.030.030.030.030.030.030.030.030.030.030.030.030.02 9.02 8.127.126.225.524.924.524.324.224.124.124.024.024.024.030.028.126.124.222.320.719.618.918.518.318.21 8.118.018.018.018.030.027.124.221.218.115.513.913.012.512.312.212.112.012.012.012.030.026.222.318.113.69.17.46.76.46.26.16.16.06.06.06.030.025.520.715.59.10.00.00.00.00.00.00.00.00.00.00.030.024.919.613.97.40.030.024.518.913.06.70.030.024.318.512.66.40.030.024.218.312.36.20.030.024.118.212.26.10.03 0.024.118.212.26.10.0弧=60.402=60.4根 据“热电模拟”实验获得的各网格节点温度（分歧点用红色标出）:30.030.030.030.030.030.030.030.030.030.030.030.030.030.030.030.030.029.028.127.12 6.225.524.924.524.324.224.124.124.024.024.024.030.028.126.12 4.222.320.719.618.918.518.318.218.11 8.118.018.018.030.027.124.221.218.115.513.913.012.512.312.212.112.112.01 2.012.030.026.222.318.113.69.17.46.76.46.26.16.16.06.06.06.030.025.520.715.59.10.00.00.00.00.00.00.00.00.00.00.030.024.919.61 3.97.40.030.024.518.913.06.70.030.024.318.512.66.40.030.024.218.312.36.20.030.2112.6.10.004.18.2230.24.1116.10.08.22.20热流量：仇=60.3502=60.35温度分布图像:对流边界条件:运算结果:计算获得的各节点温度:29.929.729.529.329.128.928.828.728.628.628.528.528.528.528.528.529.729.128.528.027.426.926.526.22 6.025.825.725.725.625.625.625.629.528.527.626.625.624.824.123.623.323.022.922.822.722.72 2.722.729.328.026.625.223.822.521.520.920.42 0.220.019.91 9.819.819.819.829.127.425.623.821.919.918.617.917.417.217.016.916.916.916.816.828.926.924.822.519.916.715.214.614.314.114.014.013.913.91 3.913.928.826.524.121.518.615.228.726.223.620.917.914.6热流量:28.626.023.32 0.517.514.322220.217.218.5.93.14.16225.2217.18.583.00.114.12 8.222.92011455.17.1.1.8弧=28.302=28.3根 据“热电模拟”实验获得的各网格节点温度（分歧点用红色标出）:29.929.729.529.329.229.028.828.728.628.628.528.528.528.528.528.529.729.128.6227.27.026.626.326225.825225.62525.685.05.9.75.6.6.02227.626.62224.1223.32222.822.822.822.7229.58.65.43.63.12.9.77.82928.02223.8221.520.2020.20.0119.919119.36.65.22.59.429.9.89.88热流量:29.127.425.623.821 .91 9.91 8.617.917.417.217.016.916.916.916.916.828.926.92 4.82 2.519.916.715.214.614.314.114.014.014.013.913.913.928.826.624.121.618.61 5.228.82 6.223.620.917.914.628,626.023.320.517.514.328.625.92 3.120.317.214.128.525.823.020.117.214.128.625.823.020.117.114.101=28.302=28.1分布图像:1012108A642028262422201816141210 x。