版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选文档传热学二维导热物体温度场的数值模拟作者:陈振兴学 号:10037005学院(系):化工学院专业:过程装备与控制工程班级:装备01指导教师:李增耀实验时间:2012-10二维导热物体温度场的数值模拟 一、物理描述有一个用砖砌成的长方形截面的冷空气通道,其截面尺寸和示意图如图1-1所示,假设在垂直纸面方向上冷空气及砖墙的温度变化很小,可以近似地予以忽略。在以下情况下试计算:(1)砖墙横截面上的温度分布;(2)垂直于纸面方向的每米长度上通过砖墙的导热量。1、内外表面均为第三类边界条件,且已知:.33 砖墙的导热系数2、内外壁分布均匀地维持在0及30; 图1-1 二、数学描述该结构的导热问题可
2、以作为二维问题处理,并且其截面如图1-1所示,由于对称性,仅研究其1/4部分即可。其网络节点划分如图2-1;上述问题为二维矩形域内的稳态、无内热源、常物性的导热问题,对于这样的物理问题,我们知道,描写其的微分方程即控制方程,就是导热微分方程: 第三类边界条件:内外表面均为第三类边界条件,且已知: 砖墙的导热系数 a f(m,n)c b = n e m d 图2-1三:方程的离散如上图2-1所示,用一系列与坐标轴平行的网络线把求解区域划分成许多子区域,以网格线的交点作为需要确定温度值的空间位置,即节点,节点的位置已该点在两个方向上的标号m、n来表示。每一个节点都可以看成是以它为中心的小区域的代表
3、,如上(m,n):对于(m,n)为内节点时:由级数展开法或热平衡法都可以得到,当=时: 对于(m,n)为边界节点时:位于平直边界上的节点:外部角点:如图2-1中a、b、d、e、f点,内部角点:如图2-1中c点,由已知条件有,当m=1或n=13时的节点的温度衡为=30,当(m=6且n<9)和(n=8且6<m<17)时的节点的温度为=10。四:编程思路及流程图 开始 输入已知参数 结束否说明边界条件取定初始试探值TA(i, j)=0计算新的内节点和边界点温度T(i, j)T=TA比较所有节点|TA(i, j)- T(i, j)|<EPSILO ?计算内外边界散热量及热平衡偏
4、差输出内外边界散热量及热平衡偏差结束图3-1五、程序及运行结果第三类边界条件:1、实验程序(C语言):/ 1.cpp : 定义控制台应用程序的入口点。/#include "stdafx.h" #include<stdio.h>#include<math.h>int _tmain(int argc, _TCHAR* argv) int i,j,l; float dt=1.0,dx=0.1,dy=0.1; float t1317,a1317; float q1=0,q2=0,q=0,e; float lmd=0.53,h1=10.33,h2=3.93,t
5、1=30,t2=10,ep=1.0e-7; /*打印出题目*/ printf("ttt二维稳态导热问题tt"); printf("ntttttt-陈振兴装备n"); printf("n题目:二维导热物体温度场的电模拟实验n"); printf("n矩形区域,l1=2.2;l2=3;l3=2;l4=1.2,假设区域内无内热源,导热系数为常熟,内外表面均为第三类边界条件且已知t1=30;t2=10;h1=10.33;h2=3.93;LMD=0.53;求该矩形区域内的温度分布及垂直于纸面方向的单位长度上通过墙体的导热量。n&quo
6、t;); /*各节点上的温度值*/ for(j=0;j<17;j+) t0j=30.0; for(i=1;i<13;i+) ti0=30.0; for(i=7;i<13;i+) ti7=10.0; for(j=8;j<17;j+) t7j=10.0; for(i=1;i<7;i+)for(j=1;j<17;j+)tij=20,aij=0; for(i=7;i<13;i+)for(j=1;j<7;j+)tij=20,aij=0; while(dt>=ep) for(i=1;i<7;i+)for(j=1;j<17;j+)aij=ti
7、j; for(i=7;i<13;i+)for(j=1;j<7;j+)aij=tij; for(i=6;i<12;i+)for(j=2;j<6;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(i=2;i<6;i+) for(j=2;j<16;j+) tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(j=2;j<6;j+)t12j=(t12j-1+t12j+1+2*t11j)/4; for(i=2;i<6;i+)ti16=(ti-116+ti+116+2*ti15)/4; for(i=2;i&l
8、t;12;i+)ti1=(dx*h1*t1/lmd+(ti+11+ti-11)/2+ti2)/(2+dx*h1/lmd); for(j=2;j<16;j+)t1j=(dy*h1*t1/lmd+(t1j+1+t1j-1)/2+t2j)/(2+dy*h1/lmd); for(i=7;i<12;i+)ti6=(dx*h2*t2/lmd+(ti+16+ti-16)/2+ti5)/(2+dx*h2/lmd); for(j=7;j<16;j+)t6j=(dy*h2*t2/lmd+(t6j+1+t6j-1)/2+t5j)/(2+dy*h2/lmd); t11=(h1*dx*t1+lmd*(
9、t21+t12)/2)/(lmd+h1*dy);t116=(h1*dx*t1+lmd*(t115+t216)/2)/(lmd+h1*dy);t616=(h2*dx*t2+lmd*(t615+t516)/2)/(lmd+h2*dy);t121=(h1*dx*t1+lmd*(t122+t111)/2)/(lmd+h1*dy);t126=(h2*dx*t2+lmd*(t125+t116)/2)/(lmd+h2*dy);t66=(h2*dy*t2+lmd*(t56+t65+t76/2+t67/2)/(3*lmd+h2*dx); for(i=1;i<7;i+)for(j=1;j<17;j+)
10、dt=dt+abs(tij-aij); dt=dt/(6*16);for(i=7;i<13;i+)for(j=1;j<7;j+)dt=dt+abs(tij-aij); dt=dt/(6*6); printf("温度分布为:ttttttttt");l=0; for(i=1;i<7;i+) for(j=1;j<17;j+) printf("%3.1f ",tij); l=l+1; if(l=16) printf("n"); l=0; l=0; for(i=7;i<13;i+) for(j=1;j<7;j
11、+)printf("%3.1f ",tij); l=l+1; if(l=6) printf("n"); l=0; for(j=2;j<16;j+)q1=q1+(30-t1j)*h1*dx; for(i=2;i<12;i+)q1=q1+(30-ti1)*h1*dy;for(j=7;j<17;j+)q2=q2+(t6j-10)*h2*dx;for(i=7;i<12;i+)q2=q2+(ti6-10)*h2*dy;q1=q1+h1*(dx/2*(30-t116)+dy/2*(30-t121)+dx*(30-t11); q2=q2+h2*
12、(dx/2*(t616-10)+dy/2*(t126-10)+dx*(t77-10);q=(q1+q2)/2;e=abs(q2-q1)/q);printf("单位长度上1/4墙体的导热量为:%4.2fW,偏差为:%3.2f",q,e);getchar();getchar();return 0;运行结果图:图3-2实验算得导热量为97.62W,与数值模拟的偏差为(26.73*4W-97.62W)/(26.73*4)W*100%=8.7%2、数值模拟程序(matlab):z=29.9 29.7 29.5 29.3 29.1 29.0 28.8 28.7 28.7 28.6 28
13、.6 28.5 28.5 28.5 28.5 28.229.7 29.1 28.5 28.0 27.4 26.9 26.5 26.2 26.0 25.9 25.8 25.7 25.6 25.6 25.6 25.529.5 28.5 27.6 26.6 25.7 24.8 24.1 23.6 23.3 23.1 22.9 22.8 22.8 22.7 22.7 22.729.3 28.0 26.6 25.2 23.8 22.5 21.6 20.9 20.5 20.2 20.0 19.9 19.9 19.8 19.8 19.829.1 27.4 25.7 23.8 21.9 19.9 18.7 1
14、7.9 17.5 17.2 17.1 17.0 16.9 16.9 16.9 16.929.0 26.9 24.8 22.5 19.9 16.7 15.2 14.6 14.3 14.2 14.1 14.0 14.0 14.0 14.0 13.928.8 26.5 24.1 21.5 18.6 15.2 0 0 0 0 0 0 0 0 0 0;28.7 26.2 23.5 20.8 17.8 14.6 0 0 0 0 0 0 0 0 0 0;28.6 25.9 23.2 20.3 17.3 14.2 0 0 0 0 0 0 0 0 0 0;28.5 25.7 22.8 20.0 17.0 14.
15、0 0 0 0 0 0 0 0 0 0 0;28.3 25.4 22.6 19.7 16.7 13.6 0 0 0 0 0 0 0 0 0 0;27.2 25.1 22.5 19.5 16.3 12.5 0 0 0 0 0 0 0 0 0 0;v=18 22 26;xx,yy=meshgrid(y,x);surf(xx,yy,z);colorbar;xlabel('x');ylabel('y');zlabel('z');az = 0;el =- 90;view(az, el);shading interp;axis tight;figure,co
16、ntour(xx,yy,z,v);grid on数值模拟图:图3-3图3-4第一类边界条件1、 实验程序(C语言):#include "stdafx.h" #include<stdio.h>#include<math.h>int _tmain(int argc, _TCHAR* argv) int i,j,l; float dt=1.0,dx=0.1,dy=0.1; float t1216,a1216; float q1=0,q2=0,q=0,e; float lmd=0.53,t1=30,t2=0,ep=1.0e-7; /*打印出题目*/ prin
17、tf("ttt二维稳态导热问题tt"); printf("ntttttt-陈振兴装备n"); printf("n题目:二维导热物体温度场的电模拟实验n"); printf("n矩形区域,l1=2.2;l2=3;l3=2;l4=1.2,假设区域内无内热源,导热系数为常熟,内外表面均为第一类边界条件且已知t1=30;t2=0;LMD=0.53;求该矩形区域内的温度分布及垂直于纸面方向的单位长度上通过墙体的导热量。n"); /*各节点上的温度值*/ for(j=0;j<16;j+) t0j=30.0; for(i=
18、1;i<12;i+) ti0=30.0; for(i=5;i<12;i+) ti5=0.0; for(j=6;j<16;j+) t5j=0.0; for(i=1;i<5;i+)for(j=1;j<15;j+)tij=20,aij=0; for(i=5;i<11;i+)for(j=1;j<5;j+)tij=20,aij=0; while(dt>=ep) for(i=1;i<5;i+)for(j=1;j<15;j+)aij=tij; for(i=5;i<11;i+)for(j=1;j<5;j+)aij=tij; for(i=5
19、;i<11;i+)for(j=1;j<5;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(i=1;i<5;i+)for(j=1;j<15;j+) tij=(ti-1j+ti+1j+tij-1+tij+1)/4; for(j=1;j<5;j+)t11j=(t11j-1+t11j+1+2*t10j)/4; for(i=1;i<5;i+)ti15=(ti-115+ti+115+2*ti14)/4; for(i=1;i<5;i+)for(j=1;j<15;j+)dt=dt+abs(tij-aij); dt=dt/(4*14
20、);for(i=5;i<11;i+)for(j=1;j<5;j+)dt=dt+abs(tij-aij); dt=dt/(4*6); printf("温度分布为:ttttttttt");l=0; for(i=0;i<6;i+) for(j=0;j<16;j+) printf("%4.1f ",tij); l=l+1; if(l=16) printf("n"); l=0; l=0; for(i=6;i<12;i+) for(j=0;j<6;j+)printf("%4.1f ",tij
21、); l=l+1; if(l=6) printf("n"); l=0; for(j=1;j<15;j+)q1=q1+(30-t1j)*lmd; for(i=1;i<11;i+)q1=q1+(30-ti1)*lmd;for(j=6;j<15;j+)q2=q2+t4j*lmd;for(i=6;i<11;i+)q2=q2+ti4*lmd;q1=q1+lmd*(30-t115)/2+(30-t111)/2); q2=q2+lmd*(t415/2+t114/2+t44);q=(q1+q2)/2;e=abs(q2-q1)/q);printf("单位长度
22、上墙体的导热量为:%4.2fW,偏差为:%3.2f",4*q,e);getchar();getchar();return 0;运行结果图:2、数值模拟程序(matlab):z=30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0;30.0 29.0 28.1 27.1 26.2 25.5 24.9 24.5 24.3 24.2 24.1 24.1 24.0 24.0 24.0 24.0;30.0 28.1 26.1 24.2 22.3 20.7 19.6 18.9 18.5 18
23、.3 18.2 18.1 18.0 18.0 18.0 18.0;30.0 27.1 24.2 21.2 18.1 15.5 13.9 13.0 12.5 12.3 12.2 12.1 12.0 12.0 12.0 12.0;30.0 26.2 22.3 18.1 13.6 9.1 7.4 6.7 6.4 6.2 6.1 6.1 6.0 6.0 6.0 6.0;30.0 25.5 20.7 15.5 9.1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.0 24.9 19.6 13.9 7.4 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.0 24.5 18.9 13.0 6.7 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.0 24.3 18.5 12.6 6.4 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;30.0 24.2 18.3 12.3 6.2 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0;3
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 共同施工合同模板
- 高端小区广告合同模板
- 银行资金借贷合同模板
- 瓦工分包合同模板
- 读后感课件教学课件
- 体育用品合同模板
- 综合工程合同模板
- 喂鸡场租赁合同模板
- 金融宣传合同模板
- 码头土方装卸合同模板
- 工业产品质量安全风险监测实施规范
- 482023年广西职业院校技能大赛中职组《职业英语技能》赛项职场应用环节样题
- 高中化学教学学生高阶思维能力培养路径分析
- 原始记录书写培训课件
- 《数学家高斯》课件
- 2023年中国石化安庆石化公司校园招聘150人历年高频难易度、易错点模拟试题(共500题)附带答案详解
- 项目质量管理与保障措施
- 2023年10月自考试题06089劳动关系与劳动法
- 岗位风险排查管理制度
- 2024野生哺乳动物及栖息地调查技术规程
- 制图-正投影基础课件
评论
0/150
提交评论