| 查看: 170 | 回复: 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 ] |
回复此楼» 猜你喜欢
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有210人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有41人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有1人回复
-
基金申请
已经有44人回复
2楼2010-05-10 13:49:10
50