| 查看: 172 | 回复: 0 | |||
| 当前主题已经存档。 | |||
hongsemenghuan木虫 (著名写手)
|
[交流]
请问关于非线性显示差分的问题
|
||
|
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 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
