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

zhanglei_1

热传导方程转化为差分方程后是:T(i,j,k,n+1)=T(i,j,k,n)+Pd(i,j,k,n)*t(n)/(p*Cp)...<BR>                                                                 V*[T(i-1,j,k,n)-T(i,j,k,n)+T(i+,j,k,n)-T(i,j,k,n)...<BR>                                                                       +T(i,j+1,k,n)-T(i,j,k,n)+T(i,j-1,k,n)-T(i,j,k,n)<BR>                                                                       +T(i,j,k+1,n)-T(i,j,k,n)+T(i,j,k-1,n)-T(i,j,k,n)<BR>                                                                     ]        %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

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

zhanglei_1

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

dundun628

呵呵，欢迎以后有问题经常讨论哦<BR><BR>关于你的这篇文章，最好发到CFD版，那里的人会更熟悉些：）（虽然大家都是每个板块转一转）