Solving heat transfer equation using finite difference method

Question

Michael Omodara 2018-4-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/396655-solving-heat-transfer-equation-using-finite-difference-method

回答： Benjamin Großmann 2018-4-23

Please, I am trying to solve for the temperature distribution in a box but I am having issues with my code error: "Index exceeds matrix dimension". This is my code

M=15.5;             
L= 1.22;                                
dx= 0.122;                              
T(1)=35;                                
Tc = 25;                                
N= L/dx;                                
St= 50400;                              
dt= 1800;                               
rho= 720;                               
t= St/dt;                               
T(:,1)=T(1);                            
K= 0.0676+(0.00054*M);                  
c= 0.334 +(0.00977*M);                  
q= 2*10^7;
h= (q*(dx)^2)/K;
for j=1:t
    for i=1:N
T(i,j+1)=F(K,dt,dx,rho)*(T(i+1,j)+T(i-1,j)+h)+(1- (2*F(K,dt,dx,rho)))*T(i,j);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Benjamin Großmann 2018-4-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/396655-solving-heat-transfer-equation-using-finite-difference-method#answer_316706

You start with i=1 and one of your indices is T(i-1), so this is addressing the 0-element of T. In Matlab we start with index 1. This could be one problem but it is not possible to debug your code as it is since there are "end"s missing and the function or Matrix "F" is not given.