Programacion termodinamica 2

...

disp('HALLAMOS VHT')

syms T

VHT=VHO+VC1*T+(VC2*T^2)/2+(VC3*T^3)/3+(VC4*T^4)/4

R=8.314

disp('~~~~~~~~~~~~~~~~~~~~~HALAMOS IR~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')

disp('UTILIZAMOS LA ECUACION DE GIBBS HELMHOLTZ')

LK=(VHT/R*T^2)

disp('INTEGRANDO LA ECUACION ANTERIOR TENEMOS:')

I=VHO/R*298-VC1*log(298)/R-VC2*298/2*R-(VC3*298^2)/6*R-(VC4*298^3)/12*R

lnk=-VHO/R*T+VC1*log(T)/R+VC2*T/2*R+(VC3*T^2)/6*R+(VC4*T^3)/12*R+I

disp('HALLAMOS VGT')

VGT=lnk*(-R*T)

disp('****************************************************************')

IR=-VG298/298+VHO/298-VC1*log(298)-VC2*298/2-(VC3*298^2)/6-(VC4*298^3)/12

disp('**************HALLAMOS LA CONSTANTE "K"*************************')

T=[298 300 350 400 450 500 550 600 650]

VHT=(VHO+VC1.*T+(VC2.*T.^2)/2+(VC3.*T.^3)/3+(VC4.*T.^4)/4)

VGT1=(VHO-VC1.*T.*log(T)-(VC2.*T.^2)./2-(VC3.*T.^3)./6-(VC4.*T.^4)./12-IR.*T)

K=log(-VGT1./R.*T)

disp('************REALIZAMOS LOS RESPECTIVOS GRAFICOS*****************')

subplot(2,2,1), plot(T,VHT,'*')

title('*****VARIACION DE LA ENTALPIA CON LA TEMPERATURA*****'), xlabel('TEMPERATURA, K'),ylabel('ENTALPIA, KJ/mol')

legend('DATOS HALLADOS')

grid

subplot(2,2,2),plot(T,VGT1,'*')

title('***VARIACION DE LA ENERGIA LIBRE CON LA TEMPERATURA***'), xlabel('TEMPERATURA, K'),ylabel('ENERGIA LIBRE, KJ/mol')

legend('DATOS HALLADOS')

grid

subplot(2,2,3),plot(T,K,'*')

title('*****VARIACION DE LA CONSTANTE DE EQUILIBRIO CON LA TEMPERATURA*****'), xlabel('TEMPERATURA, K'),ylabel('CONSTANTE DE EQUILIBRIO')

legend('DATOS HALLADOS')

grid

GRAFICO:

[pic 2]

C. ejercicio 3:

clc,clear

disp('***********EQUILIBRIO EN LAS REACCIONES QUIMICAS************')

disp('SUPONGAMOS QUE TENEMOS LA SIGUIENTE REACCION')

disp('===============================================================')

disp(' 2CO + 2H2 = CH4 + CO2 ')

disp('===============================================================')

disp('CAPACIDAD CALORIFICA DE CADA UNO DE LOS COMPONENTES DE LA REACCION, J/mol*K')

disp('ESCRIBIMOS LAS CONSTANTES DE REACCION')

a=2, b=2, c=1, d=1

disp('ESCRIBIMOS LOS DATOS DEL DATA BANK')

C1=[3.087*10 -1.285*10^-2 2.789*10^-5 -1.272*10^-8]

C2=[2.714*10 9.274*10^-3 -1.381*10^-5 7.645*10^-9]

C3=[1.925*10 5.213*10^-2 1.197*10^-5 -1.132*10^-8]

C4=[1.980*10 7.344*10^-2 -5.602*10^-5 1.715*10^-8]

disp('HALLAMOS LAS VARIACIONES')

VC1=c*C3(1,1)+d*C4(1,1)-a*C1(1,1)-b*C2(1,1)

VC2=c*C3(1,2)+d*C4(1,2)-a*C1(1,2)-b*C2(1,2)

VC3=c*C3(1,3)+d*C4(1,3)-a*C1(1,3)-b*C2(1,3)

VC4=c*C3(1,4)+d*C4(1,4)-a*C1(1,4)-b*C2(1,4)

disp('CALORES DE FORMACION Y ENERGIA LIBRE DE GIBBS')

COMP1=[-1.106*10^5 -1.374*10^5]

COMP2=[0 0]

COMP3=[-7.490*10^4 -5.087*10^4]

COMP4=[-3.938*10^5 -3.946*10^5]

disp('HALLAMOS SUS RESPECTIVAS VARIACIONES')

VH298=c*COMP3(1,1)+d*COMP4(1,1)-a*COMP1(1,1)-b*COMP2(1,1)

VG298=c*COMP3(1,2)+d*COMP4(1,2)-a*COMP1(1,2)-b*COMP2(1,2)

disp('HALLAMOS VARIACION DE HO')

VHO=VH298-VC1*298-(VC2*298^2)/2-(VC3*298^3)/3-(VC4*298^4)/4

disp('HALLAMOS VHT')

syms T

VHT=VHO+VC1*T+(VC2*T^2)/2+(VC3*T^3)/3+(VC4*T^4)/4

R=8.314

disp('~~~~~~~~~~~~~~~~~~~~~HALAMOS IR~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')

disp('UTILIZAMOS LA ECUACION DE GIBBS HELMHOLTZ')

LK=(VHT/R*T^2)

disp('INTEGRANDO LA ECUACION ANTERIOR TENEMOS:')

I=VHO/R*298-VC1*log(298)/R-VC2*298/2*R-(VC3*298^2)/6*R-(VC4*298^3)/12*R

lnk=-VHO/R*T+VC1*log(T)/R+VC2*T/2*R+(VC3*T^2)/6*R+(VC4*T^3)/12*R+I

disp('HALLAMOS VGT')

VGT=lnk*(-R*T)