zhanglei_1 发表于 2006-5-4 10:13

求助:一个热传导方程的差分解法!谢谢

热传导方程转化为差分方程后是:T(i,j,k,n+1)=T(i,j,k,n)+Pd(i,j,k,n)*t(n)/(p*Cp)...<BR>                                                               V*      %t(n)为加热时间<BR>其中T为25*21*10*10矩阵,边界条件为:T(1,:,:,:)=T(25,:,:,:)=T(:,1,:,:)=T(:,21,:,:)=T0=23...<BR>                                                            %加热四方体区域四侧面为常温<BR>                                                         T(:,:,:,1)=T0=28      %加热时间为0温度为常温<BR>          我使用了直接迭代解上面的方程,但结果是错误的,求助,指点!<BR>          谢谢!

dundun628 发表于 2006-5-4 10:24

<P>如果我没有估计错误，你的时间步长取大了，造成方程不稳定。可以去看看有关方程稳定性的东西，然后调整一下时间步长试试。另外，空间离散的步长也很关键，你这种离散的方法，会导致离散后时间步长大约为一维时的1/3，所以方程很容易不稳定！</P>

zhanglei_1 发表于 2006-5-4 10:49

好,谢谢!如果再有什么问题,希望能再次请学长指点!

dundun628 发表于 2006-5-4 11:11

呵呵，欢迎以后有问题经常讨论哦<BR><BR>关于你的这篇文章，最好发到CFD版，那里的人会更熟悉些：）（虽然大家都是每个板块转一转）
页: [1]
查看完整版本: 求助:一个热传导方程的差分解法!谢谢