传热学大作业报告 二维稳态导热

传热学大作业报告二维稳态导热二维稳态导热大作业报告导热问题是传热学中非常重要的一个研究领域。

在导热问题中，我们研究的是物体内部的温度分布、热流分布以及热传导过程。

本次大作业中，我们将研究一个二维稳态导热问题，分析材料内部的温度分布情况。

在二维稳态导热问题中，我们假设热传导发生在一个二维平面内，而且热流只在平面内的两个方向上进行。

我们的目标是研究材料内部的温度分布情况，并找到材料内各个位置的温度。

为了研究这个问题，我们首先需要建立热传导的数学模型。

根据热传导方程，在稳态下，热传导的速率是不变的。

假设材料在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 是无量纲化的热传导比。

为了求解这个二维亥姆霍兹方程，我们使用了有限差分法。

首先将平面划分成一个个小的网格单元，然后在每个网格单元中，使用二阶中央差分法对方程进行离散化。

最终得到一个线性方程组，可以通过求解该方程组，得到无量纲温度分布。

为了验证我们的计算结果，我们将研究一个简单的导热问题，即一个正方形材料中心局部加热的情况。

我们假设正方形材料的一部分区域中心加热，其余区域保持恒定温度。

我们通过计算得到了材料内部的温度分布，并且将结果与理论解进行了比较。

通过对比发现，计算结果与理论解非常吻合，验证了我们的计算方法的准确性和可靠性。

西安交⼤传热学上机实验报告传热学上机实验报告⼆维导热物体温度场的数值模拟学院：化⼯学院姓名：沈佳磊学号：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画图⼯具绘制出如下图像：。

核科学与技术学院《传热学》二维稳态导热问题的数值解法作业姓名：罗晓学号：2014151214班级：20141512任课教师：李磊，张智刚哈尔滨工程大学核科学与技术学院2016年11月28日问题重述：第一题：如图所示，一个无限长矩形柱体，其横截面的边长分别为L 1和L 2，常物性。

该问题可视为二维稳态导热问题，边界条件如图中所示，其中L 1=0.6m ，L 2=0.4m ，T w 1=60℃，T w 2=20℃，λ=200W/(mK)。

1）编写程序求解二维导热方程。

2）绘制x =L 1/2和y=L 2/2处的温度场，并与解析解进行比较。

已知矩形内的温度场的解析解为：()()()()1211w2w1sh sh sin ,L L L y L x t t y x t πππ+=。

第二题将第一题中2y L =处的边界条件变为2w t t =，其他条件不变。

1)编写程序求解二维导热方程并计算从y =0处导入的热量2Φ。

2)当21L L 时，该二维导热问题可简化为一维导热问题。

在一维的近似下，试计算从y =0处导入的热量1Φ，并比较不同L 2/L 1下21ΦΦ的比值。

由该问题的解析解可知：L 2/L 10.0070.010.050.080.121ΦΦ0.99870.99120.9560.930.912解：（第一题第一问）对于此问题，由于可以视为二维稳态导热问题，由二维稳态热传导1,1,,,1,1,22220m n m n m nm n m n m nt t t t t t xy+-+-+-+-+=∆∆基本方程：22220t tx y∂∂+=∂∂用数值法对该区域进行节点划分(如下图所示）：x 方向上一共划分M 个节点，y 方向上一共划分N 个节点。

可以将以上方程改为：如果我们取x y ∆=∆则有：1,1,1,1,,4m n m n m n m nm n t t t t t +-+-+++=考虑到节点位置的特殊性，我们在此将节点的种类进行如下划分，并给出节点的离散方程。

西安交通大学传热学大作业一、物理问题有一个用砖砌成的长方形截面的冷空气通道，其截面尺寸如下图1-1所示，假设在垂直于纸面方向上用冷空气及砖墙的温度变化很小，可以近似地予以忽略。

在下列两种情况下试计算：砖墙横截面上的温度分布；垂直于纸面方向的每米长度上通过砖墙的导热量。

第一种情况：内外壁分别均匀维持在0℃及30℃；第二种情况：内外壁均为第三类边界条件，且已知：K m W K m W h C t K m W h C t ∙=∙=︒=∙=︒=∞∞/53.0砖墙导热系数/20,10/4,30222211λ二、数学描写由对称的界面必是绝热面，可取左上方的四分之一墙角为研究对象，该问题为二维、稳态、无内热源的导热问题。

控制方程：02222=∂∂+∂∂y tx t边界条件：① 给出了边界上的温度，属于第一类边界条件：由对称性知边界1绝热： 0=w q ； 边界2、3为等温边界：t w2=0℃，t w3=30℃② 给出了边界上的边界上物体与周围流体间的表面传热系数h 及周围流体的温度t f ，属于第三类边界条件 由对称性知边界1绝热： 0=w q ；边界2为对流边界，)()(2f w w w t t h n tq -=∂∂-=λ； 边界3为对流边界，)()(3f w w w t t h n t q -=∂∂-=λ。

1-1图2-1图三、数学模型网格划分：将长方形截面离散成31×23个点，用有限个离散点的值的集合来代替整个截面上温度的分布，通过求解按傅里叶导热定律、牛顿冷却公式及热平衡法建立的代数方程，来获得整个长方形截面的温度分布，进而求出其通过壁面的冷量损失。

步长为0.1m ，记为△x=△y=0.1m 。

采用热平衡法，利用傅里叶导热定律和能量守恒定律，按照以导入元体（m,n ）方向的热流量为正，列写每个节点代表的元体的代数方程。

第一种情况：()()()()()()()()()()⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎨⎧+++==︒==︒==︒==︒==︒==︒==︒==︒=+-+-代表内部点,，点4126~6,1018,26~6,106,18~6,10,2618~6,10,631~1,3023,31~1,301,23~1,30,3123~1,30,11,1,,1,1,n m t t t t t n C m t n C m t n C n t n C n t n C m t n C m t n C n t n C n t n m n m n m n m n m 第二种情况对于外部角点（1,1）、（1,23）、（31,1）、（31,，23）有：()()02222,1,,22,,1,22=∆∆-+-∆+∆∆-+-∆±±x y t t t t x h y x t t t t yh n m n m n m f n m n m n m f λλ 得到：()()()()⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧++=++=++=++=22,3123,3023,312,311,301,3122,123,223,12,11,21,11865331400186533140018653314001865331400t t t t t t t t t t t t 同理可得：对于内部角点（6,6）（6,18）（26,6）（26,18） ，有()()()()()()()()⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧++++=++++=++++=++++=7,2618,2518,2719,2618,267,266,256,275,266,2618,717,619,618,518,67,66,75,66,56,671853359533592000718533595335920007185335953359200071853359533592000t t t t t t t t t t t t t t t t t t t t对于外部边界节点有()()()()⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=+++==+++==+++==+++=+-+-+-+-20~2,29253146537360020~2,29253146537360022~2,29253146537360022~229253146537360023,123,122,23,1,11,12,1,1,311,31,30311,11,1,21m t t t t m t t t t n t t t t n t t t t m m m m m m m m n n n n n n n n ，，， 对于内部边界节点有()()()()⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=+++==+++==+++==+++=+-+-+-+-25~7,6125330653153100025~7,6125330653153100017~7,6125330653153100017~7,6125330653153100018,118,119,18,6,16,15,6,1,261,26,27261,61,6,56n t t t t n t t t t n t t t t n t t t t m m m m m m m m n n n n n n n n ，， 对于内部节点有()1,1,,1,1,41+-+-+++=n m n m n m n m n m t t t t t传热问题的有限差分解法中主要采用迭代法。

《传热学》第四章大作业 ——二维稳态导热问题的数值解法
第一题：
如图所示，一个无限长矩形柱体，其横截面的边长分别为L 1和L 2，常物性。

该问题可视为二维稳态导热问题，边界条件如图中所示，其中L 1=0.6m ，L 2=0.4m ， T w1=60℃，T w2=20℃，λ=200W/(m·K)。

(1) 编写程序求解二维导热方程。

(2) 绘制x =L 1/2和y =L 2/2处的温度场，并与解析解进行比较。

已知矩形内
的温度场的解析解为()()()()
1211w2w1sh sh sin ,L L L y L x t t y x t πππ+=。

第二题
将第一题中y =L 2处的边界条件变为t =t w2，其他条件不变。

(1) 编写程序求解二维导热方程并计算从y =0处导入的热量Φ2。

(2) 当L 2<<L 1时，该二维导热问题可简化为一维导热问题。

在一维的近似下，试计算从y =0处导入的热量Φ1，并比较不同L 2/L 1下Φ2/Φ1的比值。

由该问。

二维稳态导热的数值计算2.1物理问题一矩形区域，其边长L=W=1，假设区域内无内热源，导热系数为常数，三个边温度为T1=0，一个边温度为T2=1，求该矩形区域内的温度分布。

2.2 数学描述 对上述问题的微分方程及其边界条件为：2222T T 0x y∂∂+=∂∂ x=0，T=T 1=0x=1，T=T 1=0y=0，T=T 1=0y=1，T=T 2=1 该问题的解析解：112121(1)sin n n n sh y T T n L x n T T n L sh W L ππππ∞=⎛⎫⋅ ⎪---⎛⎫⎝⎭=⋅ ⎪-⎛⎫⎝⎭⋅ ⎪⎝⎭∑ 2.3数值离散2.3.1区域离散区域离散x 方向总节点数为N ，y 方向总节点数为M ，区域内任一节点用I,j 表示。

2.3.2方程的离散 对于图中所有的内部节点方程可写为：2222,,0i j i jt t x y ⎛⎫⎛⎫∂∂+= ⎪ ⎪∂∂⎝⎭⎝⎭ 用I,j 节点的二阶中心差分代替上式中的二阶导数，得：+1,,-1,,+1,,-1222+2+0i j i j i ji j i j i j T T T T T T x y --+=上式整理成迭代形式：()()22,1,-1,,1,-12222+2()2()i j i j i j i j i j y x T T T T T x y x y ++=++++ (i=2,3……,N -1)，(j=2,3……,M -1)补充四个边界上的第一类边界条件得：1,1j T T = (j=1,2,3……,M),1N j T T = (j=1,2,3……,M),1i j T T = (i=1,2,3……,N),2i M T T (i=1,2,3……,N)#include<stdio.h>#include<math.h>#define N 10#define M 10main(){char s;int i,j,l;float cha,x,y;float t[N][M],a[N][M];/*打印出题目*/printf("\t\t\t 二维稳态导热问题\t\t");printf("\n\t\t\t\t\t\t----何鹏举\n");printf("\n 题目：补充材料练习题二\n");printf("\n 矩形区域，边长L=W=1，假设区域内无内热源，导热系数为常熟，三个边温度为T1=0，一个边温度为T2=1，求该矩形区域内的温度分布。

二维稳态热传导方程在热传导领域，二维稳态热传导方程是一个重要的理论基础，它描述了物体在特定条件下热传导过程的物理规律。

在拥有固定温度场的前提下，二维稳态热传导方程可以描述物体的温度分布。

它的基本的数学表达式如下：$$frac{partial^{2}T}{partialx^{2}}+frac{partial^{2}T}{partial y^{2}}=0$$其中，T代表物体的温度，x和y分别表示物体内温度场的横向和纵向方向坐标。

从它的表达式可以看出，这个方程在x和y方向上是等价的，表明温度场在物体不同位置上具有同样的温差。

二维稳态热传导方程经常用以计算物体表层温度分布，如受稳态热传导作用的墙壁、空气和水等物质的温度分布研究。

在室外热传导研究中，除了研究墙壁的温度分布外，还可以使用二维稳态热传导方程来研究空气的温度分布，以及如何影响空气中的风流和湿度等。

二维稳态热传导方程还可用于研究物体内部温度分布，如金属、砖块、石子等物质的温度场。

通过利用二维稳态热传导方程，可以计算出物体的温度场，并且可以更好地判断物体的热传导率、热容量和热导率等特性。

二维稳态热传导方程还可以用于计算液体的温度场，如汽水的温度场。

这里，液体温度的变化是由流动的温度场所决定的，可以用二维稳态热传导方程来模拟温度场变化，以及液体在不同温度场下的密度变化。

此外，二维稳态热传导方程还可以用于估算加热物体的热量传递速度，以及物体内部温度的变化趋势。

这里的热量传递速度是指，当热源传入物体时，物体内部温度所发生的变化，可以根据二维稳态热传导方程计算出来。

综上所述，二维稳态热传导方程对热传导理论和技术具有重要的指导意义。

它不仅可以用于计算墙壁、空气、水和金属等物质的温度分布，还可以用于估算物体内部温度变化和热量传递速度。

有效地利用它，可以使我们更好地把握热传导过程，从而改善热传导设计、分析和控制的能力。

总的来说，二维稳态热传导方程是热传导领域中不可缺少的理论基础，在设计和研究热传导过程中具有重要的指导意义。

传热学二维导热物体温度场的数值模拟作者：陈振兴学号：10037005学院(系)：化工学院专业：过程装备与控制工程班级：装备01指导教师：李增耀实验时间：2012-10二维导热物体温度场的数值模拟一、物理描述有一个用砖砌成的长方形截面的冷空气通道，其截面尺寸和示意图如图1-1所示，假设在垂直纸面方向上冷空气及砖墙的温度变化很小，可以近似地予以忽略。

在以下情况下试计算：（1）砖墙横截面上的温度分布；（2）垂直于纸面方向的每米长度上通过砖墙的导热量。

1、内外表面均为第三类边界条件，且已知：10,3011=︒=∞h C t .33 C m W ︒⋅2/ 93.3,1022=︒=∞h C t C m W ︒⋅2/ 砖墙的导热系数C m W ︒⋅=/3.50λ2、内外壁分布均匀地维持在0C ︒及30C ︒；11h t ，∞ 1w t 22h t ∞ 2w t图1-1二、数学描述该结构的导热问题可以作为二维问题处理，并且其截面如图1-1所示，由于对称性，仅研究其1/4部分即可。

其网络节点划分如图2-1；上述问题为二维矩形域内的稳态、无内热源、常物性的导热问题，对于这样的物理问题，我们知道，描写其的微分方程即控制方程，就是导热微分方程：02222=∂∂+∂∂ytx t 第三类边界条件：内外表面均为第三类边界条件，且已知：33.10,3011=︒=∞h C t C m W ︒⋅2/ 93.3,1022=︒=∞h C t C m W ︒⋅2/砖墙的导热系数C m W ︒⋅=/3.50λa f（m ，n ）c bx ∆=y ∆ x ∆n y ∆ e m d图2-1三：方程的离散如上图2-1所示，用一系列与坐标轴平行的网络线把求解区域划分成许多子区域，以网格线的交点作为需要确定温度值的空间位置，即节点，节点的位置已该点在两个方向上的标号m 、n 来表示。

每一个节点都可以看成是以它为中心的小区域的代表，如上（m ，n ）：对于（m ，n ）为内节点时：由级数展开法或热平衡法都可以得到，当x ∆=y ∆时：)t t t t (41t 1,1,,1,1,-+-++++=n m n m n m n m n m 对于（m ，n ）为边界节点时： 位于平直边界上的节点：)t t 2t (41t 1,,1,1,--+++=n m n m n m n m 外部角点：如图2-1中a 、b 、d 、e 、f 点，)t t (21t 1,,1,--+=n m n m n m 内部角点：如图2-1中c 点，)t t 2t 2t (61t 1,1,,1,1,-+-++++=n m n m n m n m n m 由已知条件有，当m=1或n=13时的节点的温度衡为1w t =30C ︒，当（m=6且n<9）和（n=8且6<m<17）时的节点的温度为2w t =10C ︒。

二维稳态导热实验报告【实验目的】通过进行二维稳态导热实验，探究各种材料的导热性能，并了解热传导的基本原理和计算方法。

【实验原理】热传导是指物质内部由高温区传递热量到低温区的过程。

在二维导热实验中，通过热传导实现稳态状态，即材料内部的温度分布不发生变化，形成热平衡。

根据傅里叶热传导定律，稳态条件下热流密度在导向方向的分量为常数：q_x=-k(t)*(∂T/∂x)其中，q_x为单位面积上的热流量，k(t)为温度函数，T为温度。

【实验器材】1.导热材料试样（不同材料）2.温度计3.导热实验仪器（包括将试样放在恒温环境中的装置）【实验步骤】1.准备不同材料的导热试样，并将它们放入导热实验仪器中。

2.打开实验仪器，设置恒定的环境温度，使试样与环境保持稳定的温度差。

3.通过温度计记录试样上不同位置的温度，同时记录环境的温度。

4.根据收集到的数据，计算出试样内部的热传导率。

5.将实验得到的数据整理成表格和图表，并进行分析和讨论。

【实验结果】通过实验，我们得到了不同材料的热传导率数据，并绘制了相应的图表。

根据实验数据分析，我们可以得出以下结论：不同材料的热传导性能存在较大差异，导热性能较好的材料具有较大的热传导率。

同时，我们也观察到不同位置的温度差异，并且温度差与距离的平方成正比。

【实验讨论】在实验中，我们使用了二维稳态导热实验方法，获得了不同材料的热传导率数据。

然而，还有一些因素可能会对实验结果产生影响。

首先，试样的几何形状和尺寸可能会影响热传导过程。

其次，实验环境的温度稳定性对实验结果也有较大影响。

在实验中，我们尽量保持环境温度稳定以减小误差。

另外，实验中使用的温度计也可能存在一定的误差，这也需要在分析数据时进行考虑。

【实验总结】通过二维稳态导热实验，我们了解了热传导的基本原理和计算方法，并且探究了不同材料的导热性能。

实验结果表明，材料的导热性能受其热传导率的影响，导热性能较好的材料具有较高的热传导率。

通过本次实验，我们对热传导有了更深入的理解，并且提高了实验操作的能力。

二维稳态导热实验报告实验目的：通过测量传热过程中实验样品的温度分布，了解二维稳态导热过程的基本情况，验证热传导方程推导的正确性。

实验原理：物体在稳态热传导过程中，温度对时间的偏导数为零。

因此，稳态导热过程可以用二维热传导方程表示：∂²T/∂x²+∂²T/∂y²+Q/(kA)=0其中，T为样品各点的温度，x和y为样品的两个坐标，Q为内部恒定热源产生的热量，k为样品的导热系数，A为样品的横截面积。

实验步骤：1.准备实验样品，样品形状为正方形，边长为L。

在样品的中心部位设置一个内部恒定热源，热量为Q。

2.将实验样品放在两个绝缘材料之间，并保持两侧绝缘材料的温度固定。

3.在实验样品的表面布置一系列热敏电阻温度传感器，用于测量不同位置的温度。

4.将温度传感器连接到数据采集系统，实时监测和记录不同位置的温度数据。

5.等待样品温度分布稳定后，记录各个位置的温度值。

6.根据测得的温度数据，绘制样品的温度分布图。

实验结果与分析：根据实验步骤得到的温度数据，我们可以分析样品的温度分布。

在内部恒定热源的周围，温度较高，并随着距离热源的增加而逐渐降低。

在远离热源的地方，温度趋于稳定，呈现出均匀分布的特点。

进一步分析实验结果，我们可以计算得到样品的导热系数k。

根据二维热传导方程，我们可以推导出：k=-Q/(A*∂T/∂x)其中，Q为内部恒定热源产生的热量，A为样品的横截面积，∂T/∂x为样品在x方向上的温度梯度。

通过实验数据的分析，我们可以计算出上述表达式中的各个参数，从而得到样品的导热系数。

实验误差分析：在实际实验中，由于室温的影响、传感器的精度等因素，会引入一定的误差。

1.在实验之前，将实验样品和温度传感器等设备等置于同一环境温度下，等待它们达到稳定状态。

2.在用于固定样品和温度传感器的绝缘材料上，保持温度稳定，以减小环境温度变化对实验结果的影响。

3.使用准确的温度传感器和数据采集设备，以提高数据的精确性。

传热学导热问题的数值计算实践报告姓名：学号：班级：完成日期： 2019年12月12日一、实践题目及要求例题4-5 二维肋片稳态导热问题的数值计算（1）自主编程，编程语言自定，最后提交源程序（2）提交电子报告（word格式），包括：（a）给出空间离散示意图（网格划分）（b）节点离散方程（c）图示温度等值线（可以利用origin或matlab）二、空间离散示意图三、节点离散方程由上图所示得各节点的节点离散方程节点1：T(m,n)=0.25×[T(m+1,n)+T(m-1,n)+T(m,n+1)T(m,n-1)]节点2：T(M,n)= 1/(4+2×Bi)×[T(M,n-1)+T(M,n+1)+2×T(M-1,n)] 节点3：T(m,N)=1/(4+2×Bi)×[T(m-1,N)+T(m+1,N)+2×T(m,N-1)] 节点4：T(M,N)= 1/(2+2×Bi)×[T(M-1,1)+T(M,2)]节点5：T(m,1)=0.25×[T(m-1,1)+T(m+1,1)+2×T(m,2)]节点6：T(M,1)= 1/(2+2×Bi)×[T(M,N-1)+T(M,2)]四、温度等值线1、等温线图，工况1（Bi=0.01）η= 0.96702、工况1，温度与y轴分布图3、等温线图，工况2（Bi=1）η=0.19104、温度与y轴分布图（工况2）图像分析：从四幅图的显示来看，结果是可信的。

要是网格划分过松，就会出现在肋板顶端的绝热边界上温度的分布存在问题，温度的最高值并不是在半肋板顶端边界n=1处，而是在n>1的不远处的离散点上，这与预期是相违背的，但是当网格划分到达一定的密度，就可以避免这个问题，虽然在图像上看不出来，但此问题还是存在的，不过由于网格足够小，可忽略。

五、Matlab编程源程序function exampleT0=input('T0=');Tf=input('Tf=');h=input('h=');k=input('k=');x=input('x=');H=input('H=');M=input('M=');st=H/(M-1);N=floor(x/st)+1;Bi=h*st/k;p=1;for m=1:(M)for n=1:(N)T(m,n)=0;endendfor n=1:NT(1,n)=T0-Tf;endwhile p==1;p=0;for m=1:M;n=1:N;c(m,n)=T(m,n);endfor m=2:Mfor n=1:Nif (m>=2&&m<M&&n>=2&&n<N)T(m,n)=0.25*(T(m+1,n)+T(m-1,n)+T(m,n+1)+T(m,n-1));elseif (m==M&&n>=2&&n<N)T(M,n)=1/(4+2*Bi)*(T(M,n-1)+T(M,n+1)+2*T(M-1,n));elseif (m==M&&n==N)T(M,N)=1/(2+2*Bi)*(T(M,N-1)+T(M-1,N));elseif (m>=2&&m<M&&n==N)T(m,N)=1/(4+2*Bi)*(T(m-1,N)+T(m+1,N)+2*T(m,N-1));elseif (m>=2&&m<M&&n==1)T(m,1)=0.25*(T(m-1,1)+T(m+1,1)+2*T(m,2));elseif (m==M&&n==1);T(M,1)= 1/(2+2*Bi)*(T(M-1,N)+T(M,2));endendendfor m=1:M;n=1:N;if abs(c(m,n)-T(m,n))>=1E-6;p=1;endendendT1=0;T2=0;for m=2:1:MT1=T1+T(m,N);endfor n=2:1:(N-1)T2=T2+T(M,n);endQ=(0.5*(T(1,N)+T(M,1))+T1+T2)/(((M-1)+(N-1))*80);T=rot90(T+20);disp(Q)disp(Bi)disp(N)disp(T)contour(T)end六、个人总结与心得体会在本次实践中，我取得了较大收获。