| 查看: 165 | 回复: 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 ] |
回复此楼» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有52人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有0人回复
2楼2010-05-10 13:49:10