How do i create a convection Boundary conditions
1 回表示 (過去 30 日間)
古いコメントを表示
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
0 件のコメント
回答 (0 件)
参考
カテゴリ
Help Center および File Exchange で Calculus についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!