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

..《传热学》上机大作业二维导热物体温度场的数值模拟学校：XX交通大学__张晓璐学号:10031133__能动A06问题〔4-23有一个用砖砌成的长方形截面的冷空气通道,形状和截面尺寸如下图所示,假设在垂直纸面方向冷空气和砖墙的温度变化很小,差别可以近似的予以忽略。在下列两种情况下计算：砖墙横截面上的温度分布；垂直于纸面方向上的每米长度上通过墙砖上的导热量。第一种情况：内外壁分别维持在10和30第二种情况：内外壁与流体发生对流传热,且有,,,,问题分析控制方程边界条件所研究物体关于横轴和纵轴对称,所以只研究四分之一即可,如下图：对上图所示各边界：边界1：由对称性可知：此边界绝热,。边界2：情况一：第一类边界条件情况二：第三类边界条件边界3：情况一：第一类边界条件情况二：第三类边界条件三：区域离散化及公式推导如下图所示,用一系列和坐标抽平行的相互间隔的网格线将所示区域离散化,每个交点可以看做节点,该节点的温度近似看做节点所在区域的平均温度。利用热平衡法列出各个节点温度的代数方程。第一种情况：内部角点：平直边界1：平直边界2：平直边界3：第二种情况：内部角点：平直边界1：平直边界2：平直边界3：内角点：外角点：编程计算各节点温度和冷量损失〔冷量推导在后面〔用fortran编程由以上区域离散化分析可以得到几十个方程,要求解这些方程无疑是非常繁琐的,所以采用迭代法,用计算机编程求解这些方程的解,就可以得到各点温度的数值。迭代法：在迭代法中首先要对计算的场作出假设〔设定初场,在迭代计算过程中不断予以改进,直到计算前的假定值和计算后的结果相差小于允许值为止,称为迭代计算已经收敛。这里采用高斯-赛德尔迭代法解决此问题。第一种情况,等温边界programdengwen01implicitnoneintegerm,n realt<16,12>,tt<16,12>,wc<16,12> real::wcmax=0.2 real::chuwen=20integer::t1=10integer::t2=30 realq1,q2,q3,q4,q5,tz,qq dom=7,16 t<m,7>=t1 enddo don=1,7 t<6,n>=t1 enddodon=1,12 t<1,n>=t2 enddo dom=2,16 t<m,12>=t2 enddo dom=2,5 don=1,11 t<m,n>=chuwen enddo enddo dom=6,16 don=8,11 t<m,n>=chuwen enddo enddo dowhile<wcmax>0.0000001>dom=2,16 tt<m,12>=t2 enddo don=1,12 tt<1,n>=t2 enddo dom=7,16 tt<m,7>=t1 enddo don=1,7 tt<6,n>=t1 enddo don=8,11 tt<16,n>=0.25*<2*t<16,n>+t<16,n-1>+t<16,n+1>> enddo don=8,11 dom=6,15 tt<m,n>=0.25*<t<m+1,n>+t<m-1,n>+t<m,n+1>+t<m,n-1>> enddo enddo don=2,11 dom=2,5 tt<m,n>=0.25*<t<m+1,n>+t<m-1,n>+t<m,n+1>+t<m,n-1>> enddo enddo dom=2,5 tt<m,1>=0.25*<t<m+1,1>+t<m-1,1>+2*t<m,2>> enddo dom=1,16 don=1,12 wc<m,n>=abs<t<m,n>-tt<m,n>> enddo enddo wcmax=maxval<wc<1:16,1:12>>dom=1,16 don=1,12 t<m,n>=tt<m,n> enddo enddo enddoopen<33,file='question1.txt'> dom=1,16 don=1,12 print*,"t<",m,n,">=",t<m,n> write<33,*>m,n,t<m,n> enddo enddoq1=0 q2=0don=2,6 q1=q1+<4*t<6,n>-2*t<5,n>-t<6,n+1>-t<6,n-1>>*0.53/2/0.1 enddo dom=7,15 q2=q2+<4*t<m,7>-2*t<m,8>-t<m+1,7>-t<m-1,7>>*0.53/2/0.1enddo q3=<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.1 print*,'qq=',qq write<33,*>'qq=',qq endprogramdengwen01运行结果：由于数据量较大所以各点温度见附录或者附件question1.txt第二种情况,对流边界programquestion02implicitnoneintegerm,n realt<16,12>,tt<16,12>,wc<16,12> real::wcmax=0.2 real::chuwen=10 real::drxs=0.53integer::tf1=10integer::tf2=30 realqq,tz dom=1,6 don=1,12 t<m,n>=chuwen enddo enddo dom=7,16 don=7,12 t<m,n>=chuwen enddo enddo open<01,file="question2.txt"> dowhile<wcmax>0.0001> dom=2,5 don=2,11 tt<m,n>=0.25*<t<m-1,n>+t<m+1,n>+t<m,n-1>+t<m,n+1>>!内节点 enddo enddodom=6,15 don=8,11 tt<m,n>=0.25*<t<m-1,n>+t<m+1,n>+t<m,n-1>+t<m,n+1>> enddo enddo dom=2,5 tt<m,1>=0.25*<2*t<m,2>+t<m-1,1>+t<m+1,1>>!边界1 enddo don=8,11 tt<16,n>=0.25*<2*t<15,n>+t<16,n-1>+t<16,n+1>> enddo don=2,6 tt<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> enddo dom=7,15 tt<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> enddo don=2,11 tt<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>!边界3 enddo dom=2,15 tt<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,16 don=1,12 wc<m,n>=abs<t<m,n>-tt<m,n>> enddo enddo wcmax=maxval<wc<1:16,1:12>>dom=1,16 don=1,12 t<m,n>=tt<m,n> enddo enddo enddo open<44,file="question02.txt"> dom=1,16 don=1,12 print*,"t<",m,n,">=",t<m,n> write<44,*>m,n,t<m,n> enddo enddo tz=0 don=2,7 tz=tz+t<6,n> enddo dom=7,15 tz=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=',qq endprogramquestion02运行结果：由于数据量较大所以各点温度见附录或者附件question02.txt每米长度上通过壁面的冷量的计算第一问：通过壁面的冷量既可以根据内表面计算也可以根据外表面计算,这里根据内表面计算,由于第一问是恒温边界条件,所以只需求出内边界各个节点所获得的热量,并加和即可第二问：通过壁面的冷量既可以根据内表面计算也可以根据外表面计算,这里根据内表面计算,因为这种情况下是对流换热所以只需计算出所有内表面单个单元的传热量,求和即可等温图：可以根据各点数据,用origin软件画出等温图问题一等温图：问题二等温图：附录：1.第一种情况各点温度mn温度1130.000001230.000001330.000001430.000001530.000001630.000001730.000001830.000001930.0000011030.0000011130.0000011230.000002126.068442226.080932326.122422426.205632526.353322626.597632726.971212827.481362928.0784121028.7127121129.3563421230.000003122.111913222.132863322.203113422.346793522.610013623.065983723.805843824.875843926.1195531027.4160731128.7126731230.000004118.113474218.135494318.210384418.368404518.673964619.250444720.310354822.096594924.1079141026.1193741128.0782541230.000005114.070995214.085255314.134535414.242485514.466985614.951495716.088525819.092265922.0961251024.8752451127.4809651230.000006110.000006210.000006310.000006410.000006510.000006610.000006710.000006816.087806920.3090661023.8045161126.9703761230.000007710.000007814.949897919.2478471023.0633671126.5960071230.000008710.000008814.463948918.6690381022.6050981126.3502881230.000009710.000009814.236829918.3592791022.3376891126.2000191230.0000010710.0000010814.1240810918.1935410102211210101230.0000011710.0000011814.0659511918.10447111022.10210111126.06204111230.0000012710.0000012814.0352312918.05631121022.05552121126.03395121230.0000013710.0000013814.0186813918.03000131022.02974131126.01826131230.0000014710.0000014814.0094814918.01526141022.01518141126.00934141230.0000015710.0000015814.0039715918.00640151022.00638151126.00393151230.0000016710.0000016814.0000016918.00000161022.00000161126.00000161230.00000第二种情况各点温度mn温度1126.132741226.149821326.204031426.304081526.463541626.697731727.016671827.415001927.8687411028.3518111128.8459811229.342292123.197072223.225572323.316322423.485102523.757262624.163212724.725512825.436482926.2457121027.1024321127.9740321228.845913120.204633220.239283320.350823420.562973520.917453621.472583722.285973823.359993924.5754431025.8384431127.1019731228.351454117.143144217.176364317.284994417.498774517.877304618.524014719.586134821.142434922.8579641024.5742141126.2442241227.867765114.015435214.038265314.114235414.270105514.569215615.160345716.392435818.765995921.1401251023.3565651125.4331951227.412946110.842216210.847176310.863756410.898326510.969296611.155906712.057526816.389286919.5803461022.2790661124.7193361227.012867711.154507815.153547918.5132071021.4603371124.1524971226.691178710.966758814.557358917.8588881020.8968581123.7393881226.452619710.894109814.250419917.4683791020.5291091123.4557891226.2861810710.8568110814.0819710