西安交大传热学大作业matlab,西安交通大学传热学大作业二维温度场热电比拟实验.doc...
二维导热物体温度场的数值模拟
一、物理问题
有一个用砖砌成的长方形截面的冷空气通道,其截面尺寸如下图1-1所示,假设在垂直于纸面方向上用冷空气及砖墙的温度变化很小,可以近似地予以忽略。在下列两种情况下试计算:
砖墙横截面上的温度分布;垂直于纸面方向的每米长度上通过砖墙的导热量。
第一种情况:内外壁分别均匀维持在0℃及30℃;
第二种情况:内外壁均为第三类边界条件,且已知:
二、数学描写
由对称的界面必是绝热面,可取左上方的四分之一墙角为研究对象,该问题为二维、稳态、无内热源的导热问题。
控制方程:
边界条件:
第一种情况:
由对称性知边界1绝热: ;
边界2为等温边界,满足第一类边界条件: ;
边界3为等温边界,满足第一类边界条件: 。
第一种情况:
由对称性知边界1绝热: ;
边界2为对流边界,满足第三类边界条件: ;
边界3为对流边界,满足第三类边界条件: 。
三、方程离散
用一系列与坐标轴平行的间隔0.1m的二维网格线将温度区域划分为若干子区域,如图1-3所示。
采用热平衡法,利用傅里叶导热定律和能量守恒定律,按照以导入元体(m,n)方向的热流量为正,列写每个节点代表的元体的代数方程,
第一种情况:
边界点:
边界1(绝热边界):
边界2(等温内边界):
边界3(等温外边界):
内节点:
第二种情况
边界点:
边界1(绝热边界):
边界2(内对流边界):
边界3(外对流边界):
内角点:
外角点:
内节点:
( ;)
编程思路及流程图
编程思路为设定两个二维数组t(i,j)、ta(i,j)分别表示本次迭代和上次迭代各节点的温度值,iter(实际编程时并未按照此名称来命名迭代步长)表示迭代进行的次数, 、分别表示外边界、内边界的散热量。开始时,给t(i,j)、ta(i,j) 赋相同的初始值,t(i,j)根据内节点和各边界节点的离散方程进行迭代,迭代后比较t(i,j)、ta(i,j)各个节点之间温度之差,若两个温度之差小给定的精度,则此时迭代完成,t(i,j)就是所求的温度场分布,若两温度之差不满足精度要求,则将t(i,j)的值赋给ta(i,j),t(i,j)继续迭代,直到二者各个点的温度之差满足精度要求,记下此时的迭代次数,并根据所得到的温度场分布计算内外边界上散热量以及偏差。
取定初始试探值t(i,j)=0ta
取定初始试探值
t(i,j)=0
ta(i,j)=t(i,j)
说明边界条件
输入已知参数
开始
ta(i,j)=t(
ta(i,j)=t(i,j)
计算新的内节点温度及新的边界点温度t(i,j)
iter=1
结束输出t(i,j)、iter平均导热量及偏差
结束
输出t(i,j)、iter
平均导热量及偏差
iter=iter+1
计算内外边界上散热量及其平均值、偏差
结果讨论
本次实验的实际边界条件为等温边界条件,因此以第一种情况为依据进行实验编程。程序(见附注)运行结果如图1-5所示。
使用MATLAB7.0软件将各节点温度作为原始数据做出四分之一墙角的温度场分布图
将实验测量结果和数值计算结果进行比较如表1-1
实验测量值
数值计算值
四分之一墙角外侧散热量
60.712
60.4286728
四分之一墙角内侧散热量
60.341
60.4287262
单位高度墙壁总散热量
242.106
241.715
相对平衡偏差
0.613
0.0000884
对比两种结果分析数值计算方法得到的散热量比实验测得的散热量相近,而前者的相对平衡偏差比较小。前者的误差主要来自于迭代次数和网格数的划分,但也可能存在网格独立解,可以进一步细化网格进行求解观察数值模拟结果;后者的误差主要来自于电阻的阻值误差和读数仪表的基本误差,提高标准电阻的准确性和读数仪表的精度等级可以减小实验误差。
西安交大传热学大作业matlab,西安交通大学传热学大作业二维温度场热电比拟实验.doc...相关推荐
- XJTU_ 西安交通大学2020大学计算机作业-第八周
XJTU_ 西安交通大学2020大学计算机作业-第八周 XJTU_ 西安交通大学2020大学计算机作业-第八周 注:所有题所有用例均已通过. 文章没上传到github,反正没人点星星:https:// ...
- XJTU_ 西安交通大学2020大学计算机作业-第六周
XJTU_ 西安交通大学2020大学计算机作业-第六周 XJTU_ 西安交通大学2020大学计算机作业-第六周 注:所有return 0;均可删去. 所有题均已通过编译. 第一题 请输入一个8位的十进 ...
- 温度 数值模拟 matlab,西安交通大学——温度场数值模拟matlab
<西安交通大学--温度场数值模拟matlab>由会员分享,可在线阅读,更多相关<西安交通大学--温度场数值模拟matlab(4页珍藏版)>请在读根文库上搜索. 1.温度场模拟m ...
- XJTU_ 西安交通大学2020大学计算机作业-第十二周
XJTU_ 西安交通大学2020大学计算机作业-第十二周 XJTU_ 西安交通大学2020大学计算机作业-第十二周 注:所有题所有用例均已通过. 俺又更新了,还是求点星星
- XJTU_ 西安交通大学2020大学计算机作业-第七周
XJTU_ 西安交通大学2020大学计算机作业-第七周 XJTU_ 西安交通大学2020大学计算机作业-第七周 注:所有题均已通过所有用例. 文章已上传到github:https://github.c ...
- XJTU_ 西安交通大学2020大学计算机作业-第九周
XJTU_ 西安交通大学2020大学计算机作业-第九周 XJTU_ 西安交通大学2020大学计算机作业-第九周 注:所有题所有用例均已通过. 俺又更新了,还是求点星星
- Matlab中矩阵编号方式以及一维二维三维数据间的相互转换
Matlab中矩阵编号方式以及一维二维三维数据间的相互转换 文章目录 Matlab中矩阵编号方式以及一维二维三维数据间的相互转换 一.问题的提出 二.一维数据转为二维和三维 三.三维数据转为一维和二维 ...
- 温度 数值模拟 matlab,西安交通大学——温度场数值模拟(matlab)
西安交通大学材料制备与成型实验--温度场数值模拟,matlab编程 温度场模拟matlab代码: clear,clc,clf L1=8;L2=8;N=9;M=9;% 边长为8cm的正方形划分为8*8的 ...
- 西安交通大学护理学计算机考试,西安交通大学护理专业怎么样
西安交通大学是国家教育部直属重点大学,现有教工6137人,其中专任教师3081人,现有(2018年底)全日制在校生40346人,西安交通大学护理专业为本科专业,面向全国符合条件的高考生招生. 专业招生 ...
最新文章
- Redis在windows实现将数据缓存起来定时更新读取
- wenbao与windows命令
- php oauth 扩展,PHP 安装 Oauth 扩展
- Apache Flink 零基础入门(十五)Flink DataStream编程(如何自定义DataSource)
- centons7网卡配置文件使用openvswitch bridge
- 服务器:浅谈 Nginx 性能调优,太实用了!
- js DOM——JS学习笔记2015-7-2(第73天)
- python中的yield使用详解
- 数据库设计优化经验谈(转载)
- 执行webpack发生了什么?
- 在部署服务器端的时候postgresql调试和远程连接操作需要考虑好远程连接问题。
- 程序员基础 MySQL数据库开发必备常识
- Neuralog.v2013.06 1CD (测井曲线智能矢量化软件)
- html+css+js实现小游戏flybird(完整版)
- 崂山道士“法术”已成真?还需VR+AI+5G三强紧联合
- ios适配iPhone和iPad
- 还爱着你心中曾经那朵红玫瑰吗?
- 手机 CPU 架构类型了解
- 苹果和Android用点应用对比是真的吗,苹果手机与安卓手机相比,真的是苹果手机完胜吗...
- 对淘宝双飞翼布局的的一点理解
热门文章
- [爬虫练习]爬取同程安全应急响应公开漏洞列表以及详情
- MySQL的InnoDB表如何设计主键索引-转自淘宝MySQL经典案例
- Oracle 学习笔记(四) Oracle表空间和数据文件
- Spring AOP 性能监控器
- 基于openstack的自动化测试平台设计头脑风暴
- 经典算法(1)- 求1到N的质数
- DSP及海思嵌入式板实现gb28181
- Glib2:error: Installed (but unpackaged) files found(七)
- AndroidLinux gdb用法
- Mac上搭建Go语言环境及第一个hello.go程序