传热学数值计算大作业
更新时间:2023-10-24 09:35:01 阅读量: 综合文库 文档下载
传热学数值计算大作业
一选题
《传热学》第四版P179页例题 4-3
二相关数据及计算方法
1.厚2δ=0.06m的无限大平板受对称冷却,故按一半厚度作为模型进行计算
2. δ=0.03m,初始温度t0=100℃,流体温度t∞=0℃;
λ=40W/(m.K),h=1000W/(m
2
.K),Bi=h*△x/λ=0.25;
3.设定Fo=0.25和Fo=1两种情况通过C语言编程(源程序文件见附件)进行数值分析计算;
当Fo=0.25时,Fo<1/(2*(1+Bi)),理论上出现正确的计算结果; 当Fo=1时,Fo>1/(2*(1+Bi)),Fo>0.5,理论上温度分布出现振荡,与实际情况不符。
三网格划分
将无限大平面的一半划分为6个控制体,共7个节点。 △x=0.03/N=0.03/6=0.005,即空间步长为0.005m
四节点离散方程
绝热边界节点即i=1时,tij+1=2Fo△ti+1j+(1-2Fo△)tij 内部节点即0
tij+1=tij(1-2Fo△Bo△-2Fo△)+2Fo△ti-1j+2Fo△Bo△tf
五温度分布线图(origin)
六结果分析
1 空间步长,时间步长对温度分布的影响
空间步长和时间步长决定了Bo和Fo,两者越小计算结果越精确,但同时计算所需的时间就越长。 2 Fo数的大小对计算结果的影响
编程时对Fo=1及0.25的情况分别进行了计算,发现当Fo=1时,各点温度随时间发生振荡,某点的温度高反而会使下一时刻的温度变低,违反了热力学第二定律,因此在计算中对Fo的选取有限制。为了保证各项前的系数均为正值,对于内节点,Fo>0.5;对于对流边界节点,Fo<1/(2*(1+Bi))。 3 备注
在Fo=0.25时,为了反映较长时间后温度的分布,取T=600,并选取了其中部分时刻的温度输出进行画图。图像显示,随着时间的增长,各点温度趋向一致。
而当Fo=1时由于结果会出现振荡,只取T=6观察即可。
附录1 C语言源程序 当步长为0.005,Fo=1的程序
#include
int N=6;/*网格划分后计算区域共N个控制体*/ int T=6;/*共计算0~T时刻的温度变化*/
float t[7][7],c[7][7];/*将节点和时间对应构成二维数组*/ /*参数设定*/
inti,j,tf,r,h,min; floatBi,x,Fo; Fo=1; x=0.03/N;
tf=0;/*边界对流换热温度为0摄氏度*/ r=40; h=1000; Bi=h*x/r;
/*0时刻温度均为100摄氏度,即对各点设定初值*/ for(i=0;i<=N;i++) t[i][0]=100;
/*对T>0时刻各节点的温度进行计算*/ for(j=0;j { for(i=0;i<=N;i++) c[i][j]=t[i][j]; min=1; while(min>0.00001) { for(i=0;i<=N;i++) { if(i==0) t[i][j+1]=2*Fo*t[i+1][j]+(1-2*Fo)*t[i][j];/*绝热边界离散方程*/ else if ((i!=0)&&(i!=N)) t[i][j+1]=Fo*(t[i+1][j]+t[i-1][j])+(1-2*Fo)*t[i][j];/*中间所有节点的离散方程*/ else if (i==N) t[i][j+1]=t[i][j]*(1-2*Fo*Bi-2*Fo)+2*Fo*t[i-1][j]+2*Fo*Bi*tf;/*对流换热边界节点的离散方程*/ } min=0; for(i=0;i<=N;i++) min=min+abs(t[i][j]-c[i][j]); min=min/(N+1); } } /*输出各时刻相应节点的温度值*/ printf(\当步长为0.005,Fo=1,0~6时刻,间隔1个时间步长时计算的温度分布为:\\n\int n=0; for(j=T;j>=0;j--) for(i=0;i<=N;i++) { printf(\ n=n+1; if(n==(N+1))
正在阅读:
传热学数值计算大作业10-24
亚超声遥控开关的制作及电路09-10
2022江西省鹰潭市《教育心理学知识》教师教育与答案04-13
爱特中控使用手册 - 图文10-09
《数字媒体技术导论》复习03-15
2016年北湖区中小学生田径运动会成绩册 - 图文09-20
新能源知识读本在线考试(真题修改后)05-11
河南省书画家详细联系方式通讯录12-15
- 多层物业服务方案
- (审判实务)习惯法与少数民族地区民间纠纷解决问题(孙 潋)
- 人教版新课标六年级下册语文全册教案
- 词语打卡
- photoshop实习报告
- 钢结构设计原理综合测试2
- 2014年期末练习题
- 高中数学中的逆向思维解题方法探讨
- 名师原创 全国通用2014-2015学年高二寒假作业 政治(一)Word版
- 北航《建筑结构检测鉴定与加固》在线作业三
- XX县卫生监督所工程建设项目可行性研究报告
- 小学四年级观察作文经典评语
- 浅谈110KV变电站电气一次设计-程泉焱(1)
- 安全员考试题库
- 国家电网公司变电运维管理规定(试行)
- 义务教育课程标准稿征求意见提纲
- 教学秘书面试技巧
- 钢结构工程施工组织设计
- 水利工程概论论文
- 09届九年级数学第四次模拟试卷
- 传热学
- 数值
- 作业
- 计算
- 住房补贴有哪些政策?是一次性的吗
- 职场潜规则 - 第五章 与下属相处的艺术
- 关于成立项目建设条件评审委员会的通知
- 红帽Red Hat Linux 9光盘启动安装过程图解
- 某工程弃渣场变化引起的施工成本增加变更索赔报告(含临建设施变化)
- 汽包液位知识
- 2014年跃进汽车集团公司汽车在内蒙古自治区保有量分析年报 - 图文
- 广告牌匾拆除工程施工组织设计方案
- 2007年国家公务员考试各地市面试真题
- 玉溪市药店名录109家
- 2019年中国酒店管理行业发展回顾与发展趋势分析报告
- 05-口腔颌面部解剖学复习思考题及病例分析附加
- 美丽作文之美丽的民族园作文
- 第三节 战略性人力资源管理
- 在司法巡查动员部署大会上的表态发言
- 人肉搜索:是天堂还是地狱
- 3、INTO中断实验
- 会计制度设计-万福生科
- 学习筑梦中国心得体会4000字范文
- 连续与离散信号的可视化