西安交通大学传热学大作业

1《传热学》上机大作业二维导热物体温度场的数值模拟学校：西安交通大学姓名：张晓璐学号:10031133班级：能动A062一．问题（4-23）有一个用砖砌成的长方形截面的冷空气通道，形状和截面尺寸如下图所示，假设在垂直纸面方向冷空气和砖墙的温度变化很小，差别可以近似的予以忽略。在下列两种情况下计算：砖墙横截面上的温度分布；垂直于纸面方向上的每米长度上通过墙砖上的导热量。第一种情况：内外壁分别维持在10C和30C第二种情况：内外壁与流体发生对流传热，且有Ctf101，)/(2021kmWh，Ctf302，)/(422kmWh，KmW/53.03二．问题分析1.控制方程02222ytxt2.边界条件所研究物体关于横轴和纵轴对称，所以只研究四分之一即可，如下图：对上图所示各边界：边界1：由对称性可知：此边界绝热，0wq。边界2：情况一：第一类边界条件Ctw10情况二：第三类边界条件4)()(11f边界3：情况一：第一类边界条件Ctw30情况二：第三类边界条件)()(22f三：区域离散化及公式推导如下图所示，用一系列和坐标抽平行的相互间隔cm10的网格线将所示区域离散化，每个交点可以看做节点，该节点的温度近似看做节点所在区域的平均温度。利用热平衡法列出各个节点温度的代数方程。第一种情况：内部角点：511~8,15~611~2,5~2)(411,1,,1,1,nmnmtttttnmnmnmnmnm平直边界1：11~8),2(415~2),2(411,161,16,15,161,11,12,1,nttttmttttnnnnmmmm平直边界2：7,16~7,107~1,6,10,,nmtnmtnmnm平直边界3：12,16~2,30;12~1,1,30,,nmtnmtnmnm第二种情况：内部角点：11~8,15~611~2,5~2)(411,1,,1,1,nmnmtttttnmnmnmnmnm平直边界1：11~8),2(415~2),2(411,161,16,15,161,11,12,1,nttttmttttnnnnmmmm平直边界2：67,16~7206~1,61.0,10,)2(222111111,1,,1,nmhnmmyxCtxhtxhttttffnmnmnmnm平直边界3：12,16~2411~1,11.0,30,)2(222222221,1,,1,nmhnmmyxCtxhtxhttttffnmnmnmnm内角点：20,10,)3(22)(2111116,67,78,67,57,6hCtxhtxhtttttff外角点：4,30,)1(222222211,112,212,1hCtxhtxhtttff4,30,2222222,11,21,1hCtxhtxhtttff4,30,22222212,1511,1612,16hCtxhtxhtttff20,10,2111112,61,51,6hCtxhtxhtttff20,10,2111118,167,157,16hCtxhtxhtttff7四．编程计算各节点温度和冷量损失（冷量推导在后面）（用fortran编程）由以上区域离散化分析可以得到几十个方程，要求解这些方程无疑是非常繁琐的，所以采用迭代法，用计算机编程求解这些方程的解，就可以得到各点温度的数值。迭代法：在迭代法中首先要对计算的场作出假设（设定初场），在迭代计算过程中不断予以改进，直到计算前的假定值和计算后的结果相差小于允许值为止，称为迭代计算已经收敛。这里采用高斯-赛德尔迭代法解决此问题。第一种情况，等温边界programdengwen01implicitnoneintegerm,nrealt(16,12),tt(16,12),wc(16,12)real::wcmax=0.2real::chuwen=20integer::t1=10integer::t2=30realq1,q2,q3,q4,q5,tz,qqdom=7,16t(m,7)=t1enddodon=1,7t(6,n)=t1enddodon=1,12t(1,n)=t2enddo8dom=2,16t(m,12)=t2enddodom=2,5don=1,11t(m,n)=chuwenenddoenddodom=6,16don=8,11t(m,n)=chuwenenddoenddodowhile(wcmax0.0000001)dom=2,16tt(m,12)=t2enddodon=1,12tt(1,n)=t2enddodom=7,16tt(m,7)=t1enddodon=1,7tt(6,n)=t1enddodon=8,11tt(16,n)=0.25*(2*t(16,n)+t(16,n-1)+t(16,n+1))enddodon=8,11dom=6,15tt(m,n)=0.25*(t(m+1,n)+t(m-1,n)+t(m,n+1)+t(m,n-1))enddoenddodon=2,11dom=2,5tt(m,n)=0.25*(t(m+1,n)+t(m-1,n)+t(m,n+1)+t(m,n-1))enddo9enddodom=2,5tt(m,1)=0.25*(t(m+1,1)+t(m-1,1)+2*t(m,2))enddodom=1,16don=1,12wc(m,n)=abs(t(m,n)-tt(m,n))enddoenddowcmax=maxval(wc(1:16,1:12))dom=1,16don=1,12t(m,n)=tt(m,n)enddoenddoenddoopen(33,file='question1.txt')dom=1,16don=1,12print*,t(,m,n,)=,t(m,n)write(33,*)m,n,t(m,n)enddoenddoq1=0q2=0don=2,6q1=q1+(4*t(6,n)-2*t(5,n)-t(6,n+1)-t(6,n-1))*0.53/2/0.1enddodom=7,15q2=q2+(4*t(m,7)-2*t(m,8)-t(m+1,7)-t(m-1,7))*0.53/2/0.1enddoq3=(6*t(6,7)-2*t(5,7)-2*t(6,8)-t(6,6)-t(7,7))*0.53/2/0.1!内角点q4=0.5*(4*t(6,1)-2*t(5,1)-t(6,2)-t(6,2))*0.53/2/0.1!6,1q5=0.5*(4*t(16,7)-2*t(16,8)-t(15,7)-t(15,7))*0.53/2/0.1qq=4*(q1+q2+q3+q4+q5)*0.1print*,'qq=',qqwrite(33,*)'qq=',qqendprogramdengwen01运行结果：10由于数据量较大所以各点温度见附录或者附件question1.txt第二种情况，对流边界programquestion02implicitnoneintegerm,nrealt(16,12),tt(16,12),wc(16,12)real::wcmax=0.2real::chuwen=10real::drxs=0.53integer::tf1=10integer::tf2=30realqq,tzdom=1,6don=1,12t(m,n)=chuwenenddoenddodom=7,16don=7,12t(m,n)=chuwenenddoenddoopen(01,file=question2.txt)dowhile(wcmax0.0001)11dom=2,5don=2,11tt(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n-1)+t(m,n+1))!内节点enddoenddodom=6,15don=8,11tt(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n-1)+t(m,n+1))enddoenddodom=2,5tt(m,1)=0.25*(2*t(m,2)+t(m-1,1)+t(m+1,1))!边界1enddodon=8,11tt(16,n)=0.25*(2*t(15,n)+t(16,n-1)+t(16,n+1))enddodon=2,6tt(6,n)=(2*t(5,n)+t(6,n+1)+t(6,n-1)+2*20*0.1*10/0.53)/(4+2*20*0.1/0.53)enddodom=7,15tt(m,7)=(2*t(m,8)+t(m+1,7)+t(m-1,7)+2*20*0.1*10/0.53)/(4+2*20*0.1/0.53)enddodon=2,11tt(1,n)=(2*t(2,n)+t(1,n+1)+t(1,n-1)+2*4*0.1*30/0.53)/(4+2*4*0.1/0.53)!边界3enddodom=2,15tt(m,12)=(2*t(m,11)+t(m+1,12)+t(m-1,12)+2*4*0.1*30/0.53)/(4+2*4*0.1/0.53)enddott(6,7)=(2*(t(5,7)+t(6,8))+t(7,7)+t(6,6)+2*20*0.1*10/0.53)/(6+2*20*0.1/0.53)!内角点tt(1,12)=(t(2,12)+t(1,11)+2*4*0.1*30/0.53)/(2+2*4*0.1/0.53)!外角点tt(1,1)=(t(2,1)+t(1,2)+4*0.1*30/0.53)/(2+4*0.1/0.53)tt(16,12)=(t(16,11)+t(15,12)+4*0.1*30/0.53)/(2+4*0.1/0.53)tt(6,1)=(t(5,1)+t(6,2)+20*0.1*10/0.53)/(2+20*0.1/0.53)tt(16,7)=(t(15,7)+t(16,8)+20*0.1*10/0.53)/(2+20*0.1/0.53)dom=1,16don=1,12wc(m,n)=abs(t(m,n)-tt(m,n))enddoenddowcmax=maxval(wc(1:16,1:12))dom=1,16don=1,12t(m,n)=tt(m,n)enddo12enddoenddoopen(44,file=question02.txt)dom=1,16don=1,12print*,t(,m,n,)=,t(m,n)write(44,*)m,n,t(m,n)enddoenddotz=0don=2,7tz=tz+t(6,n)enddodom=7,15tz=tz+t(m,7)enddoqq=4*(0.1*20*(tz-15*10)+0.05*20*(t(6,1)-10)+0.05*20*(t(16,7)-10))print*,'qq=',qqwrite(44,*)'qq=',qqendprogramquestion02运行结果：由于数据量较大所以各点温度见附录或者附件question02.txt五．每米长度上通过壁面的冷量的计算第一问：通过壁面的冷量既可以根据内表面计算也可以根据13外表面计算，这里根据内表面计算，由于第一问是恒温边界条件，所以只需求出内边界各个节点所获得的热量，并加和即可WAqQw1212.161*总第二问：通过壁面的冷量既可以根据内表面计算也可以根据外表面计算，这里根据内表面计算，因为这种情况下是对流换热所以只需计算出所有内表面单个单元的传热量，求和即可WtthQQifiii7063.124)(1.0111总等温图：可以