%Variáveis
k=10;
rho=100000;
cp=1.15;
Tfe=100; %temperatura de
fluido
Tfw=100; %temperatura de
fluido
L=1;
s=60; %tempo total;
he=5;
hw=5;
N=7;
tx=10;
%Passo de tempo e
critério de parada
dt=2;
PAR=(s/dt)+1;
%Malha Computacional
Uniforme
XL=L;
L1=N;
L2=L1-1;
L3=L1-2;
dx=XL/L3;
XU(2)=0;
for i=3:L1
XU(i)=XU(i-1)+dx;
end
X(1)=XU(2);
for i=2:L2
X(i)=0.5*(XU(i+1)+XU(i));
XCV(i)=XU(i+1)-XU(i);
end
X(L1)=XU(L1);
for i=2:L1
XDIF(i)=X(i)-X(i-1);
end
%Volumes de controle
internos
for t=2:tx
for i=3:(N-2)
T(i,t)=0;
%Coeficientes:
B(i,t)=k/(XDIF(i+1));
C(i,t)=k/(XDIF(i));
D(i,t)=(rho*cp*(dx/dt))*T
(i,t-1);
A(i,t)=rho*cp*(dx/dt)+B(i
,t)+C(i,t);
P(i)=B(i,t)/(A(i,t)C(i,t)*P(i-1));
Q(i)=(D(i,t)+C(i,t)*Q(i1))/(A(i,t)-C(i,t)*P(i1));
T(i,t)=P(i)*T(i+1,t)+Q(i)
;
end
end
%Fim da barra
for i=(N-1)
for t=2:tx
T(i,t)=0; %retirar qnd
colocar TDMA, apenas
teste
end
%Temperatura na parede T1
for i=1
for t=1:tx
T(i,t)=(((T(i+1)*k)/L)+hw
*Tfw)/(hw+k/L);
end
end
%Temperatura na parede TL
for i=N
for t=1:tx
T(i,t)=(((T(i1)*k)/L)+he*Tfe)/(he+k/L)
;
end
end
%Plotar Gráfico
plot (T(i,t))
title('Variação da
temperatura')
xlabel('t')
ylabel('i')