Pontificia Universidad Católica de Chile

Escuela de Ingenierı́a
Departamento de Ingenierı́a Quı́mica y Bioprocesos

Proyecto Final Fisicoquı́mica IIQ2043

Josefina Bizzarri
Constanza Farı́as
Aarón Velásquez

19 de Julio de 2021
Contenidos
1 Introducción 2

2 Equilibio Lı́quido-Vapor 3
2.1 Obtención de datos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2 Descripción del procedimiento experimental utilizado . . . . . . . . . . . . . . . . . . . . . . . 3
2.3 Descripción del modelo a utilizar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.4 Ajuste de parámetros con NRTL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.5 Análisis de consistencia termodinámica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.6 Análisis de error del modelo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.7 Utilidad de los datos experimentales seleccionados en algún proceso a escala industrial . . . . 7
2.8 Anexos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7

3 Equilibrio Sólido-Lı́quido 10
3.1 Obtención de datos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.2 Descripción del procedimiento experimental utilizado . . . . . . . . . . . . . . . . . . . . . . . 10
3.3 Descripción del modelo a utilizar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
3.4 Ajuste de parámetros con Wilson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
3.5 Análisis de error de los modelos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
3.6 Utilidad de los datos experimentales seleccionados en algún proceso a escala industrial . . . . 13
3.7 Anexos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

4 Conclusión 17

5 Bibliografı́a 18

1
1 Introducción
Para poder generar un modelamiento, o diseño de un proceso, es de suma importancia poseer datos pun-
tuales sobre las propiedades termodinámicas involucradas en el sistema, de esta forma, es posible desarrollar
modelos de equilibrio lı́quido-vapor, sólido-lı́quido, entre otros.

Los sistemas multicomponentes añaden una nueva dimensión al comportamiento de las fases de una mez-
cla. Los compuestos puros están rodeadas de moleculas iguales, pero en las mezclas pueden ser similares
como muy distintas generando distintas estructuras y propiedades fı́sicas. En estas mezclas pueden ocurrir
separaciones que permiten llevar a la mezcla a un estado donde distintas fases con distintas composiciones
coexisten. En termodinámica, el equilibrio lı́quido-vapor describe la distribución de una especie quı́mica
entre la fase de vapor y una fase lı́quida (Canales, 2021). Los sistemas de equilibrio lı́quido-vapor son de
gran utilidad en la ingenierı́a debido a su cantidad de usos en la industria para separación o mezcla de
sustancias. El equilibrio solido-lı́quido también es de gran utilidad en la Ingenierı́a porque permite generar
fases de equilibrio entre sólido y lı́quido y al mismo tiempo separar estos compuestos (Canales, 2021).

La realización de este proyecto resulta ser relevante para el estudio de equilibrio lı́quido-vapor y sólido-lı́quido.
Con la finalidad de obtener predicciones suministradas desde un modelo termodinámico que intentan acer-
carse a los resultados experimentales previamente extraı́dos. Los modelos a tratar en este proyecto son
NRTL (Non-randon Two Liquid) y Wilson, también se utilizará la ecuación de estado Peng-Robinson. Es-
tos modelos en general, requieren de una vasta cantidad de parámetros, ya que la determinación de estos,
suele ser el problema esencial para el adecuado modelamiento de los sistemas. Para la mezcla lı́quido-vapor
se eligió una mezcla de metanol-benceno que forma un azeotrópo, donde se comporta como si estuviera
formada por un solo componente. Para la mezcla de solido-lı́quido se eligieron los compuestos naftaleno
y p-diclorobenceno ya que el sólido p-diclorobenceno puro, puede coexistir en equilibrio con una solución
lı́quida de naftaleno disuelto en p-diclorobenceno. Estos compuestos no se combinan quı́micamente, no son
miscibles en estado sólido, pero son solubles al estado lı́quido (Bartolomé y Crespo, 2019), lo que permite
analizar sus propiedades fı́sicas y la importancia de esta mezcla.

Los datos del equilibrio L-V fueron correlacionados con ecuaciones basadas en el concepto de composiciones
locales, ası́ como, con la ecuación de estado adjunto con término quı́mico (EoSC) propuesto por Góral,
EoSC utiliza la ecuación de estado de Redlich-Kwong adjuntada con el término de asociación. El EoSC, que
tiene un parámetro binario ajustable y es capaz de correlacionar los sistemas con mucha presición (Góral et
al,2004). Los datos del equilibrio S-L se recolectaron utilizando calorimetrı́a diferencial de barrido (DSC).
Además, los coeficientes de actividad de los componentes se obtuvieron mediante la expresión del parámetro
de solubilidad de Scatchard-Hildebrand para investigar el rendimiento de este (Góral et al, 2004).

Los objetivos principales de este trabajo, radican en el estudio y la estimación de datos, tales como, co-
eficientes de actividad, fugacidad, presión calculada, parámetros de NRTL, y consistencia termodinámica
del equilibrio lı́quido-vapor para la mezcla de benceno y metanol. Luego para el equilibrio sólido-lı́quido de
naftaleno y p-diclorobenceno se buscó principalmente ajustar parámetros de Wilson con el fin de comprender
y analizar si este modelo se adapta o no a los datos experimentales y como lo hace.

2
2 Equilibio Lı́quido-Vapor
2.1 Obtención de datos
Los datos relacionados a la mezcla lı́quido-vapor se obtuvieron de la página de NIST Thermolit, donde se
encontró un paper relacionado a la mezcla de metanol y benceno. Este paper queda citado como:

Góral, M. Skrzecz, A. Bok, A. Maczyński,A and Oracz,

(Vapor–Liquid) Equilibrium Data. Part 3. Binary Alkanol–Aromatic Hydrocarbon Systems
Journal of Physical and Chemical Reference Data., 2004, 33, 959-962
Temperature: 338.33-353.22 K
Pressure: 101.32 kPa
DOI: https://doi.org/10.1063/1.1760761

2.2 Descripción del procedimiento experimental utilizado

El procedimiento experimental radicó en combinar las pruebas de consistencia termodinámica, correlación
de datos, comparación con entalpı́a de datos de mezcla y comparación de datos VLE para varias mezclas.
Las correlaciones de los datos se realizaron con una ecuación de estado con un término quı́mico (EoSC),
método desarrollado por Góral, (Góral et al,2004). Esta ecuación de estado, solo posee un parámetro binario
ajustable θ, el cual después se relaciona con energı́a de Gibbs de exceso a concentración equimolar:

Q = GE
0.5 /RT

Luego, Q depende de una mezcla sencilla del calor de la mezcla, la cual se describe por la ecuación de
Gibbs-Helmholtz. Lo que posteriormente sirve como predicción de ELV el cual finalmente se reduce a la
predicción de la presión experimental en base a la calculada, la cual se obtiene con N, n y m, donde N es
el número total de puntos experimentales en el conjunto de datos, n es el número de puntos de datos para
sustancias puras y m es el número de parámetros ajustables. En el caso de correlación se utilizó m=1; en el
caso de la predicción se utilizó m = 0 (Góral et al,2004).
PN −n
σ = [ k=1 ((Pexper − Pcalc )2 )/(N − n − m)]0.5

2.3 Descripción del modelo a utilizar

En esta sección se utiliza el modelo NRTL. Su principal supuesto es el concepto de composición local, que
muestra que para mezclas moleculares con fuertes interacciones no ideales, se agrupan de manera no aleato-
ria, donde la composición es diferente a la composición global debido a la naturaleza fı́sica de estas moléculas.
Este concepto es ayudado por la teorı́a de dos lı́quidos de Scott, que establece que hay dos tipos de moléculas
en una mezcla binaria, tipo 1 y tipo 2. La energı́a de Gibbs residual de tipo 1 del centro es igual a la suma
de la energı́a de Gibbs residual de las interacciones experimentadas por la molécula 1 en el centro siendo gij .
Ocurre de la misma manera con las moléculas de tipo 2 en el centro. La ecuación NRTL para el exceso de
energı́a de Gibbs se muestra a continuación,(Buentello,2007):

GE /RT = (x1 ∗ x2 ∗ [(τ12 ∗ G12 )/(x1 ∗ G12 + x2 ) + (τ21 ∗ G21 )/(x1 + x2 ∗ G21 ))] (1)

En donde:

τ12 = (g12 − g22 )/RT (2)

τ21 = (g21 − g11 )/RT (3)

G12 =Exp(-α12 *τ 12) (4)

3
G21 =Exp(-α21 *τ 21) (5)

α12 = α21 (6)

gij representa un parámetro energético de interacción entre las moléculas con subı́ndice i-j. El parámetro
αij está relacionado con las distribución no aleatorias de las moléculas en la mezcla. La ecuación de NRTL
para el caso binario tiene 3 parámetros, se ha indicado que α varı́a entre 0.2 y 0.47. Luego se deducen las
ecuaciones para los coeficientes de actividad obteniendo:

γ1 = Exp (x22 ∗ [(τ12 ∗ G12 )/(x1 ∗ G12 + x2 )2 + τ21 ∗ (G21 /(x1 + x2 ∗ G21 ))2 ]) (7)

γ2 = Exp (x21 ∗ [(τ ∗ G21 )/(x2 ∗ G21 + x1 )2 + τ12 ∗ (G12 /(x2 + x1 ∗ G12 ))2 ]) (8)

Para sistemas que se desvı́an de los sistemas ideales y los sistemas parcialmente miscibles, si los parámetros
ajustables se obtienen cuidadosamente, esta ecuación puede representar bien los datos experimentales.

2.4 Ajuste de parámetros con NRTL

Para realizar el respectivo ajuste, en primer lugar se hace una tabla con los datos experimentales extraı́dos
de la literatura que se encuentran en la ”Tabla 1” en anexos, se obtienen los valores de los parámetros de
Antoine de la literatura. Para el metanol: A = 5,20409, B = 1581,341 y C = -33,5 con los diversos valores
de temperatura involucrados.
Estos parámetros se reemplazaron en la siguiente ecuación de Antoine:

B
log10 (P sat/100)(kP a) = (A − T +C )

Posteriormente, para el benceno se obtienen los parámetros A = 14,1603, B = 2948,78 y C = -44,563

los cuales se reemplzan en la respectiva ecuación de Antoine:

B
lnP sat(kP a) = (A − T +C )

Luego se escoge α = 0,47 para ajustar los parámetros y se dejan libres τ 12 y τ 21 y se calculan los G12
y G21 con las fórmulas mencionadas anteriormente.

En consecuencia, se calculan γ1 y γ2 con toda la información obtenida anteriormente, para después calcular
las composiciones de vapor (yi ) y presiones calculadas. Por último se realiza una función objetivo con la
sumatoria (P − P calc )2 , donde P corresponde a la presión del sistema igual a 101,32 kPa y Pcal corresponde
a la presión calculada con las distintas presiones de saturación obtenidas dependientes de las temperaturas
experimentales. Finalmente a la función objetivo se le aplica solver de Excel para minimizar este valor
cambiando los parámetros τ 12 y τ 21.

Se obtienen los siguientes resultados y se compara el ajuste de los datos a los experimentales en el siguiente
gráfico. τ 12 = 1,0409, τ 21 = 1,82035, G12 = 0,6131 y G21 = 0,4250

4
Gráfico N°1: temperatura vs composiciones metanol del modelo NRTL a presión constante
101.32 KPa

Del gráfico NRTL se logra apreciar el correcto ajuste de estos datos con el modelo NRTL.

Reporte de calculos de coeficientes de de fugacidad de la mezcla en la fase vapor, coeficiente de fugacidad de

saturacion del compuesto puro y factor de Poynting:

Para realizar estos cálculos se utiliza la ecuación de Peng-Robinson, las presiones de saturación para cada
temperatura, la presión en [MPa] del sistema, las temperaturas experimentales y la constante de los gases
”R” (8,314472 (cm3MPa/molK)).

Se calcula la fugacidad en fase vapor para metanol y benceno para cada temperatura experimental con la
ayuda de la planilla ”Preos.xlsx” que se utilizó en clases de fisicoquı́mica, (Canales, 2021), para esto debemos
dejar fija la presión (en MPa), cambiar la temperatura y copiar el valor de la fugacidad. Después con la
definición se calculan los coeficientes de fugacidad.

De la tabla 2 y 3 del anexo podemos observar que para el caso en particular de la fase vapor los coeficientes
de fugacidad de ambos compuestos realmente tienden a 1 como fue supuesto.

Para calcular los coeficientes de fugacidad de saturación del compuesto puro se debe utilizar las presiones
de saturación de cada compuesto con su respectiva temperatura, se van cambiando estos dos valores y se
obtiene tanto la fugacidad de saturación como el volumen de saturación lı́quido.

Luego se calcula el coeficiente de fugacidad de saturación como la división entre la fugacidad de saturación
y la presión de saturación. De las tablas 2 y 3 del anexo se observa que estos coeficientes también tienden a 1.

5
Para calcular el factor de poynting de cada compuesto se utiliza la siguiente ecuación:

P OYi = EXP [ViL,sat ∗ (P − Pisat )/(R ∗ T )]

Para esto se utiliza el volumen de saturación calculado a cada temperatura y presión de saturación, la presión
del sistema, la presión de saturación, nuestra constante de los gases ”R” y las temperaturas. Luego de cal-
cular el factor de poynting para cada compuesto se puede observar en las tablas 2 y 3 del anexo que en todos
los casos este factor tiende a 1 por lo que queda demostrado que las suposiciones que se hicieron en un inicio
estaban en lo correcto.

2.5 Análisis de consistencia termodinámica

Para evaluar la consistencia de los datos experimentales se utiliza la siguiente fórmula:
Z 1
ln(γ1 /γ2 ) · dx1 = 0 (1)
0

Se usan los coeficientes de actividad obtenidos anteriormente y se calcula el valor de ln(γ1 /γ2 ) para todos
los datos experimentales. Una vez realizado este paso se arma un gráfico ln(γ1 /γ2 ) vs x1 para analizar su
resultado.

Gráfico N°2: consistencia termodinámica

Se ajusta una ecuación polinómica de tercer grado ya que representa de manera correcta la lı́nea del gráfico
con un R2 = 0.9996 obteniendo y = -4.5115 x3 +7.9819 x2 -7.8251 x + 2.3077. Luego se calcula la integral
sabiendo que ln(γ1 /γ2 ) = a x3 + b x2 + c x + d y al integrar entre 0 y 1 se obtiene a/4 + b/3 +c/2 + d
y al reemplazar los valores se obtiene como resultado -0.072091667 siendo un valor muy cercano a 0 por lo
que se concluye que los datos si son consistentes termodinamicamente.

2.6 Análisis de error del modelo

El análisis de error del modelo consistió en comparar las presiones experimentales con las calculadas y las
composiciones de vapor experimentales y calculadas. Luego, para esto se utilizó la fórmula entregada:

6
 Pn
(100/n) ∗ i=1 |(θiexp − θicalc )/(θicalc )|
El error obtenido para el caso de las composiciones de vapor, resultó ser de 2, 281%. Con esto es posible
concluir que en general NRTL, resulta ser un buen modelo para esta mezcla de benceno y metanol, dado
que el porcentaje de error es mı́nimo.

2.7 Utilidad de los datos experimentales seleccionados en algún proceso a escala

industrial
En la industria petroquı́mica de China ha sido de utilidad la sı́ntesis de distintos productos quı́micos, como
combustibles, explosivos, productos agrı́colas, disolventes, entre otros. Para ello ha sido requerido el uso
de tolueno el cual se ha obtenido por medio de mezclas de benceno y metanol. Por lo tanto, actualmente,
la alquilación del benceno con metanol para producir tolueno ha sido una opción muy anhelada, dada la
disminución de los recursos de petróleo crudo (Peng Dong et al, 2018). Por otra parte, las mezclas liquido
vapor de benceno y metanol han servido para dilucidar los potenciales peligros cuando se mezclan diferentes
productos quı́micos inflamables, para ello se han hecho estudios en la industria en varias proporciones de
mezcla utilizando recipientes de explosión esférico (Yi-Min Chang et al,2005). Por tanto, en general las
mezclas de benceno y metanol, ha servido para desarrollar la industria petroquı́mica, agrı́cola y además para
estudiar las posibles consecuencias relacionadas con los explosivos y los peligros que contraen los quı́micos
inflamables.

2.8 Anexos

”Tabla 1: Datos experimentales de metanol y benceno”

7
”Tabla 2: coeficientes y factor de pointing para metanol”

”Tabla 3: coeficientes y factor de pointing para benceno”

8
”Tabla 4: Planilla 1 de cálculos utilizada”

”Tabla 5: Planilla 2 de cálculos utilizada”

9
3 Equilibrio Sólido-Lı́quido
3.1 Obtención de datos
Los datos relacionados a la mezcla sólido-lı́quido se obtuvieron de la página de NIST Thermolit donde se
encontró un paper relacionado a la mezcla de naftaleno y p-diclorobenceno. Este paper queda citado como:

Wei, D.; Jin, K.

(Solid + liquid) equilibria of (naphthalene + isomeric dichlorobenzenes)
J. Chem. Thermodyn., 2009, 41, 145-149
Temperature: 305.15 K - 350.75 K
Pressure: 101 kPa
DOI: 10.1016/j.jct.2008.10.005

3.2 Descripción del procedimiento experimental utilizado

El procedimiento experimental consistió en utilizar calorimetrı́a diferencial de barrido para encontrar los
datos, mientras que los coeficientes de actividad se obtuvieron por la expresión de Scatchard-Hildebrand.
Este enfoque ofrece un procedimiento completo para estimar con buena precisión. Los productos quı́micos del
estudio fueron suministrados por Tianjin Kewei Chemical Reagent Co., con una pureza superior al 99,5%. Las
temperaturas y entalpias de fusión se midieron usando DSC un calorı́metro diferencial de barrido equipado
con la muestra y crisoles de aluminio estándar de referencia. Luego, gas nitrógeno de alta pureza se purgó a
través de la cámara de DSC. Antes del análisis, el calor del flujo se calibró utilizando el calor de fusión del
inicio. La muestra se calentó muy lentamente dentro de un tubo de vidrio a casi la temperatura de fusión
del componente principal, se solidificó con agitación continua. Se tomó una pequeña cantidad de sólido y se
selló en una muestra Mettler crisol para su posterior análisis.

3.3 Descripción del modelo a utilizar

En esta sección se utilizó el modelo de Wilson que dedujo su expresión para la energı́a Gibbs de exceso para
una mezcla binaria, en base a dos supuestos moleculares. En primer lugar supone que la energı́a libre de
mezclado tiene una relación parecida a la de la ecuación de Flory-Huggins, y en segundo lugar la distribución
de las moléculas alrededor de una molécula central siguen una relación llegando a las siguientes ecuaciones
(Buentello,2007):

GE /RT = −x1 (ln(x1 + x2 Λ12 ) − (x2 ln(x2 + x1 Λ21 )) (1)

ln(γ1 ) = − ln(x1 + x2 Λ12 ) + x2 ((Λ12 )/(x1 + x2 Λ12 ) − (Λ21 )/(x1 Λ21 + x2 )) (2)

ln(γ2 ) = − ln(x2 + x1 Λ21 ) − x1 ((Λ12 )/(x1 + x2 Λ12 ) − (Λ21 )/(x1 Λ21 + x2 )) (3)

La ecuación de Wilson para el caso binario tiene dos parámetros ajustables Λ12 y Λ21 que se relacionan
con los volúmenes molares de los componentes puros y con diferencias de energı́as formulando las siguientes
ecuaciones:

Λ12 = (v2 /v1) ∗ exp(−A12 )/RT ) (4)

Λ21 = (v1 /v2) ∗ exp(−A21 )/RT ) (5)

La ecuación de Wilson proporciona una buena representación de energı́a Gibbs exceso para muchas mezclas
miscibles, siendo muy útil para mezclas polares debido a las suposiciones moleculares que representan bien
este tipo de interacciones.

10
3.4 Ajuste de parámetros con Wilson
Para realizar el ajuste de parámetros se utilizan las siguientes fórmulas

ln(γ1 x1 ) = (−∆Hf us,1 /R) ∗ (1/T − 1/Tm,1 ) (6)

ln(γ2 x2 ) = (−∆Hf us,2 /R) ∗ (1/T − 1/Tm,2 ) (7)

Para ajustar el modelo de Wilson lo primero que se hizo fue calcular coeficientes de actividad experimentales,
para esto se obtuvieron los valores de entalpı́a de fusión y las temperaturas de fusión para ambos compuestos
desde el paper de referencia:

Tabla N°1: datos de fusión de los componentes

Teniendo en conocimiento estos ultimos datos y las composiciones experimentales se pueden despejar los
coeficientes de actividad de las ecuaciones (6) y (7) para ası́ poder calcular los coeficientes de actividad ex-
perimentales de cada composición y de esta forma poder realizar el ajuste.(T= Tf us,mezcla , R= 8.314 J/molK)

Una vez obtenidos los coeficientes de actividad experimentales para cada compuesto se utilizó la herramienta
”Solve” de excel para poder ajustar los parámetros. La función objetivo utilizada fue la sumatoria del valor
de la siguiente expresión para cada dato experimental:

(ln(γ1exp ) − ln(γ1calc.wilson ))2 + (ln(γ2exp ) − ln(γ2calc.wilson ))2

En este caso los valores del logaritmo natural de los coeficientes de actividad calculados con Wilson se dejan
expresados como la parte derecha de las ecuaciones (2) y (3) en función de los parámetros ajustables. Luego
de realizar algunas veces la función solver en excel se obtuvo que los parámetros ajustables son:

Λ12 = 0, 106326407, Λ21 = 0, 250749306

Una vez obtenidos estos parámetros se calcularon nuevos coeficientes de actividad para cada composición
y cada compuesto (con el modelos de Wilson), luego utilizando las ecuaciones (6) y (7) se despeja la com-
posición para luego obtener el nuevo valor calculado a partir del modelo ajustable.

Se sabe que no todos los datos de composiciones obtenidas con ambas ecuaciones (6) y (7), tienen sentido
fı́sico por lo que se procede a graficar todas las composiciones del compuesto 1 para ası́ encontrar el punto
eutéctico y descartar las composiciones que se encuentren por debajo de este.

11
Gráfico N°3: gráfico sólido lı́quido (P const= 101 Kpa)

Finalmente este es el gráfico que obtenemos al utilizar esta forma de ajustar y se observa que aunque se
trata de un modelo ajustable los datos no estan del todo cercanos ya que se encuentra desplazado. El punto
eutectico generado por los datos calculados es de una temperatura de aproximadamente 309 K y 0.37 de
composición de naftaleno.

Por el mal ajuste que se obtuvo, se decide cambiar de función objetivo, cambiamos los coeficientes de activi-
dad por las composiciones del compuesto 1, la nueva función objetivo queda de la siguiente forma:

F O := (X1exp − X1calc )2 + (X2exp − X2calc )2

Para este caso la expresión de las composiciones calculadas queda más compleja ya que se debe dejar expre-
sada las ecuaciones (2) y (3) para ası́ obtener un coeficiente de actividad, esta expresión debe ser utilizada
en las ecuaciones (6) y (7), despejando la composición.

Se debe tener cuidado con los datos sin sentido fı́sico, por lo que hay que seguir cierto procedimiento: hacer
una primera iteración con todos los valores obtenidos con las ecuaciones de cada compuesto, luego graficar las
temperaturas versus las composiciones y descartar las composiciones que no tienen sentido fı́sico. Luego, con
los datos que si tienen sentido fı́sico se debe utilizar nuevamente solver para volver a ajustar los parámetros.

12
Gráfico N°4: Gráfico sólido lı́quido (P const= 101Kpa)

Con esta nueva forma de ajuste se observa que los datos calculados quedan exactos con los datos experimen-
tales y se puede decir que el modelo ajusta muy bien. Los nuevos parámetros son:

Λ12 = 0, 000168811, Λ21 = 0, 00090697

El punto eutectico generado por los datos calculados es de una temperatura de aproximadamente 305K y
0.36 de composición de naftaleno.

3.5 Análisis de error de los modelos

Para analizar el error del modelo se calculó la desviación absoluta promedio con la siguiente fórmula:
Pn exp
(100/n) ∗ i=1 (θi − θicalc )/(θicalc )

Se obtuvo un error del 16,64% respecto a los datos experimentales con los calculados. Se ajustaron los datos
en base a los coeficientes de actividad y esto puede llevar a unos parámetros menos exactos como sucedió en
el informe, se generó un alto porcetaje de error con este ajuste. Luego se ajustaron los datos en base a las
composiciones calculadas y experimentales donde se obtuvo un error del 0.029% lo que deja demostrado que
este ajuste es mucho mejor para los datos ya que no se alejan de la realidad.

3.6 Utilidad de los datos experimentales seleccionados en algún proceso a escala

industrial
Una utilidad de los datos experimentales de la mezcla naftaleno y p-diclorobenceno es la investigación ex-
perimental de los coeficientes de absorción ultrasonica que se ha medido en monocristales de naftaleno y
p-diclorobenceno a frecuencias de 5 a 95 MHz para el primer cristal y de 5 a 65 MHz para el segundo.
Debido a los puntos de fusión de estos dos compuestos que son relativamente altos, se pueden obtener
estas moléculas en monocristales por lo que fue posible la aplicación de una onda acústica. También se
realizaron estudios de otros mecanismos de absorción acustica en estos sólidos orgánicos como procesos de
amortiguación acústica. El estudio de estas moléculas complicadas permite determinar la naturaleza de la
absorción acústica distinta de la generada por resonancia molecular. (Yun y Beyer, 2004) Aportando grandes
conocimientos a la industria biotecnológica.

Otra aplicación es la microextracción en fase sólida de masas de p-diclorobenceno y naftaleno en la miel.

Este proceso ayudó a ver el grado de contaminación de los panales con estas sustancias ya que se encontraron

13
cristales de naftaleno y p-diclorobenceno en ellas, lo que indica el uso de panales precontaminados. Estos
dos compuestos son absorbidos facilmente por la cera del panal y terminan en la miel. Esto es sumamente
relevante a la hora de clasificar estas mieles como aptas para el consumo humano ya que estos plaguicidas
de naftaleno y p-diclorobenceno son altamente tóxicos y no pueden ser consumidos (Harizanis et al, 2008).
Aportando grandes conocimientos a la industria alimenticia.

3.7 Anexos

x1 Tf mezcla (K)
0 326.15
0.0737 321.95
0.15 317.05
0.1922 314.15
0.237 311.65
0.2563 310.65
0.3062 307.15
0.3438 305.15
0.4228 306.25
0.4831 312.35
0.559 320.15
0.626 326.65
0.6843 331.35
0.7311 335.05
0.7958 339.65
0.8432 343.55
0.9209 347.95
0.953 350.75
1 353.25

Table 1: Datos experimentales de naftaleno y p-diclorobenceno

”Tabla 2: Planilla 1 de cálculos utilizada con ajuste de coeficientes de actividad”

14
”Tabla 3: Planilla 2 de cálculos utilizada con ajuste de coeficientes de actividad”

”Tabla 4: Planilla 1 de cálculos utilizada con ajuste de composición actividad”

15
”Tabla 5: Planilla 2 de cálculos utilizada con ajuste de composición

16
4 Conclusión
En conclusión, la elección de los compuestos para ambas mezclas, lı́quido-vapor y sólido-lı́quido fueron real-
izadas en base a la capacidad de mezcla entre ellas. Luego, se realizaron ajustes a los datos experimentales
con diferentes modelos, como Wilson, NRTL y Peng Robinson, los cuales se caracterizan por poseer una
cantidad de parámetros mayor a otros modelos, como también, son reconocidos porque suelen funcionar de
forma acertiva para distintas mezclas y procesos a estudiar. Por otra parte, la elección de los compuestos,
también se vio influenciada por la gran utilidad que poseen en las industrias, dado que es de suma relevancia
estudiar los procedimientos que pueden resultar con impactos en el desarrollo de la ingenierı́a quı́mica. Por
ende, este trabajo tuvo como objetivo estudiar los modelos ya presentados los compuetsos involucrados y los
ajustes de ellos.

Para la mezcla L-V se ajustaron los parámetros de NRTL para analizar la precisión de los datos experimen-
tales. Para un valor de α = 0.47 se encontraron los siguentes valores, τ 12 = 1,0409, τ 21 = 1,8283, G12 =
0,6131 y G21 = 0,4250. También se reportaron los calculos de coeficientes de fugacidad de la mezcla en fase
de vapor y saturación del compuesto puro y el factor de poynting. Se obtuvo que para la fase de vapor los
coeficeintes de fugacidad de ambos compuestos tienden a 1. Los coeficientes de fugacidad de saturación de
los compuestos puros tienden a 1. Los cálculos realizados demuestran también que el factor de poynting
tiende a 1 demostrando que las supociciones que se hicieron al inicio estaban en lo correcto.

Para la mezcla S-L se ajustaron los parametros de Wilson para analizar la precisión de los datos experi-
mentales. Donde se realizó una función objetivo en base a las diferencias de los coeficientes de actividad
para ajustar las composiciones nuevas. Se obtuvieron valores para Λ12 = 0, 106326407, Λ21 = 0, 250749306
Se encontró que el punto eutectico de la mezcla es aproximadamente 308 K a una composición de 0.37 de
naftaleno. Al encontrarse un error muy grande se decidió ajustar en base a las composiciones donde se ob-
tuvieron valores Λ12 = 0, 00016881, Λ21 = 0, 00090697 obteniendo un punto eutéctico de aproximadamente
305 K a una composición de 0.37 de naftaleno.

Dentro de los resultados más importantes, se pueden destacar los porcentajes de error obtenidos en ambas
mezclas. En la mezcla lı́quido-vapor se obtuvo un 2, 281% de error en las composiciones de vapor, es decir,
el modelo el NRTL, fue una buena opción para ajustar los datos experimentales. Mientras que para la
mezcla sólido-lı́quido se obtuvo un error del 16,64% respecto a los datos experimentales con los calculados.
Se ajustaron los datos en base a los coeficientes de actividad y esto puede llevar a unos parámetros menos
exactos como sucedió en el informe, luego se ajustaron los datos en base a las composiciones calculadas y
experimentales donde se obtuvo un error del 0.029% lo que que permitió concluir que este ajuste es mucho
mejor para los datos ya que no se alejan de la realidad.
De la misma manera, para la mezcla lı́quido-vapor, se obtuvo una consistencia termodinámica cercana a 0,
especı́ficamente con un valor de -0.07209, lo que indica que los datos son efectivamente consistentes termod-
inamicamente.

Resulta importante dar énfasis en el uso de los componentes en la industria, ya que quedó en evidencia en este
trabajo, que las mezclas de benceno y metanol, son de utilidad para desarrollar la industria petroquı́mica,
como también, para dilucidar los posibles peligros involucrados a la hora de hacer estudios con quı́micos
inflamables. También las mezclas de naftaleno y p-diclorobenceno tienen relevancia en los estudios de ab-
sorción ultrasónica y acústica como también para identificarlos como compuestos tóxicos en la miel ayudando
a la industria alimenticia.

17
5 Bibliografı́a
Bartolomé, M. Crespo, R. (2019). Laboratorio de fisicoquómica II. Universidad Nacional Mayor de San Mar-
cos. Recuperado de:https://bit.ly/3zbCzcd

Buntello, G (2007). Análisis de modelos termodinámicos para su aplicación en la recuperación de proteı́nas

mediante sistemas de dos fases acuosas.. Instituto tecnológico y de estudis superiores de Monterrey. Recu-
perado de: https://bit.ly/3ksSXAO

Canales, R. (2021). Clase 6 fisicoquı́mica IIQ2043. Fugacidad. Pontificia Universidad Católica de Chile

Canales, R. (2021). Clase 7 fisicoquı́mica IIQ2043. Introducción a sistemas multicomponentes. Pontificia

Universidad Católica de Chile

Canales, R. (2021). Clases 19 fisicoquı́mica IIQ2043. Equilibrio Sólido-Lı́quido. Pontificia Universidad

Católica de Chile

Góral, M. Skrzecz, A. Bok, A. Maczyński,A and Oracz, (Vapor–Liquid) Equilibrium Data. Part 3. Binary
Alkanol–Aromatic Hydrocarbon Systems. Recuperado de: https://doi.org/10.1063/1.1760761

Harizanis, P. Alissandrakis, E. Tarantilis, A. Polissious, M. (2008) Solid phase microextraction, gas chro-
matographic, mass spectrometric analysis of p dichlorobenzene and naphthalene in honey. Food Additives
and Contaminants, 25:10, 1272-1277, DOI: 10.1080/02652030802007546. Recuperado de:
https://doi-org.pucdechile.idm.oclc.org/10.1080/02652030802007546

Peng Dong, Zeyu Li, Xiaorui Wang, Hongfei Yun Guixian Li (2018) The alkylation reaction of benzene with
methanol to produce toluene: Y-C and Y-CCs catalyst, Green Chemistry Letters and Reviews, 11:2, 158-164,
DOI: 10.1080/17518253.2018.1444794. Recuperado de:https://doi.org/10.1080/17518253.2018.1444794

Yi-Min Chang,Jo-Ming Tseng,Chi-Min Shu and Kwan-Hua Hu (2005)Flammability Studies of Benzene and
Methanol with Various Vapor Mixing Ratios at 150°C (*Doctoral Candidate, Graduate School of Engineering
Science and Technology, National Yunlin University of Science and Technology (NYUST)). Recuperado de:
https://link-springer-com.pucdechile.idm.oclc.org/content/pdf/10.1007/BF02705658.pdf

Yun, S. Beyer, R. G (2004) Ultrasonic Absorption Measurements in Single Crystals of Naphthalene and
pDichlorobenzene. The Journal of chemical physics, 40, 2538, DOI: 10.1063/1.1725560. Recuperado de:
https://doi.org/10.1063/1.1725560

18