| 查看: 140 | 回复: 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 ] |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有118人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
2楼2010-05-10 13:49:10