Programacion termodinamica 2
Enviado por andy camarena • 11 de Abril de 2017 • Práctica o problema • 771 Palabras (4 Páginas) • 189 Visitas
ALUMNO: DE LA CRUZ HUAROCC JHOEL
E.A.P: INGENIERIA QUIMICA DEL GAS NATURAL Y ENERGIA
PROGRAMACION SOBRE EQUILIBRIO EN LAS REACCIONES QUIMICAS
- PROGRAMAS PARA EL CALCULO DE LA CONSTANTE DE EQUILIBRIO:
- EJERCICIO 1:
clc,clear
disp('***********EQUILIBRIO EN LAS REACCIONES QUIMICAS************')
disp('SUPONGAMOS QUE TENEMOS LA SIGUIENTE REACCION')
disp('===============================================================')
disp(' C2H5OH + 0.5O2 = CH3CHO + H2O ')
disp('===============================================================')
disp('CAPACIDAD CALORIFICA DE CADA UNO DE LOS COMPONENTES DE LA REACCION, J/mol*K')
disp('ESCRIBIMOS LAS CONSTANTES DE REACCION')
a=1, b=0.5, c=1, d=1
disp('ESCRIBIMOS LOS DATOS DEL DATA BANK')
C1=[9.014 2.141*10^-1 -8.390*10^-5 1.373*10^-9]
C2=[2.811*10 -3.680*10^-6 1.746*10^-5 -1.065*10^-8]
C3=[7.716 1.823*10^-1 -1.007*10^-4 2.380*10^-8]
C4=[3.224*10 1.924*10^-3 1.055*10^-5 -3.596*10^-9]
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=[-2.350*10^5 -1.684*10^5]
COMP2=[0 0]
COMP3=[-1.644*10^5 -1.334*10^5]
COMP4=[-2.420*10^5 -2.288*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)
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 1]
- Ejercicio 2:
clc,clear
disp('***********EQUILIBRIO EN LAS REACCIONES QUIMICAS************')
disp('SUPONGAMOS QUE TENEMOS LA SIGUIENTE REACCION')
disp('===============================================================')
disp(' CO2 + H2 = CO + H2O ')
disp('===============================================================')
disp('CAPACIDAD CALORIFICA DE CADA UNO DE LOS COMPONENTES DE LA REACCION, J/mol*K')
disp('ESCRIBIMOS LAS CONSTANTES DE REACCION')
a=1, b=1, c=1, d=1
disp('ESCRIBIMOS LOS DATOS DEL DATA BANK')
C1=[1.980*10 7.344*10^-2 -5.602*10^-5 1.715*10^-8]
C2=[2.714*10 9.274*10^-3 -1.381*10^-5 7.645*10^-9]
C3=[3.087*10 -1.285*10^-2 2.789*10^-5 -1.272*10^-8]
C4=[3.224*10 1.924*10^-3 1.055*10^-5 -3.596*10^-4]
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=[-3.938*10^5 -3.946*10^5]
COMP2=[0 0]
COMP3=[-1.106*10^5 -1.374*10^5]
COMP4=[-2.420*10^5 -2.288*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
...