传热学上机实验班级:学号：姓名：一：实验问题一个长方形截面的冷空气通道的尺寸如附图所示。

假设在垂直于纸面的方向上冷空气及通道墙壁的温度变化很小，可以忽略。

试用数值方法计算下列两种情况下通道壁面中的温度分布及每米长度上通过壁面的冷量损失：(1)内、外壁面分别维持在10℃及30℃；(2)内、外壁面与流体发生对流传热，且有λ=0.53W/（m·K），t f1=10°C、h1=20W/（m2·K）, t f2=30°C、h2=4W/（m2·K）。

二：问题分析与求解本题采用数值解法，将长方形截面离散成31×23个点，用有限个离散点的值的集合来代替整个截面上温度的分布，通过求解按傅里叶导热定律、牛顿冷却公式及热平衡法建立的代数方程，来获得整个长方形截面的温度分布，进而求出其通过壁面的冷量损失。

1. 建立控制方程及定解条件对于第一问，其给出了边界上的温度，属于第一类边界条件。

⎪⎪⎪⎩⎪⎪⎪⎨⎧︒=︒==∂∂+∂∂CC y tx t 301002222外壁温内壁温 对于第二问，其给出了边界上的边界上物体与周围流体间的表面传热系数h 及周围流体的温度 t f ，属于第三类边界条件。

()⎪⎪⎩⎪⎪⎨⎧-=⎪⎭⎫ ⎝⎛∂∂-=∂∂+∂∂f w wt t h n t yt x t λ022222. 确定节点（区域离散化）用一系列与坐标轴平行的网格线把长方形截面划分为31×23个节点。

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

3. 建立节点物理量的代数方程对于第一问有如下离散方程：()()()()()()()()()()⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎨⎧+++==︒==︒==︒==︒==︒==︒==︒==︒=+-+-代表内部点,，点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 t4. 设立温度场的迭代初值传热问题的有限差分解法中主要采用迭代法。

采用此法求解时需要对被分解的温度场预先假定一个解，称为初场。

对于本问题，本文采用内部流体温度作为初始温度t 0=10°C 。

采用高斯—赛德尔迭代法进行迭代计算。

5.求解代数方程组源程序如下：问题一：m=31;n=23;t=zeros(m,n); %将长方形截面离散化为31×23个点p=10 %赋初温t(:,:)=p;t(:,1)=30;t(:,23)=30;t(1,:)=30;t(31,:)=30; %对外边界上的点给定温度30°Cfor x=6:26for y=6:18t(x,y)=10;endend%对内边界上的点给定温度10°Cfor i=1:100000 %多次迭代保证结果准确性for n=2:22 %对内部节点进行迭代运算for m=2:5t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endfor m=27:30t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endendfor m=2:30for n=2:5t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endfor n=19:22t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endendendt' %求得温度分布矩阵contour(t',1000); %画等温线图，等温线条数1000条。

C=contour(t',10); %作等温边界条件的等温线图，等温线条数10条clabel(C,'manual')问题二：m=31;n=23;t=zeros(m,n); %将长方形截面离散化为31×23个点p=10 %赋初温t(:,:)=p;for i=1:100000 %多次迭代运算t(1,1)=400/31+53/186*(t(2,1)+t(1,2)); %外角点温度计算公式t(1,23)=400/31+53/186*(t(2,23)+t(1,22));t(31,1)=400/31+53/186*(t(30,1)+t(31,2));t(31,23)=400/31+53/186*(t(30,23)+t(31,22));t(6,6)=2000/359+53/359*(t(5,6)+t(6,5))+53/718*(t(7,6)+t(6,7)); %内角点温度计算公式t(6,18)=2000/359+53/359*(t(5,18)+t(6,19))+53/718*(t(6,17)+t(7,18));t(26,6)=2000/359+53/359*(t(26,5)+t(27,6))+53/718*(t(25,6)+t(26,7));t(26,18)=2000/359+53/359*(t(26,19)+t(27,18))+53/718*(t(25,18)+t(26,17 ));for m=2:30 %外边界温度分布t(m,1)=600/73+53/146*t(m,2)+53/292*(t(m-1,1)+t(m+1,1));t(m,23)=600/73+53/146*t(m,22)+53/292*(t(m-1,23)+t(m+1,23));endfor n=2:22t(1,n)=600/73+53/146*t(2,n)+53/292*(t(1,n-1)+t(1,n+1));t(31,n)=600/73+53/146*t(30,n)+53/292*(t(31,n-1)+t(31,n+1));endfor m=7:25 %内边界温度分布t(m,6)=1000/153+53/306*t(m,5)+53/612*(t(m-1,6)+t(m+1,6));t(m,18)=1000/153+53/306*t(m,19)+53/612*(t(m-1,18)+t(m+1,18));endfor n=7:17t(6,n)=1000/153+53/306*t(5,n)+53/612*(t(6,n-1)+t(6,n+1));t(26,n)=1000/153+53/306*t(27,n)+53/612*(t(26,n-1)+t(26,n+1));endfor m=2:30 %内部节点温度分布for n=2:5t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endfor n=19:22t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endendfor n=2:22for m=2:5t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endfor m=27:30t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endendendt' %获得对流边界条件下的温度分布矩阵contour(t',1000); %作对流边界条件的等温线图，等温线条数1000条C=contour(t',10); %作对流边界条件的等温线图，等温线条数10条clabel(C,'manual')问题一（第一类边界条件）及问题二的温度分布矩阵如下：鉴于31列，23行的矩阵在WORD 中不好排列，故在这里，本文将温度矩阵选择90°，按23列，31行排列。