TRABAJO 1ER CORTE DE CINTICA QUMICA.

INTEGRANTE: David Esteban Soler Cod: 6101433

Presentado a: Porf. Yovanny Morales Hernndez

Fundacin Universidad Amrica Facultad de Ingeniera Qumica Cintica Qumica. Bogot, Colombia 2013
1

Trabajo 1er corte de Cintica Qumica.

Enunciado:
El oxigeno y el hexafluoropropileno se hacen reaccionar en fase gaseosa para producir, Oxido de Hexafluoropropileno, sin embargo estos dos reactivos tambin reaccionan para producir fluoruro de carbonilo y fluoruro de trifluoroacetilo, los cuales son indeseados, el fluoruro de trifluoroacetilo a su vez reacciona con el oxigeno, produciendo mas fluoruro de carbonilo, con la primera partida del presupuesto de la investigacin en la universidad, se construy un reactor de laboratorio por lotes, con un sistema de calentamiento trmico y una capacidad de 905 mL, en una de las pruebas que se realizaron, se cargaron 0.051147 mol de Oxigeno y 0.052681 mol de hexafluoropropileno, mediante el sistema de calentamiento se mantuvo a 478 K el reactor, al arrancar la reaccin la presin era de 4.5 bar, el reactor se dej operar durante un tiempo, al cabo del cual la presin fue de 5.113244 bar y la cantidad de oxigeno y producto deseado fueron indetectables. a) Determine la composicin del reactor con el menor nmero de suposiciones o idealidades. b) Determine la conversin. c) Se invirti bien el dinero?

Desarrollo: I. ESTEQUIOMETRIA:

Para realizar este Ejercicio se debe hacer uso de varias herramientas que se irn nombrando en el transcurso del desarrollo, en cuanto lo primero se necesita hacer uso de la estequiometria, para esto se tienen que plantear las reacciones qumicas que se llevaron a cabo en el reactor de forma simblica.

Como la primera reaccin es reversible, entonces en realidad en el reactor suceden las siguientes 4 reacciones: A estas reacciones las llamaremos R1 a R4 respetivamente, con estas ecuaciones se realiza la matriz estequiomtrica para saber cules son las reacciones realmente independientes dentro del reactor, la matriz estequiomtrica es la siguiente y se construye segn los coeficientes estequiomtricos. Despus de construida se opera la matriz para llegar a una de ngulo inferior cero, eso se muestra en los siguientes pasos:
C3 F6 1 1 1 0 O2 C3 F6 O CF3COF COF2 1 1 0 0 2 1 1 0 0 2 1 0 1 1 1 0 1 2 2

C3 F6 1 0 1 0

O2 C3 F6 O CF3COF COF2 1 1 0 0 2 1 1 1 1 2 1 0 1 1 1 0 1 2 2

C3 F6 1 0 0 0

O2 1

C3 F6 O CF3COF COF2 2 1 1 1 0 0 1 1 1 0 1 1 2

1 2 1 2 1 2

C3 F6 1 0 0 0

O2 C3 F6 O CF3COF COF2 1 1 0 0 2 1 1 1 1 2 0 0 0 0 1 0 1 2 2

C3 F6 1 0 0 0

O2 1

C3 F6 O CF3COF COF2 2 1 1 0 0 0 1 1 0 0 1 2 0

1 2 1 2 0

C3 F6 1 0 0 0

O2 C3 F6 O CF3COF COF2 1 1 0 0 2 1 1 1 1 2 0 1 2 1 0 0 0 0

Cuando se obtiene la matriz de ngulo inferior cero se traduce para dejarla de nuevo en ecuaciones qumicas, tomando los nmeros en esta como coeficientes estequiomtricos, estas nuevas reacciones son las independientes y con las cuales se trabajar de aqu en adelante para determinar la composicin del reactor:

 Teniendo las reacciones independientes que suceden dentro del reactor R1, R4, R5 respectivamente, se pueden plantear las ecuaciones estequimtricas para el avance de reaccin adimensional, las ecuaciones generales para la estequiometria de un reactor de STR, son las siguientes: ( ) ( )[ ( ) ( ) ]

Ec. 1 Ecuacin para cada compuesto basado en la estequiometria de avance adimensional. ( ) ( )[ ]

Ec. 2 Ecuacin para moles totales basado en la estequiometria de avance adimensional. Sabiendo que Nj(t) es el nmero total de moles del compuesto j en el tiempo t, Ntot(t) es el nmero total de moles de toda la mezcla en el tiempo t, son los coeficientes estequiomtricos de cada una de las reacciones independientes dentro de la sumatoria, es el avance de reaccin m (m), divido el nmero total de moles en el estado de referencia y es el cambio de moles que causa la reaccin m. Planteando las ecuaciones estequimtricas para cada uno de los compuestos con la ecuacin 1 y para el nmero de moles totales con la ecuacin 2 quedan de la siguiente manera: ( ) ( ) ( )[ ( ) ( ) ( ( ) ( ) ( ) ( )] )]

( )[

( ) ( )

( )[

( ) ( ) ( )[ ( )[

( ) ( ) ( ) ( ) ( )

( ( ( ) ) ( (

( )] )]

)]

( )

Para hallar el nmero de moles totales son necesarios los deltas de reaccin, que se hallan de la siguiente forma: ( ) Ec. 3 Expresin para sacar los delta de reaccin. Para cada una de las reacciones se plantea un delta de reaccin:

Teniendo los deltas de las reacciones se puede establecer la ecuacin para el nmero de moles totales de la reaccin. ( ) ( )[ ]

Remplazando los datos que saben como el nmero de moles inciales y algunos nmeros de moles finales que los da el enunciado las ecuaciones de la estequiometria quedan como sigue: ( ) [ ]

( )

Ec. 4 Ecuacin que representa para el problema el nmero de moles de hexafluoropropileno en el tiempo t. [ [ ] ]

Ec. 5 Ecuacin que representa para el problema el nmero de moles de oxigeno en el tiempo t. [ ]

Ec. 6 Ecuacin que representa para el problema el nmero de moles de oxido de hexafluoropropileno en el tiempo t. ( ) [ ]

Ec. 7 Ecuacin que representa para el problema el nmero de moles de fluoruro de trifluoroacetilo en el tiempo t. ( ) [ ]

Ec. 8 Ecuacin que representa para el problema el nmero de moles de fluoruro de carbonilo en el tiempo t. ( ) [ ]

Ec. 9 Ecuacin que representa para el problema el nmero de moles totales en el tiempo t. En este punto del problema es que debemos eliminar las idealidades y suposiciones, el sistema queda con un grado de libertad ya que, hay 6 ecuaciones con 7 incgnitas, se necesita aadir una ecuacin mas para que el sistema se pueda resolver, esta ecuacin de mas ser una que satisfaga el nmero de moles totales, esto se logra hacer mediante una ecuacin de estado, en este caso es una mezcla de gases polares y apolares, segn el libro Thermodynamic Models for Industrial Applications: From Classical and Advanced Mixing Rules to Association Theories,
7

en el capitulo uno aparece en la figura 1.1 unas caractersticas para elegir el modelo termodinmico que sea capaz de modelar una mezcla, debido que esta mezcla es polar, no tiene electrolitos y esta a una presin menor de 10 bar, el libro aconseja usar cualquiera de 3 modelos termodinmicos, el de SchwartebtruberRenon, PR o RKS con la regla de mezclado de Wong-Sandler, Huron-Vidal versin 2 o PSRK. El modelo que escog fue el de PR con regla de mezclado de Wong-Sandler, ms que todo, porque es un modelo termodinmico que he usado en ocasiones anteriores y tengo el conocimiento de cmo hacerlo.

Figura 1. Thermodynamic Models for Industrial Applications: From Classical and Advanced Mixing Rules to Association Theories Capitulo 1

II.

MODELO TERMODINMICO PARA HALLAR NMERO DE MOLES.

Es necesario usar un modelo termodinmico capaz de modelar la mezcla de hexafluoropropileno, fluoruro de carbonilo, y fluoruro de trifluoroacetilo, ya que en el enunciado nos dan los valores de presin, temperatura y volumen del reactor, por medio de una EoS como la de Peng-Robinson, se debe ser capaz de hallar el nmero de moles, sin embargo como se utilizar una regla de mezclado como la WS, esta depende de las composiciones del sistema por lo cual se debe plantear un algoritmo iterativo para la solucin del problema, en base a la heurstica del estudiante se plante el siguiente algoritmo iterativo.
Calcular la composicin de la mezcla en condiciones de gas ideal.

Con estos valores calcular las constantes a y b de PR, usando la regla de mezclado de WS.

Calcular el volumen molar de la mezcla de gases y hallar con este la nueva composicin.

Si el error entre el valor del volumen molar obtenido y el anterior es de 1e-6, parar si no, continuar la iteracin.

Figura. 2 Algoritmo necesario para el problema Primero que todo hay que mostrar la ecuacin de PR para la mezcla, para entender que es lo que se necesita de esta ecuacin y como sacarlo.

Ec. 10 EoS de Peng-Robinson Donde P es la presin a la cual se quiere calcular el volumen molar, R es la constante universal de los gases, T es la temperatura a la cual se quiere calcular el volumen, V es el volumen molar de la mezcla, y las constantes a y b, son constantes de interaccin molecular, que se sacan a travs de la regla de mezclado de WongSandler.
9

La regla de mezclado de Wong-Sandler enuncia que:

Ec. 11 b de mezcla segn Wong-Sandler

Ec. 12 a de mezcla segn Wong-Sandler En donde los parmetros Q y D se calculan como sigue: Ec. 13 Parmetro para calcular el a y b de mezcla Q Ec. 14 Parmetro para calcular el a y b de mezcla D Resolviendo las sumatoria para Q: ( ) Ec. 15 Sumatoria del parmetro Q resuelto Donde Xm es la fraccin molar del compuesto m adems ( ) : ) ( ) (

Ec. 16 psi para calcular Q. Donde esto significa que b y a son las constantes de Peng-Robinson para el compuesto 1 puro adems:

) =(

*(

) +
10

Donde en estado puro a y b se calculan como sigue: [ (

)]

Ec. 17 Parmetro a para PR

Ec. 18 Parmetro b para PR Donde Tc es la temperatura critica, Pc es la presin crtica, TR es la temperatura reducida del sistema y es la funcin de factor acntrico de Ptizer que se calcula como sigue:

Ec. 19 funciones del factor acntrico Resolviendo la sumatoria para D:

Ec. 20 Sumatoria del parmetro D resuelto Donde X, a y b, son la fraccin molar y las constantes para cada compuesto, R es la constante de los gases, T es la temperatura del sistema, es la energa de gibbs

en exceso y ms adelante se explicar cmo se calcula y C es una constante y se

calcula como sigue:
( )

En cuanto la energa de Gibbs en exceso se calcula mediante la siguiente ecuacin:

Ec. 21 Gibbs en exceso en funcin del coeficiente de actividad

11

Donde , es el coeficiente de actividad de cada una de las especies dentro de la mezcla, es es necesario obtenerlo mediante UNIQUAC, un mtodo de contribucin de grupos ya que para algunos compuestos en la mezcla no se encuentran como tal constantes para otros modelos de actividad. El mtodo UNIFAC, se compone de 2 partes una combinatoria (C) debida a forma y tamaos de las molculas y otra residual (R) debida a las interacciones moleculares:

Ec. 22 Calculo coeficiente de actividad por UNIFAC Ec. 23 Parte combinatoria de UNIFAC Donde: ( ) ( )

Ec. 24 parmetro l, para ecuacin combinatoria de UNIFAC Ec. 25 parmetro , para ecuacin combinatoria de UNIFAC Ec. 26 parmetro r, para ecuacin combinatoria de UNIFAC Ec. 27 parmetro , para ecuacin combinatoria de UNIFAC Ec. 28 parmetro q, para ecuacin combinatoria de UNIFAC
12

Donde xi es la fraccin molar del componente i en la mezcla, es el nmero entero de grupos funcionales tipo k que hay en la molcula del componente i y Z=10. En Este paso combinatorial se nos presenta un impedimento enorme para proseguir con este modelo termodinmico, este impedimento es que no se tienen los parmetros R y Q, para toda la molcula de fluoruro de trifluoroacetilo, para visualizar mejor de lo que se est hablando en la siguiente imagen se visualiza la molcula:

Figura 3. Molcula de fluoruro de trifluoroacetilo. Como se puede ver se puede dividir la molcula en un grupo CF 3 y en otro COF, para el grupo CF3 se encuentran tanto parmetros R y Q de UNIFAC como aij, para la parte residual, pero la parte del COF, no se encuentran parmetros, por lo tanto se tendran que restar entre grupos, para hallar estos parmetros, lo cual traera consigo un error enorme a la hora de desarrollar el proceso de obtencin de energa de Gibbs en exceso y por lo cual tambin acarreara un error en el clculo de las constantes a y b para PR, por ende en ultimas, traera una desviacin a la hora de calcular el factor de compresibilidad y por ltimo el nmero de moles totales, en cada iteracin el error aumentara, por lo cual es necesario utilizar otro modelo que no necesite de UNIFAC para funcionar, en nuestras opciones se
13

encuentra una que es la que ms me llam la atencin se encuentra en el libro introduccin a la termodinmica de la ingeniera qumica, de autor Smith Van Ness, el tema se llama Propiedades de los fluidos a partir de las correlaciones del tipo Pitzer, en este trabajo habla de un mtodo para hallar propiedades termodinmicas con un mtodo de correlaciones como alternativa a el uso de una ecuacin cubica de estado, el trabajo de Ptizer para hallar el factor de compresibilidad, de una mezcla, a partir de este factor de compresibilidad se puede hallar el volumen molar o directamente el nmero de moles, el factor de compresibilidad esta dado mediante las siguientes ecuaciones:

Ec. 29 Factor acntrico para un gas real.

Ec. 29 Factor acntrico compuesto El y , esta dado por las presiones reducidas y las temperaturas reducidas, el libro enuncia que para las propiedades criticas para las mezclas se pueden aproximar mediante una regla lineal de mezclado dado que raras veces se conocen propiedades criticas para mezclas, se utilizan pseudoparametros crticos P pc y Tpc, en cuando al factor acntrico se calcula uno total para la mezcla, estos 3 parmetros necesarios para calcular el factor de compresibilidad se calculan como sigue: Ec. 30 Temperatura pseudocritica. Ec. 31 Presin pseudocritica.

14

Ec. 32 Factor acntrico en funcin de varios factores acntricos. Para las presiones y temperaturas pseudoreducidas, para el calcula de y :

Ec. 33 Temperatura pseudoreducida.

Ec. 34 presin pseudoreducida Con estas temperaturas reducidas segn el mtodo de Lee-Kesler, hay que hallar los y , mediante unas correlaciones con la ecuacin del virial, de esta ecuacin sale: ( )

Ec. 35 Ecuacin del virial.

Ec. 36 trmino de la ecuacin del virial simplificado. Transformando la ecuacin del virial queda:

Ec. 37 Ecuacin virial transformada Transformndola a la ecuacin 29:

Ec. 38 Trmino de la ecuacin 29

15

Ec. 39 Trmino de la ecuacin 29 Donde B0 y B1 son:

Ec. 40 Trmino

para la ecuacin 38.

Ec. 41 Trmino

para la ecuacin 38.

Para este procedimiento es necesario crear un nuevo mtodo iterativo para tener las composiciones reales aproximadas dentro del reactor:
Calcular la composicin de la mezcla en condiciones de gas ideal. Con estas composiciones calcular las temperatura y presion pseodureducida. Calcular el factor de compresibilidad. Calcular la nueva composicion.

Si el error entre el valor del volumen molar obtenido y el anterior es de 1e-6, parar si no, continuar la iteracin.

Figura 4. Nuevo algoritmo para obtener la composicin del reactor.

Como primera medida necesitamos para este procedimiento las temperaturas crticas, las presiones crticas y los factores acntricos de estas 3 molculas, para cada una de estas molculas se obtuvieron del simulador comercial Aspen
16

Properties versin 7.3, en la siguiente tabla se muestran los datos para cada una de las 3 molculas. Compuesto CF3COF COF2 C3F6 Tc (K) 297 368 Pc (bar) 57.6 29 0.3522 0.282943 0.20459

Tabla 1. Propiedades crticas y factor acntrico para compuestos del problema.

17

III.

MUESTRA DE CALCULO

Para este nuevo proceso iterativo es necesario primero calcular las moles totales con la ley de gases ideales para con estas calcular las temperaturas y presiones pseudoreducidas tanto como el factor acntrico de Ptizer para la mezcla. Como muestra de clculo se demostrar la primera iteracin. Para hallar el nmero de moles totales es necesario las condiciones a las cuales se encuentra el gas y otros datos que son los siguientes: Propiedad o dato Presin Temperatura Volumen Constante de los gases ideales Valor 5.113244 bar 478 K 0.905 L

Tabla 2. Tablas De propiedades del sistema.

Por la ley de los gases ideales que es como sigue el nmero de moles totales es:

Ec. 15 Sumatoria del parmetro Q resuelto Remplazando:

Despejando n, como nmero de moles totales: ( ) Con este el resultado se puede plantear un sistema 3x3 para remplazar los avances de reaccin adimensionales con las ecuaciones 5, 6, 9; dando as para cada Z una expresin como siguen:

Ya Teniendo Z1, se puede calcular Z2 con este como sigue.

18

Y teniendo ya estos dos se puede encontrar Z3, como sigue:

Los valores arrojados para Z1, Z2 y Z3 estn representados en la siguiente tabla y se calcularon en base a lo dicho en las 3 ecuaciones anteriores: Z1 0,371129809 Z2 0,614095791 Z3 0,242965983

Tabla 3. Avances de reaccin adimensional calculados con G.I.

Con estos avances y las ecuaciones 4, 7 y 8 se pueden hallar los nmeros de moles totales y los resultados estn reportados en la siguiente tabla: ( ) 0,013306994 ( ) 0,08898701 ( ) 0,014147332

Tabla 4. Cantidad de moles en el tiempo t para cada especie calculada con G.I.

Ya teniendo el nmero de moles de cada especie, se procede a calcular las fracciones molares de cada especie con la siguiente frmula:

Donde Ni, es el nmero de moles de la especie i en la mezcla y Ntot, es el nmero de moles totales en la mezcla segn esto, se sacan las fracciones molares y son las que se muestran en la siguiente tabla: ( ) 0,11428067

0,76422182

0,12149751

Tabla 5. Fraccin molar en el tiempo t para cada especie calculada con G.I.

Ya teniendo estos datos se pueden obtener las temperaturas y presiones pseudocriticas, el calor acntricos y las temperaturas y presiones pseudoreducidas. Con las ecuaciones 30,31 y 32 se obtienen las temperaturas y presiones pseudocriticas y el valor del factor acntrico, debido a su similitud se va ha hacer solo el clculo para la temperatura pseocudocritica, los datos necesarios para estos clculos estn en la tabla 5 y 1. La temperatura pseudocritica se calcula como sigue:
19

( Como resultado se obtiene:

Los clculos para la presin pseudocritica y el factor acntrico se hacen de manera anloga y dan los datos representados en la siguiente tabla: (K) (bar) 53,3732044 0,28133804

Tabla 6. Temperatura, presin pseudocritica y omega calculado con fracciones molares calculadas de G.I.

Con estos datos y la ecuacin 33 y 34 se calculan las temperaturas y presiones pseudoreducidas y se representa los datos en la siguiente tabla: (K) 1,47936656 (bar) 0,09580171

Tabla 7. Temperatura, presin pseudoreducidas con fracciones molares calculadas de G.I.

Teniendo las temperaturas reducidas se puede utilizar las tablas de Lee-Kesler de apndice E de Introduccin a la Termodinmica de la Ingeniera Qumica, o utilizar las correlaciones basadas en la ecuacin del virial para calcular Z, en este caso utilizaremos las correlaciones basadas en la ecuacin del Virial, en las que necesarios, B0 y B1 y se calculan con las ecuaciones 40 y 41, en este caso se mostrar el clculo de B0 ya que el de B1 es anlogo:

Los datos para B0 y B1 se muestran en la siguiente tabla:

-0,14252317

0,10579449

Tabla 8. Constantes B0 y B1 de la ecuacin del virial y dependientes de la temperatura pseudoreducida.

Con estos datos es posible ya calcular Z0 y Z1 necesarios para calcular el factor de compresibilidad total, estos factores de compresibilidad tambin se calculan de

20

forma anloga con las ecuaciones 38 y 39 por lo cual se calcular solo uno de ellos. Z0 se calcula como sigue con las ec

Los datos para Z0 y Z1 se muestran en la siguiente tabla: Z0 Z1 0,0068511

Tabla 9. Z0 y Z1 de la ecuacin del virial y dependientes de la temperatura pseudoreducida.

El factor de compresibilidad para la primera iteracin se calcula con la ecuacin 29 y es el siguiente: ( )

Teniendo el factor de compresibilidad se calcula el nmero de moles con la ecuacin 29 para un gas o una mezcla de gases reales y se despeja nmero de moles quedando: ( ) Teniendo esta ecuacin se remplazan los valores y da lo siguiente: ( )

Con este nmero de moles se vuelven a calcular los avances adimensionales de reaccin, y se vuelve a obtener al final de proceso un nmero de moles, los errores para cada iteracin se obtuvieron de la siguiente manera: ( ) Ec. 42 calculo de error porcentual. Los datos para todas las iteraciones hasta obtener un error muy pequeo se encuentran en las siguientes tablas:

21

Nmero de moles Iteracin totales al inicio del No. algoritmo

Avances de reaccin Z1

Nmero de Nmero de moles de Avances de Avances de moles de cada reaccin 2 reaccin 3 cada componente componente 2 Z2 Z3 CF3COF COF2 0,08898701 0,091556581

Nmero de moles de cada componente 3 C3F6 0,014147332 0,015003856

Gas ideal 1 2 3 4 5 6 7 8 9 10

0,116441336

0,371129809 0,61409579 0,24296598 0,013306994

0,11729785968423 0,362880361 0,62234524 0,25946488 0,010737423

0,11726373790557 0,363208998 0,6220166 0,2588076 0,010839788 0,091454216 0,014969734 0,11726507354273 0,363196135 0,62202947 0,25883333 0,010835781 0,091458222 0,01497107 0,11726502122510 0,363196638 0,62202896 0,25883232 0,010835938 0,091458065 0,014971018 0,11726502327435 0,363196619 0,62202898 0,25883236 0,010835932 0,091458072 0,01497102 0,11726502319409 0,363196619 0,62202898 0,25883236 0,010835932 0,091458071 0,01497102 0,11726502319723 0,363196619 0,62202898 0,25883236 0,010835932 0,091458071 0,01497102 0,117265023197106 0,363196619 0,62202898 0,25883236 0,010835932 0,091458071 0,01497102 0,117265023197111 0,363196619 0,62202898 0,25883236 0,010835932 0,091458071 0,01497102 0,117265023197110 0,363196619 0,62202898 0,25883236 0,010835932 0,091458071 0,01497102 Tabla 10. Contiene los avances de reaccin de las 3 reacciones y el nmero de moles de cada componente.

Fracciones Fracciones Fracciones molares de molares de Datos Iteracin molares de Datos cada cada pseudocriticos No. cada pseudocriticos componente componente 2 componente 2 3 CF3COF COF2 C3F6 Ppc 53,37320441 Tpc 323,1112655

Factor acntrico total

Gas ideal 0,114280668 0,764221821 0,12149751

0,281338042

1 0,091539801 0,780547753 0,127912446 53,33937217 320,0873732 0,279260448 2 0,092439387 0,77990193 0,127658683 53,34071051 320,2069927 0,279342634 3 0,092404165 0,779927217 0,127668619 53,34065811 320,2023091 0,279339416 4 0,092405544 0,779926226 0,127668229 53,34066016 320,2024926 0,279339542 5 0,09240549 0,779926265 0,127668245 53,34066008 320,2024854 0,279339537 6 0,092405492 0,779926264 0,127668244 53,34066008 320,2024857 0,279339537 7 0,092405492 0,779926264 0,127668244 53,34066008 320,2024857 0,279339537 8 0,092405492 0,779926264 0,127668244 53,34066008 320,2024857 0,279339537 9 0,092405492 0,779926264 0,127668244 53,34066008 320,2024857 0,279339537 10 0,092405492 0,779926264 0,127668244 53,34066008 320,2024857 0,279339537 Tabla 11. Contiene las fracciones molares, los datos pseudocriticos y el factor acntrico.

22

Datos Iteracin Datos pseudoreducidos No. pseudoreducidos 2 Ppr Gas ideal 0,095801705 Tpr 1,479366556

Constantes del virial B0 -0,142523175

Constantes del virial 2 B1 0,105794494

1 0,095862471 1,493342256 -0,139155712 0,107080274 2 0,095860065 1,492784389 -0,139288561 0,107030143 3 0,09586016 1,492806224 -0,139283359 0,107032107 4 0,095860156 1,492805369 -0,139283563 0,10703203 5 0,095860156 1,492805402 -0,139283555 0,107032033 6 0,095860156 1,492805401 -0,139283555 0,107032033 7 0,095860156 1,492805401 -0,139283555 0,107032033 8 0,095860156 1,492805401 -0,139283555 0,107032033 9 0,095860156 1,492805401 -0,139283555 0,107032033 10 0,095860156 1,492805401 -0,139283555 0,107032033 Tabla 12. Contiene los datos pseudoreducidos y las constantes del virial para el clculo de los Z.

Iteracin No. Gas ideal

Z 0,992697875

Nuevo nmero de moles 0,11729786

Error porcentual -

0,990770399

0,006851103

1 0,991067145 0,006873829 0,992986733 0,117263738 0,029098321 2 0,9910555 0,006873006 0,992975423 0,117265074 -0,00113899 3 0,991055956 0,006873039 0,992975866 0,117265021 4,46149E-05 4 0,991055938 0,006873037 0,992975849 0,117265023 -1,74754E-06 5 0,991055938 0,006873037 0,99297585 0,117265023 6,84503E-08 6 0,991055938 0,006873037 0,99297585 0,117265023 -2,68115E-09 7 0,991055938 0,006873037 0,99297585 0,117265023 1,05008E-10 8 0,991055938 0,006873037 0,99297585 0,117265023 -4,10659E-12 9 0,991055938 0,006873037 0,99297585 0,117265023 1,77518E-13 10 0,991055938 0,006873037 0,99297585 0,117265023 0 Tabla 13. Contiene los Z para el caculo del factor de compresibilidad, el nuevo nmero de moles y el error porcentual.

23

Resultados, conclusiones y aclaraciones

Aclaraciones
El modelo termodinmico usado para describir el problema, no es el ms idneo para el sistema como tal, el libro Smith van Ness, dice que las correlaciones de LeeKesler, sirven solo para mezclas de molculas apolares o poco polares, aunque las molculas con las que trabajamos no son extremadamente polares como el agua, esto puede causar errores en el clculo de las moles totales y en obviamente en la composicin del reactor al final, sin embargo, se abord este modelo, pues el modelo idneo para el desarrollo del problema que era el de PR con regla de mezclado de WS, tena el impedimento de que el coeficiente de actividad no se poda calcular con exactitud, y al repetir varias veces este mismo error podan haber desviaciones ms grandes que la tuvo hacerlo por el mtodo de Lee-Kesler. El mtodo de Lee-Kesler est restringido a bajas presiones, aunque la presin que tenemos no es alta, tipo de 10 a 50 bar, esta presin de 5 bar, crea desviaciones y errores con la prediccin de la realidad. Los clculos de las temperatura y presin critica del fluoruro de trifluoroacetilo, tienen un problema y es que no se encuentran datos experimentales para estos, por tanto en el simulador comercial Aspen Plus, tampoco estaban estos datos, el simulador los calcul por medio de contribucin de Joback y el mtodo de Ganni para obtener el factor acntrico, el clculo del factor acntrico se sabe que tiene un porcentaje bajo de error, porque para todas las molculas se hace de la misma manera, sin embargo el mtodo de Joback para el clculo de la temperatura y la presin critica, puede tener errores y diferir, en dado el caso con datos experimentales.

24

Resultados y conclusiones
Los resultados finales de la composicin del reactor en fraccin molar y el nmero de moles totales con la mayor cantidad de decimales posible estn en la siguiente tabla:
XCOF2 XC3F6 XCF3COF 0,09240549237648380 0,779926263571360 0,1276682440521560 Nmero de moles 0,117265023197110

Para poder comparar los datos con algo y sacar conclusiones en necesario ver los valores que se obtuvieron por gases ideales, eso est representado en la siguiente tabla:
XCOF2 0,114280668285240 XC3F6 XCF3COF Nmero de moles 0,116441336033489

0,7642218213064280 0,1214975104083320

Los errores porcentuales calculados con la ecuacin 42, para el error porcentual, se muestran en la siguiente tabla para cada uno de estos datos.
XCOF2 23,67302565% XC3F6 2,013580386% XCF3COF 4,833413109% Nmero de moles 0,702415043%

De esto se puede obtener que aunque el error no es muy grande para el nmero de moles totales, cambia mucho para la fraccin molar de fluoruro de carbonilo, por lo cual fue utilitario acceder a este mtodos para calcular la real composicin del reactor, por otro lado se concluye que la inversin que se hizo en el reactor se puede ver ambiguamente como mala y buena, buena desde el punto que se adquirieron unas variables capaces de dejar el sistema determinado, y mala porque no se obtuvo el compuesto de inters. Por otro lado la conversin fue del 100% ya que el oxigeno presente a principio de la reaccin se consumi completamente. Sin embargo se compar el resultado con otros modelos termodinmicos ms robustos y un modelo parecido que es el de Lee-Kesler-Plocker hecho en Aspen Properties, los resultados se recogen en la siguiente tabla comparando el factor de compresibilidad.

25

Modelo

Modelo de Lee-Kesler 0,99297585

Modelo de Lee-KeslerPlocker 0,9927552

Valor de Z

Modelo de Modelo de PR BenedictModelo de con WebbSchwartzentrubermodificacin de RubinRenon HYSYS Starling 0,9906926 0,9950574 0,9910014

En esta tabla se puede ver las variaciones en el clculo del factor de compresibilidad de la mezcla por diferentes caminos, sin embargo los que ms se ajustara realmente son el modelo de Schwartzentruber-Renon y el modelo de Peng-Robinson con la modificacin de HYSYS, ya que el modelo SR, es para mezclas de sustancias polares y no ideales, el modelo de HYSYS, tiene un impedimento y es que es funciona en las mismas condiciones que el modelo de Lee-Kesler, el de LeeKesler-Plocker y el de Benedict-Webb-Rubin-Starling, este impedimento es que funciona con mezcla poco polares o apolares, el modelo de SchwartzentruberRenon funciona con mezclas apolares pero tambin tiene un problema, es que necesita datos de interaccin binaria experimentales, por lo cual tambin hay posibilidades de error, ya que el simulador calculo el factor de compresibilidad sin estos datos.

26

Bibliografa
1. Smith Van Ness. Introduccin a la Termodinmica en Ingeniera Qumica. Quinta edicin. Mxico: McGRAW-HILL, 1996. 2. Georgios M. Kontogeorgis, Georgios K. Folas. Thermodynamic Models for Industrial Applications: From Classical and Advanced Mixing Rules to Association Theories. Primera Edicin. Gran Bretaa: WILEY, 2010. 3. Uzi Mann. Principles Of Chemical Reactor Analysis And Design: New Tools for Industrial Chemical Reactor Operations. Segunda Edicin. Estados Unidos de America: WILEY, 2009. 4. Bruce E. Poling, John M. Prausnitz, John P. OConnell. The Properties Of Gases And Liquids. Quinta Edicin. Estados Unidos de America: McGRAWHILL, 2001. 5. John M. Prausnitz, Rdinger N. Lichtenthaleer, Edmundo Gomes de Azevedo. Termodinmica Molecular de los Equilibrios de Fases. Tercera Edicin. Madrid, Espaa: PRETINCE HALL, 2000. 6. Base de datos del simulador comercial Aspen Properties. 7. Base de datos de ChemSpider, http://www.chemspider.com/.

27