传热传质上机实习题(参考资料C语言).doc_第1页
传热传质上机实习题(参考资料C语言).doc_第2页
传热传质上机实习题(参考资料C语言).doc_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

练习题三：一维非稳态导热的数值计算 非稳态导热问题由于有时间变量，其数值计算出现了一些新的特点。在非稳态导热微分方程中，与时间因素相关的非稳态项是温度对时间的一阶导数，这给差分离散带来了新的特点。由于这个特点，可以采用不同的方法构造差分方程，从而得到几种不同的差分格式，即所谓的显式、隐式和半显式。我们仍从一个具体问题出发来研究非稳态导热问题的数值计算。3.1问题一块无限大平板（如图3所示），其一半厚度为L=0.1m，初始温度T0=1000，突然将其插入温度T=20的流体介质中。平板的导热系数=34.89W/（m），密度=7800kg/m3,比热c=712J/（kg），平板与介质的对流换热系数为h=233W/（m2），求平板内各点的温度分布。程序#include#include#define N 10#define K 11main()int i,j,l;float cha;float a,x,y,Fo,Bi;float tNK,bNK;/*打印出题目*/printf(ttt一维非稳态导热问题tt);printf(ntttttt-n);printf(n题目：练习题三n);y=1;/*y代表*/x=0.05/(N-1);a=34.89/(7800*712);Fo=(a*y)/(x*x);Bi=233*x/34.89;printf(n显示格式条件:);printf(n1、Fo=%3.1f0nn,1-2*Fo*Bi-2*Fo);/*时刻为零时，赋予初场温度*/for(i=0;iN;i+)ti0=1000;/*循环开始，每次计算一个时刻*/for(j=0;jK-1;j+) for(i=0;i0.001) for(i=0;iN-1;i+) if(i=0) tij+1=Fo*(ti+1j+ti+1j)+(1-2*Fo)*tij;/*当计算t0时，要用到t-1,其中t-1=t2的(对称分布)*/ else tij+1=Fo*(ti+1j+ti-1j)+(1-2*Fo)*tij; tN-1j+1=tN-2j*(1-2*Fo*Bi-2*Fo)+2*Fo*tN-1j+2*Fo*Bi*20;/*边界点温度用热平衡法推导出公式*/ cha=0; for(i=0;i=0;j-)for(i=0;i999.99)printf(%6.1f ,tij);elseprintf(%6.2f ,tij);l=l+1;if(l=N)pr
人人文库> 全部分类> 应用文书 > 技术指导

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论