Determinación del equilibrio líquido-vapor de agua, aromáticos y sus mezclas mediante simulación molecular

Resumen

La simulación molecular presenta la ventaja de ofrecer un marco teórico importante para predecir propiedades termodinámicas y de transporte de fluidos con aplicaciones industriales. En este trabajo, se explotó está ventaja para predecir el equilibrio líquido vapor de agua, compuestos aromáticos y sus mezclas a condiciones tanto sub- como supercríticas. Se realizó una comparación de diferentes potenciales intermoleculares conocidos mediante el cálculo de propiedades termodinámicas de sistemas puros que sirvió de punto de partida para llevar a cabo una optimización de parámetros transferibles para un potencial intermolecular de agua y compuestos aromáticos. En el caso de agua, se llevo a cabo el análisis y evaluación de propiedades termodinámicas de un modelo simple de agua. En este modelo, las contribuciones electrostáticas se aproximan mediante interacciones de corto alcance en vez de las típicas fuerzas de Coulomb de largo alcance. En general, se han encontrado grandes desviaciones con respecto a los datos experimentales, tal como un valor de temperatura crítica de 360K, valor 50% alejado del valor experimental. Debido a que estos resultados nos indican la importancia de incluir las fuerzas de Coulomb en el modelo molecular empleado para reproducir correctamente las propiedades de agua, el trabajo de investigación se ha enfocado en la optimización de los parámetros de los potenciales TIP4P y SPC/E. Los resultados obtenidos muestran que es posible encontrar una mejor aproximación al punto crítico experimental a partir de la optimización del modelo SPC/E. Sin embargo, el buen acuerdo con los experimentos del modelo original a condiciones ambiente se pierde usando los parámetros del modelo optimizado. Por otro lado, la estimación de propiedades de compuestos aromáticos esta de acuerdo con los datos experimentales permitiendo la reproducción de la densidad de líquido saturado, presión de saturación y entalpía de vaporización para compuestos puros. Finalmente, en el caso de mezclas se ha aplicado el conjunto de parámetros obtenidos para aromáticos. Las propiedades termodinámicas de la mezcla binaria aromático-aromático y aromático agua son analizadas en un amplio rango de temperaturas y presiones. Las desviaciones encontradas entre los valores calculados y los experimentales sugieren aplicar un mejor método de optimización para sistemas puros o por otro lado, promover un potencial de interacción intermolecular más sofisticado. Las estimaciones a condiciones cercanas al punto crítico están en buen acuerdo con los datos experimentales. RESUM La simulació molecular presenta l'avantatge d'oferir un marc teòric important per a cercar propietats termodinàmiques i de transport de fluids amb aplicacions industrials. En aquest treball es va explotar aquesta avantatge per predir l'equilibri líquid vapor d'aigua, components aromàtics i les seves mescles, tant a condicions sub com supercrítiques. Es va realitzar una comparació de diferents potencials intermoleculars, coneguts mitjançant el càlcul de propietats termodinàmiques de sistemes purs, que ha servit de punt de sortida per portar a terme una optimització de paràmetres transferibles per a un potencial intermolecular de propietats termodinàmiques d'un model simple d'aigua. En aquest model, les contribucions electrostàtiques s'aproximen mitjançant interaccions de curt abast en lloc de les típiques forces de Coulomb de llarg abast. En general, s'ha trobat grans desviacions respecte a les dades experimentals, tal com un valor de temperatura crítica de 360K, valor 50% allunyat del valor experimental. Degut a que aquests resultats ens indiquen la importància d'incloure les forces de Coulomb en el model molecular emprat per reproduir correctament les propietats d'aigua, el treball d'investigació s'ha enfocat en l'optimització dels paràmetres dels potencials TIP4P i SPC/E. Els resultats obtinguts mostren que és possible trobar una millor aproximació al punt crític experimental a partir de l'optimització del model SPC/E. No obstant, el bon acord amb els experiments del model original a condicions ambientals es perden usant els paràmetres del model optimitzat. Per altre banda, l'estimació de propietats de compostos aromàtics esta d'acord amb les dades experimentals permetent la reproducció de la densitat de líquid saturat, pressió de saturació i entalpia de vaporització per a compostos purs mitjançant el potencial AUA-Aromátics proposat. Finalment, en el cas de mescles s'ha aplicat el conjunt de paràmetres obtinguts per aromàtics. Les propietats termodinàmiques de la mescla binària aromàtic-aromàtic i aromàtic-aigua són analitzades en un ample rang de temperatures i pressions. Les desviacions trobades entre els valors calculats i els experimentals suggereixen aplicar un millor mètode d'optimització per a sistemes purs o, per altre banda, promoure un potencial d'interacció intermolecular més sofisticat. Les estimacions en condicions properes al punt crític tenen un bon acord amb les dades experimentals.