




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
《传热学》上机实践大作业二维导热物体温度场的数值模拟能动A02赵凯2010031134一、物理问题有一个用砖砌成的长方形截面的冷空气通道,其截面尺寸如下图所示,假设在垂直于纸面方向上冷空气及砖墙的温度变化很小,可以近似地予以忽略。在下列两种情况下试计算:砖墙横截面上的温度分布;垂直于纸面方向的每米长度上通过砖墙的导热量。第一种情况:内外壁分别均匀地维持在0℃及30℃;第二种情况:内外表面均为第三类边界条件,且已知:t=30℃,/i=10.35W/m2.^col 1t=10℃,/i=3.93W/m^Koo2 2石专墙导热系猴=0.35/加二、数学描写1、控制方程该问题为无内热源的二维稳态导热问题,因此控制方程为导热微分方程:d21 d21 人——十——=0dxx2 dyy22、边界条件该问题中,导热物体在x方向上,y方向上都是对称的,因此可以只取其中的四分之一部分作为研究对象,其他部分情况完全相同,如下图所示:对于上图所示各边界:边界1:由对称性可知:其为绝热边界,即边界1:由对称性可知:其为绝热边界,即q=0。w边界2:第一种情况:其为等温边界,满足第一类边界条件。即:t=0。C即:w第二种情况:其为对流边界,满足第三类边界条件。即:边界3:边界3:第一种情况:其为等温边界,满足第一类边界条件。即:第二种情况:其为对流边界,满足第三类边界条件。第二种情况:其为对流边界,满足第三类边界条件。,at即:q=—X(一)=h(t—t)
w anwiwf即:三、方程离散如下图所示,用一系列与坐标轴平行的间隔10cm的网格线将求解区域划分成子区域。111可将上图所示各节点分成内节点与边界点两类。分别利用热平衡法列各个节点的代数方程。第一种情况(等温边界):边界点:,1边界1(绝热边界):/1— (2t)+t-+t1J,m—2~5乂,【:八、、乂,【人m,1 4m,2 m_1,1 m+1,11_ 、…,t= (21 +1 +1 ),n=8~1116,n4 15,n16,n-1 16,n+1边界2(内等温边界):t=0,m=6,n=1~7;m=7~16,n=7m,n边界3(外等温边界):t=30,m=1,n―1~12;m=2~16,n=12m,n内节点:, 1,, , , ,、.t=(t +1 +1 +1 );m,n 4m+1,n m-1,n m,n+1 m,n-1m=2~5,n=2~11;m=6~15,n=8~11第二种情况(对流边界):边界点:t边界1(绝热边界):m,11(21t边界1(绝热边界):m,11(214m,2+1m—1,1+1 ),m=2~5m+1,1t16,n1(2t4 15,n+t16,n—1+1 ),n=8~1116,n+12t+t+t+2Bit边界2(内对流边界):t6n5,n 6,n+16,n—1A边界2(内对流边界):t6n2(Bi+2)A2t+t+t +2Bitt=m,8 m+1,7 m—1,7 A11,m=7~16m,7 2(Bi+2)A12t+t+t+2Bitt= 2,n 1,n+1 1,n—1 A22.n=1~11边界3(外对流边界):1,n 2(Bi+2) ,A22t +t+t +2Bitt=m,11 m+1,12 m—1,12 A22,m=2~16m,12 2(Bi+2)A22(t+t)+t+t+2Bitt— 5,7 6,8 7,7 6,6A11内角点: 6,7 2(Bi+3)A1t+t+2Bitt―2,12 1,11 A22外角点: 1,12 2(Bi+1)A内节点:(t +t+t+t);m+1,n m—1,n m,n+1 m,n—1m—2~5,n—2~11;m—6~15,n—8~11(BiA1,t—t—10;BiA281—30)四、编程求解第一种情况(等温边界):Fortran程序代码如下所示:Programdenwengimplicitnoneinteger::t1=0integer::t2=30integerm,nreal::t(16,12),ta(16,12),et(16,12)real::epslona=1realfainei,fainei1,fainei2,fainei3,fainei4,fainei5,fainei6,fainei7realfaiwai,faiwai1,faiwai2,faiwai3,faiwai4,faiwai5,faiwai6,faiwai7realpianchadon=1,7t(6,n)=t1enddodom=7,16t(m,7)=t1enddodon=1,12t(1,n)=t2enddodom=2,16t(m,12)=t2enddodom=2,5don=1,11t(m,n)=10enddoenddodom=6,16don=8,11t(m,n)=10enddoenddoopen(01,file='dengwen.dat')dowhile(epslona>0.00000001)dom=2,5ta(m,1)=0.25*(2*t(m,2)+t(m-1,1)+t(m+1,1))enddodom=2,5don=2,11ta(m,n)=0.25*(t(m+1,n)+t(m-1,n)+t(m,n+1)+t(m,n-1))enddoenddodom=6,15don=8,11ta(m,n)=0.25*(t(m+1,n)+t(m-1,n)+t(m,n+1)+t(m,n-1))enddoenddodon=8,11ta(16,n)=0.25*(2*t(15,n)+t(16,n-1)+t(16,n+1))enddodon=1,7ta(6,n)=t1enddodom=7,16ta(m,7)=t1enddodon=1,12ta(1,n)=t2enddodom=2,16ta(m,12)=t2enddodom=1,16don=1,12et(m,n)=abs(ta(m,n)-t(m,n))enddoenddoepslona=maxval(et(1:16,1:12))dom=1,16don=1,12t(m,n)=ta(m,n)enddoenddoenddofainei1=0.5*lanbuda*t(5,1)fainei3=lanbuda*t(5,8)fainei5=0.5*lanbuda*t(16,8)fainei2=0don=2,7fainei6=lanbuda*t(5,n)fainei2=fainei2+fainei6enddofainei4=0dom=6,15fainei7=lanbuda*t(m,8)fainei4=fainei4+fainei7enddofainei=4*(fainei1+fainei2+fainei3+fainei4+fainei5)faiwai1=0.5*lanbuda*(30-t(2,1))faiwai3=lanbuda*(30-t(2,11))faiwai5=0.5*lanbuda*(30-t(16,11))faiwai2=0don=2,10faiwai6=lanbuda*(30-t(2,n))faiwai2=faiwai2+faiwai6enddofaiwai4=0dom=3,15faiwai7=lanbuda*(30-t(m,11))faiwai4=faiwai4+faiwai7enddofaiwai=4*(faiwai1+faiwai2+faiwai3+faiwai4+faiwai5)print*,' m nt'dom=1,16don=1,12print*,m,n,t(m,n)write(01,*)m,n,t(m,n)enddoenddopiancha=abs(fainei-faiwai)/((fainei+faiwai)/2)print*,'内部热流量=',faineiprint*,'外部热流量=',faiwaiprint*,'热平衡偏差=',pianchaendprogramdenweng运行结果如图所示:'C:\ProgramFiles\Mi:crosDftVisualStudio\Comm□ SDEV98\M¥PROJECTS\Debug\zui3....0,产mnt1130.09S091230.03@031330.096091430.09S091620.036031630.096091730.0BS001a20.BHGBB1930.0960911030.0BQ0011130-BBQB811230.096Q3Z1Z4.1D2662224.121402324.1S3G32424.308452524.529982G24.S9G442725.456812826.222052927.ll.7tl21028.Q690621129.03^522123fl-0BGBH|*C:\ProgramFiles\MicrasoftVisualStudio\Ccirrmon\MSDEV9B\MVPROJECT£\Debug\3.eKeb40.00O000@(E+005B.00B0QBB(E+0Q60.00O00O8(E+007B.OOOBOOBE+OO8 6.0109469 12.8176710 18.0176311 24.0108712 30.OQ00SI1O.00OO0OBIE+0020.000O000E+003B.00fl0BOB(E+QQ40.00O00OS(E1035B.EHtJEBtJBtE+En60.000000S(E+0370.00O000BIE+008 6.0091819 12.01483IQ1S.014S111 24.00914雪窗a雪窗a部部平ss内are270,6282239.GGS10.1213418tocontinue第二种情况(对流边界):Fortran程序代码如下所示:programduiliuimplicitnoneinteger::t1=10integer::t2=30integerm,nreal::t(16,12),ta(16,12),et(16,12)real::epslona=1realbi1,bi2realfainei,fainei1,fainei2,fainei3,fainei4,fainei5,fainei6,fainei7realfaiwai,faiwai1,faiwai2,faiwai3,faiwai4,faiwai5,faiwai6,faiwai7realpianchabi1=h1*detax/lanbudabi2=h2*detax/lanbudadom=1,16don=1,12t(m,n)=10enddoenddoopen(01,file='crs.dat')dowhile(epslona>0.000000001)dom=2,5ta(m,1)=0.25*(2*t(m,2)+t(m-1,1)+t(m+1,1))enddodon=8,11ta(16,n)=0.25*(2*t(15,n)+t(16,n-1)+t(16,n+1))enddota(6,n)=(2*t(5,n)+t(6,n+1)+t(6,n-1)+2*bi1*t1)/(2*bi1+4)enddodom=7,15ta(m,7)=(2*t(m,8)+t(m+1,7)+t(m-1,7)+2*bi1*t1)/(2*bi1+4)enddodon=2,11ta(1,n)=(2*t(2,n)+t(1,n+1)+t(1,n-1)+2*bi2*t2)/(2*bi2+4)enddodom=2,15ta(m,12)=(2*t(m,11)+t(m+1,12)+t(m-1,12)+2*bi2*t2)/(2*bi2+4)enddodom=2,5don=2,11ta(m,n)=0.25*(t(m+1,n)+t(m-1,n)+t(m,n+1)+t(m,n-1))enddoenddodom=6,15don=8,11ta(m,n)=0.25*(t(m+1,n)+t(m-1,n)+t(m,n+1)+t(m,n-1))enddoenddodon=2,6ta(6,7)=(2*t(5,7)+2*t(6,8)+t(7,7)+t(6,6)+2*bi1*t1)/(2*bi1+6)ta(1,12)=(t(2,12)+t(1,11)+2*bi2*t2)/(2*bi2+2)ta(6,1)=(t(5,1)+t(6,2)+bi1*t1)/(bi1+2)ta(16,7)=(t(16,8)+t(15,7)+bi1*t1)/(bi1+2)ta(16,12)=(t(16,11)+t(15,12)+bi2*t2)/(bi2+2)ta(1,1)=(t(2,1)+t(1,2)+bi2*t2)/(bi2+2)dom=1,16don=1,12et(m,n)=abs(ta(m,n)-t(m,n))enddoenddoepslona=maxval(et(1:16,1:12))dom=1,16don=1,12t(m,n)=ta(m,n)enddoenddoenddofainei1=0.05*h1*(t(6,1)-10)fainei3=0.1*h1*(t(6,7)-10)fainei5=0.05*h1*(t(16,7)-10)fainei2=0don=2,6fainei6=0.1*h1*(t(6,n)-10)fainei2=fainei2+fainei6enddofainei4=0dom=7,15fainei7=0.05*h1*(t(m,8)-10)fainei4=fainei4+fainei7enddofainei=4*(fainei1+fainei2+fainei3+fainei4+fainei5)faiwai6=0.1*h2*(30-t(1,n))faiwai2=faiwai2+faiwai6enddofaiwai4=0dom=2,15faiwai7=0.1*h2*(30-t(m,12))faiwai4=faiwai4+faiwai7enddofaiwai=4*(faiwai1+faiwai2+faiwai3+faiwai4+faiwai5)don=1,12dom=1,16print*,m,n,t(m,n)write(01,*)m,n,t(m,n)enddoenddopiancha=abs(fainei-faiwai)/((fainei+faiwai)/2)print*,'内部热流量=',faineiprint*,'外部热流量=',faiwaiprint*,'热平衡偏差=',pianchaclose(01)endprogramduiliufaiwai1=0.05*h2*(30-t(1,1))faiwai3=0.1*h2*(30-t(1,12))faiwai5=0.05*h2*(30-t(16,12))faiwai2=0don=2,11运行结果如图所示:'■'C:\Pnognam rcrasoftVr&jalStudio\Comman\NSDEV9S\MYPROJECTS\Debug\zui4....mnt1128.586092125.817023123.0O2124120.118625117.15t206114.12212?1Q.0Q0BQ00E+0QS1Q-BBBBQQ8E+BQ916.0@00603£+06101S.0S06@00E+0S111G-BBBBGBBE+BB1216.SS0iBS0dE-|-0S131B.OEJEJ0O0OE+OEJ141151S.0S0IB603E+0S1610.OS00S09E+O01228.594052225.839933223.036424220.158075217.192046214.14-464?
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 汽车资讯自媒体行业跨境出海战略研究报告
- 2025年横包卷烟机组合作协议书
- 小学语文教材选用中的问题及解决措施
- 地下管线施工管理措施计划
- 2024年度四川省护师类之主管护师押题练习试卷A卷附答案
- 2025年锰酸锂合作协议书
- 三年级上册道德与法治教学评价计划
- 国际会计学在房地产行业的风险管理范文
- 珠宝活动策划方案10篇
- 初中八年地理试卷及答案
- 《论语》全文带拼音有注释(完整版)
- 慢性心力衰竭患者液体出入量智能化管理实践
- 2024年全国职业院校技能大赛(植物病虫害防治赛项)考试题库-下(多选、判断题)
- 二年级下册口算题1000题大全-
- 幼儿园班级幼儿图书目录清单(大中小班)
- 提高急诊科高危患者转运成功率
- MOOC 英语口语进阶-南京大学 中国大学慕课答案
- 干细胞招商话术
- 招标代理服务 投标方案(技术方案)
- 中外政治思想史-形成性测试四-国开(HB)-参考资料
- 2024年医院重症专科护士培训考试题库(含答案)
评论
0/150
提交评论