传热学上机实验

传热学上机实验班级:学号：姓名：一：实验问题一个长方形截面的冷空气通道的尺寸如附图所示。假设在垂直于纸面的方向上冷空气及通道墙壁的温度变化很小，可以忽略。试用数值方法计算下列两种情况下通道壁面中的温度分布及每米长度上通过壁面的冷量损失：(1)内、外壁面分别维持在10℃及30℃；(2)内、外壁面与流体发生对流传热，且有λ=0.53W/（m·K），tf1=10°C、h1=20W/（m2·K）,tf2=30°C、h2=4W/（m2·K）。二：问题分析与求解本题采用数值解法，将长方形截面离散成31×23个点，用有限个离散点的值的集合来代替整个截面上温度的分布，通过求解按傅里叶导热定律、牛顿冷却公式及热平衡法建立的代数方程，来获得整个长方形截面的温度分布，进而求出其通过壁面的冷量损失。1.建立控制方程及定解条件对于第一问，其给出了边界上的温度，属于第一类边界条件。CCytxt301002222外壁温内壁温对于第二问，其给出了边界上的边界上物体与周围流体间的表面传热系数h及周围流体的温度tf，属于第三类边界条件。fwwtthntytxt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,nmtttttnCmtnCmtnCntnCntnCmtnCmtnCntnCntnmnmnmnmnm对于第二问有如下离散方程：对于外部角点（1,1）、（1,23）、（31,1）、（31,，23）有：02222,1,,22,,1,22xyttttxhyxttttyhnmnmnmfnmnmnmf得到：22,3123,3023,312,311,301,3122,123,223,12,11,21,11865331400186533140018653314001865331400tttttttttttt同理可得：对于内部角点（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,671853359533592000718533595335920007185335953359200071853359533592000tttttttttttttttttttt对于外部边界节点有20~2,29253146537360020~2,29253146537360022~2,29253146537360022~229253146537360023,123,122,23,1,11,12,1,1,311,31,30311,11,1,21mttttmttttnttttnttttmmmmmmmmnnnnnnnn，，，对于内部边界节点有25~7,6125330653153100025~7,6125330653153100017~7,6125330653153100017~7,6125330653153100018,118,119,18,6,16,15,6,1,261,26,27261,61,6,56nttttnttttnttttnttttmmmmmmmmnnnnnnnn，，对于内部节点有1,1,,1,1,41nmnmnmnmnmttttt4.设立温度场的迭代初值传热问题的有限差分解法中主要采用迭代法。采用此法求解时需要对被分解的温度场预先假定一个解，称为初场。对于本问题，本文采用内部流体温度作为初始温度t0=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°Cforx=6:26fory=6:18t(x,y)=10;endend%对内边界上的点给定温度10°Cfori=1:100000%多次迭代保证结果准确性forn=2:22%对内部节点进行迭代运算form=2:5t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endform=27:30t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endendform=2:30forn=2:5t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endforn=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;fori=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));form=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));endforn=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));endform=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));endforn=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));endform=2:30%内部节点温度分布forn=2:5t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endforn=19:22t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endendforn=2:22form=2:5t(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n+1)+t(m,n-1));endform=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行排列。详见EXCEL文档。问题一等温边界条件温度分布矩阵问题二对流边界条件温度分布矩阵问题一即等温边界条件下的温度分布图如下：5101520253024681012141618202212141618202224262827.625.824.122.320.618.81715.313.511.851015202530246810121416182022问题二即对流边界条件下的温度分布图如下：5101520253024681012141618202212141618202224262828.226.324.522.720.919.117.315.413.611.8510152025302468101214161820226.解的分析根据对角占优原则，迭代公式的选择应使每一个迭代变量的系数总是大于或等于该式中其他变量系数的绝对值。问题一中，1,1,,1,1,41nmnmnmnmnmttttt，满足对角占优原则。因此问题一的方程组是收敛的。问题二中，内部节点亦满足上式，对于外部角点，22,3123,3023,312,311,301,3122,123,223,12,11,21,11865331400186533140018653314001865331400tttttttttttt易得：193531865318653所以外部角点满足对角占优原则。同理易知内部角点、内部边界点、外部边界点均满足对角占优原则。因此问题二的方程组也是收敛的。综上所述：本文所得结果是合理的。7.通过壁面的冷量损失取四分之一的长方形截面进行研究，计算单位长度墙壁的导热量：等温边界条件：ytxxty11对流边界条件：txhtyh11按上式分别计算墙内外侧散热量21、。因此，整个长方形截面的单位长度墙壁总散热量为：4221对于第一问的等温边界条件按上式可得：W84.39对于第二问的对流边界条件按上式可得：W97.30三．实验总结通过本次实验，我加深了对数值模拟求解实际传热学问题的理解，对于工程中的传热