How do i create a convection Boundary conditions

1 次查看（过去 30 天）

george mobol 2016-8-8

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/298836-how-do-i-create-a-convection-boundary-conditions

I am trying to solve a heat transfer problem with convection boundaries at the bottom, fixed temperature boundaries at the top, left and write side. these are the codes, but the answers differ from analytic solutions. N=5; % Number of columns M=5; % Number of rows T_amb = 290;% initial temp T = ones(5,5)*T_amb; T(:,1)=500;% boundary condition at the left T(:,5)=500;% boundary condition at the right T(1,:)=500;% boundary condition at the top % convection equation boundary at bottom,m=5, is T(m,n)= 2*T(m-1,n)+T(m,n-1)+T(m,n+1)+1500-9*T(m,n);

codes

N=5; % Number of columns
M=5; % Number of rows
T_amb = 290;% initial temp
T = ones(5,5)*T_amb;
Tol = 10^-6;
pmax = 10*10^6;
for p=1:pmax;
    T_old=T;
    for m=2:M-1;
        for n=2:N-1;
            if m==5;
                T(m,n)= 2*T(m-1,n)+T(m,n-1)+T(m,n+1)+1500-9*T(m,n);
                %T(5,n)= 2*T(4,n)+T(5,n-1)+T(5,n+1)+1500-9*T(5,n);% convection equation boundary condition at m=5
            end
          T(m,n) = 1/4*(T(m+1,n)+T(m-1,n)+T(m,n+1)+T(m,n-1)); 
          end
      end
      diff = max(max(abs(T - Tol)));
      if (diff < Tol)
          break
      end
  end