24小时热门版块排行榜

返回列表

当前主题已经存档。

hongsemenghuan

木虫 (著名写手)

应助: 0 (幼儿园)
金币: 1433.5
散金: 84
红花: 2
帖子: 1223
在线: 132小时
虫号: 253430
注册: 2006-05-21
性别: MM
专业: 光学

[交流] 请问关于非线性显示差分的问题

s=0；T0=60000；

初始时刻：
for i=1:numx+2
if (x(i)-x(1))^2<=r0^2
      T(i,1)=T0;
else
      T(i,1)=300;
end
end

  差分方程如下：
for t=1:numt
for i=2:numx+1
当温度T(i,t)>=427时
T(i,t+1)=T(i,t)+(dt*k/(dens*cp*dx^2))*(T(i+1,t)-2*T(i,t)+T(i-1,t))+(Qal*A*dt/
(dens*cp))*exp(-E./(RR*T(i,t))) (非线性的)
每执行一次上述情况的计算，s=s+1；每一点只进行一次这样的计算。

温度达不到427的点
T(i,t+1)=T(i,t)+(dt*k/(dens*cp*dx^2))*(T(i+1,t)-2*T(i,t)+T(i-1,t))  （线性的）

end

边界：
T(1,t)=T(2,t);T(numx+2,t)=300;
end

取dx=0.5*10^(-9)m;dt=1*10^(-14)s 满足稳定性条件k*dt/(dens*cp*dx^2)<=1/2

计算结果 s 不随初始温度T0增大而增大。  这回是什么原因造成的呢？对于非线性方程来
说，稳定性条件对吗？谢谢帮助！
Sample Text

回复此楼

» 猜你喜欢

1楼 2006-06-11 16:34:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 hongsemenghuan 的主题更新

返回列表