| 查看: 185 | 回复: 1 | |||
[交流]
【求助】麻烦大家帮忙看下，为何温度会升那么快。 已有1人参与
|
|
clc,clear all; length=12e-2; width=12e-2; height=18e-2; delta_x=0.4e-2; IE=round(length/delta_x); JE=round(width/delta_x); KE=round(height/delta_x); te=300*ones(IE+1,JE+1,KE+1); lamda=0.6;%water ro=1000;% c=4200; a=lamda/ro/c; delta_time=0.1*delta_x^2/a; power=646;%加热器功率，单位瓦 Qv=power/(11*delta_x)^3;%功率密度，即功率除以热源所占体积 qv=zeros(IE+1,JE+1,KE+1); qv(11:21,11:21,11:21)=Qv;% for n=1:100; %main domain for i=2:IE; for j=2:JE; for k=2:KE; te(i,j,k)=te(i,j,k)+(delta_time*lamda/ro/c)*((te(i-1,j,k)-2*te(i,j,k)+te(i+1,j,k))/delta_x^2+(te(i,j-1,k)-2*te(i,j,k)+te(i,j+1,k))/delta_x^2+(te(i,j,k-1)-2*te(i,j,k)+te(i,j,k+1))/delta_x^2)+0.001*delta_time*qv(i,j,k)/ro/c; end end end mesh(te(:,:,15)) drawnow; end [ Last edited by 木易寒 on 2010-5-13 at 09:44 ] |
回复此楼» 猜你喜欢
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有227人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有107人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有62人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
2楼2010-05-10 13:49:10