西安交通大学传热学大作业二维温度场热电比拟实验1
传热学大作业
传热学大作业——二维物体热传导问题的数值解法1.二维热传导问题的物理描述：本次需要解决的问题是结合给定的边界条件，通过二维导热物体的数值解法，求解出某建筑物墙角稳态下的温度分布t以及单位长度壁面上的热流量φ。
1.1关于边界条件和研究对象选取的物理描述：如图所示为本次作业需要求解的建筑物墙壁的截面。
尺寸如图中所标注。
1.2由于墙角的对称性，A-A，B-B截面都是绝热面，并且由于对称性，我们只需要研究墙角的1/4即可(图中阴影部分)。
假设在垂直纸面方向上不存在热量的传递，我们只需要对墙角进行二维问题的研究即可。
1.3 关于导热量计算截面的物理描述：本次大作业需要解决对流边界条件和等温边界条件下两类边界条件的问题。
由于对称性，我们只需研究1/4墙角外表面和内表面的导热量再乘4，即是墙壁的总导热量。
2.二维热传导问题的数学描写：本次实验的墙角满足二维，稳态无内热源的条件，因此：壁面内满足导热微分方程：∂2t ∂x2+∂2t∂x2=0。
在绝热面处，满足边界条件：−λ(∂t∂n)=0。
在对流边界处满足边界条件：−λ(∂t∂n )x=x(x x−x x)3.二维热传导问题离散方程的建立：本次作业中墙角的温度场是一个稳态的连续的场。
本次作业中将1/4墙角的温度场离散化，划分成若干小的网格，每个网格的节点看成以它为中心的一个小区域的代表。
通过这些节点，采用“热平衡法”，建立起相应的离散方程，通过高斯-赛德尔迭代法，得到最终收敛的温度场，从而完成对墙角温度场的数值解。
对1/4墙角的网格划分如下：选取步长Δx=Δy=0.1m，为了方便研究，对导热物体的网格节点进行编码，编码规则如下：x,y坐标轴的方向如图所示，x,y轴的单位长度为步长Δx, 取左下角点为(1,1)点，其他点的标号为其在x,y轴上的坐标。
以此进行编码，进行离散方程的建立。
建立离散方程，要对导热物体中的节点根据其边界条件进行分类（特殊节点用阴影标出）：首先以对流边界条件下的墙角为例1.外壁面上，平直边界节点：建立离散方程：λΔy x x+1,x−x x,xΔx+λΔx2x x,x+1−x x,xΔy+λΔx2x x,x−1−x x,xΔy+hoΔx(x xx−x x,x)=0以(i,j)为中心节点，进一步整理得：x x,x=x2·(x x,x−1+x x,x+1)+x·x x+1,x+xx·Δx·x xx2x+xx·Δx2.外部角点：建立离散方程：ho·Δx(x xx−x x,x)+λΔy2x x,x+1−x x,xΔx+xΔx2·x x,x−1−x x,xΔx=0以(i,j)为中心节点，进一步整理得：x x,x=x2·(x x+1,x+x x,x−1)+xx·Δx·x xxx+xx·Δx3.绝热+对流边界角点：建立离散方程：ho·Δx2·(x xx−x x,x)+xΔx2·x x,x+1−x x,xΔx+xΔx2·x x+1,x−x x,xΔx=0以(i,j)为中心节点，进一步整理得：x x,x=x2·(x x,x+1+x x+1,x)+xx·Δx2·x xxx+xx·Δx24.内部角点：建立离散方程：hi·Δx·(x xx−x x,x)+x·Δx·x x,x+1−x x,xΔx+xΔx·x x−1,x−x x,xΔx+xΔx2·x x+1,x−x x,xΔx+xΔx2·x x,x−1−x x,xΔx=0以(i,j)为中心节点，进一步整理得：x x,x=x2·(x x+1,x+x x,x−1)+x(x x,x+1+x x−1,x)+xx·Δx·x xx3x+xx·Δx5.绝热平直边界节点：建立离散方程：x Δx2·x x,x+1−x x,xΔx+xΔx2·x x,x−1−x x,xΔx+xΔx·x x−1,x−x x,xΔx =0以(i,j)为中心节点，进一步整理得：x x,x=x2·(x x,x−1+x x,x+1)+x·x x−1,x2x6.对于普通内部节点：建立离散方程：xΔx·x x,x+1−x x,xΔx+xΔx·x x,x−1−x x,xΔx+xΔx·x x−1,x−x x,xΔx +xΔxx x+1,x−x x,xΔx=0以(i,j)为中心节点，进一步整理得：x x,x=x·(x x,x−1+x x,x+1+x x−1,x+x x+1,x)4x等温边界条件下：等温边界下内部节点和绝热边界下的节点离散方程与上述5，6式形式相同，在等温壁面处，节点方程只需写成x x,x=x x即可4.方程的求解：由上图可知，本题中有16*12=192个节点，相应地，就会有192个待求解的离散方程。
传热学实验指导书
差和热系统中温度差的比例尺度，C= e1 − e2 ( v / °c )。当两个表面均为对流边界条件时， t1 − t2
C=
e∞1 t ∞1
− e∞2 − t∞2
(v/ °c)
，其中
e1, e2
分别为相应于外墙和内墙壁温的电势值，
e∞1, e∞ 2
为相应流体温度的电动势。也就是图中节点上的电动势。在先定比例系数后就可先定加 在电模型最外层两边界上电动势差值。利用系数可以从测得的电动势值换算相应的温度 值。
图 1-2
模拟墙角的几何尺寸为 L1 = 2.2m, L2 = 3.0m, L3 = 2.0m, L4 = 1.2m ;材料的导热系 数为 λ =0.53W\(m.K,).. 等温边界条件时，墙角外壁面温度 t1 = 30°c ，内壁面温度 t2 =0°c ；模拟墙角两端应
维持 2V 的电压差，电压、温度比例系数 c1
相应的电网络节点上的电动势方程为
ei−1, j − ei, j + ei, j−1 − ei, j + ei+1, j − ei, j + ei, j+1 − ei, j = 0 (1-4)
R1
R2
R3
R4
图 1—1
只要满足 R1 = R2 = R3 = R4 = R
的条件，式（ 1-3）和式（1-4 ）完全类似。
导热现象和导电现象之间的相似之处可以从他们的数学描写式可以看出。 在导热系统中，二维稳定导热微分方程为：
∂ 2t + ∂ 2t = 0 (1-1) ∂x 2 ∂y 2 在导电系统中，二维稳定导电微分方程为：
∂ 2e + ∂ 2e = 0 (1-2) ∂x 2 ∂y 2
西安交通大学传热学大作业
《传热学》上机大作业二维导热物体温度场的数值模拟学校：西安交通大学姓名：张晓璐学号:10031133班级：能动A06一．问题（4-23）有一个用砖砌成的长方形截面的冷空气通道，形状和截面尺寸如下图所示，假设在垂直纸面方向冷空气和砖墙的温度变化很小，差别可以近似的予以忽略。
在下列两种情况下计算：砖墙横截面上的温度分布；垂直于纸面方向上的每米长度上通过墙砖上的导热量。
第一种情况：内外壁分别维持在10C ︒和30C ︒第二种情况：内外壁与流体发生对流传热，且有C t f ︒=101，)/(2021k m W h ⋅=，C t f ︒=302，)/(422k m W h ⋅=，K m W ⋅=/53.0λ二．问题分析 1.控制方程02222=∂∂+∂∂ytx t 2.边界条件所研究物体关于横轴和纵轴对称，所以只研究四分之一即可，如下图：对上图所示各边界：边界1：由对称性可知：此边界绝热，0=w q 。
边界2：情况一：第一类边界条件C t w ︒=10情况二：第三类边界条件)()(11f w w w t t h ntq -=∂∂-=λ 边界3：情况一：第一类边界条件C t w ︒=30情况二：第三类边界条件)()(22f w w w t t h ntq -=∂∂-=λ 三：区域离散化及公式推导如下图所示，用一系列和坐标抽平行的相互间隔cm 10的网格线将所示区域离散化，每个交点可以看做节点，该节点的温度近似看做节点所在区域的平均温度。
利用热平衡法列出各个节点温度的代数方程。
第一种情况： 内部角点：11~8,15~611~2,5~2)(411,1,,1,1,====++++=+-+-n m n m t t t t t n m n m n m n m n m 平直边界1：11~8),2(415~2),2(411,161,16,15,161,11,12,1,=++==++=+-+-n t t t t m t t t t n n n nm m m m平直边界2：7,16~7,107~1,6,10,,======n m t n m t n m n m平直边界3：12,16~2,30;12~1,1,30,,======n m t n m t n m n m第二种情况： 内部角点：11~8,15~611~2,5~2)(411,1,,1,1,====++++=+-+-n m n m t t t t t n m n m n m n m n m 平直边界1：11~8),2(415~2),2(411,161,16,15,161,11,12,1,=++==++=+-+-n t t t t m t t t t n n n nm m m m平直边界2：7,16~7206~1,61.0,10,)2(222111111,1,,1,======∆=∆︒=+∆∆+++=-+-n m h n m m y x C t xh t xh t t t t f f n m n m n m n m λλ平直边界3：12,16~2411~1,11.0,30,)2(222222221,1,,1,======∆=∆︒=+∆∆+++=-+-n m h n m m y x C t xh t xh t t t t f f n m n m n m n m λλ内角点：20,10,)3(22)(2111116,67,78,67,57,6=︒=+∆∆++++=h C t xh t xh t t t t t f f λλ外角点：4,30,)1(222222211,112,212,1=︒=+∆∆++=h C t xh t x h t t t f f λλ4,30,2222222,11,21,1=︒=+∆∆++=h C t xh t xh t t t f f λλ4,30,22222212,1511,1612,16=︒=+∆∆++=h C t xh t xh t t t f f λλ20,10,2111112,61,51,6=︒=+∆∆++=h C t xh t xh t t t f f λλ20,10,2111118,167,157,16=︒=+∆∆++=h C t xh t xh t t t f f λλ四．编程计算各节点温度和冷量损失（冷量推导在后面）（用fortran编程）由以上区域离散化分析可以得到几十个方程，要求解这些方程无疑是非常繁琐的，所以采用迭代法，用计算机编程求解这些方程的解，就可以得到各点温度的数值。
西安交通大学传热学上机实验报告
φ1 − φ2 E= （φ1 + φ2）2
三、计算过程
用 MATLAB 编写计算程序，取网格步长 ∆x = ∆y = 0.1m 。 1、第一类边界条件 （1）运行程序 1（见附录 1） ，得到等温边界条件下计算墙角温度分布图：
图 4 等温边界条件下计算等温线分布（左图中每两条线间隔为三摄氏度） 运行程序 2（见附录 2） ，得到等温边界条件下实测墙角温度分布图：
s1=0; for i=2:11 s1=s1+(30-T(i,2))*0.53; end for j=2:15 s1=s1+(30-T(11,j))*0.53; end s1=s1+(30-T(1,2))*0.53/2+(30-T(11,16))*0.53/2
%墙角外侧换热量
s2=0; for i=2:6 s2=s2+T(i,5)*0.53; end for j=7:15 s2=s2+T(8,j)*0.53; end s2=s2+T(1,5)*0.53/2+T(8,16)*0.53/2+T(7,5)*0.53/2+T(8,6)*0.53/2 %墙角内侧换热量 s=2*(s1+s2) %单位长度墙壁的总换热量 e=abs(s1-s2)/((s1+s2)/2)
图3
内节点和绝热边界
图 3 所示的内节点和绝热边界节点方程如下： 内节点：
⎡(t −t )∆x (t −t )∆x (t −t )∆y (t −t )∆y⎤ ΦN +ΦS +ΦE +ΦW = λ⋅1⋅ ⎢ i, j+1 i, j + i, j−1 i, j + i+1, j i, j + i−1, j i, j ⎥ = 0 ∆y ∆y ∆x ∆x ⎣ ⎦
传热学大作业报告二维稳态导热
传热学大作业报告二维稳态导热二维稳态导热大作业报告导热问题是传热学中非常重要的一个研究领域。
在导热问题中，我们研究的是物体内部的温度分布、热流分布以及热传导过程。
本次大作业中，我们将研究一个二维稳态导热问题，分析材料内部的温度分布情况。
在二维稳态导热问题中，我们假设热传导发生在一个二维平面内，而且热流只在平面内的两个方向上进行。
我们的目标是研究材料内部的温度分布情况，并找到材料内各个位置的温度。
为了研究这个问题，我们首先需要建立热传导的数学模型。
根据热传导方程，在稳态下，热传导的速率是不变的。
假设材料在x和y两个方向上的热传导系数分别为kx和ky，温度分布函数为T(x, y)，则可以得到以下的二维热传导方程：kx * d^2T/dx^2 + ky * d^2T/dy^2 = 0这是一个二维的亥姆霍兹方程，我们可以通过求解它来得到材料内部的温度分布。
为了进一步分析问题，我们对热传导方程进行了无量纲化处理。
使用无量纲化可以简化计算，并且使得结果更加清晰。
我们引入了一个无量纲化的温度变量θ，通过以下公式进行计算：θ=(T-T0)/(T1-T0)其中T是位置(x,y)处的温度，T0是最低温度，T1是最高温度。
这样处理之后，热传导方程可以写成：d^2θ/dx^2 + σ * d^2θ/dy^2 = 0其中σ = ky / kx 是无量纲化的热传导比。
为了求解这个二维亥姆霍兹方程，我们使用了有限差分法。
首先将平面划分成一个个小的网格单元，然后在每个网格单元中，使用二阶中央差分法对方程进行离散化。
最终得到一个线性方程组，可以通过求解该方程组，得到无量纲温度分布。
为了验证我们的计算结果，我们将研究一个简单的导热问题，即一个正方形材料中心局部加热的情况。
我们假设正方形材料的一部分区域中心加热，其余区域保持恒定温度。
我们通过计算得到了材料内部的温度分布，并且将结果与理论解进行了比较。
通过对比发现，计算结果与理论解非常吻合，验证了我们的计算方法的准确性和可靠性。
传热学实验指导书---实验一(本部)
实验一 非稳态法测量材料的导热性能实验一、实验目的1. 快速测量绝热材料的导热系数和比热。
2. 掌握使用热电偶测量温差的方法。
二、实验原理X图1 第二类边界条件无限大平板导热的物理模型本实验是根据第二类边界条件，无限大平板的导热问题来设计的。
设平板厚度为2δ。
初始温度为t 0，平板两面受恒定的热流密度q c 均匀加热（见图1）。
求任何瞬间沿平板厚度方向的温度分布t(x,τ)。
导热微分方程式、初始条件和第二类边界条件如下：22),(),(x x t a x t ∂∂=∂∂τττ初始条件 0)0,(t x t =边界条件x=0,0),0(=∂∂xt τX=δ,0),(=+∂∂λτδcq x t 方程的解为：⎥⎦⎤⎢⎣⎡--+--=-∑∞=+1221220)exp(cos(2)1(63),(n o n n n n c F x x a q t x t μδμμδδδδλττq c式中： t —温度； τ—时间； t 0 — 初始温度；ɑ — 平板的导温系数； μn — n π n=1,2,3,……2δτa Fo =— 傅立叶准则； q c— 沿方向从端面向平板加热的恒热流密度；随着时间t 的延长，Fo 数变大，上式中级数和项愈小。
当Fo>0.5时，级数和项变得很小，可以忽略，上式变成：)612(),(220-+-=-δτδτλδτa a q t x t c 由此可见，当Fo>0.5后，平板各处温度和时间成线性关系，温度随时间变化的速率是常数，并且到处相同。
这种状态称为准稳态。
在准稳态时，平板中心面x=0处的温度为：)61(),0(20-=-δτλδτa q t t c 平板加热面X=δ处为：)31(),(20+-=-δτλδτδa q t t c 此两面的温差为：λδττδcq t t t 21),0(),(=-=∆如已知q c 和δ，再测出t ∆，就可以由上式求出导热系数：tq c∆=2δλ式中，λ—平板的导热系数，oW /(m C)⋅ cq —沿x 方向给平板加热的恒定热流密度，2W /mδ—平板的厚度，mt ∆—平板中心面x=0处和平板加热面x=δ处两面的温差，o C又，根据热平衡原理，在准稳态有下列关系：式中，F —平板的横截面积ρ—试件材料的密度C —试件材料的比热—准稳态时的温升速率由上式可求得比热为：实验时， 以试件中心处为准。
西安交通大学数值传热学大作业
5
数值传热学论文
能量方程： u 边界条件：
T T 2T 2T v ( 2 2 ) x y x y
T ( x, Tp ) T ( x,0)
(4)
(5) (6) (7) (8) (9) (10)
u( x, Tp ) u( x,0)
v( x, T v( x, 0 ) p )
图 2 网格划分
百叶窗翅片通道内周期性充分发展流动与换热的控制方程如下： 连续性方程： 动量方程： u
u v 0 x y
(1) (2) (3)
u u 1 p 2 u 2u v v( 2 2 ) x y x x y u v v 1 p 2v 2v v v( 2 2 ) x y y x y
2
数值传热学论文
主要符号表
f
Nu m
摩擦因数 平均 Nusselt 数 Prandtl（普朗特）数 雷诺数 竖直平板和封闭方腔壁面间的距离,热扩散系数（定义 u r ） 表面换热系数 导热系数 温度 平均温度 内部翅片的温度 W/(m2℃) W/(m℃) ℃ ℃ ℃ m/s m m m Pa W/m2
数值传热学论文
数值传热学大作业
1
数值传热学论文
百叶窗翅片流动换热的数值模拟
（西安交通大学能源与动力工程学院，710049，西安） 摘要： 针对具有一定倾斜角度的流动和换热已经进入周期性充分发展的百叶窗换 热问题，在稳态、层流、常物性和翅片温度恒定的条件下，采用 SIMPLER 算法， 对百叶窗的一个翅片单元进行了数值模拟计算。在翅片倾角θ =25°，雷诺数 Re 在 10 到 500 范围内变化时，得到了平均 Nusselt 数与阻力系数 f 的计算结果。计 算结果表明：随着 Re 的增大，平均 Nusselt 数逐渐增大，f 却随之逐渐减少。 关键词：百叶窗；周期性发展；数值模拟；SIMPLER 算法
西安交大传热学上机实验报告
西安交⼤传热学上机实验报告传热学上机实验报告⼆维导热物体温度场的数值模拟学院：化⼯学院姓名：沈佳磊学号：2110307016班级：装备11⼀、物理问题有⼀个⽤砖砌成的长⽅形截⾯的冷空⽓空道，其截⾯尺⼨如下图所⽰，假设在垂直于纸⾯⽅向上冷空⽓及砖墙的温度变化很⼩，可以近似地予以忽略。
在下列两种情况下试计算：（1）砖墙横截⾯上的温度分布；（2）垂直于纸⾯⽅向的每⽶长度上通过砖墙的导热量。
外矩形长为3.0m，宽为2.2m；内矩形长为2.0m，宽为1.2m。
第⼀种情况：内外壁分别均匀地维持在0℃及30℃；第⼆种情况：内外表⾯均为第三类边界条件，且已知：外壁：30℃，h1=10W/m2·℃,内壁：10℃，h2= 4 W/m2·℃砖墙的导热系数λ=0.53 W/m·℃由于对称性，仅研究1/4部分即可。
⼆、数学描写对于⼆维稳态导热问题，描写物体温度分布的微分⽅程为拉普拉斯⽅程22220t t x x ??+=??这是描写实验情景的控制⽅程。
三、⽅程离散⽤⼀系列与坐标轴平⾏的⽹格线把求解区域划分成许多⼦区域，以⽹格线的交点作为确定温度值的空间位置，即节点。
每⼀个节点都可以看成是以它为中⼼的⼀个⼩区域的代表。
由于对称性，仅研究1/4部分即可。
依照实验时得点划分⽹格。
建⽴节点物理量的代数⽅程对于内部节点，由?x=?y ，有,1,1,,1,11()4m n m n m n m n m n t t t t t +-+-=+++由于本实验为恒壁温，不涉及对流，故内⾓点，边界点代数⽅程与该式相同。
设⽴迭代初场，求解代数⽅程组图中，除边界上各节点温度为已知且不变外，其余各节点均需建⽴类似3中的离散⽅程，构成⼀个封闭的代数⽅程组。
以t ? =0°C 为场的初始温度，代⼊⽅程组迭代，直⾄相邻两次内外传热值之差⼩于0.01，认为已达到迭代收敛。
四、编程及结果program mainimplicit nonereal ,dimension(1:16,1:12)::treal ,dimension(1:16,1:12)::t1real q,q1,q2,q3,q4,q5,q6,q7,q8,q9,q10,q11,a integer m,n,z logical::converged=.false.z=1t=0a=0.53do n=1,12t(1,n)=30end dodo m=2,16t(m,12)=30end dodo n=1,7t(6,n)=0end dodo m=7,16t(m,7)=0end dodo while(.not.converged.and.z<10000)t1=tdo m=2,5do n=1,11if( n==1 )thent(m,n)=0.25*(t(m-1,n)+t(m+1,n)+2*t(m,n+1))elset(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n-1)+t(m,n+1)) end if end doend dodo n=8,11do m=6,16if (m==16) thent(m,n)=0.25*(t(m,n-1)+t(m,n+1)+2*t(m-1,n)) elset(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n-1)+t(m,n+1)) end if end doend doz=z+1do m=1,16do n=1,12if(abs(t(m,n)-t1(m,n))>0.000001) thenconverged=.false.exitelseconverged=.true.end ifend doend doend dowrite(*,'(16f5.1)',advance='no')((t(m,n),m=1,16),n=12,7,-1) write(*,*) write(*,'(6f5.1)',advance='no')((t(m,n),m=1,6),n=6,1,-1)do n=2,11q1=(t(1,n)-t(2,n))*a+q1end dodo m=2,15q2=(t(m,12)-t(m,11))*a+q2end doq3=(t(1,1)-t(2,1))*a*0.5q4=(t(16,12)-t(16,11))*a*0.5q10=q1+q2+q3+q4write(*,*)do n=2,6q5=(t(5,n)-t(6,n))*a+q5end dodo m=7,15q6=(t(m,8)-t(m,7))*a+q6end doq7=(t(5,1)-t(6,1))*a*0.5q8=(t(16,8)-t(16,7))*a*0.5q9=(t(5,7)-t(6,7))*a*2q11=q5+q6+q7+q8+q9q=(q10+q11)*0.5*4print*,"外表⾯导量=",q10,"内表⾯导热量",q11,"每⽶⾼砖墙导热量",q end结果截图：将以上结果⽤matlab画图⼯具绘制出如下图像：。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
二维导热物体温度场的数值模拟
一、物理问题
有一个用砖砌成的长方形截面的冷空气通道， 于纸面方向上用冷空气及砖墙的温度变化很小， 可以近似地予以忽略。
在下列两种情况下试计算： 砖墙横截面上的温度分布；垂直于纸面方向的每 米长度上通过砖墙的导热量。
第一种情况：内外壁分别均匀维持在 0℃及
30℃； 第二种情况：内外壁均为第三类边界条
件， 且已知：
t 1 30 C,h 1 10.35W / m 2
K
2
t 2 10 C, h 2 3.93W / m 2
K
砖墙导热系数 0.35/ m K
二、数学描写
由对称的界面必是绝热面， 态、无内热源的导热问题。
控制方程：
22
tt 22 xy
边界条件： 第一种情况：
由对称性知边界 1 绝热：
边界 2 为等温边界，满足第一类边界条件： t w 0 C ； 边界 3 为等温边界，满足第一类边界条件： t w 30 C 。
第一种情况：
由对称性知边界 1 绝热： q w 0；
边界 2 为对流边界，满足第三类边界条件： q w
( t
)w h 2(t w
可取左上方的四分之一墙角为研究对象， 该问题为二维、 稳
图1-
t f )；
n
t 边界3 为对流边界，满足第三类边界条件：q w
( )
w h
2
(t
w
t
f
)。
w
n
w 2 w f
0,m 6,n 1~ 7;m 7 ~ 16,n 7
30,m 1,n 1~12;m 2 ~ 16,n 12
三、方程离散
用一系列与坐标轴平行的间隔 0.1m 的二维网格线 将温度区域划分为若干子区域，如图 1-3 所示。
采用热平衡法， 利用傅里叶导热定律和能量守恒定 律，按照以导入元体（ m,n ）方向的热流量为正，列写 每个节点代表的元体的代数方程， 第一种情况： 边界点： 1
边界 绝热边界） ： 边界
图1-3
t
m ,1
t
16,n
等温内边界） ： 14
(2t m,2 1 4
(2t 15,n
t m 1,1 t m 1,1)，m 2 ~ 5
t 16,n 1 t 16,n 1)，
n 8 ~ 11
边界
等温外边界） ： 内节
点：
1
(t t t t )
4
m 1,n m 1,n m ,n 1 m,n 1
m 2 ~ 5,n 2 ~11;m 6 ~ 15,n 8 ~ 11
t
m,n
第二种情况 边界点： 边界 1（绝热边界） ：
t
m ,1
1
4
(2t m,2 t m 1,1 t m 1,1)，m 2 ~ 5
t
16,n
1
4
(2t 15,n t 16,n 1 t 16,n 1)，
n 8 ~11 4
边界 2（内对流边界） ：
t6,n
2t 5,n t 6,n 1 t 6,n 1 2Bi 1t 1 ，n 1~ 6
6,n
2(Bi 2)
t m,n
t
m,n
2t t t 2Bi t
t m,72t m,8 t m 1,7 t m 1,7 2Bi1t1，
m 7~16
边界3
(
外对流边界) ：
t 2t2,n t1,n 1 t1,n 1 2Bi 2t2 ，
n 1~11
t1,n ，n 1~11
1,n 2(Bi
2 2)
2t m,11 t m 1,12 t m 1,12 2Bi2t2
t m,12 ，m 2~16 内角
点：
2(t5,7 t6,8) t7,7 t6,6 2Bi 1t1 t6,7
6,7 2(Bi
1 3)
外角
点：
t
t2,12 t1,11 2Bi 2t2 t1,122(Bi 1)
内节
点：
t 1(t t t t );
t m,n (t m 1,n t m 1,n t m,n 1 t m,n 1);
4
m 2 ~5,n 2~11;m 6~15,n 8~11 ( Bi 12h2 x,t1
t 210 ；Bi 22h1 x,t2t 130 )
四、编程思路及流程图
编程思路为设定两个二维数组t(i,j)、ta(i,j) 分别表示本次迭代和上次迭代各节点的温度值，iter(实际编程时并未按照此名称来命名迭代步长)表示迭代进行的次数, Q1 、Q2 分
别表示外边界、内边界的散热量。
开始时，给t(i,j)、ta(i,j) 赋相同的初始值，t(i,j) 根据内节
点和各边界节点的离散方程进行迭代，迭代后比较t(i,j)、ta(i,j) 各个节点之间温度之差，若两个温度之差小给定的精度，则此时迭代完成，t(i,j) 就是所求的温度场分布，若两温度之差
不满足精度要求，则将t(i,j )的值赋给ta(i,j)，t(i,j)继续迭代，直到二者各个点的温度之差满足精度要求，记下此时的迭代次数，并根据所得到的温度场分布计算内外边界上散热量以及偏
差。
图1-4
五、结果讨论本次实验的实际边界条件为等温边界条件，因此以第一种情况为依据进行实验编程。
程序（见附注）运行结果如图1-5 所示。
图1-5
使用MATLAB7.0 软件将各节点温度作为原始数据做出四分之一墙角的温度场分布图
图1-6
图1-7
将实验测量结果和数值计算结果进行比较如表1-1
对比两种结果分析数值计算方法得到的散热量比实验测得的散热量相近，而前者的相对平衡偏差比较小。
前者的误差主要来自于迭代次数和网格数的划分，但也可能存在网格独立解，可以进一步细化网格进行求解观察数值模拟结果；后者的误差主要来自于电阻的阻值误差和读数仪表的基本误差，提高标准电阻的准确性和读数仪表的精度等级可以减小实验误差。