Infinitely Long Fin (Temperature distribution)

MATLAB Script

%% Incropera Ex 3.8

% Arif Hossain-06/05/2016

close all, clear all; clc;

Tb = 373;

T_inf = 298;

d = 5e-3;

P = pi*d;

Ac = 0.25*pi*d.^2;

h = 100;

k_cu = 401;

k_al = 177;

k_st = 13.4;

k = [k_cu k_al k_st];

m = sqrt((h*P)./(k*Ac));

x = [0:0.0001:0.299];

int = 0;

for i = 1:length(m)

    int = int+1;

    T(int,:) = (T_inf + (Tb-T_inf)*exp(-(m(i)*x)))-273;

end

%%

figure()

plot(x,T(1,:),'-r','linewidth',3)

hold on

plot(x,T(2,:),'-b','linewidth',3)

plot(x,T(3,:),'-k','linewidth',3)

set(gca,'fontsize',14)

xlabel('length [ m ]');

ylabel('Temperature [ ^0C ]')

legend('Pure Cu','2024 Al alloy', '316 stainless Steel')

grid on

xlim([0 0.3]);