A.
MATLAB (explicit method)
x = 1:10; % Geometric Parameters
dx = abs( x(1) - x(2) );
nx = length(x);
y = 1:10;
dy = abs( y(1) - y(2) );
ny = length (y);
T_left = 25; % Boundary conditions
T_top = 0;
T_right = 25;
T_bottom = 50;
T = 10*ones(nx,ny); % Initial condition
T(:,1) = T_left; % Assigning IC & BC
T(1,:) = T_top;
T(:,nx) = T_right;
T(ny,:) = T_bottom;
T(1,1) = (T_left + T_bottom)/2;
T(1,ny) = (T_top + T_left)/2;
T(nx,ny) = (T_top + T_right)/2;
T(nx,1) = (T_right + T_bottom)/2;
T_initial = T; % Initializing
T_old = T;
alpha = 1.11; % Thermal diffusivity – Copper (cm^2/s)
dt = 0.01 ; % Time steps
k1 = alpha*(dt/(dx^2));
k2 = alpha*(dt/(dy^2));
conv = 1e-4; % Convergence criteria
error = 1;
iteration = 0;
for nt = 1: 1000 % Starting the time loop.
while error > conv % Until error is within convergence criteria
for j = 2:(ny-1) % Spatial loop
for i = 2:(nx-1)
node1 = (T_old(i-1,j) - 2*T_old(i,j) + T_old(i+1,j));
node2 = (T_old(i,j-1) - 2*T_old(i,j) + T_old(i,j+1));
T(i,j) = T_old(i,j) + (node1*k1) + (node2*k2);
end
end
error = max(max(abs(T_old - T))); % Calculating error
T_old = T;
iteration = iteration + 1;
end
T_initial = T;
end
[c,h] = contourf(x,y,T); % Plotting each iteration results
clabel(c,h);
colorbar
colormap(jet)
set(gca,'YDIR','reverse');
xlabel('X Axis');
ylabel('Y Axis');
title({['2D Transient Heat Conduction.'],['Number of iterations = ',num2str(iteration)]});
pause(0.01);
B. MATLAB (implicit method)
x = 1:10; % Geometric Parameters
dx = abs( x(1) - x(2) );
nx = length(x);
y = 1:10;
dy = abs( y(1) - y(2) );
ny = length (y);
T_left = 25; % Boundary conditions
T_top = 0;
T_right = 25;
T_bottom = 50;
T = 10*ones(nx,ny); % Initial condition
T(:,1) = T_left; % Assigning IC & BC
T(1,:) = T_top;
T(:,nx) = T_right;
T(ny,:) = T_bottom;
T(1,1) = (T_left + T_bottom)/2;
T(1,ny) = (T_top + T_left)/2;
T(nx,ny) = (T_top + T_right)/2;
T(nx,1) = (T_right + T_bottom)/2;
T_initial = T; % Initializing
T_old = T;
alpha = 1.17; % Thermal diffusivity – Copper (cm^2/s)
dt = 10; % Time steps
k1 = alpha*(dt/(dx^2));
k2 = alpha*(dt/(dy^2));
conv = 1e-4; % Convergence criteria
error = 1;
iteration = 0;
% Time loop
for nt = 1: 500 % Starting the time loop.
While error>conv % Until error is within convergence criteria
for j = 2(ny-1) % Spatial loop
for i = 2:(nx-1)
node1 = 1/(1+(2*k1)+(2*k2)); % Jacobi method
node2 = k1*node1;
node3 = k2*node1;
a = (T_old((i-1), j) + T_old((i+1), j));
b = (T_old(i, (j-1)) + T_old(i, (j+1)));
T(i,j) =(T_initial(i,j) *
node1)+(a*node2)+(b*node3);
end
end
error = max(max(abs(T_old - T))); % Calculating error
T_old = T;
iteration = iteration + 1;
end
T_initial = T;
end
figure(1) % Plotting
[c,h] = contourf(x,y,T);
clabel(c,h);
colorbar
colormap(jet)
set(gca,'YDIR','reverse');
xlabel('X Axis');
ylabel('Y Axis');
title({['2D Heat Conduction in Transient State - Implicit method using Jacobi
method.'],['Number of iterations = ',num2str(iteration)]});
error = 1; % Reinitialising the error value and the iteration
count and repeat the above steps for Gauss Seidel method.
Iteration = 0;
for nt = 1: 500
while error > conv
for j = 2:(ny-1)
for i = 2:(nx-1)
node1 = 1/(1+(2*k1)+(2*k2)); % Gauss Seidel method
node2 = k1*node1;
node3 = k2*node1;
a = (T((i-1), j) + T_old((i+1), j));
b = (T(i, (j-1)) + T_old(i, (j+1)));
T(i,j) =(T_initial(i,j) * node1)+(a*node2)+(b*node3);
end
end
error = max(max(abs(T_old - T))); % Calculating error
T_old = T;
iteration = iteration + 1;
end
T_initial = T;
end
figure(2)
[c,h] = contourf(x,y,T);
clabel(c,h);
colorbar
colormap(jet)
set(gca,'YDIR','reverse');
xlabel('X Axis');
ylabel('Y Axis');
title({['2D Transient Heat Conduction - Implicit method using Gauss Seidel
method.'],['Number of iterations = ',num2str(iteration)]});
error = 1; %reinitialising the error value and the iteration count and
repeat the above steps for SOR method.
iteration = 0;
for nt = 1: 500
while error > conv
for j = 2:(ny-1)
for i = 2:(nx-1)
node1 = 1/(1+(2*k1)+(2*k2));
node2 = k1*node1;
node3 = k2*node1;
a = (T((i-1), j) + T_old((i+1), j));
b = (T(i, (j-1)) + T_old(i, (j+1)));
T(i,j) =(T_initial(i,j) *node1)+(a*node2)+(b*node3);
T(i,j) = T_old(i,j)*(1-alpha)+alpha*T(i,j);
end
end
error = max(max(abs(T_old - T)));
T_old = T;
iteration = iteration + 1;
end
T_initial = T;
end
figure(3)
[c,h] = contourf(x,y,T);
clabel(c,h);
colorbar
colormap(jet)
set(gca,'YDIR','reverse');
xlabel('X Axis');
ylabel('Y Axis');
title({['2D Transient Heat Conduction - Implicit method using SOR method.'],['Number of
iterations = ',num2str(iteration)]});
