% Introduccion a Octave - 2023 - FCAI-UNCuyo - MoCCAI
% 
% Trabajo Practico 4 - Ejercicio 6
%
% El fluido A entra al mezclador a T_{A1}=400C y 
% el B entra a una temperatura T_{B1} entre 500C y 800C. 
% Al mezclador entra el doble de fluido A que B. 
close all
clear all

%================================================
% 1) Cree funciones annimas de c_P para el fluido A y B. 
%    Grafique el comportamiento de los calores especificos 
%    entre 400C y 800C.
cp1=@(T) 3.381 + 1.804*1e-2*T - 4.300*1e-6*T.^2;
cp2=@(T) 8.592 + 1.290*1e-1*T - 4.078*1e-5*T.^2;

TT = 273+linspace(400,800,1000)';
TA1= 273+400;
TB1= 273+linspace(500,800,500)';


figure(1);clf
hold on
plot(TT-273,cp1(TT),'r-','linewidth',4)
plot(TT-273,cp2(TT),'b-','linewidth',4)
hold off
xlabel('T [C]','fontsize',24)
ylabel('c_P [cal/mol K]','fontsize',24)
h=legend('c_{P1}(T)','c_{P2}(T)');
set(h,'fontsize',24)
set(gca,'fontsize',24)

%================================================
% 2)A que temperatura T_2 salen los dos fluidos del mezclador
%   en el rango analizado? Grafique la respuesta.

% Diferencia de entalpias de cada fluido
aA=1; %reemplazar por el valor correcto
bA=1; %reemplazar por el valor correcto
cA=1; %reemplazar por el valor correcto
DHA=@(T1,T) aA*(T-T1) + bA*(T.^2-T1.^2) + cA*(T.^3-T1.^3);
            
aB=1; %reemplazar por el valor correcto
bB=1; %reemplazar por el valor correcto
cB=1; %reemplazar por el valor correcto
DHB=@(T1,T) aB*(T-T1) + bB*(T.^2-T1.^2) + cB*(T.^3-T1.^3);
Teq = zeros(500,1);

%Determinamos mediante fsolve el valor de la Temperatura de equilibrio
for i=1:500
  DH=@(T) 2.0*DHA(TA1,T) + DHB(TB1(i),T);
  T0 = 0.5*(TA1+TB1(i));
  %Teq(i) = fsolve(COMPLETAR);
end

%Graficamos resultados en el equilibrio
figure(2);clf
plot(TB1-273,Teq-273,'r-','linewidth',4)
xlabel('T_B[C]','fontsize',24)
ylabel('T_{eq}[C]','fontsize',24)
set(gca,'fontsize',24)

%================================================
% 3) Muestre grafiquemente el incremento de entalpia \Delta H
%    en kcal/mol ganado o perdido en funcion de T_{B1} por el
%    fluido A y B, respectivamente.

% Recordamos que: 2DHA = DHB
figure(3);clf
hold on
%plot(TB1-273, COMPLETAR,'k-','linewidth',10)
%plot(TB1-273, COMPLETAR,'r.-','linewidth',3)
hold off
xlabel('T_B [C]','fontsize',24)
ylabel('\Delta H[kcal/mol]','fontsize',24)
%h=legend('2\Delta H_A','-\Delta H_B');
%set(h,'fontsize',24)
set(gca,'fontsize',24)
